Evolutionary Analysis of Transcriptional Regulation Mediated by Cdx2 in Rodents

Abstract

Differences in gene expression, which arise from divergence in cis-regulatory elements or alterations in transcription factors (TFs) binding specificity, are one of the most important causes of phenotypic diversity during evolution. On one hand, changes in the cis-elements located in the vicinity of target genes affect TF binding and/or local chromatin environment, thereby modulating gene expression in cis. On the other hand, alterations in trans-factors influence the expression of their target genes in a more pleiotropic fashion. Although the evolution of amino acid sequences is much slower than that of non-coding regulatory elements, particularly for the TF DNA binding domains (DBDs), it is still possible that changes in TF-DBD might have the potential to drive large phenotypic changes if the resulting effects have a net positive effect on the organism's fitness. If so, species-specific changes in TF-DBD might be positively selected. So far, however, this possibility has been largely unexplored. By protein sequence analysis, we observed high sequence conservation in the DBD of the TF caudal-type homeobox 2 across many vertebrates, whereas three amino acid changes were exclusively found in mouse Cdx2 (mCdx2), suggesting potential positive selection in the mouse lineage. Multi-omics analyses were then carried out to investigate the effects of these changes. Surprisingly, there were no significant functional differences between mCdx2 and its rat homologue (rCdx2), and none of the three amino acid changes had any impact on its function. Finally, we used rat-mouse allodiploid embryonic stem cells to study the cis effects of Cdx2-mediated gene regulation between the two rodents. Interestingly, whereas Cdx2 binding is largely divergent between mouse and rat, the transcriptional effect induced by Cdx2 is conserved to a much larger extent. There were no significant functional differences between mCdx2 and its rat homologue (rCdx2), and none of the three amino acid changes had any impact on its function. Moreover, Cdx2 binding is largely divergent between mouse and rat; the transcriptional effect induced by Cdx2 is conserved to a much larger extent.

Keywords

Cdx2 / cis and trans / evolution / gene regulation / mouse-rat allodiploid fusion cells

References

[1] V. Hinman and G. Cary, “The Evolution of Gene Regulation,” eLife 6 (2017): e27291.

[2] M. C. King and A. C. Wilson, “Evolution at Two Levels in Humans and Chimpanzees,” Science 188, no. 4184 (1975): 107–116.

[3] P. J. Wittkopp, B. K. Haerum, and A. G. Clark, “Evolutionary Changes in Cis and Trans Gene Regulation,” Nature 430, no. 6995 (2004): 85–88.

[4] S. A. Signor and S. V. Nuzhdin, “The Evolution of Gene Expression in Cis and Trans,” Trends in Genetics 34, no. 7 (2018): 532–544.

[5] A. Goncalves, S. Leigh-Brown, D. Thybert, et al., “Extensive Compensatory Cis-Trans Regulation in the Evolution of Mouse Gene Expression,” Genome Research 22, no. 12 (2012): 2376–2384.

[6] C. D. Meiklejohn, J. D. Coolon, D. L. Hartl, and P. J. Wittkopp, “The Roles of Cis- and Trans- Regulation in the Evolution of Regulatory Incompatibilities and Sexually Dimorphic Gene Expression,” Genome Research 24, no. 1 (2014): 84–95.

[7] G. A. Wray, M. W. Hahn, E. Abouheif, et al., “The Evolution of Transcriptional Regulation in Eukaryotes,” Molecular Biology and Evolution 20, no. 9 (2003): 1377–1419.

[8] G. Yvert, R. B. Brem, J. Whittle, et al., “Trans-Acting Regulatory Variation in Saccharomyces cerevisiae and the Role of Transcription Factors,” Nature Genetics 35, no. 1 (2003): 57–64.

[9] J. K. Pickrell, J. C. Marioni, A. A. Pai, et al., “Understanding Mechanisms Underlying Human Gene Expression Variation With RNA Sequencing,” Nature 464, no. 7289 (2010): 768–772.

[10] Q. Gao, W. Sun, M. Ballegeer, C. Libert, and W. Chen, “Predominant Contribution of Cis-Regulatory Divergence in the Evolution of Mouse Alternative Splicing,” Molecular Systems Biology 11, no. 7 (2015): 816.

[11] J. Hou, X. Wang, E. McShane, et al., “Extensive Allele-Specific Translational Regulation in Hybrid Mice,” Molecular Systems Biology 11, no. 8 (2015): 825.

[12] P. J. Wittkopp, B. K. Haerum, and A. G. Clark, “Regulatory Changes Underlying Expression Differences Within and Between Drosophila Species,” Nature Genetics 40, no. 3 (2008): 346–350.

[13] E. S. Wong, B. M. Schmitt, A. Kazachenka, et al., “Interplay of Cis and Trans Mechanisms Driving Transcription Factor Binding and Gene Expression Evolution,” Nature Communications 8 (2017): 1092.

[14] G. D. Bell, N. C. Kane, L. H. Rieseberg, and K. L. Adams, “RNA-seq Analysis of Allele-Specific Expression, Hybrid Effects, and Regulatory Divergence in Hybrids Compared With Their Parents From Natural Populations,” Genome Biology and Evolution 5, no. 7 (2013): 1309–1323.

[15] J. Chen, V. Nolte, and C. Schlötterer, “Temperature Stress Mediates Decanalization and Dominance of Gene Expression in Drosophila melanogaster,” PLoS Genetics 11, no. 2 (2015): e1004883.

[16] J. D. Coolon, C. J. McManus, K. R. Stevenson, B. R. Graveley, and P. J. Wittkopp, “Tempo and Mode of Regulatory Evolution in Drosophila,” Genome Research 24, no. 5 (2014): 797–808.

[17] J. J. Emerson, L. C. Hsieh, H. M. Sung, et al., “Natural Selection on Cis and Trans Regulation in Yeasts,” Genome Research 20, no. 6 (2010): 826–836.

[18] J. D. Gruber, K. Vogel, G. Kalay, and P. J. Wittkopp, “Contrasting Properties of Gene-Specific Regulatory, Coding, and Copy Number Mutations in Saccharomyces cerevisiae: Frequency, Effects, and Dominance,” PLoS Genetics 8, no. 2 (2012): e1002497.

[19] B. P. Metzger, F. Duveau, D. C. Yuan, S. Tryban, B. Yang, and P. J. Wittkopp, “Contrasting Frequencies and Effects of Cis- and Trans-Regulatory Mutations Affecting Gene Expression,” Molecular Biology and Evolution 33, no. 5 (2016): 1131–1146.

[20] B. Schaefke, J. J. Emerson, T. Y. Wang, M. Y. Lu, L. C. Hsieh, and W. H. Li, “Inheritance of Gene Expression Level and Selective Constraints on Trans- and Cis-Regulatory Changes in Yeast,” Molecular Biology and Evolution 30, no. 9 (2013): 2121–2133.

[21] H. M. Sung, T. Y. Wang, D. Wang, et al., “Roles of Trans and Cis Variation in Yeast Intraspecies Evolution of Gene Expression,” Molecular Biology and Evolution 26, no. 11 (2009): 2533–2538.

[22] A. Suvorov, V. Nolte, R. V. Pandey, S. U. Franssen, A. Futschik, and C. Schlötterer, “Intra-Specific Regulatory Variation in Drosophila pseudoobscura,” PLoS One 8, no. 12 (2013): e83547.

[23] D. Wang, H. M. Sung, T. Y. Wang, et al., “Expression Evolution in Yeast Genes of Single-Input Modules is Mainly due to Changes in Trans-Acting Factors,” Genome Research 17, no. 8 (2007): 1161–1169.

[24] B. P. H. Metzger, P. J. Wittkopp, and J. D. Coolon, “Evolutionary Dynamics of Regulatory Changes Underlying Gene Expression Divergence Among Saccharomyces Species,” Genome Biology and Evolution 9, no. 4 (2017): 843–854.

[25] J. J. Crowley, V. Zhabotynsky, W. Sun, et al., “Analyses of Allele-Specific Gene Expression in Highly Divergent Mouse Crosses Identifies Pervasive Allelic Imbalance,” Nature Genetics 47, no. 4 (2015): 353–360.

[26] V. J. Lynch, G. May, and G. P. Wagner, “Regulatory Evolution Through Divergence of a Phosphoswitch in the Transcription Factor CEBPB,” Nature 480, no. 7377 (2011): 383–386.

[27] D. Villar, P. Flicek, and D. T. Odom, “Evolution of Transcription Factor Binding in Metazoans–Mechanisms and Functional Implications,” Nature Reviews. Genetics 15, no. 4 (2014): 221–233.

[28] K. R. Nitta, A. Jolma, Y. Yin, et al., “Conservation of Transcription Factor Binding Specificities Across 600 Million Years of Bilateria Evolution,” eLife 4 (2015): e04837.

[29] G. Slodkowicz and N. Goldman, “Integrated Structural and Evolutionary Analysis Reveals Common Mechanisms Underlying Adaptive Evolution in Mammals,” Proceedings of the National Academy of Sciences of the United States of America 117, no. 11 (2020): 5977–5986.

[30] C. F. Olson-Manning, M. R. Wagner, and T. Mitchell-Olds, “Adaptive Evolution: Evaluating Empirical Support for Theoretical Predictions,” Nature Reviews. Genetics 13, no. 12 (2012): 867–877.

[31] M. A. S. Crissey, R. J. Guo, S. Funakoshi, J. Kong, J. Liu, and J. P. Lynch, “Cdx2 Levels Modulate Intestinal Epithelium Maturity and Paneth Cell Development,” Gastroenterology 140, no. 2 (2011): 517–528.e518.

[32] T. Ezaki, R.-J. Guo, H. Li, A. B. Reynolds, and J. P. Lynch, “The Homeodomain Transcription Factors Cdx1 and Cdx2 Induce E-Cadherin Adhesion Activity by Reducing β-and p120-Catenin Tyrosine Phosphorylation,” American Journal of Physiology. Gastrointestinal and Liver Physiology 293, no. 1 (2007): G54–G65.

[33] K. Aoki, Y. Tamai, S. Horiike, M. Oshima, and M. M. Taketo, “Colonic Polyposis Caused by mTOR-Mediated Chromosomal Instability in Apc+/Δ716 Cdx2+/− Compound Mutant Mice,” Nature Genetics 35, no. 4 (2003): 323–330.

[34] C. Bonhomme, I. Duluc, E. Martin, et al., “The Cdx2 Homeobox Gene has a Tumour Suppressor Function in the Distal Colon in Addition to a Homeotic Role During Gut Development,” Gut 52, no. 10 (2003): 1465–1471.

[35] A. Nishiyama, L. Xin, A. A. Sharov, et al., “Uncovering Early Response of Gene Regulatory Networks in ESCs by Systematic Induction of Transcription Factors,” Cell Stem Cell 5, no. 4 (2009): 420–433.

[36] D. Jana, H. T. Kale, and P. C. Shekar, “Generation of Cdx2-mCherry Knock-in Murine ES Cell Line to Model Trophectoderm and Intestinal Lineage Differentiation,” Stem Cell Research 39 (2019): 101521.

[37] S. Amin, R. Neijts, S. Simmini, et al., “Cdx and T Brachyury Co-Activate Growth Signaling in the Embryonic Axial Progenitor Niche,” Cell Reports 17, no. 12 (2016): 3165–3177.

[38] K. Peng, X. Li, C. Wu, et al., “Derivation of Haploid Trophoblast Stem Cells via Conversion In Vitro,” IScience 11 (2019): 508–518.

[39] X. Li, X. L. Cui, J. Q. Wang, et al., “Generation and Application of Mouse-Rat Allodiploid Embryonic Stem Cells,” Cell 164, no. 1–2 (2016): 279–292.

[40] D. Thybert, M. Roller, F. C. P. Navarro, et al., “Repeat Associated Mechanisms of Genome Evolution and Function Revealed by the Mus caroli and Mus pahari Genomes,” Genome Research 28, no. 4 (2018): 448–459.

[41] N. Kumar, Y. H. Tsai, L. Chen, et al., “The Lineage-Specific Transcription Factor CDX2 Navigates Dynamic Chromatin to Control Distinct Stages of Intestine Development,” Development 146, no. 5 (2019): dev172189.

[42] T. Cui, L. Jiang, T. Li, et al., “Derivation of Mouse Haploid Trophoblast Stem Cells,” Cell Reports 26, no. 2 (2019): 407–414.e405.

[43] M. Spivakov, “Spurious Transcription Factor Binding: Non-Functional or Genetically Redundant?,” BioEssays 36, no. 8 (2014): 798–806.

[44] D. T. Odom, R. D. Dowell, E. S. Jacobsen, et al., “Tissue-Specific Transcriptional Regulation has Diverged Significantly Between Human and Mouse,” Nature Genetics 39, no. 6 (2007): 730–732.

[45] E. T. Dermitzakis and A. G. Clark, “Evolution of Transcription Factor Binding Sites in Mammalian Gene Regulatory Regions: Conservation and Turnover,” Molecular Biology and Evolution 19, no. 7 (2002): 1114–1121.

[46] M. Ke, J. Liu, W. Chen, et al., “Integrated and Quantitative Proteomic Approach for Charting Temporal and Endogenous Protein Complexes,” Analytical Chemistry 90, no. 21 (2018): 12574–12583.

[47] R. C. Edgar, “MUSCLE: A Multiple Sequence Alignment Method With Reduced Time and Space Complexity,” BMC Bioinformatics 5 (2004): 113.

[48] N. P. Brown, C. Leroy, and C. Sander, “MView: A Web-Compatible Database Search or Multiple Alignment Viewer,” Bioinformatics 14, no. 4 (1998): 380–381.

[49] A. Dobin, C. A. Davis, F. Schlesinger, et al., “STAR: Ultrafast Universal RNA-seq Aligner,” Bioinformatics 29, no. 1 (2013): 15–21.

[50] Y. Liao, G. K. Smyth, and W. Shi, “featureCounts: An Efficient General Purpose Program for Assigning Sequence Reads to Genomic Features,” Bioinformatics 30, no. 7 (2014): 923–930.

[51] J. Feng, C. A. Meyer, Q. Wang, J. S. Liu, X. Shirley Liu, and Y. Zhang, “GFOLD: A Generalized Fold Change for Ranking Differentially Expressed Genes From RNA-seq Data,” Bioinformatics 28, no. 21 (2012): 2782–2788.

[52] M. I. Love, W. Huber, and S. Anders, “Moderated Estimation of Fold Change and Dispersion for RNA-seq Data With DESeq2,” Genome Biology 15, no. 12 (2014): 550.

[53] S. Chen, Y. Zhou, Y. Chen, and J. Gu, “fastp: An Ultra-Fast All-in-One FASTQ Preprocessor,” Bioinformatics 34, no. 17 (2018): i884–i890.

[54] B. Langmead and S. L. Salzberg, “Fast Gapped-Read Alignment With Bowtie 2,” Nature Methods 9, no. 4 (2012): 357–359.

[55] A. Tarasov, A. J. Vilella, E. Cuppen, I. J. Nijman, and P. Prins, “Sambamba: Fast Processing of NGS Alignment Formats,” Bioinformatics 31, no. 12 (2015): 2032–2034.

[56] Y. Zhang, T. Liu, C. A. Meyer, et al., “Model-Based Analysis of ChIP-Seq (MACS),” Genome Biology 9, no. 9 (2008): R137.

[57] S. Heinz, C. Benner, N. Spann, et al., “Simple Combinations of Lineage-Determining Transcription Factors Prime Cis-Regulatory Elements Required for Macrophage and B Cell Identities,” Molecular Cell 38, no. 4 (2010): 576–589.

[58] F. Ramírez, F. Dündar, S. Diehl, B. A. Grüning, and T. Manke, “deepTools: A Flexible Platform for Exploring Deep-Sequencing Data,” Nucleic Acids Research 42, no. Web Server issue (2014): W187–W191.

[59] M. Martin, “Cutadapt Removes Adapter Sequences From High-Throughput Sequencing Reads,” EMBnetjournal 17, no. 1 (2011): 10–12, https://doi.org/10.14806/ej171200.

[60] G. Yu, L. G. Wang, and Q. Y. He, “ChIPseeker: An R/Bioconductor Package for ChIP Peak Annotation, Comparison and Visualization,” Bioinformatics 31, no. 14 (2015): 2382–2383.

[61] A. R. Quinlan and I. M. Hall, “BEDTools: A Flexible Suite of Utilities for Comparing Genomic Features,” Bioinformatics 26, no. 6 (2010): 841–842.

[62] M. Wang and L. Kong, “pblat: A Multithread Blat Algorithm Speeding up Aligning Sequences to Genomes,” BMC Bioinformatics 20, no. 1 (2019): 28.

RIGHTS & PERMISSIONS

2025 The Author(s). Cell Proliferation published by Beijing Institute for Stem Cell and Regenerative Medicine and John Wiley & Sons Ltd.