Comparative analysis of regulatory information and circuits across distant species.

Despite the large evolutionary distances between metazoan species, they can show remarkable commonalities in their biology, and this has helped to establish fly and worm as model organisms for human biology. Although studies of individual elements and factors have explored similarities in gene regulation, a large-scale comparative analysis of basic principles of transcriptional regulatory features is lacking. Here we map the genome-wide binding locations of 165 human, 93 worm and 52 fly transcription regulatory factors, generating a total of 1,019 data sets from diverse cell types, developmental stages, or conditions in the three species, of which 498 (48.9%) are presented here for the first time. We find that structural properties of regulatory networks are remarkably conserved and that orthologous regulatory factor families recognize similar binding motifs in vivo and show some similar co-associations. Our results suggest that gene-regulatory properties previously observed for individual factors are general principles of metazoan regulation that are remarkably well-preserved despite extensive functional divergence of individual network connections. The comparative maps of regulatory circuitry provided here will drive an improved understanding of the regulatory underpinnings of model organism biology and how these relate to human biology, development and disease.

Evolution of H3K27me3-marked chromatin is linked to gene expression evolution and to patterns of gene duplication and diversification.

Histone modifications are critical for the regulation of gene expression, cell type specification, and differentiation. However, evolutionary patterns of key modifications that regulate gene expression in differentiating organisms have not been examined. Here we mapped the genomic locations of the repressive mark histone 3 lysine 27 trimethylation (H3K27me3) in four species of Drosophila, and compared these patterns to those in C. elegans. We found that patterns of H3K27me3 are highly conserved across species, but conservation is substantially weaker among duplicated genes. We further discovered that retropositions are associated with greater evolutionary changes in H3K27me3 and gene expression than tandem duplications, indicating that local chromatin constraints influence duplicated gene evolution. These changes are also associated with concomitant evolution of gene expression. Our findings reveal the strong conservation of genomic architecture governed by an epigenetic mark across distantly related species and the importance of gene duplication in generating novel H3K27me3 profiles.

Diverse patterns of genomic targeting by transcriptional regulators in Drosophila melanogaster.

Annotation of regulatory elements and identification of the transcription-related factors (TRFs) targeting these elements are key steps in understanding how cells interpret their genetic blueprint and their environment during development, and how that process goes awry in the case of disease. One goal of the modENCODE (model organism ENCyclopedia of DNA Elements) Project is to survey a diverse sampling of TRFs, both DNA-binding and non-DNA-binding factors, to provide a framework for the subsequent study of the mechanisms by which transcriptional regulators target the genome. Here we provide an updated map of the Drosophila melanogaster regulatory genome based on the location of 84 TRFs at various stages of development. This regulatory map reveals a variety of genomic targeting patterns, including factors with strong preferences toward proximal promoter binding, factors that target intergenic and intronic DNA, and factors with distinct chromatin state preferences. The data also highlight the stringency of the Polycomb regulatory network, and show association of the Trithorax-like (Trl) protein with hotspots of DNA binding throughout development. Furthermore, the data identify more than 5800 instances in which TRFs target DNA regions with demonstrated enhancer activity. Regions of high TRF co-occupancy are more likely to be associated with open enhancers used across cell types, while lower TRF occupancy regions are associated with complex enhancers that are also regulated at the epigenetic level. Together these data serve as a resource for the research community in the continued effort to dissect transcriptional regulatory mechanisms directing Drosophila development.

A cis-regulatory map of the Drosophila genome.

Systematic annotation of gene regulatory elements is a major challenge in genome science. Direct mapping of chromatin modification marks and transcriptional factor binding sites genome-wide has successfully identified specific subtypes of regulatory elements. In Drosophila several pioneering studies have provided genome-wide identification of Polycomb response elements, chromatin states, transcription factor binding sites, RNA polymerase II regulation and insulator elements; however, comprehensive annotation of the regulatory genome remains a significant challenge. Here we describe results from the modENCODE cis-regulatory annotation project. We produced a map of the Drosophila melanogaster regulatory genome on the basis of more than 300 chromatin immunoprecipitation data sets for eight chromatin features, five histone deacetylases and thirty-eight site-specific transcription factors at different stages of development. Using these data we inferred more than 20,000 candidate regulatory elements and validated a subset of predictions for promoters, enhancers and insulators in vivo. We identified also nearly 2,000 genomic regions of dense transcription factor binding associated with chromatin activity and accessibility. We discovered hundreds of new transcription factor co-binding relationships and defined a transcription factor network with over 800 potential regulatory relationships.

Evidence of abundant stop codon readthrough in Drosophila and other metazoa.

While translational stop codon readthrough is often used by viral genomes, it has been observed for only a handful of eukaryotic genes. We previously used comparative genomics evidence to recognize protein-coding regions in 12 species of Drosophila and showed that for 149 genes, the open reading frame following the stop codon has a protein-coding conservation signature, hinting that stop codon readthrough might be common in Drosophila. We return to this observation armed with deep RNA sequence data from the modENCODE project, an improved higher-resolution comparative genomics metric for detecting protein-coding regions, comparative sequence information from additional species, and directed experimental evidence. We report an expanded set of 283 readthrough candidates, including 16 double-readthrough candidates; these were manually curated to rule out alternatives such as A-to-I editing, alternative splicing, dicistronic translation, and selenocysteine incorporation. We report experimental evidence of translation using GFP tagging and mass spectrometry for several readthrough regions. We find that the set of readthrough candidates differs from other genes in length, composition, conservation, stop codon context, and in some cases, conserved stem-loops, providing clues about readthrough regulation and potential mechanisms. Lastly, we expand our studies beyond Drosophila and find evidence of abundant readthrough in several other insect species and one crustacean, and several readthrough candidates in nematode and human, suggesting that functionally important translational stop codon readthrough is significantly more prevalent in Metazoa than previously recognized.

Effect of D. melanogaster larval density on pupal size.

Many genetic pathways and environmental factors have been shown to affect Drosophila melanogaster adult body size. Larval density often varies considerably between vials, even when the same number of females of the same genotype are allowed to lay eggs in the vials for the same amount of time. To more accurately quantify the effects that larval population density has on pupal size, we established cultures of 1, 2, 10, 25, 50, 75 or 100 first instar larvae into vials and measured pupal length. We collected Oregon-R eggs on apple juice plates in six different cages and generated replicate cultures. We found that pupal size decreases as larval density in the culture increases by 25 individuals. The difference between male and female length remained relatively constant at each density (0.2 mm), but overall size decreased. The mean size differences between vials with 1 larvae and 100 larvae is 0.1(+/-0.02) mm in females and 0.11(+/-0.02) mm in males. These results suggest that fecundity and sex ratio could complicate results in Drosophila size studies.

Versatile P[acman] BAC libraries for transgenesis studies in Drosophila melanogaster.

We constructed Drosophila melanogaster bacterial artificial chromosome libraries with 21-kilobase and 83-kilobase inserts in the P[acman] system. We mapped clones representing 12-fold coverage and encompassing more than 95% of annotated genes onto the reference genome. These clones can be integrated into predetermined attP sites in the genome using UC31 integrase to rescue mutations. They can be modified through recombineering, for example, to incorporate protein tags and assess expression patterns.

Centrosomin represses dendrite branching by orienting microtubule nucleation.

Neuronal dendrite branching is fundamental for building nervous systems. Branch formation is genetically encoded by transcriptional programs to create dendrite arbor morphological diversity for complex neuronal functions. In Drosophila sensory neurons, the transcription factor Abrupt represses branching via an unknown effector pathway. Targeted screening for branching-control effectors identified Centrosomin, the primary centrosome-associated protein for mitotic spindle maturation. Centrosomin repressed dendrite branch formation and was used by Abrupt to simplify arbor branching. Live imaging revealed that Centrosomin localized to the Golgi cis face and that it recruited microtubule nucleation to Golgi outposts for net retrograde microtubule polymerization away from nascent dendrite branches. Removal of Centrosomin enabled the engagement of wee Augmin activity to promote anterograde microtubule growth into the nascent branches, leading to increased branching. The findings reveal that polarized targeting of Centrosomin to Golgi outposts during elaboration of the dendrite arbor creates a local system for guiding microtubule polymerization.

Pri peptides are mediators of ecdysone for the temporal control of development.

Animal development fundamentally relies on the precise control, in space and time, of genome expression. Whereas we have a wealth of information about spatial patterning, the mechanisms underlying temporal control remain poorly understood. Here we show that Pri peptides, encoded by small open reading frames, are direct mediators of the steroid hormone ecdysone for the timing of developmental programs in Drosophila. We identify a previously uncharacterized enzyme of ecdysone biosynthesis, GstE14, and find that ecdysone triggers pri expression to define the onset of epidermal trichome development, through post-translational control of the Shavenbaby transcription factor. We show that manipulating pri expression is sufficient to either put on hold or induce premature differentiation of trichomes. Furthermore, we find that ecdysone-dependent regulation of pri is not restricted to epidermis and occurs over various tissues and times. Together, these findings provide a molecular framework to explain how systemic hormonal control coordinates specific programs of differentiation with developmental timing.

Myc-dependent genome instability and lifespan in Drosophila.

The Myc family of transcription factors are key regulators of cell growth and proliferation that are dysregulated in a large number of human cancers. When overexpressed, Myc family proteins also cause genomic instability, a hallmark of both transformed and aging cells. Using an in vivo lacZ mutation reporter, we show that overexpression of Myc in Drosophila increases the frequency of large genome rearrangements associated with erroneous repair of DNA double-strand breaks (DSBs). In addition, we find that overexpression of Myc shortens adult lifespan and, conversely, that Myc haploinsufficiency reduces mutation load and extends lifespan. Our data provide the first evidence that Myc may act as a pro-aging factor, possibly through its ability to greatly increase genome instability.

Genome-wide analyses of Shavenbaby target genes reveals distinct features of enhancer organization.

BACKGROUND: Developmental programs are implemented by regulatory interactions between Transcription Factors (TFs) and their target genes, which remain poorly understood. While recent studies have focused on regulatory cascades of TFs that govern early development, little is known about how the ultimate effectors of cell differentiation are selected and controlled. We addressed this question during late Drosophila embryogenesis, when the finely tuned expression of the TF Ovo/Shavenbaby (Svb) triggers the morphological differentiation of epidermal trichomes. RESULTS: We defined a sizeable set of genes downstream of Svb and used in vivo assays to delineate 14 enhancers driving their specific expression in trichome cells. Coupling computational modeling to functional dissection, we investigated the regulatory logic of these enhancers. Extending the repertoire of epidermal effectors using genome-wide approaches showed that the regulatory models learned from this first sample are representative of the whole set of trichome enhancers. These enhancers harbor remarkable features with respect to their functional architectures, including a weak or non-existent clustering of Svb binding sites. The in vivo function of each site relies on its intimate context, notably the flanking nucleotides. Two additional cis-regulatory motifs, present in a broad diversity of composition and positioning among trichome enhancers, critically contribute to enhancer activity. CONCLUSIONS: Our results show that Svb directly regulates a large set of terminal effectors of the remodeling of epidermal cells. Further, these data reveal that trichome formation is underpinned by unexpectedly diverse modes of regulation, providing fresh insights into the functional architecture of enhancers governing a terminal differentiation program.

Broad Complex isoforms have unique distributions during central nervous system metamorphosis in Drosophila melanogaster.

Broad Complex (BRC) is a highly conserved, ecdysone-pathway gene essential for metamorphosis in Drosophila melanogaster, and possibly all holometabolous insects. Alternative splicing among duplicated exons produces several BRC isoforms, each with one zinc-finger DNA-binding domain (Z1, Z2, Z3, or Z4), highly expressed at the onset of metamorphosis. BRC-Z1, BRC-Z2, and BRC-Z3 represent distinct genetic functions (BRC complementation groups rbp, br, and 2Bc, respectively) and are required at discrete stages spanning final-instar larva through very young pupa. We showed previously that morphogenetic movements necessary for adult CNS maturation require BRC-Z1, -Z2, and -Z3, but not at the same time: BRC-Z1 is required in the mid-prepupa, BRC-Z2 and -Z3 are required earlier, at the larval-prepupal transition. To explore how BRC isoforms controlling the same morphogenesis events do so at different times, we examined their central nervous system (CNS) expression patterns during the approximately 16 hours bracketing the hormone-regulated start of metamorphosis. Each isoform had a unique pattern, with BRC-Z3 being the most distinctive. There was some colocalization of isoform pairs, but no three-way overlap of BRC-Z1, -Z2, and -Z3. Instead, their most prominent expression was in glia (BRC-Z1), neuroblasts (BRC-Z2), or neurons (BRC-Z3). Despite sequence similarity to BRC-Z1, BRC-Z4 was expressed in a unique subset of neurons. These data suggest a switch in BRC isoform choice, from BRC-Z2 in proliferating cells to BRC-Z1, BRC-Z3, or BRC-Z4 in differentiating cells. Together with isoform-selective temporal requirements and phenotype considerations, this cell-type-selective expression suggests a model of BRC-dependent CNS morphogenesis resulting from intercellular interactions, culminating in BRC-Z1-controlled, glia-mediated CNS movements in late prepupa.

Anciently duplicated Broad Complex exons have distinct temporal functions during tissue morphogenesis.

Broad Complex (BRC) is an essential ecdysone-pathway gene required for entry into and progression through metamorphosis in Drosophila melanogaster. Mutations of three BRC complementation groups cause numerous phenotypes, including a common suite of morphogenesis defects involving central nervous system (CNS), adult salivary glands (aSG), and male genitalia. These defects are phenocopied by the juvenile hormone mimic methoprene. Four BRC isoforms are produced by alternative splicing of a protein-binding BTB/POZ-encoding exon (BTBBRC) to one of four tandemly duplicated, DNA-binding zinc-finger-encoding exons (Z1BRC, Z2BRC, Z3BRC, Z4BRC). Highly conserved orthologs of BTBBRC and all four ZBRC were found among published cDNA sequences or genome databases from Diptera, Lepidoptera, Hymenoptera, and Coleoptera, indicating that BRC arose and underwent internal exon duplication before the split of holometabolous orders. Tramtrack subfamily members, abrupt, tramtrack, fruitless, longitudinals lacking (lola), and CG31666 were characterized throughout Holometabola and used to root phylogenetic analyses of ZBRC exons, which revealed that the ZBRC clade includes Zabrupt. All four ZBRC domains, including Z4BRC, which has no known essential function, are evolving in a manner consistent with selective constraint. We used transgenic rescue to explore how different BRC isoforms contribute to shared tissue-morphogenesis functions. As predicted from earlier studies, the common CNS and aSG phenotypes were rescued by BRC-Z1 in rbp mutants, BRC-Z2 in br mutants, and BRC-Z3 in 2Bc mutants. However, the isoforms are required at two different developmental stages, with BRC-Z2 and -Z3 required earlier than BRC-Z1. The sequential action of BRC isoforms indicates subfunctionalization of duplicated ZBRC exons even when they contribute to common developmental processes.

The transcription factor Sox9 is required for cranial neural crest development in Xenopus.

The SOX family of transcription factors has been implicated in cell fate specification during embryogenesis. One member of this family, Sox9, has been shown to regulate both chondrogenesis and sex determination in the mouse embryo. Heterozygous mutations in Sox9 result in Campomelic Dysplasia (CD), a lethal human disorder characterized by autosomal XY sex reversal, severe skeletal malformations and several craniofacial defects. Sox9 is also expressed in neural crest progenitors but very little is known about the function of Sox9 in the neural crest. We have cloned the Xenopus homolog of the Sox9 gene. It is expressed maternally and accumulates shortly after gastrulation at the lateral edges of the neural plate, in the neural crest-forming region. As development proceeds, Sox9 expression persists in migrating cranial crest cells as they populate the pharyngeal arches. Depletion of Sox9 protein in developing embryos, using morpholino antisense oligos, causes a dramatic loss of neural crest progenitors and an expansion of the neural plate. Later during embryogenesis, morpholino-treated embryos have a specific loss or reduction of neural crest-derived skeletal elements, mimicking one aspect of the craniofacial defects observed in CD patients. We propose that Sox9 is an essential component of the regulatory pathway that leads to cranial neural crest formation.