The kinetics of pre-mRNA splicing in the Drosophila genome and the influence of gene architecture

  1. Athma A Pai
  2. Telmo Henriques
  3. Kayla McCue
  4. Adam Burkholder
  5. Karen Adelman
  6. Christopher B Burge  Is a corresponding author
  1. Massachusetts Institute of Technology, United States
  2. National Institute for Environmental Health Sciences, United States
  3. Harvard Medical School, United States

Abstract

Production of most eukaryotic mRNAs requires splicing of introns from pre-mRNA. The splicing reaction requires definition of splice sites, which are initially recognized in either intron-spanning ('intron definition') or exon-spanning ('exon definition') pairs. To understand how exon and intron length and splice site recognition mode impact splicing, we measured splicing rates genome-wide in Drosophila, using metabolic labeling/RNA sequencing and new mathematical models to estimate rates. We found that the modal intron length range of 60-70 nt represents a local maximum of splicing rates, but that much longer exon-defined introns are spliced even faster and more accurately. Surprisingly, we observed low variation in splicing rates across introns in the same gene, suggesting the presence of gene-level influences, and we identified multiple gene level variables associated with splicing rate. Together our data suggest that developmental and stress response genes may have preferentially evolved exon definition in order to enhance rates of splicing.

Data availability

The following data sets were generated

Article and author information

Author details

  1. Athma A Pai

    Department of Biology, Massachusetts Institute of Technology, Cambridge, United States
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-7995-9948
  2. Telmo Henriques

    Epigenetics and Stem Cell Biology Laboratory, National Institute for Environmental Health Sciences, Durham, United States
    Competing interests
    No competing interests declared.
  3. Kayla McCue

    Program in Computational and Systems Biology, Massachusetts Institute of Technology, Cambridge, United States
    Competing interests
    No competing interests declared.
  4. Adam Burkholder

    Center for Integrative Bioinformatics, National Institute for Environmental Health Sciences, Durham, United States
    Competing interests
    No competing interests declared.
  5. Karen Adelman

    Department of Biological Chemistry and Molecular Pharmacology, Harvard Medical School, Boston, United States
    Competing interests
    Karen Adelman, Reviewing editor, eLife.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-5364-334X
  6. Christopher B Burge

    Department of Biology, Massachusetts Institute of Technology, Cambridge, United States
    For correspondence
    cburge@mit.edu
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-9047-5648

Funding

National Institutes of Health (Z01-ES101987)

  • Telmo Henriques
  • Adam Burkholder
  • Karen Adelman

National Institutes of Health (R01-GM085319)

  • Athma A Pai
  • Christopher B Burge

Jane Coffin Childs Memorial Fund for Medical Research

  • Athma A Pai

U.S. Department of Energy (FG02-97ER25308)

  • Kayla McCue

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Reviewing Editor

  1. Timothy W Nilsen, Case Western Reserve University, United States

Publication history

  1. Received: October 5, 2017
  2. Accepted: December 22, 2017
  3. Accepted Manuscript published: December 27, 2017 (version 1)
  4. Version of Record published: January 10, 2018 (version 2)

Copyright

This is an open-access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.

Metrics

  • 5,621
    Page views
  • 783
    Downloads
  • 37
    Citations

Article citation count generated by polling the highest count across the following sources: Crossref, Scopus, PubMed Central.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Athma A Pai
  2. Telmo Henriques
  3. Kayla McCue
  4. Adam Burkholder
  5. Karen Adelman
  6. Christopher B Burge
(2017)
The kinetics of pre-mRNA splicing in the Drosophila genome and the influence of gene architecture
eLife 6:e32537.
https://doi.org/10.7554/eLife.32537

Further reading

    1. Computational and Systems Biology
    David Elkind, Hannah Hochgerner ... Amit Zeisel
    Research Article Updated

    The mouse brain is by far the most intensively studied among mammalian brains, yet basic measures of its cytoarchitecture remain obscure. For example, quantifying cell numbers, and the interplay of sex, strain, and individual variability in cell density and volume is out of reach for many regions. The Allen Mouse Brain Connectivity project produces high-resolution full brain images of hundreds of brains. Although these were created for a different purpose, they reveal details of neuroanatomy and cytoarchitecture. Here, we used this population to systematically characterize cell density and volume for each anatomical unit in the mouse brain. We developed a DNN-based segmentation pipeline that uses the autofluorescence intensities of images to segment cell nuclei even within the densest regions, such as the dentate gyrus. We applied our pipeline to 507 brains of males and females from C57BL/6J and FVB.CD1 strains. Globally, we found that increased overall brain volume does not result in uniform expansion across all regions. Moreover, region-specific density changes are often negatively correlated with the volume of the region; therefore, cell count does not scale linearly with volume. Many regions, including layer 2/3 across several cortical areas, showed distinct lateral bias. We identified strain-specific or sex-specific differences. For example, males tended to have more cells in extended amygdala and hypothalamic regions (MEA, BST, BLA, BMA, and LPO, AHN) while females had more cells in the orbital cortex (ORB). Yet, inter-individual variability was always greater than the effect size of a single qualifier. We provide the results of this analysis as an accessible resource for the community.

    1. Computational and Systems Biology
    2. Immunology and Inflammation
    Magdalena L Russell, Noah Simon ... Frederick A Matsen IV
    Research Article

    To appropriately defend against a wide array of pathogens, humans somatically generate highly diverse repertoires of B cell and T cell receptors (BCRs and TCRs) through a random process called V(D)J recombination. Receptor diversity is achieved during this process through both the combinatorial assembly of V(D)J-genes and the junctional deletion and insertion of nucleotides. While the Artemis protein is often regarded as the main nuclease involved in V(D)J recombination, the exact mechanism of nucleotide trimming is not understood. Using a previously published TCRβ repertoire sequencing data set, we have designed a flexible probabilistic model of nucleotide trimming that allows us to explore various mechanistically interpretable sequence-level features. We show that local sequence context, length, and GC nucleotide content in both directions of the wider sequence, together, can most accurately predict the trimming probabilities of a given V-gene sequence. Because GC nucleotide content is predictive of sequence-breathing, this model provides quantitative statistical evidence regarding the extent to which double-stranded DNA may need to be able to breathe for trimming to occur. We also see evidence of a sequence motif that appears to get preferentially trimmed, independent of GC-content-related effects. Further, we find that the inferred coefficients from this model provide accurate prediction for V- and J-gene sequences from other adaptive immune receptor loci. These results refine our understanding of how the Artemis nuclease may function to trim nucleotides during V(D)J recombination and provide another step toward understanding how V(D)J recombination generates diverse receptors and supports a powerful, unique immune response in healthy humans.