Brenton R. Graveley

Alternative splicing: increasing diversity in the proteomic world.

How can the genome of Drosophila melanogaster contain fewer genes than the undoubtedly simpler organism Caenorhabditis elegans? The answer must lie within their proteomes. It is becoming clear that alternative splicing has an extremely important role in expanding protein diversity and might therefore partially underlie the apparent discrepancy between gene number and organismal complexity. Alternative splicing can generate more transcripts from a single gene than the number of genes in an entire genome. However, for the vast majority of alternative splicing events, the functional significance is unknown. Developing a full catalog of alternatively spliced transcripts and determining each of their functions will be a major challenge of the upcoming proteomic era.

Sorting out the complexity of SR protein functions.

been clear whether all of these activities occur during the removal of each intron+ Recent studies now suggest that all of the proposed SR protein functions are carried out during each round of splicing, and at least some of these functions are performed by independent SR protein molecules+ This review discusses recent advances in understanding the diverse functions of SR proteins in metazoan pre-mRNA splicing and presents a model that takes these new findings into account+ Although the reader should keep in mind that the activity of SR proteins in vivo can be influenced by modulating their subcellular localization, this review focuses on biochemical experiments that specifically address the mechanisms by which SR proteins directly function in the splicing reaction+ (The reader may also refer to related reviews recently published elsewhere (Tacke & Manley, 1999; Blencowe, 2000)+)

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.

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.

Genome-wide Analysis of Drosophila Circular RNAs Reveals Their Structural and Sequence Properties and Age-Dependent Neural Accumulation

Circularization was recently recognized to broadly expand transcriptome complexity. Here, we exploit massive Drosophila total RNA-sequencing data, >5 billion paired-end reads from >100 libraries covering diverse developmental stages, tissues, and cultured cells, to rigorously annotate >2,500 fruit fly circular RNAs. These mostly derive from back-splicing of protein-coding genes and lack poly(A) tails, and the circularization of hundreds of genes is conserved across multiple Drosophila species. We elucidate structural and sequence properties of Drosophila circular RNAs, which exhibit commonalities and distinctions from mammalian circles. Notably, Drosophila circular RNAs harbor >1,000 well-conserved canonical miRNA seed matches, especially within coding regions, and coding conserved miRNA sites reside preferentially within circularized exons. Finally, we analyze the developmental and tissue specificity of circular RNAs and note their preferred derivation from neural genes and enhanced accumulation in neural tissues. Interestingly, circular isoforms increase substantially relative to linear isoforms during CNS aging and constitute an aging biomarker.

Dynamic Integration of Splicing within Gene Regulatory Pathways

Precursor mRNA splicing is one of the most highly regulated processes in metazoan species. In addition to generating vast repertoires of RNAs and proteins, splicing has a profound impact on other gene regulatory layers, including mRNA transcription, turnover, transport, and translation. Conversely, factors regulating chromatin and transcription complexes impact the splicing process. This extensive crosstalk between gene regulatory layers takes advantage of dynamic spatial, physical, and temporal organizational properties of the cell nucleus, and further emphasizes the importance of developing a multidimensional understanding of splicing control.

Regulatory divergence in Drosophila revealed by mRNA-seq

The regulation of gene expression is critical for organismal function and is an important source of phenotypic diversity between species. Understanding the genetic and molecular mechanisms responsible for regulatory divergence is therefore expected to provide insight into evolutionary change. Using deep sequencing, we quantified total and allele-specific mRNA expression levels genome-wide in two closely related Drosophila species (D. melanogaster and D. sechellia) and their F(1) hybrids. We show that 78% of expressed genes have divergent expression between species, and that cis- and trans-regulatory divergence affects 51% and 66% of expressed genes, respectively, with 35% of genes showing evidence of both. This is a relatively larger contribution of trans-regulatory divergence than was expected based on prior studies, and may result from the unique demographic history of D. sechellia. Genes with antagonistic cis- and trans-regulatory changes were more likely to be misexpressed in hybrids, consistent with the idea that such regulatory changes contribute to hybrid incompatibilities. In addition, cis-regulatory differences contributed more to divergent expression of genes that showed additive rather than nonadditive inheritance. A correlation between sequence similarity and the conservation of cis-regulatory activity was also observed that appears to be a general feature of regulatory evolution. Finally, we examined regulatory divergence that may have contributed to the evolution of a specific trait--divergent feeding behavior in D. sechellia. Overall, this study illustrates the power of mRNA sequencing for investigating regulatory evolution, provides novel insight into the evolution of gene expression in Drosophila, and reveals general trends that are likely to extend to other species.

Codon optimality is a major determinant of mRNA stability

mRNA degradation represents a critical regulated step in gene expression. Although the major pathways in turnover have been identified, accounting for disparate half-lives has been elusive. We show that codon optimality is one feature that contributes greatly to mRNA stability. Genome-wide RNA decay analysis revealed that stable mRNAs are enriched in codons designated optimal, whereas unstable mRNAs contain predominately non-optimal codons. Substitution of optimal codons with synonymous, non-optimal codons results in dramatic mRNA destabilization, whereas the converse substitution significantly increases stability. Further, we demonstrate that codon optimality impacts ribosome translocation, connecting the processes of translation elongation and decay through codon optimality. Finally, we show that optimal codon content accounts for the similar stabilities observed in mRNAs encoding proteins with coordinated physiological function. This work demonstrates that codon optimization exists as a mechanism to finely tune levels of mRNAs and, ultimately, proteins.

Mutually Exclusive Splicing of the Insect Dscam Pre-mRNA Directed by Competing Intronic RNA Secondary Structures

Drosophila Dscam encodes 38,016 distinct axon guidance receptors through the mutually exclusive alternative splicing of 95 variable exons. Importantly, known mechanisms that ensure the mutually exclusive splicing of pairs of exons cannot explain this phenomenon in Dscam. I have identified two classes of conserved elements in the Dscam exon 6 cluster, which contains 48 alternative exons--the docking site, located in the intron downstream of constitutive exon 5, and the selector sequences, which are located upstream of each exon 6 variant. Strikingly, each selector sequence is complementary to a portion of the docking site, and this pairing juxtaposes one, and only one, alternative exon to the upstream constitutive exon. The mutually exclusive nature of the docking site:selector sequence interactions suggests that the formation of these competing RNA structures is a central component of the mechanism guaranteeing that only one exon 6 variant is included in each Dscam mRNA.

Arginine/Serine-Rich Domains of SR Proteins Can Function as Activators of Pre-mRNA Splicing

Serine/arginine (SR)-rich splicing factors contain an RNA binding domain and an arginine/serine (RS)-rich domain required for protein-protein interactions. In addition to their roles in the basic splicing reaction, SR proteins function as components of splicing enhancer complexes. Here, we investigate the role of RS domains in splicing enhancer function. Hybrid proteins containing RS domains fused to the MS2 RNA binding protein were tested in vitro with RNA substrates bearing an MS2 recognition sequence. These hybrid proteins activated splicing in nuclear extracts, but not in S100 extracts lacking SR proteins. However, intact recombinant SR proteins could complement the activity of the hybrid proteins in S100 extracts. These data demonstrate that RS domains function as splicing activators and suggest that the general and enhancer-dependent functions of SR proteins can be uncoupled.