David Sturgill

The developmental transcriptome of Drosophila melanogaster

Drosophila melanogaster is one of the most well studied genetic model organisms, nonetheless its genome still contains unannotated coding and non-coding genes, transcripts, exons, and RNA editing sites. Full discovery and annotation are prerequisites for understanding how the regulation of transcription, splicing, and RNA editing directs development of this complex organism. We used RNA-Seq, tiling microarrays, and cDNA sequencing to explore the transcriptome in 30 distinct developmental stages. We identified 111,195 new elements, including thousands of genes, coding and non-coding transcripts, exons, splicing and editing events and inferred protein isoforms that previously eluded discovery using established experimental, prediction and conservation-based approaches. Together, these data substantially expand the number of known transcribed elements in the Drosophila genome and provide a high-resolution view of transcriptome dynamics throughout development.

Identification of functional elements and regulatory circuits by Drosophila modENCODE

From Genome to Regulatory Networks For biologists, having a genome in hand is only the beginning—much more investigation is still needed to characterize how the genome is used to help to produce a functional organism (see the Perspective by Blaxter). In this vein, Gerstein et al. (p. 1775) summarize for the Caenorhabditis elegans genome, and The modENCODE Consortium (p. 1787) summarize for the Drosophila melanogaster genome, full transcriptome analyses over developmental stages, genome-wide identification of transcription factor binding sites, and high-resolution maps of chromatin organization. Both studies identified regions of the nematode and fly genomes that show highly occupied targets (or HOT) regions where DNA was bound by more than 15 of the transcription factors analyzed and the expression of related genes were characterized. Overall, the studies provide insights into the organization, structure, and function of the two genomes and provide basic information needed to guide and correlate both focused and genome-wide studies. The Drosophila modENCODE project demonstrates the functional regulatory network of flies. To gain insight into how genomic information is translated into cellular and developmental programs, the Drosophila model organism Encyclopedia of DNA Elements (modENCODE) project is comprehensively mapping transcripts, histone modifications, chromosomal proteins, transcription factors, replication proteins and intermediates, and nucleosome properties across a developmental time course and in multiple cell lines. We have generated more than 700 data sets and discovered protein-coding, noncoding, RNA regulatory, replication, and chromatin elements, more than tripling the annotated portion of the Drosophila genome. Correlated activity patterns of these elements reveal a functional regulatory network, which predicts putative new functions for genes, reveals stage- and tissue-specific regulators, and enables gene-expression prediction. Our results provide a foundation for directed experimental and computational studies in Drosophila and related species and also a model for systematic data integration toward comprehensive genomic and functional annotation.

Constraint and turnover in sex-biased gene expression in the genus Drosophila

Both genome content and deployment contribute to phenotypic differences between species. Sex is the most important difference between individuals in a species and has long been posited to be rapidly evolving. Indeed, in the Drosophila genus, traits such as sperm length, genitalia, and gonad size are the most obvious differences between species. Comparative analysis of sex-biased expression should deepen our understanding of the relationship between genome content and deployment during evolution. Using existing and newly assembled genomes, we designed species-specific microarrays to examine sex-biased expression of orthologues and species-restricted genes in D. melanogaster, D. simulans, D. yakuba, D. ananassae, D. pseudoobscura, D. virilis and D. mojavensis. We show that averaged sex-biased expression changes accumulate monotonically over time within the genus. However, different genes contribute to expression variance within species groups compared to between groups. We observed greater turnover of species-restricted genes with male-biased expression, indicating that gene formation and extinction may play a significant part in species differences. Genes with male-biased expression also show the greatest expression and DNA sequence divergence. This higher divergence and turnover of genes with male-biased expression may be due to high transcription rates in the male germline, greater functional pleiotropy of genes expressed in females, and/or sexual competition.

Diversity and dynamics of the Drosophila transcriptome

Animal transcriptomes are dynamic, with each cell type, tissue and organ system expressing an ensemble of transcript isoforms that give rise to substantial diversity. Here we have identified new genes, transcripts and proteins using poly(A)+ RNA sequencing from Drosophila melanogaster in cultured cell lines, dissected organ systems and under environmental perturbations. We found that a small set of mostly neural-specific genes has the potential to encode thousands of transcripts each through extensive alternative promoter usage and RNA splicing. The magnitudes of splicing changes are larger between tissues than between developmental stages, and most sex-specific splicing is gonad-specific. Gonads express hundreds of previously unknown coding and long non-coding RNAs (lncRNAs), some of which are antisense to protein-coding genes and produce short regulatory RNAs. Furthermore, previously identified pervasive intergenic transcription occurs primarily within newly identified introns. The fly transcriptome is substantially more complex than previously recognized, with this complexity arising from combinatorial usage of promoters, splice sites and polyadenylation sites.

Evolution of protein-coding genes in Drosophila

Several contributing factors have been implicated in evolutionary rate heterogeneity among proteins, but their evolutionary mechanisms remain poorly characterized. The recently sequenced 12 Drosophila genomes provide a unique opportunity to shed light on these unresolved issues. Here, we focus on the role of natural selection in shaping evolutionary rates. We use the Drosophila genomic data to distinguish between factors that increase the strength of purifying selection on proteins and factors that affect the amount of positive selection experienced by proteins. We confirm the importance of translational selection in shaping protein evolution in Drosophila and show that factors such as tissue bias in expression, gene essentiality, intron number, and recombination rate also contribute to evolutionary rate variation among proteins.

Demasculinization of X chromosomes in the Drosophila genus

X chromosomes evolve differently from autosomes, but general governing principles have not emerged. For example, genes with male-biased expression are under-represented on the X chromosome of D. melanogaster, but are randomly distributed in the genome of Anopheles gambiae. In direct global profiling experiments using species-specific microarrays, we find a nearly identical paucity of genes with male-biased expression on D. melanogaster, D. simulans, D. yakuba, D. ananassae, D. virilis and D. mojavensis X chromosomes. We observe the same under-representation on the neo-X of D. pseudoobscura. It has been suggested that precocious meiotic silencing of the X chromosome accounts for reduced X chromosome male-biased expression in nematodes, mammals and Drosophila. We show that X chromosome genes with male-biased expression are under-represented in somatic cells and in mitotic male germ cells. These data are incompatible with simple X chromosome inactivation models. Using expression profiling and comparative sequence analysis, we show that selective gene extinction on the X chromosome, creation of new genes on autosomes and changed genomic location of existing genes contribute to the unusual X chromosome gene content.

Comparative genomics of Drosophila and human core promoters

BackgroundThe core promoter region plays a critical role in the regulation of eukaryotic gene expression. We have determined the non-random distribution of DNA sequences relative to the transcriptional start site in Drosophila melanogaster promoters to identify sequences that may be biologically significant. We compare these results with those obtained for human promoters.ResultsWe determined the distribution of all 65,536 octamer (8-mers) DNA sequences in 10,914 Drosophila promoters and two sets of human promoters aligned relative to the transcriptional start site. In Drosophila, 298 8-mers have highly significant (p ≤ 1 × 10-16) non-random distributions peaking within 100 base-pairs of the transcriptional start site. These sequences were grouped into 15 DNA motifs. Ten motifs, termed directional motifs, occur only on the positive strand while the remaining five motifs, termed non-directional motifs, occur on both strands. The only directional motifs to localize in human promoters are TATA, INR, and DPE. The directional motifs were further subdivided into those precisely positioned relative to the transcriptional start site and those that are positioned more loosely relative to the transcriptional start site. Similar numbers of non-directional motifs were identified in both species and most are different. The genes associated with all 15 DNA motifs, when they occur in the peak, are enriched in specific Gene Ontology categories and show a distinct mRNA expression pattern, suggesting that there is a core promoter code in Drosophila.ConclusionDrosophila and human promoters use different DNA sequences to regulate gene expression, supporting the idea that evolution occurs by the modulation of gene regulation.

Comparative validation of the D. melanogaster modENCODE transcriptome annotation

Accurate gene model annotation of reference genomes is critical for making them useful. The modENCODE project has improved the D. melanogaster genome annotation by using deep and diverse high-throughput data. Since transcriptional activity that has been evolutionarily conserved is likely to have an advantageous function, we have performed large-scale interspecific comparisons to increase confidence in predicted annotations. To support comparative genomics, we filled in divergence gaps in the Drosophila phylogeny by generating draft genomes for eight new species. For comparative transcriptome analysis, we generated mRNA expression profiles on 81 samples from multiple tissues and developmental stages of 15 Drosophila species, and we performed cap analysis of gene expression in D. melanogaster and D. pseudoobscura. We also describe conservation of four distinct core promoter structures composed of combinations of elements at three positions. Overall, each type of genomic feature shows a characteristic divergence rate relative to neutral models, highlighting the value of multispecies alignment in annotating a target genome that should prove useful in the annotation of other high priority genomes, especially human and other mammalian genomes that are rich in noncoding sequences. We report that the vast majority of elements in the annotation are evolutionarily conserved, indicating that the annotation will be an important springboard for functional genetic testing by the Drosophila community.

Cajal bodies are linked to genome conformation

The mechanisms underlying nuclear body (NB) formation and their contribution to genome function are unknown. Here we examined the non-random positioning of Cajal bodies (CBs), major NBs involved in spliceosomal snRNP assembly and their role in genome organization. CBs are predominantly located at the periphery of chromosome territories at a multi-chromosome interface. Genome-wide chromosome conformation capture analysis (4C-seq) using CB-interacting loci revealed that CB-associated regions are enriched with highly expressed histone genes and U small nuclear or nucleolar RNA (sn/snoRNA) loci that form intra- and inter-chromosomal clusters. In particular, we observed a number of CB-dependent gene-positioning events on chromosome 1. RNAi-mediated disassembly of CBs disrupts the CB-targeting gene clusters and suppresses the expression of U sn/snoRNA and histone genes. This loss of spliceosomal snRNP production results in increased splicing noise, even in CB-distal regions. Therefore, we conclude that CBs contribute to genome organization with global effects on gene expression and RNA splicing fidelity.

TET-catalyzed oxidation of intragenic 5-methylcytosine regulates CTCF-dependent alternative splicing

Intragenic 5‐methylcytosine and CTCF mediate opposing effects on pre‐mRNA splicing: CTCF promotes inclusion of weak upstream exons through RNA polymerase II pausing, whereas 5‐methylcytosine evicts CTCF, leading to exon exclusion. However, the mechanisms governing dynamic DNA methylation at CTCF‐binding sites were unclear. Here, we reveal the methylcytosine dioxygenases TET1 and TET2 as active regulators of CTCF‐mediated alternative splicing through conversion of 5‐methylcytosine to its oxidation derivatives. 5‐hydroxymethylcytosine and 5‐carboxylcytosine are enriched at an intragenic CTCF‐binding sites in the CD45 model gene and are associated with alternative exon inclusion. Reduced TET levels culminate in increased 5‐methylcytosine, resulting in CTCF eviction and exon exclusion. In vitro analyses establish the oxidation derivatives are not sufficient to stimulate splicing, but efficiently promote CTCF association. We further show genomewide that reciprocal exchange of 5‐hydroxymethylcytosine and 5‐methylcytosine at downstream CTCF‐binding sites is a general feature of alternative splicing in naïve and activated CD4+ T cells. These findings significantly expand our current concept of the pre‐mRNA “splicing code” to include dynamic intragenic DNA methylation catalyzed by the TET proteins.