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.

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.