Extensive Divergence of Transcription Factor Binding in Embryos with Highly Conserved Gene Expression

English version České info

To better characterize how variation in regulatory sequences drives divergence in gene expression, we undertook a systematic study of transcription factor binding and gene expression in blastoderm embryos of four species, which sample much of the diversity in the 40 million-year old genus Drosophila: D. melanogaster, D. yakuba, D. pseudoobscura and D. virilis. We compared gene expression, measured by mRNA-seq, to the genome-wide binding, measured by ChIP-seq, of four transcription factors involved in early anterior-posterior patterning. We found that mRNA levels are much better conserved than individual transcription factor binding events, and that changes in a gene's expression were poorly explained by changes in adjacent transcription factor binding. However, highly bound sites, sites in regions bound by multiple factors and sites near genes are conserved more frequently than other binding, suggesting that a considerable amount of transcription factor binding is weakly or non-functional and not subject to purifying selection.

Vyšlo v časopise: Extensive Divergence of Transcription Factor Binding in Embryos with Highly Conserved Gene Expression. PLoS Genet 9(9): e32767. doi:10.1371/journal.pgen.1003748
Kategorie: Research Article
prolekare.web.journal.doi_sk: https://doi.org/10.1371/journal.pgen.1003748

Souhrn

Zdroje

1. KingMC, WilsonAC (1975) Evolution at two levels in humans and chimpanzees. Science 188: 107–116.

2. CarrollSB (2008) Evo-devo and an expanding evolutionary synthesis: a genetic theory of morphological evolution. Cell 134: 25–36 doi:10.1016/j.cell.2008.06.030

3. LudwigMZ, KreitmanM (1995) Evolutionary dynamics of the enhancer region of even-skipped in Drosophila. Mol Biol Evol 12: 1002–1011.

4. MosesAM, PollardDA, NixDA, IyerVN, LiX-Y, et al. (2006) Large-scale turnover of functional transcription factor binding sites in Drosophila. PLoS Comput Biol 2: e130 doi:10.1371/journal.pcbi.0020130

5. RichardsS, LiuY, BettencourtBR, HradeckyP, LetovskyS, et al. (2005) Comparative genome sequencing of Drosophila pseudoobscura: chromosomal, gene, and cis-element evolution. Genome research 15: 1–18 doi:10.1101/gr.3059305

6. BalhoffJP, WrayGA (2005) Evolutionary analysis of the well characterized endo16 promoter reveals substantial variation within functional sites. Proc Natl Acad Sci USA 102: 8591–8596 doi:10.1073/pnas.0409638102

7. SchmidtD, WilsonMD, BallesterB, SchwaliePC, BrownGD, et al. (2010) Five-vertebrate ChIP-seq reveals the evolutionary dynamics of transcription factor binding. Science 328: 1036–1040 doi:10.1126/science.1186176

8. OdomDT, DowellRD, JacobsenES, GordonW, DanfordTW, et al. (2007) Tissue-specific transcriptional regulation has diverged significantly between human and mouse. Nat Genet 39: 730–732 doi:10.1038/ng2047

9. BornemanAR, GianoulisTA, ZhangZD, YuH, RozowskyJ, et al. (2007) Divergence of transcription factor binding sites across related yeast species. Science 317: 815–819 doi:10.1126/science.1140748

10. NiX, ZhangYE, NègreN, ChenS, LongM, et al. (2012) Adaptive evolution and the birth of CTCF binding sites in the Drosophila genome. PLoS Biol 10: e1001420 doi:10.1371/journal.pbio.1001420

11. ZhengW, ZhaoH, ManceraE, SteinmetzLM, SnyderM (2010) Genetic analysis of variation in transcription factor binding in yeast. Nature 464: 1187–1191 doi:10.1038/nature08934

12. KasowskiM, GrubertF, HeffelfingerC, HariharanM, AsabereA, et al. (2010) Variation in transcription factor binding among humans. Science 328: 232–235 doi:10.1126/science.1183621

13. BradleyRK, LiX-Y, TrapnellC, DavidsonS, PachterL, et al. (2010) Binding Site Turnover Produces Pervasive Quantitative Changes in Transcription Factor Binding between Closely Related Drosophila Species. PLoS Biol 8: e1000343 doi:10.1371/journal.pbio.1000343

14. WilsonMD, Barbosa-MoraisNL, SchmidtD, ConboyCM, VanesL, et al. (2008) Species-specific transcription in mice carrying human chromosome 21. Science 322: 434–438 doi:10.1126/science.1160930

15. McmanusCJ, CoolonJD, DuffMO, Eipper-MainsJ, GraveleyBR, et al. (2010) Regulatory divergence in Drosophila revealed by mRNA-seq. Genome research 20: 816–825 doi:10.1101/gr.102491.109

16. BrawandD, SoumillonM, NecsuleaA, JulienP, CsárdiG, et al. (2011) The evolution of gene expression levels in mammalian organs. Nature 478: 343–348 doi:10.1038/nature10532

17. KalinkaAT, VargaKM, GerrardDT, PreibischS, CorcoranDL, et al. (2010) Gene expression divergence recapitulates the developmental hourglass model. Nature 468: 811–814 doi:10.1038/nature09634

18. TiroshI, WeinbergerA, BezalelD, KaganovichM, BarkaiN (2008) On the relation between promoter divergence and gene expression evolution. Mol Syst Biol 4: 159 doi:10.1038/msb4100198

19. LinMF, DeorasAN, RasmussenMD, KellisM (2008) Performance and scalability of discriminative metrics for comparative gene identification in 12 Drosophila genomes. PLoS Comput Biol 4: e1000067 doi:10.1371/journal.pcbi.1000067

20. RussoCA, TakezakiN, NeiM (1995) Molecular phylogeny and divergence times of drosophilid species. Mol Biol Evol 12: 391–404.

21. LiX-Y, ThomasS, SaboPJ, EisenMB, StamatoyannopoulosJA, et al. (2011) The role of chromatin accessibility in directing the widespread, overlapping patterns of Drosophila transcription factor binding. Genome Biol 12: R34 doi:10.1186/gb-2011-12-4-r34

22. LiX-Y, MacArthurS, BourgonR, NixD, PollardDA, et al. (2008) Transcription factors bind thousands of active and inactive regions in the Drosophila blastoderm. PLoS Biol 6: e27 doi:10.1371/journal.pbio.0060027

23. LangmeadB, TrapnellC, PopM, SalzbergSL (2009) Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol 10: R25 doi:10.1186/gb-2009-10-3-r25

24. CapaldiAP, KaplanT, LiuY, HabibN, RegevA, et al. (2008) Structure and function of a transcriptional network activated by the MAPK Hog1. Nat Genet 40: 1300–1306 doi:10.1038/ng.235

25. HarrisonMM, LiX-Y, KaplanT, BotchanMR, EisenMB (2011) Zelda binding in the early Drosophila melanogaster embryo marks regions subsequently activated at the maternal-to-zygotic transition. PLoS Genet 7: e1002266 doi:10.1371/journal.pgen.1002266

26. ZhangY, LiuT, MeyerCA, EeckhouteJ, JohnsonDS, et al. (2008) Model-based analysis of ChIP-Seq (MACS). Genome Biol 9: R137 doi:10.1186/gb-2008-9-9-r137

27. Dewey CN (2006) Whole-genome alignments and polytopes for comparative genomics. GRADUATE DIVISION of the UNIVERSITY OF CALIFORNIA, BERKELEY.

28. PatenB, HerreroJ, BealK, FitzgeraldS, BirneyE (2008) Enredo and Pecan: genome-wide mammalian consistency-based multiple alignment with paralogs. Genome research 18: 1814–1828 doi:10.1101/gr.076554.108

29. FelsensteinJ (1973) Maximum-likelihood estimation of evolutionary trees from continuous characters. Am J Hum Genet 25: 471–492.

30. FelsensteinJ (1985) Phylogenies and the comparative method. American Naturalist 125: 1–15.

31. MartinsEP, GarlandTJr (1991) Phylogenetic analyses of the correlated evolution of continuous characters: a simulation study. Evolution 45: 534–557.

32. MattilaTM, BokmaF (2008) Extant mammal body masses suggest punctuated equilibrium. Proc Biol Sci 275: 2195–2199 doi:10.1098/rspb.2008.0354

33. ChaixR, SomelM, KreilDP, KhaitovichP, LunterGA (2008) Evolution of primate gene expression: drift and corrective sweeps? Genetics 180: 1379–1389 doi:10.1534/genetics.108.089623

34. BedfordT, HartlDL (2009) Optimization of gene expression by natural selection. Proc Natl Acad Sci USA 106: 1133–1138 doi:10.1073/pnas.0812009106

35. KaplanT, LiX-Y, SaboPJ, ThomasS, StamatoyannopoulosJA, et al. (2011) Quantitative models of the mechanisms that control genome-wide patterns of transcription factor binding during early Drosophila development. PLoS Genet 7: e1001290 doi:10.1371/journal.pgen.1001290

36. MacArthurS, LiX-Y, LiJ, BrownJB, ChuHC, et al. (2009) Developmental roles of 21 Drosophila transcription factors are determined by quantitative differences in binding to an overlapping set of thousands of genomic regions. Genome Biol 10: R80 doi:10.1186/gb-2009-10-7-r80

37. HarrisonPW, WrightAE, MankJE (2012) The evolution of gene expression and the transcriptome-phenotype relationship. Semin Cell Dev Biol 23: 222–229 doi:10.1016/j.semcdb.2011.12.004

38. KvonEZ, StampfelG, Yáñez-CunaJO, DicksonBJ, StarkA (2012) HOT regions function as patterned developmental enhancers and have a distinct cis-regulatory signature. Genes Dev 26: 908–913 doi:10.1101/gad.188052.112

39. GalloSM, GerrardDT, MinerD, SimichM, Soye DesB, et al. (2011) REDfly v3.0: toward a comprehensive database of transcriptional regulatory elements in Drosophila. Nucleic Acids Res 39: D118–D123 doi:10.1093/nar/gkq999

40. modENCODE Consortium (2010) RoyS, ErnstJ, KharchenkoPV, KheradpourP, et al. (2010) Identification of functional elements and regulatory circuits by Drosophila modENCODE. Science 330: 1787–1797 doi:10.1126/science.1198374

41. FisherWW, LiJJ, HammondsAS, BrownJB, PfeifferBD, et al. (2012) DNA regions bound at low occupancy by transcription factors do not drive patterned reporter gene expression in Drosophila. Proc Natl Acad Sci USA 109: 21330–21335 doi:10.1073/pnas.1209589110

42. RobertsA, PimentelH, TrapnellC, PachterL (2011) Identification of novel transcripts in annotated genomes using RNA-Seq. Bioinformatics 27: 2325–9 doi:10.1093/bioinformatics/btr355

43. LottSE, villaltaJE, SchrothGP, LuoS, TonkinLA, et al. (2011) Noncanonical Compensation of Zygotic X Transcription in Early Drosophila melanogaster Development Revealed through Single-Embryo RNA-Seq. PLoS Biol 9: e1000590 doi:10.1371/journal.pbio.1000590

44. HeQ, BardetAF, PattonB, PurvisJ, JohnstonJ, et al. (2011) High conservation of transcription factor binding and evidence for combinatorial regulation across six Drosophila species. Nat Genet 43: 414–420 doi:doi:10.1038/ng.808

45. SpivakovM, AkhtarJ, KheradpourP, BealK, GirardotC, et al. (2012) Analysis of variation at transcription factor binding sites in Drosophila and humans. Genome Biol 13: R49 doi:10.1186/gb-2012-13-9-r49

46. ConboyCM, SpyrouC, ThorneNP, WadeEJ, Barbosa-MoraisNL, et al. (2007) Cell cycle genes are the evolutionarily conserved targets of the E2F4 transcription factor. PLoS ONE 2: e1061 doi:10.1371/journal.pone.0001061

47. HeQ, BardetAF, PattonB, PurvisJ, JohnstonJ, et al. (2011) High conservation of transcription factor binding and evidence for combinatorial regulation across six Drosophila species. Nat Genet 43: 414–420 doi:10.1038/ng.808

48. HarrowJ, FrankishA, GonzalezJM, TapanariE, DiekhansM, et al. (2012) GENCODE: the reference human genome annotation for The ENCODE Project. Genome research 22: 1760–1774 doi:10.1101/gr.135350.111

49. LudwigMZ, PatelNH, KreitmanM (1998) Functional analysis of eve stripe 2 enhancer evolution in Drosophila: rules governing conservation and change. Development 125: 949–958.

50. LudwigMZ, PalssonA, AlekseevaE, BergmanCM, NathanJ, et al. (2005) Functional evolution of a cis-regulatory module. PLoS Biol 3: e93 doi:10.1371/journal.pbio.0030093

51. LudwigMZ, BergmanC, PatelNH, KreitmanM (2000) Evidence for stabilizing selection in a eukaryotic enhancer element. Nature 403: 564–567 doi:10.1038/35000615

52. SwansonCI, SchwimmerDB, BaroloS (2011) Rapid evolutionary rewiring of a structurally constrained eye enhancer. Curr Biol 21: 1186–1196 doi:10.1016/j.cub.2011.05.056

53. HareEE, PetersonBK, IyerVN, MeierR, EisenMB (2008) Sepsid even-skipped enhancers are functionally conserved in Drosophila despite lack of sequence conservation. PLoS Genet 4: e1000106 doi:10.1371/journal.pgen.1000106

54. CrockerJ, TamoriY, ErivesA (2008) Evolution acts on enhancer organization to fine-tune gradient threshold readouts. PLoS Biol 6: e263 doi:10.1371/journal.pbio.0060263

55. RheeHS, PughBF (2011) Comprehensive genome-wide protein-DNA interactions detected at single-nucleotide resolution. Cell 147: 1408–1419 doi:10.1016/j.cell.2011.11.013

56. SchaefferSW, BhutkarA, McAllisterBF, MatsudaM, MatzkinLM, et al. (2008) Polytene chromosomal maps of 11 Drosophila species: the order of genomic scaffolds inferred from genetic and physical maps. Genetics 179: 1601–1655 doi:10.1534/genetics.107.086074

57. DeweyCN (2012) Whole-genome alignment. Methods Mol Biol 855: 237–257 doi:_10.1007/978-1-61779-582-4_8

58. GuindonS, GascuelO (2003) A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol 52: 696–704.

59. TrapnellC, PachterL, SalzbergSL (2009) TopHat: discovering splice junctions with RNA-Seq. Bioinformatics 25: 1105–1111 doi:10.1093/bioinformatics/btp120

60. TrapnellC, WilliamsBA, PerteaG, MortazaviA, KwanG, et al. (2010) Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol 28: 511–515 doi:10.1038/nbt.1621

61. Team RC (n.d.) R: A Language and Environment for Statistical Computing. Available: http://www.R-project.org/.

62. Luengo HendriksCL, KeränenSVE, FowlkesCC, SimirenkoL, WeberGH, et al. (2006) Three-dimensional morphology and gene expression in the Drosophila blastoderm at cellular resolution I: data acquisition pipeline. Genome Biol 7: R123 doi:10.1186/gb-2006-7-12-r123

63. FowlkesCC, EckenrodeKB, BragdonMD, MeyerM, WunderlichZ, et al. (2011) A conserved developmental patterning network produces quantitatively different output in multiple species of Drosophila. PLoS Genet 7: e1002346 doi:10.1371/journal.pgen.1002346