Rate-limiting steps in yeast protein translation.

Deep sequencing now provides detailed snapshots of ribosome occupancy on mRNAs. We leverage these data to parameterize a computational model of translation, keeping track of every ribosome, tRNA, and mRNA molecule in a yeast cell. We determine the parameter regimes in which fast initiation or high codon bias in a transgene increases protein yield and infer the initiation rates of endogenous Saccharomyces cerevisiae genes, which vary by several orders of magnitude and correlate with 5' mRNA folding energies. Our model recapitulates the previously reported 5'-to-3' ramp of decreasing ribosome densities, although our analysis shows that this ramp is caused by rapid initiation of short genes rather than slow codons at the start of transcripts. We conclude that protein production in healthy yeast cells is typically limited by the availability of free ribosomes, whereas protein production under periods of stress can sometimes be rescued by reducing initiation or elongation rates.

XACT-Seq Comprehensively Defines the Promoter-Position and Promoter-Sequence Determinants for Initial-Transcription Pausing.

Pausing by RNA polymerase (RNAP) during transcription elongation, in which a translocating RNAP uses a "stepping" mechanism, has been studied extensively, but pausing by RNAP during initial transcription, in which a promoter-anchored RNAP uses a "scrunching" mechanism, has not. We report a method that directly defines the RNAP-active-center position relative to DNA with single-nucleotide resolution (XACT-seq; "crosslink-between-active-center-and-template sequencing"). We apply this method to detect and quantify pausing in initial transcription at 411 (∼4,000,000) promoter sequences in vivo in Escherichia coli. The results show initial-transcription pausing can occur in each nucleotide addition during initial transcription, particularly the first 4 to 5 nucleotide additions. The results further show initial-transcription pausing occurs at sequences that resemble the consensus sequence element for transcription-elongation pausing. Our findings define the positional and sequence determinants for initial-transcription pausing and establish initial-transcription pausing is hard coded by sequence elements similar to those for transcription-elongation pausing.

The LIN28B-IMP1 post-transcriptional regulon has opposing effects on oncogenic signaling in the intestine.

RNA-binding proteins (RBPs) are expressed broadly during both development and malignant transformation, yet their mechanistic roles in epithelial homeostasis or as drivers of tumor initiation and progression are incompletely understood. Here we describe a novel interplay between RBPs LIN28B and IMP1 in intestinal epithelial cells. Ribosome profiling and RNA sequencing identified IMP1 as a principle node for gene expression regulation downstream from LIN28B In vitro and in vivo data demonstrate that epithelial IMP1 loss increases expression of WNT target genes and enhances LIN28B-mediated intestinal tumorigenesis, which was reversed when we overexpressed IMP1 independently in vivo. Furthermore, IMP1 loss in wild-type or LIN28B-overexpressing mice enhances the regenerative response to irradiation. Together, our data provide new evidence for the opposing effects of the LIN28B-IMP1 axis on post-transcriptional regulation of canonical WNT signaling, with implications in intestinal homeostasis, regeneration and tumorigenesis.

Intragenomic variation in non-adaptive nucleotide biases causes underestimation of selection on synonymous codon usage.

Patterns of non-uniform usage of synonymous codons vary across genes in an organism and between species across all domains of life. This codon usage bias (CUB) is due to a combination of non-adaptive (e.g. mutation biases) and adaptive (e.g. natural selection for translation efficiency/accuracy) evolutionary forces. Most models quantify the effects of mutation bias and selection on CUB assuming uniform mutational and other non-adaptive forces across the genome. However, non-adaptive nucleotide biases can vary within a genome due to processes such as biased gene conversion (BGC), potentially obfuscating signals of selection on codon usage. Moreover, genome-wide estimates of non-adaptive nucleotide biases are lacking for non-model organisms. We combine an unsupervised learning method with a population genetics model of synonymous coding sequence evolution to assess the impact of intragenomic variation in non-adaptive nucleotide bias on quantification of natural selection on synonymous codon usage across 49 Saccharomycotina yeasts. We find that in the absence of a priori information, unsupervised learning can be used to identify genes evolving under different non-adaptive nucleotide biases. We find that the impact of intragenomic variation in non-adaptive nucleotide bias varies widely, even among closely-related species. We show that the overall strength and direction of translational selection can be underestimated by failing to account for intragenomic variation in non-adaptive nucleotide biases. Interestingly, genes falling into clusters identified by machine learning are also physically clustered across chromosomes. Our results indicate the need for more nuanced models of sequence evolution that systematically incorporate the effects of variable non-adaptive nucleotide biases on codon frequencies.

Promoter-sequence determinants and structural basis of primer-dependent transcription initiation in Escherichia coli.

Chemical modifications of RNA 5'-ends enable "epitranscriptomic" regulation, influencing multiple aspects of RNA fate. In transcription initiation, a large inventory of substrates compete with nucleoside triphosphates for use as initiating entities, providing an ab initio mechanism for altering the RNA 5'-end. In Escherichia coli cells, RNAs with a 5'-end hydroxyl are generated by use of dinucleotide RNAs as primers for transcription initiation, "primer-dependent initiation." Here, we use massively systematic transcript end readout (MASTER) to detect and quantify RNA 5'-ends generated by primer-dependent initiation for ∼410 (∼1,000,000) promoter sequences in E. coli The results show primer-dependent initiation in E. coli involves any of the 16 possible dinucleotide primers and depends on promoter sequences in, upstream, and downstream of the primer binding site. The results yield a consensus sequence for primer-dependent initiation, YTSS-2NTSS-1NTSSWTSS+1, where TSS is the transcription start site, NTSS-1NTSS is the primer binding site, Y is pyrimidine, and W is A or T. Biochemical and structure-determination studies show that the base pair (nontemplate-strand base:template-strand base) immediately upstream of the primer binding site (Y:RTSS-2, where R is purine) exerts its effect through the base on the DNA template strand (RTSS-2) through interchain base stacking with the RNA primer. Results from analysis of a large set of natural, chromosomally encoded Ecoli promoters support the conclusions from MASTER. Our findings provide a mechanistic and structural description of how TSS-region sequence hard-codes not only the TSS position but also the potential for epitranscriptomic regulation through primer-dependent transcription initiation.

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.

Implementation of Nanopore sequencing as a pragmatic workflow for copy number variant confirmation in the clinic.

BACKGROUND: Diagnosis of rare genetic diseases can be a long, expensive and complex process, involving an array of tests in the hope of obtaining an actionable result. Long-read sequencing platforms offer the opportunity to make definitive molecular diagnoses using a single assay capable of detecting variants, characterizing methylation patterns, resolving complex rearrangements, and assigning findings to long-range haplotypes. Here, we demonstrate the clinical utility of Nanopore long-read sequencing by validating a confirmatory test for copy number variants (CNVs) in neurodevelopmental disorders and illustrate the broader applications of this platform to assess genomic features with significant clinical implications. METHODS: We used adaptive sampling on the Oxford Nanopore platform to sequence 25 genomic DNA samples and 5 blood samples collected from patients with known or false-positive copy number changes originally detected using short-read sequencing. Across the 30 samples (a total of 50 with replicates), we assayed 35 known unique CNVs (a total of 55 with replicates) and one false-positive CNV, ranging in size from 40 kb to 155 Mb, and assessed the presence or absence of suspected CNVs using normalized read depth. RESULTS: Across 50 samples (including replicates) sequenced on individual MinION flow cells, we achieved an average on-target mean depth of 9.5X and an average on-target read length of 4805 bp. Using a custom read depth-based analysis, we successfully confirmed the presence of all 55 known CNVs (including replicates) and the absence of one false-positive CNV. Using the same CNV-targeted data, we compared genotypes of single nucleotide variant loci to verify that no sample mix-ups occurred between assays. For one case, we also used methylation detection and phasing to investigate the parental origin of a 15q11.2-q13 duplication with implications for clinical prognosis. CONCLUSIONS: We present an assay that efficiently targets genomic regions to confirm clinically relevant CNVs with a concordance rate of 100%. Furthermore, we demonstrate how integration of genotype, methylation, and phasing data from the Nanopore sequencing platform can potentially simplify and shorten the diagnostic odyssey.

IMP1/IGF2BP1 in human colorectal cancer extracellular vesicles.

Colorectal cancer (CRC) is a leading cause of cancer-related death. There is an urgent need for new methods of early CRC detection and monitoring to improve patient outcomes. Extracellular vesicles (EVs) are secreted, lipid-bilayer bound, nanoparticles that carry biological cargo throughout the body and in turn exhibit cancer-related biomarker potential. RNA binding proteins (RBPs) are posttranscriptional regulators of gene expression that may provide a link between host cell gene expression and EV phenotypes. Insulin-like growth factor 2 RNA binding protein 1 (IGF2BP1/IMP1) is an RBP that is highly expressed in CRC with higher levels of expression correlating with poor prognosis. IMP1 binds and potently regulates tumor-associated transcripts that may impact CRC EV phenotypes. Our objective was to test whether IMP1 expression levels impact EV secretion and/or cargo. We used RNA sequencing, in vitro CRC cell lines, ex vivo colonoid models, and xenograft mice to test the hypothesis that IMP1 influences EV secretion and/or cargo in human CRC. Our data demonstrate that IMP1 modulates the RNA expression of transcripts associated with extracellular vesicle pathway regulation, but it has no effect on EV secretion levels in vitro or in vivo. Rather, IMP1 appears to affect EV regulation by directly entering EVs in a transformation-dependent manner. These findings suggest that IMP1 has the ability to shape EV cargo in human CRC, which could serve as a diagnostic/prognostic circulating tumor biomarker.NEW & NOTEWORTHY This work demonstrates that the RNA binding protein IGF2BP1/IMP1 alters the transcript profile of colorectal cancer cell (CRC) mRNAs from extracellular vesicle (EV) pathways. IMP1 does not alter EV production or secretion in vitro or in vivo, but rather enters CRC cells where it may further impact EV cargo. Our work shows that IMP1 has the ability to shape EV cargo in human CRC, which could serve as a diagnostic/prognostic circulating tumor biomarker.

Neurodevelopmental disorders and neurodiversity: definition of terms from Scotland's National Autism Implementation Team.

Adults with neurodevelopmental disorders frequently present to, but fit uneasily into, adult mental health services. We offer definitions of important terms related to neurodevelopmental disorders through unifying research data, medical and other viewpoints. This may improve understanding, clinical practice and development of neurodevelopmental disorder pathways within adult mental health services.

Ribosomal stalk proteins RPLP1 and RPLP2 promote biogenesis of flaviviral and cellular multi-pass transmembrane proteins.

The ribosomal stalk proteins, RPLP1 and RPLP2 (RPLP1/2), which form the ancient ribosomal stalk, were discovered decades ago but their functions remain mysterious. We had previously shown that RPLP1/2 are exquisitely required for replication of dengue virus (DENV) and other mosquito-borne flaviviruses. Here, we show that RPLP1/2 function to relieve ribosome pausing within the DENV envelope coding sequence, leading to enhanced protein stability. We evaluated viral and cellular translation in RPLP1/2-depleted cells using ribosome profiling and found that ribosomes pause in the sequence coding for the N-terminus of the envelope protein, immediately downstream of sequences encoding two adjacent transmembrane domains (TMDs). We also find that RPLP1/2 depletion impacts a ribosome density for a small subset of cellular mRNAs. Importantly, the polarity of ribosomes on mRNAs encoding multiple TMDs was disproportionately affected by RPLP1/2 knockdown, implying a role for RPLP1/2 in multi-pass transmembrane protein biogenesis. These analyses of viral and host RNAs converge to implicate RPLP1/2 as functionally important for ribosomes to elongate through ORFs encoding multiple TMDs. We suggest that the effect of RPLP1/2 at TMD associated pauses is mediated by improving the efficiency of co-translational folding and subsequent protein stability.

Posttranscriptional regulation of colonic epithelial repair by RNA binding protein IMP1/IGF2BP1.

RNA binding proteins, including IMP1/IGF2BP1, are essential regulators of intestinal development and cancer. Imp1 hypomorphic mice exhibit gastrointestinal growth defects, yet the specific role for IMP1 in colon epithelial repair is unclear. Our prior work revealed that intestinal epithelial cell-specific Imp1 deletion (Imp1ΔIEC ) was associated with better regeneration in mice after irradiation. Here, we report increased IMP1 expression in patients with Crohn's disease and ulcerative colitis. We demonstrate that Imp1ΔIEC mice exhibit enhanced recovery following dextran sodium sulfate (DSS)-mediated colonic injury. Imp1ΔIEC mice exhibit Paneth cell granule changes, increased autophagy flux, and upregulation of Atg5. In silico and biochemical analyses revealed direct binding of IMP1 to MAP1LC3B, ATG3, and ATG5 transcripts. Genetic deletion of essential autophagy gene Atg7 in Imp1ΔIEC mice revealed increased sensitivity of double-mutant mice to colonic injury compared to control or Atg7 single mutant mice, suggesting a compensatory relationship between Imp1 and the autophagy pathway. The present study defines a novel interplay between IMP1 and autophagy, where IMP1 may be transiently induced during damage to modulate colonic epithelial cell responses to damage.

O-GlcNAcylation alters the selection of mRNAs for translation and promotes 4E-BP1-dependent mitochondrial dysfunction in the retina.

Diabetes promotes the posttranslational modification of proteins by O-linked addition of GlcNAc (O-GlcNAcylation) to Ser/Thr residues of proteins and thereby contributes to diabetic complications. In the retina of diabetic mice, the repressor of mRNA translation, eIF4E-binding protein 1 (4E-BP1), is O-GlcNAcylated, and sequestration of the cap-binding protein eukaryotic translation initiation factor (eIF4E) is enhanced. O-GlcNAcylation has also been detected on several eukaryotic translation initiation factors and ribosomal proteins. However, the functional consequence of this modification is unknown. Here, using ribosome profiling, we evaluated the effect of enhanced O-GlcNAcylation on retinal gene expression. Mice receiving thiamet G (TMG), an inhibitor of the O-GlcNAc hydrolase O-GlcNAcase, exhibited enhanced retinal protein O-GlcNAcylation. The principal effect of TMG on retinal gene expression was observed in ribosome-associated mRNAs (i.e. mRNAs undergoing translation), as less than 1% of mRNAs exhibited changes in abundance. Remarkably, ∼19% of the transcriptome exhibited TMG-induced changes in ribosome occupancy, with 1912 mRNAs having reduced and 1683 mRNAs having increased translational rates. In the retina, the effect of O-GlcNAcase inhibition on translation of specific mitochondrial proteins, including superoxide dismutase 2 (SOD2), depended on 4E-BP1/2. O-GlcNAcylation enhanced cellular respiration and promoted mitochondrial superoxide levels in WT cells, and 4E-BP1/2 deletion prevented O-GlcNAcylation-induced mitochondrial superoxide in cells in culture and in the retina. The retina of diabetic WT mice exhibited increased reactive oxygen species levels, an effect not observed in diabetic 4E-BP1/2-deficient mice. These findings provide evidence for a mechanism whereby diabetes-induced O-GlcNAcylation promotes oxidative stress in the retina by altering the selection of mRNAs for translation.

Non-optimal codon usage is a mechanism to achieve circadian clock conditionality.

Circadian rhythms are oscillations in biological processes that function as a key adaptation to the daily rhythms of most environments. In the model cyanobacterial circadian clock system, the core oscillator proteins are encoded by the gene cluster kaiABC. Genes with high expression and functional importance, such as the kai genes, are usually encoded by optimal codons, yet the codon-usage bias of the kaiBC genes is not optimized for translational efficiency. We discovered a relationship between codon usage and a general property of circadian rhythms called conditionality; namely, that endogenous rhythmicity is robustly expressed under some environmental conditions but not others. Despite the generality of circadian conditionality, however, its molecular basis is unknown for any system. Here we show that in the cyanobacterium Synechococcus elongate, non-optimal codon usage was selected as a post-transcriptional mechanism to switch between circadian and non-circadian regulation of gene expression as an adaptive response to environmental conditions. When the kaiBC sequence was experimentally optimized to enhance expression of the KaiB and KaiC proteins, intrinsic rhythmicity was enhanced at cool temperatures that are experienced by this organism in its natural habitat. However, fitness at those temperatures was highest in cells in which the endogenous rhythms were suppressed at cool temperatures as compared with cells exhibiting high-amplitude rhythmicity. These results indicate natural selection against circadian systems in cyanobacteria that are intrinsically robust at cool temperatures. Modulation of circadian amplitude is therefore crucial to its adaptive significance. Moreover, these results show the direct effects of codon usage on a complex phenotype and organismal fitness. Our work also challenges the long-standing view of directional selection towards optimal codons, and provides a key example of natural selection against optimal codons to achieve adaptive responses to environmental changes.

Contingency and entrenchment in protein evolution under purifying selection.

The phenotypic effect of an allele at one genetic site may depend on alleles at other sites, a phenomenon known as epistasis. Epistasis can profoundly influence the process of evolution in populations and shape the patterns of protein divergence across species. Whereas epistasis between adaptive substitutions has been studied extensively, relatively little is known about epistasis under purifying selection. Here we use computational models of thermodynamic stability in a ligand-binding protein to explore the structure of epistasis in simulations of protein sequence evolution. Even though the predicted effects on stability of random mutations are almost completely additive, the mutations that fix under purifying selection are enriched for epistasis. In particular, the mutations that fix are contingent on previous substitutions: Although nearly neutral at their time of fixation, these mutations would be deleterious in the absence of preceding substitutions. Conversely, substitutions under purifying selection are subsequently entrenched by epistasis with later substitutions: They become increasingly deleterious to revert over time. Our results imply that, even under purifying selection, protein sequence evolution is often contingent on history and so it cannot be predicted by the phenotypic effects of mutations assayed in the ancestral background.

riboviz: analysis and visualization of ribosome profiling datasets.

BACKGROUND: Using high-throughput sequencing to monitor translation in vivo, ribosome profiling can provide critical insights into the dynamics and regulation of protein synthesis in a cell. Since its introduction in 2009, this technique has played a key role in driving biological discovery, and yet it requires a rigorous computational toolkit for widespread adoption. DESCRIPTION: We have developed a database and a browser-based visualization tool, riboviz, that enables exploration and analysis of riboseq datasets. In implementation, riboviz consists of a comprehensive and flexible computational pipeline that allows the user to analyze private, unpublished datasets, along with a web application for comparison with published yeast datasets. Source code and detailed documentation are freely available from https://github.com/shahpr/RiboViz . The web-application is live at www.riboviz.org. CONCLUSIONS: riboviz provides a comprehensive database and analysis and visualization tool to enable comparative analyses of ribosome-profiling datasets. This toolkit will enable both the community of systems biologists who study genome-wide ribosome profiling data and also research groups focused on individual genes to identify patterns of transcriptional and translational regulation across different organisms and conditions.

A codon model of nucleotide substitution with selection on synonymous codon usage.

The quality of phylogenetic inference made from protein-coding genes depends, in part, on the realism with which the codon substitution process is modeled. Here we propose a new mechanistic model that combines the standard M0 substitution model of Yang (1997) with a simplified model from Gilchrist (2007) that includes selection on synonymous substitutions as a function of codon-specific nonsense error rates. We tested the newly proposed model by applying it to 104 protein-coding genes in brewer's yeast, and compared the fit of the new model to the standard M0 model and to the mutation-selection model of Yang and Nielsen (2008) using the AIC. Our new model provided significantly better fit in approximately 85% of the cases considered for the basic M0 model and in approximately 25% of the cases for the M0 model with estimated codon frequencies, but only in a few cases when the mutation-selection model was considered. However, our model includes a parameter that can be interpreted as a measure of the rate of protein production, and the estimates of this parameter were highly correlated with an independent measure of protein production for the yeast genes considered here. Finally, we found that in some cases the new model led to the preference of a different phylogeny for a subset of the genes considered, indicating that substitution model choice may have an impact on the estimated phylogeny.

Explaining complex codon usage patterns with selection for translational efficiency, mutation bias, and genetic drift.

The genetic code is redundant with most amino acids using multiple codons. In many organisms, codon usage is biased toward particular codons. Understanding the adaptive and nonadaptive forces driving the evolution of codon usage bias (CUB) has been an area of intense focus and debate in the fields of molecular and evolutionary biology. However, their relative importance in shaping genomic patterns of CUB remains unsolved. Using a nested model of protein translation and population genetics, we show that observed gene level variation of CUB in Saccharomyces cerevisiae can be explained almost entirely by selection for efficient ribosomal usage, genetic drift, and biased mutation. The correlation between observed codon counts within individual genes and our model predictions is 0.96. Although a variety of factors shape patterns of CUB at the level of individual sites within genes, our results suggest that selection for efficient ribosome usage is a central force in shaping codon usage at the genomic scale. In addition, our model allows direct estimation of codon-specific mutation rates and elongation times and can be readily applied to any organism with high-throughput expression datasets. More generally, we have developed a natural framework for integrating models of molecular processes to population genetics models to quantitatively estimate parameters underlying fundamental biological processes, such a protein translation.

Adults with attention-deficit hyperactivity disorder--diagnosis or normality?

Attention-deficit hyperactivity disorder in adults evokes extreme responses within British psychiatrists, because its diagnostic validity and pharmacological treatments are heavily contested. We propose a model that accommodates apparently divergent evidence, and provides a clinical framework for clinicians and patients, allowing safe, responsible and ethically balanced clinical practice.