A Machine Learning Approach for Identifying Novel Cell Type–Specific Transcriptional Regulators of Myogenesis

English version České info

Transcriptional enhancers integrate the contributions of multiple classes of transcription factors (TFs) to orchestrate the myriad spatio-temporal gene expression programs that occur during development. A molecular understanding of enhancers with similar activities requires the identification of both their unique and their shared sequence features. To address this problem, we combined phylogenetic profiling with a DNA–based enhancer sequence classifier that analyzes the TF binding sites (TFBSs) governing the transcription of a co-expressed gene set. We first assembled a small number of enhancers that are active in Drosophila melanogaster muscle founder cells (FCs) and other mesodermal cell types. Using phylogenetic profiling, we increased the number of enhancers by incorporating orthologous but divergent sequences from other Drosophila species. Functional assays revealed that the diverged enhancer orthologs were active in largely similar patterns as their D. melanogaster counterparts, although there was extensive evolutionary shuffling of known TFBSs. We then built and trained a classifier using this enhancer set and identified additional related enhancers based on the presence or absence of known and putative TFBSs. Predicted FC enhancers were over-represented in proximity to known FC genes; and many of the TFBSs learned by the classifier were found to be critical for enhancer activity, including POU homeodomain, Myb, Ets, Forkhead, and T-box motifs. Empirical testing also revealed that the T-box TF encoded by org-1 is a previously uncharacterized regulator of muscle cell identity. Finally, we found extensive diversity in the composition of TFBSs within known FC enhancers, suggesting that motif combinatorics plays an essential role in the cellular specificity exhibited by such enhancers. In summary, machine learning combined with evolutionary sequence analysis is useful for recognizing novel TFBSs and for facilitating the identification of cognate TFs that coordinate cell type–specific developmental gene expression patterns.

Vyšlo v časopise: A Machine Learning Approach for Identifying Novel Cell Type–Specific Transcriptional Regulators of Myogenesis. PLoS Genet 8(3): e32767. doi:10.1371/journal.pgen.1002531
Kategorie: Research Article
prolekare.web.journal.doi_sk: https://doi.org/10.1371/journal.pgen.1002531

Souhrn

Zdroje

1. CarrollSBGrenierJKWeatherbeeSD 2005 From DNA to Diversity. Molecular Genetics and the Evolution of Animal Design Malden, Massachusetts Blackwell Publishing

2. DavidsonE 2006 The Regulatory Genome: Gene Regulatory Networks In Development And Evolution: Academic Press 304

3. BermanBPNibuYPfeifferBDTomancakPCelnikerSE 2002 Exploiting transcription factor binding site clustering to identify cis-regulatory modules involved in pattern formation in the Drosophila genome. Proc Natl Acad Sci U S A 99 757 762

4. MarksteinMMarksteinPMarksteinVLevineM 2002 Genome-wide analysis of clustered Dorsal binding sites identifies putative target genes in the Drosophila embryo. Proc Natl Acad Sci U S A 99 763 768

5. PhilippakisAABusserBWGisselbrechtSSHeFSEstradaB 2006 Expression-guided in silico evaluation of candidate cis regulatory codes for Drosophila muscle founder cells. PLoS Comput Biol 2 e53 doi:10.1371/journal.pcbi.0020053

6. MahaffeyJW 2005 Assisting Hox proteins in controlling body form: are there new lessons from flies (and mammals)? Curr Opin Genet Dev 15 422 429

7. NarlikarLSakabeNJBlanskiAAArimuraFEWestlundJM 2010 Genome-wide discovery of human heart enhancers. Genome Res 20 381 392

8. MacArthurSLiXYLiJBrownJBChuHC 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

9. BigginMDMcGinnisW 1997 Regulation of segmentation and segmental identity by Drosophila homeoproteins: the role of DNA binding in functional activity and specificity. Development 124 4425 4433

10. IvanAHalfonMSSinhaS 2008 Computational discovery of cis-regulatory modules in Drosophila without prior knowledge of motifs. Genome Biol 9 R22

11. LeungGEisenMB 2009 Identifying cis-regulatory sequences by word profile similarity. PLoS ONE 4 e6901 doi:10.1371/journal.pone.0006901

12. NazinaAGPapatsenkoDA 2003 Statistical extraction of Drosophila cis-regulatory modules using exhaustive assessment of local word frequency. BMC Bioinformatics 4 65

13. KantorovitzMRKazemianMKinstonSMiranda-SaavedraDZhuQ 2009 Motif-blind, genome-wide discovery of cis-regulatory modules in Drosophila and mouse. Dev Cell 17 568 579

14. AbnizovaIte BoekhorstRWalterKGilksWR 2005 Some statistical properties of regulatory DNA sequences, and their use in predicting regulatory regions in the Drosophila genome: the fluffy-tail test. BMC Bioinformatics 6 109

15. NazinaAGPapatsenkoDA 2003 Statistical extraction of Drosophila cis-regulatory modules using exhaustive assessment of local word frequency. BMC Bioinformatics 4 65

16. LiLZhuQHeXSinhaSHalfonMS 2007 Large-scale analysis of transcriptional cis-regulatory modules reveals both common features and distinct subclasses. Genome Biol 8 R101

17. BusserBWBulykMLMichelsonAM 2008 Toward a systems-level understanding of developmental regulatory networks. Curr Opin Genet Dev 18 521 529

18. Ruiz-GomezM 1998 Muscle patterning and specification in Drosophila. Int J Dev Biol 42 283 290

19. BayliesMKBateMRuiz GomezM 1998 Myogenesis: a view from Drosophila. Cell 93 921 927

20. TixierVBatailleLJaglaK 2010 Diversification of muscle types: recent insights from Drosophila. Exp Cell Res 316 3019 3027

21. BerkesCATapscottSJ 2005 MyoD and the transcriptional control of myogenesis. Semin Cell Dev Biol 16 585 595

22. EstradaBChoeSEGisselbrechtSSMichaudSRajL 2006 An Integrated Strategy for Analyzing the Unique Developmental Programs of Different Myoblast Subtypes. PLoS Genet 2 e16 doi:10.1371/journal.pgen.0020016

23. CrockerJPotterNErivesA 2010 Dynamic evolution of precise regulatory encodings creates the clustered site signature of enhancers. Nature communications 1 99

24. CrockerJTamoriYErivesA 2008 Evolution acts on enhancer organization to fine-tune gradient threshold readouts. PLoS Biol 6 e263 doi:10.1371/journal.pbio.0060263

25. ErivesALevineM 2004 Coordinate enhancers share common organizational features in the Drosophila genome. Proc Natl Acad Sci USA 101 3851 3856

26. LudwigMBergmanCPatelNKreitmanM 2000 Evidence for stabilizing selection in a eukaryotic enhancer element. Nature 403 564 567

27. SwansonCIEvansNCBaroloS 2010 Structural rules and complex regulatory circuitry constrain expression of a Notch- and EGFR-regulated eye enhancer. Dev Cell 18 359 370

28. HalfonMSCarmenaAGisselbrechtSSackersonCMJiménezF 2000 Ras pathway specificity is determined by the integration of multiple signal-activated and tissue-restricted transcription factors. Cell 103 63 74

29. HalfonMSGradYChurchGMMichelsonAM 2002 Computation-based discovery of related transcriptional regulatory modules and motifs using an experimentally validated combinatorial model. Genome Res 12 1019 1028

30. HanZFujiokaMSuMLiuMJaynesJB 2002 Transcriptional integration of competence modulated by mutual repression generates cell-type specificity within the cardiogenic mesoderm. Dev Biol 252 225 240

31. KnirrSFraschM 2001 Molecular integration of inductive and mesoderm-intrinsic inputs governs even-skipped enhancer activity in a subset of pericardial and dorsal muscle progenitors. Developmental Biology 238 13 26

32. GalloSMGerrardDTMinerDSimichMDes SoyeB 2011 REDfly v3.0: toward a comprehensive database of transcriptional regulatory elements in Drosophila. Nucleic Acids Res 39 D118 123

33. ViselAMinovitskySDubchakIPennacchioLA 2007 VISTA Enhancer Browser–a database of tissue-specific human enhancers. Nucleic acids research 35 D88 92

34. WoolfeAGoodsonMGoodeDKSnellPMcEwenGK 2005 Highly conserved non-coding sequences are associated with vertebrate development. PLoS Biol 3 e7 doi:10.1371/journal.pbio.0030007

35. CapovillaMKambrisZBotasJ 2001 Direct regulation of the muscle-identity gene apterous by a Hox protein in the somatic mesoderm. Development 128 1221 1230

36. DuboisLEnriquezJDaburonVCrozetFLebretonG 2007 Collier transcription in a single Drosophila muscle lineage: the combinatorial control of muscle identity. Development 134 4347 4355

37. HollowayAKBegunDJSiepelAPollardKS 2008 Accelerated sequence divergence of conserved genomic elements in Drosophila melanogaster. Genome Res 18 1592 1601

38. LootsGGLocksleyRMBlankespoorCMWangZEMillerW 2000 Identification of a coordinate regulator of interleukins 4, 13, and 5 by cross-species sequence comparisons. Science 288 136 140

39. OvcharenkoILootsGGHardisonRCMillerWStubbsL 2004 zPicture: dynamic alignment and visualization tool for analyzing conservation profiles. Genome research 14 472 477

40. HareEEPetersonBKIyerVNMeierREisenMB 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

41. ArnostiDNKulkarniMM 2005 Transcriptional enhancers: Intelligent enhanceosomes or flexible billboards? J Cell Biochem 94 890 898

42. Bryson-RichardsonRJCurriePD 2008 The genetics of vertebrate myogenesis. Nat Rev Genet 9 632 646

43. KimSY 2009 Effects of sample size on robustness and prediction accuracy of a prognostic gene signature. BMC bioinformatics 10 147

44. HanczarBHuaJSimaCWeinsteinJBittnerM 2010 Small-sample precision of ROC-related estimates. Bioinformatics 26 822 830

45. RaudysSJJainAK 1991 Small sample size effects in statistical pattern recognition: recommendations for practitioners. IEEE Transactions on Pattern Analysis and Machine Intelligence 13 252 264

46. GrothACFishMNusseRCalosMP 2004 Construction of transgenic Drosophila by using the site-specific integrase from phage phiC31. Genetics 166 1775 1782

47. MarksteinMPitsouliCVillaltaCCelnikerSEPerrimonN 2008 Exploiting position effects and the gypsy retrovirus insulator to engineer precisely expressed transgenes. Nat Genet 40 476 483

48. FraschMHoeyTRushlowCDoyleHLevineM 1987 Characterization and localization of the even-skipped protein of Drosophila. EMBO J 6 749 759

49. CarmenaAGisselbrechtSHarrisonJJiménezFMichelsonAM 1998 Combinatorial signaling codes for the progressive determination of cell fates in the Drosophila embryonic mesoderm. Genes Dev 12 3910 3922

50. SengerKArmstrongGWRowellWJKwanJMMarksteinM 2004 Immunity regulatory DNAs share common organizational features in Drosophila. Mol Cell 13 19 32

51. NarlikarLGordanROhlerUHarteminkAJ 2006 Informative priors based on transcription factor structural class improve de novo motif discovery. Bioinformatics 22 e384 392

52. BeerMATavazoieS 2004 Predicting gene expression from sequence. Cell 117 185 198

53. BlowMJMcCulleyDJLiZZhangTAkiyamaJA 2010 ChIP-Seq identification of weakly conserved heart enhancers. Nature genetics 42 806 810

54. KnirrSAzpiazuNFraschM 1999 The role of the NK-homeobox gene slouch (S59) in somatic muscle patterning. Development 126 4525 4535

55. ArnoneMIDavidsonEH 1997 The hardwiring of development: organization and function of genomic regulatory systems. Development 124 1851 1864

56. GrossniklausUPearsonRKGehringWJ 1992 The Drosophila sloppy paired locus encodes two proteins involved in segmentation that show homology to mammalian transcription factors. Genes Dev 6 1030 1051

57. LiuYHJakobsenJSValentinGAmarantosIGilmourDT 2009 A systematic analysis of Tinman function reveals Eya and JAK-STAT signaling as essential regulators of muscle development. Dev Cell 16 280 291

58. JaglaKBellardMFraschM 2001 A cluster of Drosophila homeobox genes involved in mesoderm differentiation programs. Bioessays 23 125 133

59. MichelsonAM 1994 Muscle pattern diversification in Drosophila is determined by the autonomous function of homeotic genes in the embryonic mesoderm. Development 120 755 768

60. MahonySBenosPV 2007 STAMP: a web tool for exploring DNA-binding motif similarities. Nucleic Acids Res 35 W253 258

61. BodmerR 1993 The gene tinman is required for specification of the heart and visceral muscles in Drosophila. Development 118 719 729

62. ReimIFraschM 2005 The Dorsocross T-box genes are key components of the regulatory network controlling early cardiogenesis in Drosophila. Development 132 4911 4925

63. ReimIMohlerJPFraschM 2005 Tbx20-related genes, mid and H15, are required for tinman expression, proper patterning, and normal differentiation of cardioblasts in Drosophila. Mech Dev 122 1056 1069

64. OlsonEN 2006 Gene regulatory networks in the evolution and development of the heart. Science 313 1922 1927

65. PorschMSauerMSchulzeSBahloARothM 2005 The relative role of the T-domain and flanking sequences for developmental control and transcriptional regulation in protein chimeras of Drosophila OMB and ORG-1. Mech Dev 122 81 96

66. BlackBLOlsonEN 1998 Transcriptional control of muscle development by myocyte enhancer factor-2 MEF-2 proteins. Annu Rev Cell Dev Biol 14 167 196

67. EnriquezJBoukhatmiHDuboisLPhilippakisAABulykML 2010 Multi-step control of muscle diversity by Hox proteins in the Drosophila embryo. Development 137 457 466

68. BourgouinCLundgrenSEThomasJB 1992 apterous is a Drosophila LIM domain gene required for the development of a subset of embryonic muscles. Neuron 9 549 561

69. NoseAIsshikiTTakeichiM 1998 Regional specification of muscle progenitors in Drosophila: the role of the msh homeobox gene. Development 125 215 223

70. CarmenaABuffEHalfonMSGisselbrechtSJimenezF 2002 Reciprocal regulatory interactions between the Notch and Ras signaling pathways in the Drosophila embryonic mesoderm. Dev Biol 244 226 242

71. ViselABlowMJLiZZhangTAkiyamaJA 2009 ChIP-seq accurately predicts tissue-specific activity of enhancers. Nature 457 854 858

72. ErnstJKheradpourPMikkelsenTSShoreshNWardLD 2011 Mapping and analysis of chromatin state dynamics in nine human cell types. Nature 473 43 49

73. NegreNBrownCDMaLBristowCAMillerSW 2011 A cis-regulatory map of the Drosophila genome. Nature 471 527 531

74. HonGRenBWangW 2008 ChromaSig: a probabilistic approach to finding common chromatin signatures in the human genome. PLoS Comput Biol 4 e1000201 doi:10.1371/journal.pcbi.1000201

75. NarlikarLSakabeNJBlanskiAAArimuraFEWestlundJM 2010 Genome-wide discovery of human heart enhancers. Genome research 20 381 392

76. RouaultHMazouniKCouturierLHakimVSchweisguthF 2010 Genome-wide identification of cis-regulatory motifs and modules underlying gene coregulation using statistics and phylogeny. Proc Natl Acad Sci U S A 107 14615 14620

77. GriffithOLMontgomerySBBernierBChuBKasaianK 2008 ORegAnno: an open-access community-driven resource for regulatory annotation. Nucleic acids research 36 D107 113

78. JiangCXuanZZhaoFZhangMQ 2007 TRED: a transcriptional regulatory element database, new entries and other development. Nucleic acids research 35 D137 140

79. CrawfordGEHoltIEWhittleJWebbBDTaiD 2006 Genome-wide mapping of DNase hypersensitive sites using massively parallel signature sequencing (MPSS). Genome research 16 123 131

80. PedersonJALaFolletteJWGrossCVeraksaAMcGinnisW 2000 Regulation by homeoproteins: a comparison of deformed-responsive elements. Genetics 156 677 686

81. TaylorMV 2006 Comparison of muscle development in Drosophila and vertebrates. HS Muscle development in Drosophila: Landes BIoscience

82. ZinzenRPGirardotCGagneurJBraunMFurlongEE 2009 Combinatorial binding predicts spatio-temporal cis-regulatory activity. Nature 462 65 70

83. NoyesMBChristensenRGWakabayashiAStormoGDBrodskyMH 2008 Analysis of homeodomain specificities allows the family-wide prediction of preferred recognition sites. Cell 133 1277 1289

84. ReevesNPosakonyJW 2005 Genetic programs activated by proneural proteins in the developing Drosophila PNS. Dev Cell 8 413 425

85. CarmenaABateMJiménezF 1995 lethal of scute, a proneural gene, participates in the specification of muscle progenitors during Drosophila embryogenesis. Genes Dev 9 2373 2383

86. JakobsenJSBraunMAstorgaJGustafsonEHSandmannT 2007 Temporal ChIP-on-chip reveals Biniou as a universal regulator of the visceral muscle transcriptional network. Genes Dev 21 2448 2460

87. JunionGJaglaTDuplantSTapinRDa PonteJP 2005 Mapping Dmef2-binding regulatory modules by using a ChIP-enriched in silico targets approach. Proc Natl Acad Sci U S A 102 18479 18484

88. OdomDTZizlspergerNGordonDBBellGWRinaldiNJ 2004 Control of pancreas and liver gene expression by HNF transcription factors. Science 303 1378 1381

89. SandmannTGirardotCBrehmeMTongprasitWStolcV 2007 A core transcriptional network for early mesoderm development in Drosophila melanogaster. Genes Dev 21 436 449

90. SandmannTJensenLJJakobsenJSKarzynskiMMEichenlaubMP 2006 A temporal map of transcription factor activity: mef2 directly regulates target genes at all stages of muscle development. Dev Cell 10 797 807

91. ZeitlingerJZinzenRPStarkAKellisMZhangH 2007 Whole-genome ChIP-chip analysis of Dorsal, Twist, and Snail suggests integration of diverse patterning processes in the Drosophila embryo. Genes Dev 21 385 390

92. OngCTCorcesVG 2011 Enhancer function: new insights into the regulation of tissue-specific gene expression. Nat Rev Genet 12 283 293

93. KlämbtC 1993 The Drosophila gene pointed encodes two ETS-like proteins which are involved in the development of midline glial cells. Development 117 163 176

94. BayliesMKBateM 1996 Twist: A myogenic switch in Drosophila. Science 272 1481 1484

95. MichelsonAMAbmayrSMBateMMartinez AriasAManiatisT 1990 Expression of a MyoD family member prefigures muscle pattern in Drosophila embryos. Genes Dev 4 2086 2097

96. BayliesMKAriasAMBateM 1995 wingless is required for the formation of a subset of muscle founder cells during Drosophila embryogenesis. Development 121 3829 3837

97. RanganayakuluGSchulzRAOlsonEN 1996 Wingless signaling induces nautilus expression in the ventral mesoderm of the Drosophila embryo. Dev Biol 176 143 148

98. TomancakPBermanBPBeatonAWeiszmannRKwanE 2007 Global analysis of patterns of gene expression during Drosophila embryogenesis. Genome Biol 8 R145

99. RobaskyKBulykML 2011 UniPROBE, update 2011: expanded content and search tools in the online database of protein-binding microarray data on protein-DNA interactions. Nucleic Acids Res 39 D124 128

100. RamsayRG 2005 c-Myb a stem-progenitor cell regulator in multiple tissue compartments. Growth Factors 23 253 261

101. PaululatABreuerSRenkawitz-PohlR 1999 Determination and development of the larval muscle pattern in Drosophila melanogaster. Cell Tissue Res 296 151 160

102. ClaycombJMOrr-WeaverTL 2005 Developmental gene amplification: insights into DNA replication and gene expression. Trends Genet 21 149 162

103. SchwaigerMSchubelerD 2006 A question of timing: emerging links between transcription and replication. Curr Opin Genet Dev 16 177 183

104. NessSA 2003 Myb protein specificity: evidence of a context-specific transcription factor code. Blood Cells Mol Dis 31 192 200

105. ShihHPGrossMKKioussiC 2008 Muscle development: forming the head and trunk muscles. Acta Histochem 110 97 108

106. WardleFCPapaioannouVE 2008 Teasing out T-box targets in early mesoderm. Curr Opin Genet Dev 18 418 425

107. CostelloIPimeislIMDragerSBikoffEKRobertsonEJ 2011 The T-box transcription factor Eomesodermin acts upstream of Mesp1 to specify cardiac mesoderm during mouse gastrulation. Nature cell biology 13 1084 1091

108. JaglaTBellardFLutzYDretzenGBellardM 1998 ladybird determines cell fate decisions during diversification of Drosophila somatic muscles. Development 125 3699 3708

109. DastjerdiARobsonLWalkerRHadleyJZhangZ 2007 Tbx1 regulation of myogenic differentiation in the limb and cranial mesoderm. Dev Dyn 236 353 363

110. BrohmannHJaglaKBirchmeierC 2000 The role of Lbx1 in migration of muscle precursor cells. Development 127 437 445

111. GrossMKMoran-RivardLVelasquezTNakatsuMNJaglaK 2000 Lbx1 is required for muscle precursor migration along a lateral pathway into the limb. Development 127 413 424

112. PedersonJDKiehartDPMahaffeyJW 1996 The role of HOM-C genes in segmental transformations: reexamination of the Drosophila Sex combs reduced embryonic phenotype. Dev Biol 180 131 142

113. WeatherbeeSDNijhoutHFGrunertLWHalderGGalantR 1999 Ultrabithorax function in butterfly wings and the evolution of insect wing patterns. Curr Biol 9 109 115

114. GalantRCarrollSB 2002 Evolution of a transcriptional repression domain in an insect Hox protein. Nature 415 910 913

115. BischofJMaedaRKHedigerMKarchFBaslerK 2007 An optimized transgenesis system for Drosophila using germ-line-specific phiC31 integrases. Proc Natl Acad Sci U S A 104 3312 3317

116. ClynePJCertelSJde BruyneMZaslavskyLJohnsonWA 1999 The odor specificities of a subset of olfactory receptor neurons are governed by Acj6, a POU-domain transcription factor. Neuron 22 339 347

117. BaroloSCarverLAPosakonyJW 2000 GFP and beta-galactosidase transformation vectors for promoter/enhancer analysis in Drosophila. Biotechniques 29 726 728730732

118. BaileyTLElkanC 1994 Fitting a mixture model by expectation maximization to discover motifs in biopolymers. Proc Int Conf Intell Syst Mol Biol 2 28 36

119. MatysVKel-MargoulisOVFrickeELiebichILandS 2006 TRANSFAC and its module TRANSCompel: transcriptional gene regulation in eukaryotes. Nucleic Acids Res 34 D108 110

120. BenjaminiYHochbergY 1995 Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Roy Statist Soc Ser B 57 289 300

121. ChangCLinC 2011 LIBSVM: a library for support vector machines. ACM Transactions on Intelligent Systems and Technology 2 27:1 27:27

122. Shawe-TaylorJCristianiniN 2002 On the Generalization of Soft Margin Algorithms. IEEE Transactions on Information Theory 48 2721 2735

123. GuyonIWestonJBarnhillSVapnikV 2002 Gene Selection for Cancer Classification using Support Vector Machines. Machine Learning 46 389 422

124. SubramanianATamayoPMoothaVKMukherjeeSEbertBL 2005 Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A 102 15545 15550

125. SiepelABejeranoGPedersenJSHinrichsASHouM 2005 Evolutionarily conserved elements in vertebrate, insect, worm, and yeast genomes. Genome research 15 1034 1050

126. BergerMFBadisGGehrkeARTalukderSPhilippakisAA 2008 Variation in homeodomain DNA binding revealed by high-resolution analysis of sequence preferences. Cell 133 1266 1276