Comparative motif discovery combined with comparative transcriptomics yields accurate targetome and enhancer predictions

Synergistic activation by Glass and Pointed promotes neuronal identity in the Drosophila eye disc

The integration of extrinsic signaling with cell-intrinsic transcription factors can direct progenitor cells to differentiate into distinct cell fates. In the developing Drosophila eye, differentiation of photoreceptors R1-R7 requires EGFR signaling mediated by the transcription factor Pointed, and our single-cell RNA-Seq analysis shows that the same photoreceptors require the eye-specific transcription factor Glass. We find that ectopic expression of Glass and activation of EGFR signaling synergistically induce neuronal gene expression in the wing disc in a Pointed-dependent manner. Targeted DamID reveals that Glass and Pointed share many binding sites in the genome of developing photoreceptors. Comparison with transcriptomic data shows that Pointed and Glass induce photoreceptor differentiation through intermediate transcription factors, including the redundant homologs Scratch and Scrape, as well as directly activating neuronal effector genes. Our data reveal synergistic activation of a multi-layered transcriptional network as the mechanism by which EGFR signaling induces neuronal identity in Glass-expressing cells.

Intra-Exon Motif Correlations as a Proxy Measure for Mean Per-Tile Sequence Quality Data in RNA-Seq

Given the wide variability in the quality of next-generation sequencing data submitted to public repositories, it is essential to identify methods that can perform quality control on these data sets when additional quality control data, such as mean tile data, are missing from public repositories. In this study, we present evidence that correlating counts of reads corresponding to pairs of motifs separated over specific distances on individual exons can be used as a proxy mean tile data in the data sets we analyzed and hence could be used when mean tile data are not available. As test data sets we use the Homo sapiens in vitro transcribed (IVT) data set, and a Drosophila melanogaster data set comprising wild and mutant types. We find that a FastQC analysis of the available parts of these data sets demonstrates that the per-tile sequencing quality is good for all the data sets apart from the mutant-type data where the mutant-r3 data are worse than the mutant-r2 data. Correspondingly, intra-exon motif correlations are reasonably large for all data sets except this latter case where the mutant-r2 correlations are low and the mutant-r3 correlations close to zero. We propose that these extremely low correlations are indicative of bias of technical origin, such as flowcell errors. In addition to this, the intra-exon motif correlations as a function of both guanosine-cytosine (GC) content parameters are somewhat higher and less dependent on the GC content parameters in the IVT-Plasmids messenger RNA (mRNA) selection free RNA-Seq sample (control) than in the other RNA-Seq samples that did undergo mRNA selection: both ribosomal depletion (IVT-Only) and PolyA selection (IVT-PolyA, wild type, and mutant).

A TRiP RNAi screen to identify molecules necessary for Drosophila photoreceptor differentiation

Drosophila rhabdomeric terminal photoreceptor differentiation is an extended process taking several days to complete. Following ommatidial patterning by the morphogenetic furrow, photoreceptors are sequentially recruited and specified, and terminal differentiation begins. Key events of terminal differentiation include the establishment of apical and basolateral domains, rhabdomere and stalk formation, inter-rhabdomeral space formation, and expression of phototransduction machinery. While many key regulators of these processes have been identified, the complete network of transcription factors to downstream effector molecules necessary for regulating each of these major events remains incomplete. Here, we report an RNAi screen to identify additional molecules and cellular pathways required for photoreceptor terminal differentiation. First, we tested several eye-specific GAL4 drivers for correct spatial and temporal specificity and identified Pph13-GAL4 as the most appropriate GAL4 line for our screen. We screened lines available through the Transgenic RNAi Project and isolated lines that when combined with Pph13-GAL4 resulted in the loss of the deep pseudopupil, as a readout for abnormal differentiation. In the end, we screened 6,189 lines, representing 3,971 genes, and have identified 64 genes, illuminating potential new regulatory molecules and cellular pathways for the differentiation and organization of Drosophila rhabdomeric photoreceptors.

R7 photoreceptor axon targeting depends on the relative levels of lost and found expression in R7 and its synaptic partners

As neural circuits form, growing processes select the correct synaptic partners through interactions between cell surface proteins. The presence of such proteins on two neuronal processes may lead to either adhesion or repulsion; however, the consequences of mismatched expression have rarely been explored. Here, we show that the Drosophila CUB-LDL protein Lost and found (Loaf) is required in the UV-sensitive R7 photoreceptor for normal axon targeting only when Loaf is also present in its synaptic partners. Although targeting occurs normally in loaf mutant animals, removing loaf from photoreceptors or expressing it in their postsynaptic neurons Tm5a/b or Dm9 in a loaf mutant causes mistargeting of R7 axons. Loaf localizes primarily to intracellular vesicles including endosomes. We propose that Loaf regulates the trafficking or function of one or more cell surface proteins, and an excess of these proteins on the synaptic partners of R7 prevents the formation of stable connections.

New regulators of Drosophila eye development identified from temporal transcriptome changes

In the last larval instar, uncommitted progenitor cells in the Drosophila eye primordium start to adopt individual retinal cell fates, arrest their growth and proliferation, and initiate terminal differentiation into photoreceptor neurons and other retinal cell types. To explore the regulation of these processes, we have performed mRNA-Seq studies of the larval eye and antennal primordial at multiple developmental stages. A total of 10,893 fly genes were expressed during these stages and could be adaptively clustered into gene groups, some of whose expression increases or decreases in parallel with the cessation of proliferation and onset of differentiation. Using in situ hybridization of a sample of 98 genes to verify spatial and temporal expression patterns, we estimate that 534 genes or more are transcriptionally upregulated during retinal differentiation, and 1367 or more downregulated as progenitor cells differentiate. Each group of co-expressed genes is enriched for regulatory motifs recognized by co-expressed transcription factors, suggesting that they represent coherent transcriptional regulatory programs. Using available mutant strains, we describe novel roles for the transcription factors SoxNeuro (SoxN), H6-like homeobox (Hmx), CG10253, without children (woc), Structure specific recognition protein (Ssrp), and multisex combs (mxc).

Integrative genomic analysis reveals novel regulatory mechanisms of eyeless during Drosophila eye development

Eyeless (ey) is one of the most critical transcription factors for initiating the entire eye development in Drosophila. However, the molecular mechanisms through which Ey regulates target genes and pathways have not been characterized at the genomic level. Using ChIP-Seq, we generated an endogenous Ey-binding profile in Drosophila developing eyes. We found that Ey binding occurred more frequently at promoter compared to non-promoter regions. Ey promoter binding was correlated with the active transcription of genes involved in development and transcription regulation. An integrative analysis revealed that Ey directly regulated a broad and highly connected genetic network, including many essential patterning pathways, and known and novel eye genes. Interestingly, we observed that Ey could target multiple components of the same pathway, which might enhance its control of these pathways during eye development. In addition to protein-coding genes, we discovered Ey also targeted non-coding RNAs, which represents a new regulatory mechanism employed by Ey. These findings suggest that Ey could use multiple molecular mechanisms to regulate target gene expression and pathway function, which might enable Ey to exhibit a greater flexibility in controlling different processes during eye development.

Multilevel regulation of the glass locus during Drosophila eye development

Development of eye tissue is initiated by a conserved set of transcription factors termed retinal determination network (RDN). In the fruit fly Drosophila melanogaster, the zinc-finger transcription factor Glass acts directly downstream of the RDN to control identity of photoreceptor as well as non-photoreceptor cells. Tight control of spatial and temporal gene expression is a critical feature during development, cell-fate determination as well as maintenance of differentiated tissues. The molecular mechanisms that control expression of glass, however, remain largely unknown. We here identify complex regulatory mechanisms controlling expression of the glass locus. All information to recapitulate glass expression are contained in a compact 5.2 kb cis-acting genomic element by combining different cell-type specific and general enhancers with repressor elements. Moreover, the immature RNA of the locus contains an alternative small open reading frame (smORF) upstream of the actual glass translation start, resulting in a small peptide instead of the three possible Glass protein isoforms. CRISPR/Cas9-based mutagenesis shows that the smORF is not required for the formation of functioning photoreceptors, but is able to attenuate effects of glass misexpression. Furthermore, editing the genome to generate glass loci eliminating either one or two isoforms shows that only one of the three proteins is critical for formation of functioning photoreceptors, while removing the two other isoforms did not cause defects in developmental or photoreceptor function. Our results show that eye development and function is largely unaffected by targeted manipulations of critical features of the glass transcript, suggesting a strong selection pressure to allow the formation of a functioning eye.

winged eye Induces Transdetermination of Drosophila Imaginal Disc by Acting in Concert with a Histone Methyltransferase, Su(var)3-9

Drosophila imaginal disc cells exhibit a remarkable ability to convert cell fates in response to various perturbations, a phenomenon called transdetermination (TD). We previously identified winged eye (wge) as a factor that induces eye-to-wing TD upon overexpression in eye imaginal discs, but the molecular mechanisms underlying TD have remained largely unclear. Here, we found that wge induces various histone modifications and enhances the methylation of Lys9 on histone H3 (H3K9), a feature of heterochromatin. A histone methyltransferase, Su(var)3-9, is required for wge-mediated H3K9 methylation and eye-to-wing TD. Su(var)3-9 is also required for classical wound-induced TD but not for normal development, suggesting its involvement in several types of imaginal disc TDs. Transcriptome analysis revealed that wge represses eye identity genes independently of Su(var)3-9 and activates TD-related genes by acting together with Su(var)3-9. These findings provide new insights into diverse types of chromatin regulation at progressive steps of cell-fate conversions.

Systems Biology Reveals NR2F6 and TGFB1 as Key Regulators of Feed Efficiency in Beef Cattle

Systems biology approaches are used as strategy to uncover tissue-specific perturbations and regulatory genes related to complex phenotypes. We applied this approach to study feed efficiency (FE) in beef cattle, an important trait both economically and environmentally. Poly-A selected RNA of five tissues (adrenal gland, hypothalamus, liver, skeletal muscle and pituitary) of eighteen young bulls, selected for high and low FE, were sequenced (Illumina HiSeq 2500, 100 bp, pared-end). From the 17,354 expressed genes considering all tissues, 1,335 were prioritized by five selection categories (differentially expressed, harboring SNPs associated with FE, tissue-specific, secreted in plasma and key regulators) and used for network construction. NR2F6 and TGFB1 were identified and validated by motif discovery as key regulators of hepatic inflammatory response and muscle tissue development, respectively, two biological processes demonstrated to be associated with FE. Moreover, we indicated potential biomarkers of FE, which are related to hormonal control of metabolism and sexual maturity. By using robust methodologies and validation strategies, we confirmed the main biological processes related to FE in Bos indicus and indicated candidate genes as regulators or biomarkers of superior animals.

Glass confers rhabdomeric photoreceptor identity in Drosophila, but not across all metazoans

Across metazoans, visual systems employ different types of photoreceptor neurons (PRs) to detect light. These include rhabdomeric PRs, which exist in distantly related phyla and possess an evolutionarily conserved phototransduction cascade. While the development of rhabdomeric PRs has been thoroughly studied in the fruit fly Drosophila melanogaster, we still know very little about how they form in other species. To investigate this question, we tested whether the transcription factor Glass, which is crucial for instructing rhabdomeric PR formation in Drosophila, may play a similar role in other metazoans. Glass homologues exist throughout the animal kingdom, indicating that this protein evolved prior to the metazoan radiation. Interestingly, our work indicates that glass is not expressed in rhabdomeric photoreceptors in the planarian Schmidtea mediterranea nor in the annelid Platynereis dumerilii. Combined with a comparative analysis of the Glass DNA-binding domain, our data suggest that the fate of rhabdomeric PRs is controlled by Glass-dependent and Glass-independent mechanisms in different animal clades.

A switch in transcription and cell fate governs the onset of an epigenetically-deregulated tumor in Drosophila

Tumor initiation is often linked to a loss of cellular identity. Transcriptional programs determining cellular identity are preserved by epigenetically-acting chromatin factors. Although such regulators are among the most frequently mutated genes in cancer, it is not well understood how an abnormal epigenetic condition contributes to tumor onset. In this work, we investigated the gene signature of tumors caused by disruption of the Drosophila epigenetic regulator, polyhomeotic (ph). In larval tissue ph mutant cells show a shift towards an embryonic-like signature. Using loss- and gain-of-function experiments we uncovered the embryonic transcription factor knirps (kni) as a new oncogene. The oncogenic potential of kni lies in its ability to activate JAK/STAT signaling and block differentiation. Conversely, tumor growth in ph mutant cells can be substantially reduced by overexpressing a differentiation factor. This demonstrates that epigenetically derailed tumor conditions can be reversed when targeting key players in the transcriptional network.

When an animal is developing as an embryo, different cells start to specialize into the specific cell types needed to form the tissues and organs of the body. How an individual cell commits to become a certain type of cell is mostly determined by which of the genes in its DNA are active. In animal cells, DNA is wrapped around proteins called histones, and one way that cells can maintain their distinct pattern of gene activity is via chemical tags on the histones. These tags can switch nearby genes on or off, and are added or removed by other proteins called epigenetic regulators. The epigenetic tags are also stably inherited when the cell divides, meaning that a cell’s identity can be maintained over many cell generations.

If epigenetic regulators fail to work properly or get disrupted, the pattern of gene activity in a cell becomes altered. As a consequence, that cell can lose its identity and will often turn into a cancer cell. In fact, mutations in epigenetic regulators are found in several human cancers. It is not yet understood how these changes in gene expression lead cells to become cancerous.

Torres et al. have now analyzed an epigenetic regulator called Polyhomeotic in developing larvae of the fruit fly, Drosophila melanogaster. The results show that when Polyhomeotic is not produced the fly larvae develop tumors. Moreover, the mutant cells without Polyhomeotic had different gene expression profiles compared to normal cells. This in turn caused the mutant cells, which had previously committed to a certain fate, to become more like the unspecialized cells found in early embryos.

Torres et al. next showed that, among the genes that were incorrectly regulated when Polyhomeotic’s activity was compromised, one gene called knirps was switched on by mistake, which led the mutant cells to become tumor cells. When the activity of knirps was reduced instead, almost no tumors grew. Additionally, Torres et al. found that the protein encoded by knirps activates a signaling pathway that keeps tumor cells unspecialized by blocking their normal progress to a more mature and specialized state – a process known as differentiation. Experimentally raising the levels of a different molecule that ultimately promotes differentiation caused the tumor cells to grow less.

These findings suggest that tumors caused when epigenetic regulation goes awry may be reversed by targeting key genes such as knirps. Further work is now needed to test whether these findings will also extend to humans. Forcing cancer cells from a highly dividing, non-specialized state into a dead-end, mature state may lead to new ways to treat cancer.

Dynamic genome wide expression profiling of Drosophila head development reveals a novel role of Hunchback in retinal glia cell development and blood-brain barrier integrity

Drosophila melanogaster head development represents a valuable process to study the developmental control of various organs, such as the antennae, the dorsal ocelli and the compound eyes from a common precursor, the eye-antennal imaginal disc. While the gene regulatory network underlying compound eye development has been extensively studied, the key transcription factors regulating the formation of other head structures from the same imaginal disc are largely unknown. We obtained the developmental transcriptome of the eye-antennal discs covering late patterning processes at the late 2^nd larval instar stage to the onset and progression of differentiation at the end of larval development. We revealed the expression profiles of all genes expressed during eye-antennal disc development and we determined temporally co-expressed genes by hierarchical clustering. Since co-expressed genes may be regulated by common transcriptional regulators, we combined our transcriptome dataset with publicly available ChIP-seq data to identify central transcription factors that co-regulate genes during head development. Besides the identification of already known and well-described transcription factors, we show that the transcription factor Hunchback (Hb) regulates a significant number of genes that are expressed during late differentiation stages. We confirm that hb is expressed in two polyploid subperineurial glia cells (carpet cells) and a thorough functional analysis shows that loss of Hb function results in a loss of carpet cells in the eye-antennal disc. Additionally, we provide for the first time functional data indicating that carpet cells are an integral part of the blood-brain barrier. Eventually, we combined our expression data with a de novo Hb motif search to reveal stage specific putative target genes of which we find a significant number indeed expressed in carpet cells.

The development of different cell types must be tightly coordinated, and the eye-antennal imaginal discs of Drosophila melanogaster represent an excellent model to study the molecular mechanisms underlying this coordination. These imaginal discs contain the anlagen of nearly all adult head structures, such as the antennae, the head cuticle, the ocelli and the compound eyes. While large scale screens have been performed to unravel the gene regulatory network underlying compound eye development, a comprehensive understanding of genome wide expression dynamics throughout head development is still missing to date. We studied the genome wide gene expression dynamics during eye-antennal disc development in D. melanogaster to identify new central regulators of the underlying gene regulatory network. Expression based gene clustering and transcription factor motif enrichment analyses revealed a central regulatory role of the transcription factor Hunchback (Hb). We confirmed that hb is expressed in two polyploid retinal subperineurial glia cells (carpet cells). Our functional analysis shows that Hb is necessary for carpet cell development and we show for the first time that the carpet cells are an integral part of the blood-brain barrier.

Glass promotes the differentiation of neuronal and non-neuronal cell types in the Drosophila eye

Transcriptional regulators can specify different cell types from a pool of equivalent progenitors by activating distinct developmental programs. The Glass transcription factor is expressed in all progenitors in the developing Drosophila eye, and is maintained in both neuronal and non-neuronal cell types. Glass is required for neuronal progenitors to differentiate as photoreceptors, but its role in non-neuronal cone and pigment cells is unknown. To determine whether Glass activity is limited to neuronal lineages, we compared the effects of misexpressing it in neuroblasts of the larval brain and in epithelial cells of the wing disc. Glass activated overlapping but distinct sets of genes in these neuronal and non-neuronal contexts, including markers of photoreceptors, cone cells and pigment cells. Coexpression of other transcription factors such as Pax2, Eyes absent, Lozenge and Escargot enabled Glass to induce additional genes characteristic of the non-neuronal cell types. Cell type-specific glass mutations generated in cone or pigment cells using somatic CRISPR revealed autonomous developmental defects, and expressing Glass specifically in these cells partially rescued glass mutant phenotypes. These results indicate that Glass is a determinant of organ identity that acts in both neuronal and non-neuronal cells to promote their differentiation into functional components of the eye.

A Novel Method to Detect Bias in Short Read NGS Data

Detecting sources of bias in transcriptomic data is essential to determine signals of Biological significance. We outline a novel method to detect sequence specific bias in short read Next Generation Sequencing data. This is based on determining intra-exon correlations between specific motifs. This requires a mild assumption that short reads sampled from specific regions from the same exon will be correlated with each other. This has been implemented on Apache Spark and used to analyse two D. melanogaster eye-antennal disc data sets generated at the same laboratory. The wild type data set in drosophila indicates a variation due to motif GC content that is more significant than that found due to exon GC content. The software is available online and could be applied for cross-experiment transcriptome data analysis in eukaryotes.

Successive requirement of Glass and Hazy for photoreceptor specification and maintenance in Drosophila

Development of the insect compound eye requires a highly controlled interplay between transcription factors. However, the genetic mechanisms that link early eye field specification to photoreceptor terminal differentiation and fate maintenance remain largely unknown. Here, we decipher the function of 2 transcription factors, Glass and Hazy, which play a central role during photoreceptor development. The regulatory interactions between Glass and Hazy suggest that they function together in a coherent feed-forward loop in all types of Drosophila photoreceptors. While the glass mutant eye lacks the expression of virtually all photoreceptor genes, young hazy mutants correctly express most phototransduction genes. Interestingly, the expression of these genes is drastically reduced in old hazy mutants. This age-dependent loss of the phototransduction cascade correlates with a loss of phototaxis in old hazy mutant flies. We conclude that Glass can either directly or indirectly initiate the expression of most phototransduction proteins in a Hazy-independent manner, and that Hazy is mainly required for the maintenance of functional photoreceptors in adult flies.

Two temporal functions of Glass: Ommatidium patterning and photoreceptor differentiation

Much progress has been made in elucidating the molecular networks required for specifying retinal cells, including photoreceptors, but the downstream mechanisms that maintain identity and regulate differentiation remain poorly understood. Here, we report that the transcription factor Glass has a dual role in establishing a functional Drosophila eye. Utilizing conditional rescue approaches, we confirm that persistent defects in ommatidium patterning combined with cell death correlate with the overall disruption of eye morphology in glass mutants. In addition, we reveal that Glass exhibits a separable role in regulating photoreceptor differentiation. In particular, we demonstrate the apparent loss of glass mutant photoreceptors is not only due to cell death but also a failure of the surviving photoreceptors to complete differentiation. Moreover, the late reintroduction of Glass in these developmentally stalled photoreceptors is capable of restoring differentiation in the absence of correct ommatidium patterning. Mechanistically, transcription profiling at the time of differentiation reveals that Glass is necessary for the expression of many genes implicated in differentiation, i.e. rhabdomere morphogenesis, phototransduction, and synaptogenesis. Specifically, we show Glass directly regulates the expression of Pph13, which encodes a transcription factor necessary for opsin expression and rhabdomere morphogenesis. Finally, we demonstrate the ability of Glass to choreograph photoreceptor differentiation is conserved between Drosophila and Tribolium, two holometabolous insects. Altogether, our work identifies a fundamental regulatory mechanism to generate the full complement of cells required for a functional rhabdomeric visual system and provides a critical framework to investigate the basis of differentiation and maintenance of photoreceptor identity.

The transcription factor Glass links eye field specification with photoreceptor differentiation in Drosophila

Eye development requires an evolutionarily conserved group of transcription factors, termed the retinal determination network (RDN). However, little is known about the molecular mechanism by which the RDN instructs cells to differentiate into photoreceptors. We show that photoreceptor cell identity in Drosophila is critically regulated by the transcription factor Glass, which is primarily expressed in photoreceptors and whose role in this process was previously unknown. Glass is both required and sufficient for the expression of phototransduction proteins. Our results demonstrate that the RDN member Sine oculis directly activates glass expression, and that Glass activates the expression of the transcription factors Hazy and Otd. We identified hazy as a direct target of Glass. Induced expression of Hazy in the retina partially rescues the glass mutant phenotype. Together, our results provide a transcriptional link between eye field specification and photoreceptor differentiation in Drosophila, placing Glass at a central position in this developmental process.

Ancient default activators of terminal photoreceptor differentiation in the pancrustacean compound eye: the homeodomain transcription factors Otd and Pph13

The origin of the Drosophila compound eye predates the ancestor of Pancrustacea, the arthropod clade that includes insects and Crustaceans. Recent studies in emerging model systems for pancrustacean development-the red flour beetle Tribolium castaneum and water flea Daphnia pulex-have begun to shed light on the evolutionary conservation of transcriptional mechanisms found for the Drosophila compound eye. Here, we discuss the conserved roles of the transcription factors Otd and Pph13, which complement each other in two terminal events of photoreceptor differentiation: rhabdomere morphogenesis and transcriptional default activation of opsin gene expression. The synthesis of these data allows us to frame an evolutionary developmental model of the earliest events that generated the wavelength-specific photoreceptor subtypes of pancrustacean compound eyes.

cis-regulatory analysis of the Drosophila pdm locus reveals a diversity of neural enhancers

Background

One of the major challenges in developmental biology is to understand the regulatory events that generate neuronal diversity. During Drosophila embryonic neural lineage development, cellular temporal identity is established in part by a transcription factor (TF) regulatory network that mediates a cascade of cellular identity decisions. Two of the regulators essential to this network are the POU-domain TFs Nubbin and Pdm-2, encoded by adjacent genes collectively known as pdm. The focus of this study is the discovery and characterization of cis-regulatory DNA that governs their expression.

Results

Phylogenetic footprinting analysis of a 125 kb genomic region that spans the pdm locus identified 116 conserved sequence clusters. To determine which of these regions function as cis-regulatory enhancers that regulate the dynamics of pdm gene expression, we tested each for in vivo enhancer activity during embryonic development and postembryonic neurogenesis. Our screen revealed 77 unique enhancers positioned throughout the noncoding region of the pdm locus. Many of these activated neural-specific gene expression during different developmental stages and many drove expression in overlapping patterns. Sequence comparisons of functionally related enhancers that activate overlapping expression patterns revealed that they share conserved elements that can be predictive of enhancer behavior. To facilitate data accessibility, the results of our analysis are catalogued in cisPatterns, an online database of the structure and function of these and other Drosophila enhancers.

Conclusions

These studies reveal a diversity of modular enhancers that most likely regulate pdm gene expression during embryonic and adult development, highlighting a high level of temporal and spatial expression specificity. In addition, we discovered clusters of functionally related enhancers throughout the pdm locus. A subset of these enhancers share conserved elements including sequences that correspond to known TF DNA binding sites. Although comparative analysis of the nubbin and pdm-2 encoding sequences indicate that these two genes most likely arose from a duplication event, we found only partial evidence of sequence duplication between their enhancers, suggesting that after the putative duplication their cis-regulatory DNA diverged at a higher rate than their coding sequences.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1897-2) contains supplementary material, which is available to authorized users.

Identification of Lineage-Specific Cis-Regulatory Modules Associated with Variation in Transcription Factor Binding and Chromatin Activity Using Ornstein-Uhlenbeck Models

Scoring the impact of noncoding variation on the function of cis-regulatory regions, on their chromatin state, and on the qualitative and quantitative expression levels of target genes is a fundamental problem in evolutionary genomics. A particular challenge is how to model the divergence of quantitative traits and to identify relationships between the changes across the different levels of the genome, the chromatin activity landscape, and the transcriptome. Here, we examine the use of the Ornstein-Uhlenbeck (OU) model to infer selection at the level of predicted cis-regulatory modules (CRMs), and link these with changes in transcription factor binding and chromatin activity. Using publicly available cross-species ChIP-Seq and STARR-Seq data we show how OU can be applied genome-wide to identify candidate transcription factors for which binding site and CRM turnover is correlated with changes in regulatory activity. Next, we profile open chromatin in the developing eye across three Drosophila species. We identify the recognition motifs of the chromatin remodelers, Trithorax-like and Grainyhead as mostly correlating with species-specific changes in open chromatin. In conclusion, we show in this study that CRM scores can be used as quantitative traits and that motif discovery approaches can be extended towards more complex models of divergence.

iRegulon: from a gene list to a gene regulatory network using large motif and track collections

Identifying master regulators of biological processes and mapping their downstream gene networks are key challenges in systems biology. We developed a computational method, called iRegulon, to reverse-engineer the transcriptional regulatory network underlying a co-expressed gene set using cis-regulatory sequence analysis. iRegulon implements a genome-wide ranking-and-recovery approach to detect enriched transcription factor motifs and their optimal sets of direct targets. We increase the accuracy of network inference by using very large motif collections of up to ten thousand position weight matrices collected from various species, and linking these to candidate human TFs via a motif2TF procedure. We validate iRegulon on gene sets derived from ENCODE ChIP-seq data with increasing levels of noise, and we compare iRegulon with existing motif discovery methods. Next, we use iRegulon on more challenging types of gene lists, including microRNA target sets, protein-protein interaction networks, and genetic perturbation data. In particular, we over-activate p53 in breast cancer cells, followed by RNA-seq and ChIP-seq, and could identify an extensive up-regulated network controlled directly by p53. Similarly we map a repressive network with no indication of direct p53 regulation but rather an indirect effect via E2F and NFY. Finally, we generalize our computational framework to include regulatory tracks such as ChIP-seq data and show how motif and track discovery can be combined to map functional regulatory interactions among co-expressed genes. iRegulon is available as a Cytoscape plugin from http://iregulon.aertslab.org.

Gene regulatory networks control developmental, homeostatic, and disease processes by governing precise levels and spatio-temporal patterns of gene expression. Determining their topology can provide mechanistic insight into these processes. Gene regulatory networks consist of interactions between transcription factors and their direct target genes. Each regulatory interaction represents the binding of the transcription factor to a specific DNA binding site near its target gene. Here we present a computational method, called iRegulon, to identify master regulators and direct target genes in a human gene signature, i.e. a set of co-expressed genes. iRegulon relies on the analysis of the regulatory sequences around each gene in the gene set to detect enriched TF motifs or ChIP-seq peaks, using databases of nearly 10.000 TF motifs and 1000 ChIP-seq data sets or “tracks”. Next, it associates enriched motifs and tracks with candidate transcription factors and determines the optimal subset of direct target genes. We validate iRegulon on ENCODE data, and use it in combination with RNA-seq and ChIP-seq data to map a p53 downstream network with new predicted co-factors and targets. iRegulon is available as a Cytoscape plugin, supporting human, mouse, and Drosophila genes, and provides access to hundreds of cancer-related TF-target subnetworks or “regulons”.

Ndae1 expression and regulation in Drosophila embryos

The construction and prediction of cell fate maps at the whole embryo level require the establishment of an accurate atlas of gene expression patterns throughout development and the identification of the corresponding cis-regulatory sequences. However, while the expression and regulation of genes encoding upstream developmental regulators such as transcription factors or signaling pathway components have been analyzed in detail, up to date the number of cis-regulatory sequences identified for downstream effector genes, like ion channels, pumps and exchangers, is very low. The control and regulation of ion homeostasis in each cell, including at blastoderm stages, are essential for normal embryonic development. In this study, we analyzed in detail the embryonic expression pattern and cis-regulatory modules of the Drosophila Na+-driven anion exchanger 1 (Ndae1) gene, involved in the regulation of pH homeostasis. We show that Ndae1 is expressed in a tight and complex spatial-temporal pattern. In particular, we report that this downstream effector gene is under the control of the canonical dorsal-ventral patterning cascade through dorsal, Toll, twist and snail at early embryogenesis. Moreover, we identify several cis-regulatory modules, some of which control discrete and non-overlapping aspects of endogenous gene expression throughout development.