Codon optimality modulates protein output by tuning translation initiation.

The impact of synonymous codon choice on protein output has important implications for understanding endogenous gene expression and design of synthetic mRNAs. Previously, we used a neural network model to design a series of synonymous fluorescent reporters whose protein output in yeast spanned a seven-fold range corresponding to their predicted translation speed. Here, we show that this effect is not due primarily to the established impact of slow elongation on mRNA stability, but rather, that an active mechanism further decreases the number of proteins made per mRNA. We combine simulations and careful experiments on fluorescent reporters to argue that translation initiation is limited on non-optimally encoded transcripts. Using a genome-wide CRISPRi screen to discover factors modulating the output from non-optimal transcripts, we identify a set of translation initiation factors including multiple subunits of eIF3 whose depletion restored protein output of a non-optimal reporter. Our results show that codon usage can directly limit protein production, across the full range of endogenous variability in codon usage, by limiting translation initiation.

Streamlined and sensitive mono- and di-ribosome profiling in yeast and human cells.

Ribosome profiling has unveiled diverse regulation and perturbations of translation through a transcriptome-wide survey of ribosome occupancy, read out by sequencing of ribosome-protected messenger RNA fragments. Generation of ribosome footprints and their conversion into sequencing libraries is technically demanding and sensitive to biases that distort the representation of physiological ribosome occupancy. We address these challenges by producing ribosome footprints with P1 nuclease rather than RNase I and replacing RNA ligation with ordered two-template relay, a single-tube protocol for sequencing library preparation that incorporates adaptors by reverse transcription. Our streamlined approach reduced sequence bias and enhanced enrichment of ribosome footprints relative to ribosomal RNA. Furthermore, P1 nuclease preserved distinct juxtaposed ribosome complexes informative about yeast and human ribosome fates during translation initiation, stalling and termination. Our optimized methods for mRNA footprint generation and capture provide a richer translatome profile with low input and fewer technical challenges.

choros: correction of sequence-based biases for accurate quantification of ribosome profiling data.

Ribosome profiling quantifies translation genome-wide by sequencing ribosome-protected fragments, or footprints. Its single-codon resolution allows identification of translation regulation, such as ribosome stalls or pauses, on individual genes. However, enzyme preferences during library preparation lead to pervasive sequence artifacts that obscure translation dynamics. Widespread over- and under-representation of ribosome footprints can dominate local footprint densities and skew estimates of elongation rates by up to five fold. To address these biases and uncover true patterns of translation, we present choros, a computational method that models ribosome footprint distributions to provide bias-corrected footprint counts. choros uses negative binomial regression to accurately estimate two sets of parameters: (i) biological contributions from codon-specific translation elongation rates; and (ii) technical contributions from nuclease digestion and ligation efficiencies. We use these parameter estimates to generate bias correction factors that eliminate sequence artifacts. Applying choros to multiple ribosome profiling datasets, we are able to accurately quantify and attenuate ligation biases to provide more faithful measurements of ribosome distribution. We show that a pattern interpreted as pervasive ribosome pausing near the beginning of coding regions is likely to arise from technical biases. Incorporating choros into standard analysis pipelines will improve biological discovery from measurements of translation.

Rapid detection of SARS-CoV-2 RNA in saliva via Cas13.

Rapid nucleic acid testing is central to infectious disease surveillance. Here, we report an assay for rapid COVID-19 testing and its implementation in a prototype microfluidic device. The assay, which we named DISCoVER (for diagnostics with coronavirus enzymatic reporting), involves extraction-free sample lysis via shelf-stable and low-cost reagents, multiplexed isothermal RNA amplification followed by T7 transcription, and Cas13-mediated cleavage of a quenched fluorophore. The device consists of a single-use gravity-driven microfluidic cartridge inserted into a compact instrument for automated running of the assay and readout of fluorescence within 60 min. DISCoVER can detect severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in saliva with a sensitivity of 40 copies µl-1, and was 94% sensitive and 100% specific when validated (against quantitative PCR) using total RNA extracted from 63 nasal-swab samples (33 SARS-CoV-2-positive, with cycle-threshold values of 13-35). The device correctly identified all tested clinical saliva samples (10 SARS-CoV-2-positive out of 13, with cycle-threshold values of 23-31). Rapid point-of-care nucleic acid testing may broaden the use of molecular diagnostics.

Identifying cell state-associated alternative splicing events and their coregulation.

Alternative splicing shapes the transcriptome and contributes to each cell's unique identity, but single-cell RNA sequencing (scRNA-seq) has struggled to capture the impact of alternative splicing. We previously showed that low recovery of mRNAs from single cells led to erroneous conclusions about the cell-to-cell variability of alternative splicing. Here, we present a method, Psix, to confidently identify splicing that changes across a landscape of single cells, using a probabilistic model that is robust against the data limitations of scRNA-seq. Its autocorrelation-inspired approach finds patterns of alternative splicing that correspond to patterns of cell identity, such as cell type or developmental stage, without the need for explicit cell clustering, labeling, or trajectory inference. Applying Psix to data that follow the trajectory of mouse brain development, we identify exons whose alternative splicing patterns cluster into modules of coregulation. We show that the exons in these modules are enriched for binding by distinct neuronal splicing factors and that their changes in splicing correspond to changes in expression of these splicing factors. Thus, Psix reveals cell type-dependent splicing patterns and the wiring of the splicing regulatory networks that control them. Our new method will enable scRNA-seq analysis to go beyond transcription to understand the roles of post-transcriptional regulation in determining cell identity.

riboviz 2: a flexible and robust ribosome profiling data analysis and visualization workflow.

MOTIVATION: Ribosome profiling, or Ribo-seq, is the state-of-the-art method for quantifying protein synthesis in living cells. Computational analysis of Ribo-seq data remains challenging due to the complexity of the procedure, as well as variations introduced for specific organisms or specialized analyses. RESULTS: We present riboviz 2, an updated riboviz package, for the comprehensive transcript-centric analysis and visualization of Ribo-seq data. riboviz 2 includes an analysis workflow built on the Nextflow workflow management system for end-to-end processing of Ribo-seq data. riboviz 2 has been extensively tested on diverse species and library preparation strategies, including multiplexed samples. riboviz 2 is flexible and uses open, documented file formats, allowing users to integrate new analyses with the pipeline. AVAILABILITY AND IMPLEMENTATION: riboviz 2 is freely available at github.com/riboviz/riboviz.

Accelerated RNA detection using tandem CRISPR nucleases.

Direct, amplification-free detection of RNA has the potential to transform molecular diagnostics by enabling simple on-site analysis of human or environmental samples. CRISPR-Cas nucleases offer programmable RNA-guided RNA recognition that triggers cleavage and release of a fluorescent reporter molecule, but long reaction times hamper their detection sensitivity and speed. Here, we show that unrelated CRISPR nucleases can be deployed in tandem to provide both direct RNA sensing and rapid signal generation, thus enabling robust detection of ~30 molecules per µl of RNA in 20 min. Combining RNA-guided Cas13 and Csm6 with a chemically stabilized activator creates a one-step assay that can detect severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA extracted from respiratory swab samples with quantitative reverse transcriptase PCR (qRT-PCR)-derived cycle threshold (Ct) values up to 33, using a compact detector. This Fast Integrated Nuclease Detection In Tandem (FIND-IT) approach enables sensitive, direct RNA detection in a format that is amenable to point-of-care infection diagnosis as well as to a wide range of other diagnostic or research applications.

Transmission, infectivity, and neutralization of a spike L452R SARS-CoV-2 variant.

We identified an emerging severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variant by viral whole-genome sequencing of 2,172 nasal/nasopharyngeal swab samples from 44 counties in California, a state in the western United States. Named B.1.427/B.1.429 to denote its two lineages, the variant emerged in May 2020 and increased from 0% to >50% of sequenced cases from September 2020 to January 2021, showing 18.6%-24% increased transmissibility relative to wild-type circulating strains. The variant carries three mutations in the spike protein, including an L452R substitution. We found 2-fold increased B.1.427/B.1.429 viral shedding in vivo and increased L452R pseudovirus infection of cell cultures and lung organoids, albeit decreased relative to pseudoviruses carrying the N501Y mutation common to variants B.1.1.7, B.1.351, and P.1. Antibody neutralization assays revealed 4.0- to 6.7-fold and 2.0-fold decreases in neutralizing titers from convalescent patients and vaccine recipients, respectively. The increased prevalence of a more transmissible variant in California exhibiting decreased antibody neutralization warrants further investigation.

Accelerated RNA detection using tandem CRISPR nucleases.

Direct, amplification-free detection of RNA has the potential to transform molecular diagnostics by enabling simple on-site analysis of human or environmental samples. CRISPR-Cas nucleases offer programmable RNA-guided recognition of RNA that triggers cleavage and release of a fluorescent reporter molecule1,2, but long reaction times hamper sensitivity and speed when applied to point-of-care testing. Here we show that unrelated CRISPR nucleases can be deployed in tandem to provide both direct RNA sensing and rapid signal generation, thus enabling robust detection of ~30 RNA copies/microliter in 20 minutes. Combining RNA-guided Cas13 and Csm6 with a chemically stabilized activator creates a one-step assay that detected SARS-CoV-2 RNA from nasopharyngeal samples with PCR-derived Ct values up to 29 in microfluidic chips, using a compact imaging system. This Fast Integrated Nuclease Detection In Tandem (FIND-IT) approach enables direct RNA detection in a format amenable to point-of-care infection diagnosis, as well as to a wide range of other diagnostic or research applications.

Transmission, infectivity, and antibody neutralization of an emerging SARS-CoV-2 variant in California carrying a L452R spike protein mutation.

We identified a novel SARS-CoV-2 variant by viral whole-genome sequencing of 2,172 nasal/nasopharyngeal swab samples from 44 counties in California. Named B.1.427/B.1.429 to denote its 2 lineages, the variant emerged around May 2020 and increased from 0% to >50% of sequenced cases from September 1, 2020 to January 29, 2021, exhibiting an 18.6-24% increase in transmissibility relative to wild-type circulating strains. The variant carries 3 mutations in the spike protein, including an L452R substitution. Our analyses revealed 2-fold increased B.1.427/B.1.429 viral shedding in vivo and increased L452R pseudovirus infection of cell cultures and lung organoids, albeit decreased relative to pseudoviruses carrying the N501Y mutation found in the B.1.1.7, B.1.351, and P.1 variants. Antibody neutralization assays showed 4.0 to 6.7-fold and 2.0-fold decreases in neutralizing titers from convalescent patients and vaccine recipients, respectively. The increased prevalence of a more transmissible variant in California associated with decreased antibody neutralization warrants further investigation.

Rapid, point-of-care molecular diagnostics with Cas13.

Rapid nucleic acid testing is a critical component of a robust infrastructure for increased disease surveillance. Here, we report a microfluidic platform for point-of-care, CRISPR-based molecular diagnostics. We first developed a nucleic acid test which pairs distinct mechanisms of DNA and RNA amplification optimized for high sensitivity and rapid kinetics, linked to Cas13 detection for specificity. We combined this workflow with an extraction-free sample lysis protocol using shelf-stable reagents that are widely available at low cost, and a multiplexed human gene control for calling negative test results. As a proof-of-concept, we demonstrate sensitivity down to 40 copies/µL of SARS-CoV-2 in unextracted saliva within 35 minutes, and validated the test on total RNA extracted from patient nasal swabs with a range of qPCR Ct values from 13-35. To enable sample-to-answer testing, we integrated this diagnostic reaction with a single-use, gravity-driven microfluidic cartridge followed by real-time fluorescent detection in a compact companion instrument. We envision this approach for Diagnostics with Coronavirus Enzymatic Reporting (DISCoVER) will incentivize frequent, fast, and easy testing.

Coverage-dependent bias creates the appearance of binary splicing in single cells.

Single-cell RNA sequencing provides powerful insight into the factors that determine each cell's unique identity. Previous studies led to the surprising observation that alternative splicing among single cells is highly variable and follows a bimodal pattern: a given cell consistently produces either one or the other isoform for a particular splicing choice, with few cells producing both isoforms. Here, we show that this pattern arises almost entirely from technical limitations. We analyze alternative splicing in human and mouse single-cell RNA-seq datasets, and model them with a probabilistic simulator. Our simulations show that low gene expression and low capture efficiency distort the observed distribution of isoforms. This gives the appearance of binary splicing outcomes, even when the underlying reality is consistent with more than one isoform per cell. We show that accounting for the true amount of information recovered can produce biologically meaningful measurements of splicing in single cells.

Accurate design of translational output by a neural network model of ribosome distribution.

Synonymous codon choice can have dramatic effects on ribosome speed and protein expression. Ribosome profiling experiments have underscored that ribosomes do not move uniformly along mRNAs. Here, we have modeled this variation in translation elongation by using a feed-forward neural network to predict the ribosome density at each codon as a function of its sequence neighborhood. Our approach revealed sequence features affecting translation elongation and characterized large technical biases in ribosome profiling. We applied our model to design synonymous variants of a fluorescent protein spanning the range of translation speeds predicted with our model. Levels of the fluorescent protein in budding yeast closely tracked the predicted translation speeds across their full range. We therefore demonstrate that our model captures information determining translation dynamics in vivo; that this information can be harnessed to design coding sequences; and that control of translation elongation alone is sufficient to produce large quantitative differences in protein output.

Improved genome assembly of American alligator genome reveals conserved architecture of estrogen signaling.

The American alligator, Alligator mississippiensis, like all crocodilians, has temperature-dependent sex determination, in which the sex of an embryo is determined by the incubation temperature of the egg during a critical period of development. The lack of genetic differences between male and female alligators leaves open the question of how the genes responsible for sex determination and differentiation are regulated. Insight into this question comes from the fact that exposing an embryo incubated at male-producing temperature to estrogen causes it to develop ovaries. Because estrogen response elements are known to regulate genes over long distances, a contiguous genome assembly is crucial for predicting and understanding their impact. We present an improved assembly of the American alligator genome, scaffolded with in vitro proximity ligation (Chicago) data. We use this assembly to scaffold two other crocodilian genomes based on synteny. We perform RNA sequencing of tissues from American alligator embryos to find genes that are differentially expressed between embryos incubated at male- versus female-producing temperature. Finally, we use the improved contiguity of our assembly along with the current model of CTCF-mediated chromatin looping to predict regions of the genome likely to contain estrogen-responsive genes. We find that these regions are significantly enriched for genes with female-biased expression in developing gonads after the critical period during which sex is determined by incubation temperature. We thus conclude that estrogen signaling is a major driver of female-biased gene expression in the post-temperature sensitive period gonads.

Regulation of splicing factors by alternative splicing and NMD is conserved between kingdoms yet evolutionarily flexible.

Ultraconserved elements, unusually long regions of perfect sequence identity, are found in genes encoding numerous RNA-binding proteins including arginine-serine rich (SR) splicing factors. Expression of these genes is regulated via alternative splicing of the ultraconserved regions to yield mRNAs that are degraded by nonsense-mediated mRNA decay (NMD), a process termed unproductive splicing (Lareau et al. 2007; Ni et al. 2007). As all human SR genes are affected by alternative splicing and NMD, one might expect this regulation to have originated in an early SR gene and persisted as duplications expanded the SR family. But in fact, unproductive splicing of most human SR genes arose independently (Lareau et al. 2007). This paradox led us to investigate the origin and proliferation of unproductive splicing in SR genes. We demonstrate that unproductive splicing of the splicing factor SRSF5 (SRp40) is conserved among all animals and even observed in fungi; this is a rare example of alternative splicing conserved between kingdoms, yet its effect is to trigger mRNA degradation. As the gene duplicated, the ancient unproductive splicing was lost in paralogs, and distinct unproductive splicing evolved rapidly and repeatedly to take its place. SR genes have consistently employed unproductive splicing, and while it is exceptionally conserved in some of these genes, turnover in specific events among paralogs shows flexible means to the same regulatory end.

Distinct stages of the translation elongation cycle revealed by sequencing ribosome-protected mRNA fragments.

During translation elongation, the ribosome ratchets along its mRNA template, incorporating each new amino acid and translocating from one codon to the next. The elongation cycle requires dramatic structural rearrangements of the ribosome. We show here that deep sequencing of ribosome-protected mRNA fragments reveals not only the position of each ribosome but also, unexpectedly, its particular stage of the elongation cycle. Sequencing reveals two distinct populations of ribosome footprints, 28-30 nucleotides and 20-22 nucleotides long, representing translating ribosomes in distinct states, differentially stabilized by specific elongation inhibitors. We find that the balance of small and large footprints varies by codon and is correlated with translation speed. The ability to visualize conformational changes in the ribosome during elongation, at single-codon resolution, provides a new way to study the detailed kinetics of translation and a new probe with which to identify the factors that affect each step in the elongation cycle.DOI: http://dx.doi.org/10.7554/eLife.01257.001.

Ribosome profiling of mouse embryonic stem cells reveals the complexity and dynamics of mammalian proteomes.

The ability to sequence genomes has far outstripped approaches for deciphering the information they encode. Here we present a suite of techniques, based on ribosome profiling (the deep sequencing of ribosome-protected mRNA fragments), to provide genome-wide maps of protein synthesis as well as a pulse-chase strategy for determining rates of translation elongation. We exploit the propensity of harringtonine to cause ribosomes to accumulate at sites of translation initiation together with a machine learning algorithm to define protein products systematically. Analysis of translation in mouse embryonic stem cells reveals thousands of strong pause sites and unannotated translation products. These include amino-terminal extensions and truncations and upstream open reading frames with regulatory potential, initiated at both AUG and non-AUG codons, whose translation changes after differentiation. We also define a class of short, polycistronic ribosome-associated coding RNAs (sprcRNAs) that encode small proteins. Our studies reveal an unanticipated complexity to mammalian proteomes.

Genome-wide identification of alternative splice forms down-regulated by nonsense-mediated mRNA decay in Drosophila.

Alternative mRNA splicing adds a layer of regulation to the expression of thousands of genes in Drosophila melanogaster. Not all alternative splicing results in functional protein; it can also yield mRNA isoforms with premature stop codons that are degraded by the nonsense-mediated mRNA decay (NMD) pathway. This coupling of alternative splicing and NMD provides a mechanism for gene regulation that is highly conserved in mammals. NMD is also active in Drosophila, but its effect on the repertoire of alternative splice forms has been unknown, as has the mechanism by which it recognizes targets. Here, we have employed a custom splicing-sensitive microarray to globally measure the effect of alternative mRNA processing and NMD on Drosophila gene expression. We have developed a new algorithm to infer the expression change of each mRNA isoform of a gene based on the microarray measurements. This method is of general utility for interpreting splicing-sensitive microarrays and high-throughput sequence data. Using this approach, we have identified a high-confidence set of 45 genes where NMD has a differential effect on distinct alternative isoforms, including numerous RNA-binding and ribosomal proteins. Coupled alternative splicing and NMD decrease expression of these genes, which may in turn have a downstream effect on expression of other genes. The NMD-affected genes are enriched for roles in translation and mitosis, perhaps underlying the previously observed role of NMD factors in cell cycle progression. Our results have general implications for understanding the NMD mechanism in fly. Most notably, we found that the NMD-target mRNAs had significantly longer 3' untranslated regions (UTRs) than the nontarget isoforms of the same genes, supporting a role for 3' UTR length in the recognition of NMD targets in fly.

Unproductive splicing of SR genes associated with highly conserved and ultraconserved DNA elements.

The human and mouse genomes share a number of long, perfectly conserved nucleotide sequences, termed ultraconserved elements. Whereas these regions can act as transcriptional enhancers when upstream of genes, those within genes are less well understood. In particular, the function of ultraconserved elements that overlap alternatively spliced exons of genes encoding RNA-binding proteins is unknown. Here we report that in every member of the human SR family of splicing regulators, highly or ultraconserved elements are alternatively spliced, either as alternative 'poison cassette exons' containing early in-frame stop codons, or as alternative introns in the 3' untranslated region. These alternative splicing events target the resulting messenger RNAs for degradation by means of an RNA surveillance pathway called nonsense-mediated mRNA decay. Mouse orthologues of the human SR proteins exhibit the same unproductive splicing patterns. Three SR proteins have been previously shown to direct splicing of their own transcripts, and one of these is known to autoregulate its expression by coupling alternative splicing with decay; our results suggest that unproductive splicing is important for regulation of the entire SR family. We find that unproductive splicing associated with conserved regions has arisen independently in different SR genes, suggesting that splicing factors may readily acquire this form of regulation.