C. Joel McManus

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.

RNA structure and the mechanisms of alternative splicing

Alternative splicing is a widespread means of increasing protein diversity and regulating gene expression in eukaryotes. Much progress has been made in understanding the proteins involved in regulating alternative splicing, the sequences they bind to, and how these interactions lead to changes in splicing patterns. However, several recent studies have identified other players involved in regulating alternative splicing. A major theme emerging from these studies is that RNA secondary structures play an under appreciated role in the regulation of alternative splicing. This review provides an overview of the basic aspects of splicing regulation and highlights recent progress in understanding the role of RNA secondary structure in this process.

Mod-seq: high-throughput sequencing for chemical probing of RNA structure

The functions of RNA molecules are intimately linked to their ability to fold into complex secondary and tertiary structures. Thus, understanding how these molecules fold is essential to determining how they function. Current methods for investigating RNA structure often use small molecules, enzymes, or ions that cleave or modify the RNA in a solvent-accessible manner. While these methods have been invaluable to understanding RNA structure, they can be fairly labor intensive and often focus on short regions of single RNAs. Here we present a new method (Mod-seq) and data analysis pipeline (Mod-seeker) for assaying the structure of RNAs by high-throughput sequencing. This technique can be utilized both in vivo and in vitro, with any small molecule that modifies RNA and consequently impedes reverse transcriptase. As proof-of-principle, we used dimethyl sulfate (DMS) to probe the in vivo structure of total cellular RNAs in Saccharomyces cerevisiae. Mod-seq analysis simultaneously revealed secondary structural information for all four ribosomal RNAs and 32 additional noncoding RNAs. We further show that Mod-seq can be used to detect structural changes in 5.8S and 25S rRNAs in the absence of ribosomal protein L26, correctly identifying its binding site on the ribosome. While this method is applicable to RNAs of any length, its high-throughput nature makes Mod-seq ideal for studying long RNAs and complex RNA mixtures.

Ribosome profiling reveals post-transcriptional buffering of divergent gene expression in yeast

Understanding the patterns and causes of phenotypic divergence is a central goal in evolutionary biology. Much work has shown that mRNA abundance is highly variable between closely related species. However, the extent and mechanisms of post-transcriptional gene regulatory evolution are largely unknown. Here we used ribosome profiling to compare transcript abundance and translation efficiency in two closely related yeast species (S. cerevisiae and S. paradoxus). By comparing translation regulatory divergence to interspecies differences in mRNA sequence features, we show that differences in transcript leaders and codon bias substantially contribute to divergent translation. Globally, we find that translation regulatory divergence often buffers species differences in mRNA abundance, such that ribosome occupancy is more conserved than transcript abundance. We used allele-specific ribosome profiling in interspecies hybrids to compare the relative contributions of cis- and trans-regulatory divergence to species differences in mRNA abundance and translation efficiency. The mode of gene regulatory divergence differs for these processes, as trans-regulatory changes play a greater role in divergent mRNA abundance than in divergent translation efficiency. Strikingly, most genes with aberrant transcript abundance in F1 hybrids (either over- or underexpressed compared to both parent species) did not exhibit aberrant ribosome occupancy. Our results show that interspecies differences in translation contribute substantially to the evolution of gene expression. Compensatory differences in transcript abundance and translation efficiency may increase the robustness of gene regulation.

Global analysis of trans-splicing in Drosophila

Precursor mRNA (pre-mRNA) splicing can join exons contained on either a single pre-mRNA (cis) or on separate pre-mRNAs (trans). It is exceedingly rare to have trans-splicing between protein-coding exons and has been demonstrated for only two Drosophila genes: mod(mdg4) and lola. It has also been suggested that trans-splicing is a mechanism for the generation of chimeric RNA products containing sequence from multiple distant genomic sites. Because most high-throughput approaches cannot distinguish cis- and trans-splicing events, the extent to which trans-splicing occurs between protein-coding exons in any organism is unknown. Here, we used paired-end deep sequencing of mRNA to identify genes that undergo trans-splicing in Drosophila interspecies hybrids. We did not observe credible evidence for the existence of chimeric RNAs generated by trans-splicing of RNAs transcribed from distant genomic loci. Rather, our data suggest that experimental artifacts are the source of most, if not all, apparent chimeric RNA products. We did, however, identify 80 genes that appear to undergo trans-splicing between homologous alleles and can be classified into three categories based on their organization: (i) genes with multiple 3′ terminal exons, (ii) genes with multiple first exons, and (iii) genes with very large introns, often containing other genes. Our results suggest that trans-splicing between homologous alleles occurs more commonly in Drosophila than previously believed and may facilitate expression of architecturally complex genes.

Tempo and mode of regulatory evolution in Drosophila

Genetic changes affecting gene expression contribute to phenotypic divergence; thus, understanding how regulatory networks controlling gene expression change over time is critical for understanding evolution. Prior studies of expression differences within and between species have identified properties of regulatory divergence, but technical and biological differences among these studies make it difficult to assess the generality of these properties or to understand how regulatory changes accumulate with divergence time. Here, we address these issues by comparing gene expression among strains and species of Drosophila with a range of divergence times and use F1 hybrids to examine inheritance patterns and disentangle cis- and trans-regulatory changes. We find that the fixation of compensatory changes has caused the regulation of gene expression to diverge more rapidly than gene expression itself. Specifically, we observed that the proportion of genes with evidence of cis-regulatory divergence has increased more rapidly with divergence time than the proportion of genes with evidence of expression differences. Surprisingly, the amount of expression divergence explained by cis-regulatory changes did not increase steadily with divergence time, as was previously proposed. Rather, one species (Drosophila sechellia) showed an excess of cis-regulatory divergence that we argue most likely resulted from positive selection in this lineage. Taken together, this work reveals not only the rate at which gene expression evolves, but also the molecular and evolutionary mechanisms responsible for this evolution.

DNA copy number evolution in Drosophila cell lines

BackgroundStructural rearrangements of the genome resulting in genic imbalance due to copy number change are often deleterious at the organismal level, but are common in immortalized cell lines and tumors, where they may be an advantage to cells. In order to explore the biological consequences of copy number changes in the Drosophila genome, we resequenced the genomes of 19 tissue-culture cell lines and generated RNA-Seq profiles.ResultsOur work revealed dramatic duplications and deletions in all cell lines. We found three lines of evidence indicating that copy number changes were due to selection during tissue culture. First, we found that copy numbers correlated to maintain stoichiometric balance in protein complexes and biochemical pathways, consistent with the gene balance hypothesis. Second, while most copy number changes were cell line-specific, we identified some copy number changes shared by many of the independent cell lines. These included dramatic recurrence of increased copy number of the PDGF/VEGF receptor, which is also over-expressed in many cancer cells, and of bantam, an anti-apoptosis miRNA. Third, even when copy number changes seemed distinct between lines, there was strong evidence that they supported a common phenotypic outcome. For example, we found that proto-oncogenes were over-represented in one cell line (S2-DRSC), whereas tumor suppressor genes were under-represented in another (Kc167).ConclusionOur study illustrates how genome structure changes may contribute to selection of cell lines in vitro. This has implications for other cell-level natural selection progressions, including tumorigenesis.

Competing RNA secondary structures are required for mutually exclusive splicing of the Dscam exon 6 cluster

Alternative splicing of eukaryotic pre-mRNAs is an important mechanism for generating proteome diversity and regulating gene expression. The Drosophila melanogaster Down Syndrome Cell Adhesion Molecule (Dscam) gene is an extreme example of mutually exclusive splicing. Dscam contains 95 alternatively spliced exons that potentially encode 38,016 distinct mRNA and protein isoforms. We previously identified two sets of conserved sequence elements, the docking site and selector sequences in the Dscam exon 6 cluster, which contains 48 mutually exclusive exons. These elements were proposed to engage in competing RNA secondary structures required for mutually exclusive splicing, though this model has not yet been experimentally tested. Here we describe a new system that allowed us to demonstrate that the docking site and selector sequences are indeed required for exon 6 mutually exclusive splicing and that the strength of these RNA structures determines the frequency of exon 6 inclusion. We also show that the function of the docking site has been conserved for ~500 million years of evolution. This work demonstrates that conserved intronic sequences play a functional role in mutually exclusive splicing of the Dscam exon 6 cluster.

Genome-wide analysis of the skeletogenic gene regulatory network of sea urchins

A central challenge of developmental and evolutionary biology is to understand the transformation of genetic information into morphology. Elucidating the connections between genes and anatomy will require model morphogenetic processes that are amenable to detailed analysis of cell/tissue behaviors and to systems-level approaches to gene regulation. The formation of the calcified endoskeleton of the sea urchin embryo is a valuable experimental system for developing such an integrated view of the genomic regulatory control of morphogenesis. A transcriptional gene regulatory network (GRN) that underlies the specification of skeletogenic cells (primary mesenchyme cells, or PMCs) has recently been elucidated. In this study, we carried out a genome-wide analysis of mRNAs encoded by effector genes in the network and uncovered transcriptional inputs into many of these genes. We used RNA-seq to identify >400 transcripts differentially expressed by PMCs during gastrulation, when these cells undergo a striking sequence of behaviors that drives skeletal morphogenesis. Our analysis expanded by almost an order of magnitude the number of known (and candidate) downstream effectors that directly mediate skeletal morphogenesis. We carried out genome-wide analysis of (1) functional targets of Ets1 and Alx1, two pivotal, early transcription factors in the PMC GRN, and (2) functional targets of MAPK signaling, a pathway that plays an essential role in PMC specification. These studies identified transcriptional inputs into >200 PMC effector genes. Our work establishes a framework for understanding the genomic regulatory control of a major morphogenetic process and has important implications for reconstructing the evolution of biomineralization in metazoans.

Genomic imprinting absent in Drosophila melanogaster adult females.

Genomic imprinting occurs when expression of an allele differs based on the sex of the parent that transmitted the allele. In D. melanogaster, imprinting can occur, but its impact on allelic expression genome-wide is unclear. Here, we search for imprinted genes in D. melanogaster using RNA-seq to compare allele-specific expression between pools of 7- to 10-day-old adult female progeny from reciprocal crosses. We identified 119 genes with allelic expression consistent with imprinting, and these genes showed significant clustering within the genome. Surprisingly, additional analysis of several of these genes showed that either genomic heterogeneity or high levels of intrinsic noise caused imprinting-like allelic expression. Consequently, our data provide no convincing evidence of imprinting for D. melanogaster genes in their native genomic context. Elucidating sources of false-positive signals for imprinting in allele-specific RNA-seq data, as done here, is critical given the growing popularity of this method for identifying imprinted genes.