Developmental mechanisms and cis-regulatory codes

An Intronic cis-Regulatory Element Is Crucial for the Alpha Tubulin Pl-Tuba1a Gene Activation in the Ciliary Band and Animal Pole Neurogenic Domains during Sea Urchin Development

In sea urchin development, structures derived from neurogenic territory control the swimming and feeding responses of the pluteus as well as the process of metamorphosis. We have previously isolated an alpha tubulin family member of Paracentrotus lividus (Pl-Tuba1a, formerly known as Pl-Talpha2) that is specifically expressed in the ciliary band and animal pole neurogenic domains of the sea urchin embryo. In order to identify cis-regulatory elements controlling its spatio-temporal expression, we conducted gene transfer experiments, transgene deletions and site specific mutagenesis. Thus, a genomic region of about 2.6 Kb of Pl-Tuba1a, containing four Interspecifically Conserved Regions (ICRs), was identified as responsible for proper gene expression. An enhancer role was ascribed to ICR1 and ICR2, while ICR3 exerted a pivotal role in basal expression, restricting Tuba1a expression to the proper territories of the embryo. Additionally, the mutation of the forkhead box consensus sequence binding site in ICR3 prevented Pl-Tuba1a expression.

The core promoter composition establishes a new dimension in developmental gene networks

Developmental processes are highly dependent on transcriptional regulation by RNA polymerase II, which initiates transcription at the core promoter. The dorsal-ventral gene regulatory network (GRN) includes multiple genes that are activated by different nuclear concentrations of the Dorsal transcription factor along the dorsal-ventral axis. Downstream core promoter element (DPE)-containing genes are conserved and highly prevalent among Dorsal target genes. Moreover, the DPE motif is functional in multiple Dorsal target genes, as mutation of the DPE results in the loss of transcriptional activity. Furthermore, analysis of hybrid enhancer-promoter constructs reveals that the core promoter composition plays a pivotal role in the transcriptional output. Importantly, we provide in vivo evidence that expression driven by the homeotic Antennapedia P2 promoter during Drosophila embryogenesis is dependent on the DPE. Taken together, we propose that transcriptional regulation results from the interplay between enhancers and core promoter composition, thus establishing a novel dimension in developmental GRNs.

Promoter or enhancer, what's the difference? Deconstruction of established distinctions and presentation of a unifying model

Gene transcription is strictly controlled by the interplay of regulatory events at gene promoters and gene-distal regulatory elements called enhancers. Despite extensive studies of enhancers, we still have a very limited understanding of their mechanisms of action and their restricted spatio-temporal activities. A better understanding would ultimately lead to fundamental insights into the control of gene transcription and the action of regulatory genetic variants involved in disease. Here, I review and discuss pros and cons of state-of-the-art genomics methods to localize and infer the activity of enhancers. Among the different approaches, profiling of enhancer RNAs yields the highest specificity and may be superior in detecting in vivo activity. I discuss their apparent similarities to promoters, which challenge the established view of enhancers and promoters as distinct entities, and present a unifying model of regulatory elements in transcriptional regulation, in which activity, transcriptional output and regulatory function is context specific.

Core promoter functions in the regulation of gene expression of Drosophila dorsal target genes

Developmental processes are highly dependent on transcriptional regulation by RNA polymerase II. The RNA polymerase II core promoter is the ultimate target of a multitude of transcription factors that control transcription initiation. Core promoters consist of core promoter motifs, e.g. the initiator, TATA box, and the downstream core promoter element (DPE), which confer specific properties to the core promoter. Here, we explored the importance of core promoter functions in the dorsal-ventral developmental gene regulatory network. This network includes multiple genes that are activated by different nuclear concentrations of Dorsal, an NFκB homolog transcription factor, along the dorsal-ventral axis. We show that over two-thirds of Dorsal target genes contain DPE sequence motifs, which is significantly higher than the proportion of DPE-containing promoters in Drosophila genes. We demonstrate that multiple Dorsal target genes are evolutionarily conserved and functionally dependent on the DPE. Furthermore, we have analyzed the activation of key Dorsal target genes by Dorsal, as well as by another Rel family transcription factor, Relish, and the dependence of their activation on the DPE motif. Using hybrid enhancer-promoter constructs in Drosophila cells and embryo extracts, we have demonstrated that the core promoter composition is an important determinant of transcriptional activity of Dorsal target genes. Taken together, our results provide evidence for the importance of core promoter composition in the regulation of Dorsal target genes.

Quantitative modeling of a gene's expression from its intergenic sequence

Modeling a gene's expression from its intergenic locus and trans-regulatory context is a fundamental goal in computational biology. Owing to the distributed nature of cis-regulatory information and the poorly understood mechanisms that integrate such information, gene locus modeling is a more challenging task than modeling individual enhancers. Here we report the first quantitative model of a gene's expression pattern as a function of its locus. We model the expression readout of a locus in two tiers: 1) combinatorial regulation by transcription factors bound to each enhancer is predicted by a thermodynamics-based model and 2) independent contributions from multiple enhancers are linearly combined to fit the gene expression pattern. The model does not require any prior knowledge about enhancers contributing toward a gene's expression. We demonstrate that the model captures the complex multi-domain expression patterns of anterior-posterior patterning genes in the early Drosophila embryo. Altogether, we model the expression patterns of 27 genes; these include several gap genes, pair-rule genes, and anterior, posterior, trunk, and terminal genes. We find that the model-selected enhancers for each gene overlap strongly with its experimentally characterized enhancers. Our findings also suggest the presence of sequence-segments in the locus that would contribute ectopic expression patterns and hence were "shut down" by the model. We applied our model to identify the transcription factors responsible for forming the stripe boundaries of the studied genes. The resulting network of regulatory interactions exhibits a high level of agreement with known regulatory influences on the target genes. Finally, we analyzed whether and why our assumption of enhancer independence was necessary for the genes we studied. We found a deterioration of expression when binding sites in one enhancer were allowed to influence the readout of another enhancer. Thus, interference between enhancer activities was a possible factor necessitating enhancer independence in our model.

Evolutionary psychology and design reincarnation: A reply to Peters

In this riposte it will be argued that the critique by Peters (2013) of the theoretical foundation of evolutionary psychology misses the mark, and, in the process, unfortunately repeats many common and egregious misunderstandings. This reply will attempt to outline the real position of evolutionary psychologists with respect to modularity and the development, flexibility, and learning capacities of cognitive adaptations. In particular, evolutionary psychologists’ concept of the developmental target of naturally selected design will be made salient. I also aim to provide a more accurate treatment of the neurobiological and genomic issues at stake. In sum, it will be shown that Peters’ rendering of the theoretical foundation of evolutionary psychology is a straw man representation and that the real position of evolutionary psychologists is far more interesting once some of the nuances of their theoretical foundation are brought to light.

Huckebein is part of a combinatorial repression code in the anterior blastoderm

The hierarchy of the segmentation cascade responsible for establishing the Drosophila body plan is composed by gap, pair-rule and segment polarity genes. However, no pair-rule stripes are formed in the anterior regions of the embryo. This lack of stripe formation, as well as other evidence from the literature that is further investigated here, led us to the hypothesis that anterior gap genes might be involved in a combinatorial mechanism responsible for repressing the cis-regulatory modules (CRMs) of hairy (h), even-skipped (eve), runt (run), and fushi-tarazu (ftz) anterior-most stripes. In this study, we investigated huckebein (hkb), which has a gap expression domain at the anterior tip of the embryo. Using genetic methods we were able to detect deviations from the wild-type patterns of the anterior-most pair-rule stripes in different genetic backgrounds, which were consistent with Hkb-mediated repression. Moreover, we developed an image processing tool that, for the most part, confirmed our assumptions. Using an hkb misexpression system, we further detected specific repression on anterior stripes. Furthermore, bioinformatics analysis predicted an increased significance of binding site clusters in the CRMs of h 1, eve 1, run 1 and ftz 1when Hkb was incorporated in the analysis, indicating that Hkb plays a direct role in these CRMs. We further discuss that Hkb and Slp1, which is the other previously identified common repressor of anterior stripes, might participate in a combinatorial repression mechanism controlling stripe CRMs in the anterior parts of the embryo and define the borders of these anterior stripes.

Decoding the genome with an integrative analysis tool: combinatorial CRM Decoder

The identification of genome-wide cis-regulatory modules (CRMs) and characterization of their associated epigenetic features are fundamental steps toward the understanding of gene regulatory networks. Although integrative analysis of available genome-wide information can provide new biological insights, the lack of novel methodologies has become a major bottleneck. Here, we present a comprehensive analysis tool called combinatorial CRM decoder (CCD), which utilizes the publicly available information to identify and characterize genome-wide CRMs in a species of interest. CCD first defines a set of the epigenetic features which is significantly associated with a set of known CRMs as a code called ‘trace code’, and subsequently uses the trace code to pinpoint putative CRMs throughout the genome. Using 61 genome-wide data sets obtained from 17 independent mouse studies, CCD successfully catalogued ∼12 600 CRMs (five distinct classes) including polycomb repressive complex 2 target sites as well as imprinting control regions. Interestingly, we discovered that ∼4% of the identified CRMs belong to at least two different classes named ‘multi-functional CRM’, suggesting their functional importance for regulating spatiotemporal gene expression. From these examples, we show that CCD can be applied to any potential genome-wide datasets and therefore will shed light on unveiling genome-wide CRMs in various species.

A unique H3K4me2 profile marks tissue-specific gene regulation

Characterization of the epigenetic landscape fundamentally contributes toward deciphering the regulatory mechanisms that govern gene expression. However, despite an increasing flow of newly generated data, no clear pattern of chromatin modifications has so far been linked to specific modes of transcriptional regulation. Here, we used high-throughput genomic data from CD4(+) T lymphocytes to provide a comprehensive analysis of histone H3 lysine 4 dimethylation (H3K4me2) enrichment in genomic regions surrounding transcriptional start sites (TSSs). We discovered that a subgroup of genes linked to T cell functions displayed high levels of H3K4me2 within their gene body, in sharp contrast to the TSS-centered profile typical of housekeeping genes. Analysis of additional chromatin modifications and DNase I hypersensitive sites (DHSS) revealed a combinatorial chromatin signature characteristic of this subgroup. We propose that this epigenetic feature reflects the activity of an as yet unrecognized, intragenic cis-regulatory platform dedicated to refining tissue-specificity in gene expression.

Deciphering the genome's regulatory code: the many languages of DNA

The generation of patterns and the diversity of cell types in a multicellular organism require differential gene regulation. At the heart of this process are enhancers or cis-regulatory modules (CRMs), genomic regions that are bound by transcription factors (TFs) that control spatio-temporal gene expression in developmental networks. To date, only a few CRMs have been studied in detail and the underlying cis-regulatory code is not well understood. Here, we review recent progress on the genome-wide identification of CRMs with chromatin immunoprecipitation of TF-DNA complexes followed by microarrays (ChIP-on-chip). We focus on two computational approaches that have succeeded in predicting the expression pattern driven by a CRM either based on TF binding site preferences and their expression levels, or quantitative analysis of CRM occupancy by key TFs. We also discuss the current limits of these methods and highlight some of the key problems that have to be solved to gain a more complete understanding of the structure and function of CRMs.

The developmental expression dynamics of Drosophila melanogaster transcription factors

BACKGROUND

Site-specific transcription factors (TFs) are coordinators of developmental and physiological gene expression programs. Their binding to cis-regulatory modules of target genes mediates the precise cell- and context-specific activation and repression of genes. The expression of TFs should therefore reflect the core expression program of each cell.

RESULTS

We studied the expression dynamics of about 750 TFs using the available genomics resources in Drosophila melanogaster. We find that 95% of these TFs are expressed at some point during embryonic development, with a peak roughly between 10 and 12 hours after egg laying, the core stages of organogenesis. We address the differential utilization of DNA-binding domains in different developmental programs systematically in a spatio-temporal context, and show that the zinc finger class of TFs is predominantly early expressed, while Homeobox TFs exhibit later expression in embryogenesis.

CONCLUSIONS

Previous work, dissecting cis-regulatory modules during Drosophila development, suggests that TFs are deployed in groups acting in a cooperative manner. In contrast, we find that there is rapid exchange of co-expressed partners amongst the fly TFs, at rates similar to the genome-wide dynamics of co-expression clusters. This suggests there may also be a high level of combinatorial complexity of TFs at cis-regulatory modules.

Combinatorial binding predicts spatio-temporal cis-regulatory activity

Development requires the establishment of precise patterns of gene expression, which are primarily controlled by transcription factors binding to cis-regulatory modules. Although transcription factor occupancy can now be identified at genome-wide scales, decoding this regulatory landscape remains a daunting challenge. Here we used a novel approach to predict spatio-temporal cis-regulatory activity based only on in vivo transcription factor binding and enhancer activity data. We generated a high-resolution atlas of cis-regulatory modules describing their temporal and combinatorial occupancy during Drosophila mesoderm development. The binding profiles of cis-regulatory modules with characterized expression were used to train support vector machines to predict five spatio-temporal expression patterns. In vivo transgenic reporter assays demonstrate the high accuracy of these predictions and reveal an unanticipated plasticity in transcription factor binding leading to similar expression. This data-driven approach does not require previous knowledge of transcription factor sequence affinity, function or expression, making it widely applicable.

Design flexibility in cis-regulatory control of gene expression: synthetic and comparative evidence

In early Drosophila embryos, the transcription factor Dorsal regulates patterns of gene expression and cell fate specification along the dorsal-ventral axis. How gene expression is produced within the broad lateral domain of the presumptive neurogenic ectoderm is not understood. To investigate transcriptional control during neurogenic ectoderm specification, we examined divergence and function of an embryonic cis-regulatory element controlling the gene short gastrulation (sog). While transcription factor binding sites are not completely conserved, we demonstrate that these sequences are bona fide regulatory elements, despite variable regulatory architecture. Mutation of conserved sequences revealed that putative transcription factor binding sites for Dorsal and Zelda, a ubiquitous maternal transcription factor, are required for proper sog expression. When Zelda and Dorsal sites are paired in a synthetic regulatory element, broad lateral expression results. However, synthetic regulatory elements that contain Dorsal and an additional activator also drive expression throughout the neurogenic ectoderm. Our results suggest that interaction between Dorsal and Zelda drives expression within the presumptive neurogenic ectoderm, but they also demonstrate that regulatory architecture directing expression in this domain is flexible. We propose a model for neurogenic ectoderm specification in which gene regulation occurs at the intersection of temporal and spatial transcription factor inputs.

Caudal, a key developmental regulator, is a DPE-specific transcriptional factor

The regulation of gene transcription is critical for the proper development and growth of an organism. The transcription of protein-coding genes initiates at the RNA polymerase II core promoter, which is a diverse module that can be controlled by many different elements such as the TATA box and downstream core promoter element (DPE). To understand the basis for core promoter diversity, we explored potential biological functions of the DPE. We found that nearly all of the Drosophila homeotic (Hox) gene promoters, which lack TATA-box elements, contain functionally important DPE motifs that are conserved from Drosophila melanogaster to Drosophila virilis. We then discovered that Caudal, a sequence-specific transcription factor and key regulator of the Hox gene network, activates transcription with a distinct preference for the DPE relative to the TATA box. The specificity of Caudal activation for the DPE is particularly striking when a BRE(u) core promoter motif is associated with the TATA box. These findings show that Caudal is a DPE-specific activator and exemplify how core promoter diversity can be used to establish complex regulatory networks.

Pattern formation in the vertebrate neural tube: a sonic hedgehog morphogen-regulated transcriptional network

Neuronal subtype specification in the vertebrate neural tube is one of the best-studied examples of embryonic pattern formation. Distinct neuronal subtypes are generated in a precise spatial order from progenitor cells according to their location along the anterior-posterior and dorsal-ventral axes. Underpinning this organization is a complex network of multiple extrinsic and intrinsic factors. This review focuses on the molecular mechanisms and general strategies at play in ventral regions of the forming spinal cord, where sonic hedgehog-based morphogen signaling is a key determinant. We discuss recent advances in our understanding of these events and highlight unresolved questions.

Genome-wide analysis of mRNA decay rates and their determinants in Arabidopsis thaliana

To gain a global view of mRNA decay in Arabidopsis thaliana, suspension cell cultures were treated with a transcriptional inhibitor, and microarrays were used to measure transcript abundance over time. The deduced mRNA half-lives varied widely, from minutes to >24 h. Three features of the transcript displayed a correlation with decay rates: (1) genes possessing at least one intron produce mRNA transcripts significantly more stable than those of intronless genes, and this was not related to overall length, sequence composition, or number of introns; (2) various sequence elements in the 3' untranslated region are enriched among short- and long-lived transcripts, and their multiple occurrence suggests combinatorial control of transcript decay; and (3) transcripts that are microRNA targets generally have short half-lives. The decay rate of transcripts correlated with subcellular localization and function of the encoded proteins. Analysis of transcript decay rates for genes encoding orthologous proteins between Arabidopsis, yeast, and humans indicated that yeast and humans had a higher percentage of transcripts with shorter half-lives and that the relative stability of transcripts from genes encoding proteins involved in cell cycle, transcription, translation, and energy metabolism is conserved. Comparison of decay rates with changes in transcript abundance under a variety of abiotic stresses reveal that a set of transcription factors are downregulated with similar kinetics to decay rates, suggesting that inhibition of their transcription is an important early response to abiotic stress.

Quantitative expression profiling of identified neurons reveals cell-specific constraints on highly variable levels of gene expression

The postdevelopmental basis of cellular identity and the unique cellular output of a particular neuron type are of particular interest in the nervous system because a detailed understanding of circuits responsible for complex processes in the brain is impeded by the often ambiguous classification of neurons in these circuits. Neurons have been classified by morphological, electrophysiological, and neurochemical techniques. More recently, molecular approaches, particularly microarray, have been applied to the question of neuronal identity. With the realization that proteins expressed exclusively in only one type of neuron are rare, expression profiles obtained from neuronal subtypes are analyzed to search for diagnostic patterns of gene expression. However, this expression profiling hinges on one critical and implicit assumption: that neurons of the same type in different animals achieve their conserved functional output via conserved levels and quantitative relationships of gene expression. Here we exploit the unambiguously identifiable neurons in the crab stomatogastric ganglion to investigate the precise quantitative expression profiling of neurons at the level of single-cell ion channel expression. By measuring absolute mRNA levels of six different channels in the same individually identified neurons, we demonstrate that not only do individual cell types possess highly variable levels of channel expression but that this variability is constrained by unique patterns of correlated channel expression.

Clustering of phosphorylation site recognition motifs can be exploited to predict the targets of cyclin-dependent kinase

Protein kinases are critical to cellular signalling and post-translational gene regulation, but their biological substrates are difficult to identify. We show that cyclin-dependent kinase (CDK) consensus motifs are frequently clustered in CDK substrate proteins. Based on this, we introduce a new computational strategy to predict the targets of CDKs and use it to identify new biologically interesting candidates. Our data suggest that regulatory modules may exist in protein sequence as clusters of short sequence motifs.

A multilayered defense against infection: combinatorial control of insect immune genes

The innate immune defense system involves the activity of endogenous antimicrobial peptides (AMPs), which inhibit the growth of most microbes. In insects, genes encoding AMPs are expressed at basal levels in barrier epithelia and are upregulated systemically in response to infection. To achieve this differentiated immune defense, Drosophila immune gene promoters combine tissue-specific enhancers and signal-dependent response elements. Transcription factors of the Hox, POU and GATA families control tissue-specific expression of AMP genes, either constitutively or in combination with NF-kappaB/Rel family factors that function as 'on-off switches' during infection. Here, we review these different modes of AMP expression and provide a model for transcriptional regulation of AMP genes.

Identification of cis-regulatory elements from the C. elegans Hox gene lin-39 required for embryonic expression and for regulation by the transcription factors LIN-1, LIN-31 and LIN-39

Expression of the Caenorhabditis elegans Hox gene lin-39 begins in the embryo and continues in multiple larval cells, including the P cell lineages that generate ventral cord neurons (VCNs) and vulval precursor cells (VPCs). lin-39 is regulated by several factors and by Wnt and Ras signaling pathways; however, no cis-acting sites mediating lin-39 regulation have been identified. Here, we describe three elements controlling lin-39 expression: a 338-bp upstream fragment that directs embryonic expression in P5-P8 and their descendants in the larva, a 247-bp intronic region sufficient for VCN expression, and a 1.3-kb upstream cis-regulatory module that drives expression in the VPC P6.p in a Ras-dependent manner. Three trans-acting factors regulate expression via the 1.3-kb element. A single binding site for the ETS factor LIN-1 mediates repression in VPCs other than P6.p; however, loss of LIN-1 decreases expression in P6.p. Therefore, LIN-1 acts both negatively and positively on lin-39 in different VPCs. The Forkhead domain protein LIN-31 also acts positively on lin-39 in P6.p via this module. Finally, LIN-39 itself binds to this element, suggesting that LIN-39 autoregulates its expression in P6.p. Therefore, we have begun to unravel the cis-acting sites regulating lin-39 Hox gene expression and have shown that lin-39 is a direct target of the Ras pathway acting via LIN-1 and LIN-31.