Robert P. Zinzen

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.

Tissue-specific analysis of chromatin state identifies temporal signatures of enhancer activity during embryonic development

Chromatin modifications are associated with many aspects of gene expression, yet their role in cellular transitions during development remains elusive. Here, we use a new approach to obtain cell type–specific information on chromatin state and RNA polymerase II (Pol II) occupancy within the multicellular Drosophila melanogaster embryo. We directly assessed the relationship between chromatin modifications and the spatio-temporal activity of enhancers. Rather than having a unique chromatin state, active developmental enhancers show heterogeneous histone modifications and Pol II occupancy. Despite this complexity, combined chromatin signatures and Pol II presence are sufficient to predict enhancer activity de novo. Pol II recruitment is highly predictive of the timing of enhancer activity and seems dependent on the timing and location of transcription factor binding. Chromatin modifications typically demarcate large regulatory regions encompassing multiple enhancers, whereas local changes in nucleosome positioning and Pol II occupancy delineate single active enhancers. This cell type–specific view identifies dynamic enhancer usage, an essential step in deciphering developmental networks.

Computational Models for Neurogenic Gene Expression in the Drosophila Embryo

The early Drosophila embryo is emerging as a premiere model system for the computational analysis of gene regulation in development because most of the genes, and many of the associated regulatory DNAs, that control segmentation and gastrulation are known. The comprehensive elucidation of Drosophila gene networks provides an unprecedented opportunity to apply quantitative models to metazoan enhancers that govern complex patterns of gene expression during development. Models based on the fractional occupancy of defined DNA binding sites have been used to describe the regulation of the lac operon in E. coli and the lysis/lysogeny switch of phage lambda. Here, we apply similar models to enhancers regulated by the Dorsal gradient in the ventral neurogenic ectoderm (vNE) of the early Drosophila embryo. Quantitative models based on the fractional occupancy of Dorsal, Twist, and Snail binding sites raise the possibility that cooperative interactions among these regulatory proteins mediate subtle differences in the vNE expression patterns. Variations in cooperativity may be attributed to differences in the detailed linkage of Dorsal, Twist, and Snail binding sites in vNE enhancers. We propose that binding site occupancy is the key rate-limiting step for establishing localized patterns of gene expression in the early Drosophila embryo.

A regulatory code for neurogenic gene expression in the Drosophila embryo

Bioinformatics methods have identified enhancers that mediate restricted expression in the Drosophila embryo. However, only a small fraction of the predicted enhancers actually work when tested in vivo. In the present study, co-regulated neurogenic enhancers that are activated by intermediate levels of the Dorsal regulatory gradient are shown to contain several shared sequence motifs. These motifs permitted the identification of new neurogenic enhancers with high precision: five out of seven predicted enhancers direct restricted expression within ventral regions of the neurogenic ectoderm. Mutations in some of the shared motifs disrupt enhancer function, and evidence is presented that the Twist and Su(H) regulatory proteins are essential for the specification of the ventral neurogenic ectoderm prior to gastrulation. The regulatory model of neurogenic gene expression defined in this study permitted the identification of a neurogenic enhancer in the distant Anopheles genome. We discuss the prospects for deciphering regulatory codes that link primary DNA sequence information with predicted patterns of gene expression.

Evolution of the dorsal-ventral patterning network in the mosquito, Anopheles gambiae.

The dorsal-ventral patterning of the Drosophila embryo is controlled by a well-defined gene regulation network. We wish to understand how changes in this network produce evolutionary diversity in insect gastrulation. The present study focuses on the dorsal ectoderm in two highly divergent dipterans, the fruitfly Drosophila melanogaster and the mosquito Anopheles gambiae. In D. melanogaster, the dorsal midline of the dorsal ectoderm forms a single extra-embryonic membrane, the amnioserosa. In A. gambiae, an expanded domain forms two distinct extra-embryonic tissues, the amnion and serosa. The analysis of approximately 20 different dorsal-ventral patterning genes suggests that the initial specification of the mesoderm and ventral neurogenic ectoderm is highly conserved in flies and mosquitoes. By contrast, there are numerous differences in the expression profiles of genes active in the dorsal ectoderm. Most notably, the subdivision of the extra-embryonic domain into separate amnion and serosa lineages in A. gambiae correlates with novel patterns of gene expression for several segmentation repressors. Moreover, the expanded amnion and serosa anlage correlates with a broader domain of Dpp signaling as compared with the D. melanogaster embryo. Evidence is presented that this expanded signaling is due to altered expression of the sog gene.

Cell type-specific chromatin immunoprecipitation from multicellular complex samples using BiTS-ChIP

This protocol describes the batch isolation of tissue-specific chromatin for immunoprecipitation (BiTS-ChIP) for analysis of histone modifications, transcription factor binding, or polymerase occupancy within the context of a multicellular organism or tissue. Embryos expressing a cell type–specific nuclear marker are formaldehyde cross-linked and then subjected to dissociation. Fixed nuclei are isolated and sorted using FACS on the basis of the cell type–specific nuclear marker. Tissue-specific chromatin is extracted, sheared by sonication and used for ChIP-seq or other analyses. The key advantages of this method are the covalent cross-linking before embryo dissociation, which preserves the transcriptional context, and the use of FACS of nuclei, yielding very high purity. The protocol has been optimized for Drosophila, but with minor modifications should be applicable to any model system. The full protocol, including sorting, immunoprecipitation and generation of sequencing libraries, can be completed within 5 d.

The Drosophila embryo at single-cell transcriptome resolution

3D gene expression blueprint of the fly When looking at populations of cells, features such as cell heterogeneity and localization are masked. However, single-cell sequencing reveals cellular heterogeneity and rare cell types. At the onset of gastrulation, the fly embryo consists of about 6000 cells with distinct gene expression profiles. Karaiskos et al. developed an algorithm to generate an interactive three-dimensional (3D) “virtual embryo,” with the expression of more than 8000 genes per cell measured for most cells (see the Perspective by Stadler and Eisen). The virtual embryo offers insights into developmental mechanisms—from local expression of regulators such as transcription factors and long noncoding RNAs to spatial modulation of signaling pathways. Science, this issue p. 194; see also p. 172 Single-cell sequencing and an interactive database enable the generation of a 3D virtual fly embryo and a gene expression blueprint. By the onset of morphogenesis, Drosophila embryos consist of about 6000 cells that express distinct gene combinations. Here, we used single-cell sequencing of precisely staged embryos and devised DistMap, a computational mapping strategy to reconstruct the embryo and to predict spatial gene expression approaching single-cell resolution. We produced a virtual embryo with about 8000 expressed genes per cell. Our interactive Drosophila Virtual Expression eXplorer (DVEX) database generates three-dimensional virtual in situ hybridizations and computes gene expression gradients. We used DVEX to uncover patterned expression of transcription factors and long noncoding RNAs, as well as signaling pathway components. Spatial regulation of Hippo signaling during early embryogenesis suggests a mechanism for establishing asynchronous cell proliferation. Our approach is suitable to generate transcriptomic blueprints for other complex tissues.

Enhancer responses to similarly distributed antagonistic gradients in development.

Formation of spatial gene expression patterns in development depends on transcriptional responses mediated by gene control regions, enhancers. Here, we explore possible responses of enhancers to overlapping gradients of antagonistic transcriptional regulators in the Drosophila embryo. Using quantitative models based on enhancer structure, we demonstrate how a pair of antagonistic transcription factor gradients with similar or even identical spatial distributions can lead to the formation of distinct gene expression domains along the embryo axes. The described mechanisms are sufficient to explain the formation of the anterior and the posterior knirps expression, the posterior hunchback expression domain, and the lateral stripes of rhomboid expression and of other ventral neurogenic ectodermal genes. The considered principles of interaction between antagonistic gradients at the enhancer level can also be applied to diverse developmental processes, such as domain specification in imaginal discs, or even eyespot pattern formation in the butterfly wing.

MiR-184 regulates pancreatic β-cell function according to glucose metabolism.

Background: Upon entering the pancreatic β-cell, glucose is metabolized to ultimately induce both proliferation and the release of insulin. Results: miR-184 targets Argonaute2 to impact the microRNA pathway according to glucose metabolism. Conclusion: miR-184 is a highly regulated microRNA impacting the growth and function of the β-cell. Significance: These results highlight the adaptive role of the microRNA pathway based on metabolic state. In response to fasting or hyperglycemia, the pancreatic β-cell alters its output of secreted insulin; however, the pathways governing this adaptive response are not entirely established. Although the precise role of microRNAs (miRNAs) is also unclear, a recurring theme emphasizes their function in cellular stress responses. We recently showed that miR-184, an abundant miRNA in the β-cell, regulates compensatory proliferation and secretion during insulin resistance. Consistent with previous studies showing miR-184 suppresses insulin release, expression of this miRNA was increased in islets after fasting, demonstrating an active role in the β-cell as glucose levels lower and the insulin demand ceases. Additionally, miR-184 was negatively regulated upon the administration of a sucrose-rich diet in Drosophila, demonstrating strong conservation of this pathway through evolution. Furthermore, miR-184 and its target Argonaute2 remained inversely correlated as concentrations of extracellular glucose increased, underlining a functional relationship between this miRNA and its targets. Lastly, restoration of Argonaute2 in the presence of miR-184 rescued suppression of miR-375-targeted genes, suggesting these genes act in a coordinated manner during changes in the metabolic context. Together, these results highlight the adaptive role of miR-184 according to glucose metabolism and suggest the regulatory role of this miRNA in energy homeostasis is highly conserved.

Template Requirements for Telomerase Translocation in Kluyveromyces lactis

ABSTRACT Telomeres are synthesized by telomerase, a specialized reverse transcriptase, which contains a template in its intrinsic RNA component. In Kluyveromyces lactis, the repeats synthesized by the wild-type telomerase are 25 nucleotides (nt) in length and uniform in sequence. To determine the role of the 5-nt repeats defining the ends of the K. lactis telomerase RNA template in telomerase translocation, we have made mutations in and around them and observed their effects on telomere length and the sequence of newly made telomeric repeats. These template mutations typically result in telomeres that are shorter than those of wild-type cells. The mismatches between the telomerase template and the telomeric tip that occur after telomerase-mediated incorporation of the mutations are normally not removed. Instead, the mutations lead to the synthesis of aberrant repeats that range in size from 31 to 13 bp. Therefore, the specificity with which the telomeric tip aligns with the telomere is critical for the production of the uniform repeats seen in K. lactis. In addition, the region immediately 3′ of the template may play an important role in translocation of the enzyme.