Codon stabilization coefficient as a metric to gain insights into mRNA stability and codon bias and their relationships with translation

Assessment of mRNA Decay and Calculation of Codon Occurrence to mRNA Stability Correlation Coefficients after 5-EU Metabolic Labeling

mRNA translation and decay are tightly connected. This chapter describes a method to assess the influence of each codon identity on mRNA stability in cultured cells. The technique involves metabolic labeling of the nascent mRNAs by addition of the nucleoside analog 5-ethynyluridine (5-EU), purification of the RNA at different time-points after chase of the 5-EU, then biotinylation with Click chemistry, pull-down, and sequencing. The transcripts' half-lives are calculated from the expression level of each mRNA at the different time-points. Finally, the method describes the calculation of the Codon occurrence to mRNA Stability correlation Coefficient, or CSC, as a correlation between the codon occurrence in a transcript and the transcript half-life, for each codon.

Implications of tRNA abundance on translation elongation across bovine tissues

Introduction: Translation is a crucial stage of gene expression. It may also act as an additional layer of regulation that plays an important role in gene expression and function. Highly expressed genes are believed to be codon-biased to support increased protein production, in which quickly translated codons correspond to highly abundant tRNAs. Synonymous SNPs, considered to be silent due to the degeneracy of the genetic code, may shift protein abundance and function through alterations in translational efficiency and suboptimal pairing to lowly abundant tRNAs. Methods: Here, we applied Quantitative Mature tRNA sequencing (QuantM-tRNAseq) and ribosome profiling across bovine tissues in order to investigate the relationship between tRNA expression and slowed translation. Results: Moreover, we have identified genes modulated at transcriptional and/or translational levels underlying tissue-specific biological processes. We have also successfully defined pausing sites that depict the regulatory information encoded within the open reading frame of transcripts, which could be related to translation rate and facilitate proper protein folding. This work offers an atlas of distinctive pausing sites across three bovine tissues, which provides an opportunity to predict codon optimality and understand tissue-specific mechanisms of regulating protein synthesis.

Polysome propensity and tunable thresholds in coding sequence length enable differential mRNA stability

The half-life of mRNAs, as well as their translation, increases in proportion to the optimal codons, indicating a tight coupling of codon-dependent differential translation and degradation. Little is known about the regulation of this coupling. We found that the mRNA stability gain in yeast depends on the mRNA coding sequence length. Below a critical length, codon optimality fails to affect the stability of mRNAs although they can be efficiently translated into short peptides and proteins. Above this threshold length, codon optimality-dependent differential mRNA stability emerges in a switch-like fashion, which coincides with a similar increase in the polysome propensity of the mRNAs. This threshold length can be tuned by the untranslated regions (UTR). Some of these UTRs can destabilize mRNAs without reducing translation, which plays a role in controlling the amplitude of the oscillatory expression of cell cycle genes. Our findings help understand the translation of short peptides from noncoding RNAs and the translation by localized monosomes in neurons.

Differential regulation of mRNA stability modulates transcriptional memory and facilitates environmental adaptation

Transcriptional memory, by which cells respond faster to repeated stimuli, is key for cellular adaptation and organism survival. Chromatin organization has been shown to play a role in the faster response of primed cells. However, the contribution of post-transcriptional regulation is not yet explored. Here we perform a genome-wide screen to identify novel factors modulating transcriptional memory in S. cerevisiae in response to galactose. We find that depletion of the nuclear RNA exosome increases GAL1 expression in primed cells. Our work shows that gene-specific differences in intrinsic nuclear surveillance factor association can enhance both gene induction and repression in primed cells. Finally, we show that primed cells present altered levels of RNA degradation machinery and that both nuclear and cytoplasmic mRNA decay modulate transcriptional memory. Our results demonstrate that mRNA post-transcriptional regulation, and not only transcription regulation, should be considered when investigating gene expression memory.

Relative specificity as an important consideration in the big data era

Technological breakthroughs such as high-throughput methods, genomics, single-cell studies, and machine learning have fundamentally transformed research and ushered in the big data era of biology. Nevertheless, current data collections, analyses, and modeling frequently overlook relative specificity, a crucial property of molecular interactions in biochemical systems. Relative specificity describes how, for example, an enzyme reacts with its many substrates at different rates, and how this discriminatory action alone is sufficient to modulate the substrates and downstream events. As a corollary, it is not only important to comprehensively identify an enzyme’s substrates, but also critical to quantitatively determine how the enzyme interacts with the substrates and to evaluate how it shapes subsequent biological outcomes. Genomics and high-throughput techniques have greatly facilitated the studies of relative specificity in the 21st century, and its functional significance has been demonstrated in complex biochemical systems including transcription, translation, protein kinases, RNA-binding proteins, and animal microRNAs (miRNAs), although it remains ignored in most work. Here we analyze recent findings in big data and relative specificity studies and explain how the incorporation of relative specificity concept might enhance our mechanistic understanding of gene functions, biological phenomena, and human diseases.

Oligodendrocyte differentiation alters tRNA modifications and codon optimality-mediated mRNA decay

Oligodendrocytes are specialized cells that confer neuronal myelination in the central nervous system. Leukodystrophies associated with oligodendrocyte deficits and hypomyelination are known to result when a number of tRNA metabolism genes are mutated. Thus, for unknown reasons, oligodendrocytes may be hypersensitive to perturbations in tRNA biology. In this study, we survey the tRNA transcriptome in the murine oligodendrocyte cell lineage and find that specific tRNAs are hypomodified in oligodendrocytes within or near the anticodon compared to oligodendrocyte progenitor cells (OPCs). This hypomodified state may be the result of differential expression of key modification enzymes during oligodendrocyte differentiation. Moreover, we observe a concomitant relationship between tRNA hypomodification and tRNA decoding potential; observing oligodendrocyte specific alterations in codon optimality-mediated mRNA decay and ribosome transit. Our results reveal that oligodendrocytes naturally maintain a delicate, hypersensitized tRNA/mRNA axis. We suggest this axis is a potential mediator of pathology in leukodystrophies and white matter disease when further insult to tRNA metabolism is introduced.

Mitotic checkpoint gene expression is tuned by codon usage bias

The mitotic checkpoint (also called spindle assembly checkpoint, SAC) is a signaling pathway that safeguards proper chromosome segregation. Correct functioning of the SAC depends on adequate protein concentrations and appropriate stoichiometries between SAC proteins. Yet very little is known about the regulation of SAC gene expression. Here, we show in the fission yeast Schizosaccharomyces pombe that a combination of short mRNA half-lives and long protein half-lives supports stable SAC protein levels. For the SAC genes mad2+ and mad3+ , their short mRNA half-lives are caused, in part, by a high frequency of nonoptimal codons. In contrast, mad1+ mRNA has a short half-life despite a higher frequency of optimal codons, and despite the lack of known RNA-destabilizing motifs. Hence, different SAC genes employ different strategies of expression. We further show that Mad1 homodimers form co-translationally, which may necessitate a certain codon usage pattern. Taken together, we propose that the codon usage of SAC genes is fine-tuned to ensure proper SAC function. Our work shines light on gene expression features that promote spindle assembly checkpoint function and suggests that synonymous mutations may weaken the checkpoint.

Codon optimality-mediated mRNA degradation: Linking translational elongation to mRNA stability

Messenger RNA (mRNA) translation by the ribosome represents the final step of a complicated molecular dance from DNA to protein. Although classically considered a decipherer that translates a 64-word genetic code into a proteome of astonishing complexity, the ribosome can also shape the transcriptome by controlling mRNA stability. Recent work has discovered that the ribosome is an arbiter of the general mRNA degradation pathway, wherein the ribosome transit rate serves as a major determinant of transcript half-lives. Specifically, members of the degradation complex sense ribosome translocation rates as a function of ribosome elongation rates. Central to this notion is the concept of codon optimality: although all codons impact translation rates, some are deciphered quickly, whereas others cause ribosome hesitation as a consequence of relative cognate tRNA concentration. These transient pauses induce a unique ribosome conformational state that is probed by the deadenylase complex, thereby inducing an orchestrated set of events that enhance both poly(A) shortening and cap removal. Together, these data imply that the coding region of an mRNA not only encodes for protein content but also impacts protein levels through determining the transcript's fate.

Determinants of the temperature adaptation of mRNA degradation

The rate of chemical reactions increases proportionally with temperature, but the interplay of biochemical reactions permits deviations from this relation and adaptation. The degradation of individual mRNAs in yeast increased to varying degrees with temperature. We examined how these variations are influenced by the translation and codon composition of mRNAs. We developed a method that revealed the existence of a neutral half-life above which mRNAs are stabilized by translation but below which they are destabilized. The proportion of these two mRNA subpopulations remained relatively constant under different conditions, even with slow cell growth due to nutrient limitation, but heat shock reduced the proportion of translationally stabilized mRNAs. At the same time, the degradation of these mRNAs was partially temperature-compensated through Upf1, the mediator of nonsense-mediated decay. Compensation was also promoted by some asparagine and serine codons, whereas tyrosine codons promote temperature sensitization. These codons play an important role in the degradation of mRNAs encoding key cell membrane and cell wall proteins, which promote cell integrity.

Downregulation of METTL6 mitigates cell progression, migration, invasion and adhesion in hepatocellular carcinoma by inhibiting cell adhesion molecules

RNA modifications have attracted increasing interest in recent years because they have been frequently implicated in various human diseases, including cancer, highlighting the importance of dynamic post‑transcriptional modifications. Methyltransferase‑like 6 (METTL6) is a member of the RNA methyltransferase family that has been identified in many cancers; however, little is known about its specific role or mechanism of action. In the present study, we aimed to study the expression levels and functional role of METTL6 in hepatocellular carcinoma (HCC), and further investigate the relevant pathways. To this end, we systematically conducted bioinformatics analysis of METTL6 in HCC using gene expression data and clinical information from a publicly available dataset. The mRNA expression levels of METTL6 were significantly upregulated in HCC tumor tissues compared to that in adjacent non‑tumor tissues and strongly associated with poorer survival outcomes in patients with HCC. CRISPR/Cas9‑mediated knockout of METTL6 in HCC cell lines remarkably inhibited colony formation, cell proliferation, cell migration, cell invasion and cell attachment ability. RNA sequencing analysis demonstrated that knockout of METTL6 significantly suppressed the expression of cell adhesion‑related genes. However, chromatin immunoprecipitation sequencing results revealed no significant differences in enhancer activities between cells, which suggests that METTL6 may regulate genes of interest post‑transcriptionally. In addition, it was demonstrated for the first time that METTL6 was localized in the cytosol as detected by immunofluorescence analysis, which indicates the plausible location of RNA modification mediated by METTL6. Our findings provide further insight into the function of RNA modifications in cancer and suggest a possible role of METTL6 as a therapeutic target in HCC.

Codon optimality in cancer

A key characteristic of cancer cells is their increased proliferative capacity, which requires elevated levels of protein synthesis. The process of protein synthesis involves the translation of codons within the mRNA coding sequence into a string of amino acids to form a polypeptide chain. As most amino acids are encoded by multiple codons, the nucleotide sequence of a coding region can vary dramatically without altering the polypeptide sequence of the encoded protein. Although mutations that do not alter the final amino acid sequence are often thought of as silent/synonymous, these can still have dramatic effects on protein output. Because each codon has a distinct translation elongation rate and can differentially impact mRNA stability, each codon has a different degree of 'optimality' for protein synthesis. Recent data demonstrates that the codon preference of a transcriptome matches the abundance of tRNAs within the cell and that this supply and demand between tRNAs and mRNAs varies between different cell types. The largest observed distinction is between mRNAs encoding proteins associated with proliferation or differentiation. Nevertheless, precisely how codon optimality and tRNA expression levels regulate cell fate decisions and their role in malignancy is not fully understood. This review describes the current mechanistic understanding on codon optimality, its role in malignancy and discusses the potential to target codon optimality therapeutically in the context of cancer.

tRNA copy number and codon usage in the sea cucumber genome provide insights into adaptive translation for saponin biosynthesis

Genomic tRNA copy numbers determine cytoplasmic tRNA abundances, which in turn influence translation efficiency, but the underlying mechanism is not well understood. Using the sea cucumber Apostichopus japonicus as a model, we combined genomic sequence, transcriptome expression and ecological food resource data to study its codon usage adaptation. The results showed that, unlike intragenic non-coding RNAs, transfer RNAs (tRNAs) tended to be transcribed independently. This may be attributed to their specific Pol III promoters that lack transcriptional regulation, which may underlie the correlation between genomic copy number and cytoplasmic abundance of tRNAs. Moreover, codon usage optimization was mostly restrained by a gene's amino acid sequence, which might be a compromise between functionality and translation efficiency for stress responses were highly optimized for most echinoderms, while enzymes for saponin biosynthesis (LAS, CYPs and UGTs) were especially optimized in sea cucumbers, which might promote saponin synthesis as a defence strategy. The genomic tRNA content of A. japonicus was positively correlated with amino acid content in its natural food particles, which should promote its efficiency in protein synthesis. We propose that coevolution between genomic tRNA content and codon usage of sea cucumbers facilitates their saponin synthesis and survival using food resources with low nutrient content.

A multispecies BCO2 beak color polymorphism in the Darwin's finch radiation

Carotenoid-based polymorphisms are widespread in populations of birds, fish, and reptiles,¹ but generally little is known about the factors affecting their maintenance in populations.² We report a combined field and molecular-genetic investigation of a nestling beak color polymorphism in Darwin's finches. Beaks are pink or yellow, and yellow is recessive.³ Here we show that the polymorphism arose in the Galápagos half a million years ago through a mutation associated with regulatory change in the BCO2 gene and is shared by 14 descendant species. The polymorphism is probably a balanced polymorphism, maintained by ecological selection associated with survival and diet. In cactus finches, the frequency of the yellow genotype is correlated with cactus fruit abundance and greater hatching success and may be altered by introgressive hybridization. Polymorphisms that are hidden as adults, as here, may be far more common than is currently recognized, and contribute to diversification in ways that are yet to be discovered.

Differential regulation of mRNA fate by the human Ccr4-Not complex is driven by coding sequence composition and mRNA localization

BACKGROUND

Regulation of protein output at the level of translation allows for a rapid adaptation to dynamic changes to the cell's requirements. This precise control of gene expression is achieved by complex and interlinked biochemical processes that modulate both the protein synthesis rate and stability of each individual mRNA. A major factor coordinating this regulation is the Ccr4-Not complex. Despite playing a role in most stages of the mRNA life cycle, no attempt has been made to take a global integrated view of how the Ccr4-Not complex affects gene expression.

RESULTS

This study has taken a comprehensive approach to investigate post-transcriptional regulation mediated by the Ccr4-Not complex assessing steady-state mRNA levels, ribosome position, mRNA stability, and protein production transcriptome-wide. Depletion of the scaffold protein CNOT1 results in a global upregulation of mRNA stability and the preferential stabilization of mRNAs enriched for G/C-ending codons. We also uncover that mRNAs targeted to the ER for their translation have reduced translational efficiency when CNOT1 is depleted, specifically downstream of the signal sequence cleavage site. In contrast, translationally upregulated mRNAs are normally localized in p-bodies, contain disorder-promoting amino acids, and encode nuclear localized proteins. Finally, we identify ribosome pause sites that are resolved or induced by the depletion of CNOT1.

CONCLUSIONS

We define the key mRNA features that determine how the human Ccr4-Not complex differentially regulates mRNA fate and protein synthesis through a mechanism linked to codon composition, amino acid usage, and mRNA localization.

Integrated transcriptomics and proteomics revealed the distinct toxicological effects of multi-metal contamination on oysters

The Pearl River Estuary (PRE) is the largest estuary in southern China and under high metal stress. In the present study, we employed an integrated method of transcriptomics and proteomics to investigate the ecotoxicological effects of trace metals on the Hong Kong oyster Crassostrea hongkongensis. Three oyster populations with distinct spatial distributions of metals were sampled, including the Control (Station QA, the lowest metal levels), the High Cd (Station JZ, the highest Cd), and the High Zn-Cu-Cr-Ni (Station LFS, with the highest levels of zinc, copper, chromium, and nickel). Dominant metals in oysters were differentiated by principal component analysis (PCA), and theirgene and protein profiles were studied using RNA-seq and iTRAQ techniques. Of the 2250 proteins identified at both protein and RNA levels, 70 proteins exhibited differential expressions in response to metal stress in oysters from the two contaminated stations. There were 8 proteins altered at both stations, with the potential effects on mitochondria and endoplasmic reticulum by Ag. The genotoxicity, including impaired DNA replication and transcription, was specifically observed in the High Cd oysters with the dominating influence of Cd. The structural components (cytoskeleton and chromosome-associated proteins) were impaired by the over-accumulated Cu, Zn, Cr, and Ni at Station LFS. However, enhanced tRNA biogenesis and exosome activity might help the oysters to alleviate the toxicities resulting from their exposure to these metals. Our study provided comprehensive information on the molecular changes in oysters at both protein and RNA levels in responding to multi-levels of trace metal stress.

High-throughput 5'P sequencing enables the study of degradation-associated ribosome stalls

RNA degradation is critical for gene expression and mRNA quality control. mRNA degradation is connected to the translation process up to the degree that 5'-3' mRNA degradation follows the last translating ribosome. Here, we present an improved high-throughput 5'P degradome RNA-sequencing method (HT-5Pseq). HT-5Pseq is easy, scalable, and uses affordable duplex-specific nuclease-based rRNA depletion. We investigate in vivo ribosome stalls focusing on translation termination. By comparing ribosome stalls identified by ribosome profiling, disome-seq and HT-5Pseq, we find that degradation-associated ribosome stalls are often enriched in Arg preceding the stop codon. On the contrary, mRNAs depleted for those stalls use more frequently a TAA stop codon preceded by hydrophobic amino acids. Finally, we show that termination stalls found by HT-5Pseq, and not by other approaches, are associated with decreased mRNA stability. Our work suggests that ribosome stalls associated with mRNA decay can be easily captured by investigating the 5'P degradome.

Genome-wide analysis of lncRNA stability in human

Transcript stability is associated with many biological processes, and the factors affecting mRNA stability have been extensively studied. However, little is known about the features related to human long noncoding RNA (lncRNA) stability. By inhibiting transcription and collecting samples in 10 time points, genome-wide RNA-seq studies was performed in human lung adenocarcinoma cells (A549) and RNA half-life datasets were constructed. The following observations were obtained. First, the half-life distributions of both lncRNAs and messanger RNAs (mRNAs) with one exon (lnc-human1 and m-human1) were significantly different from those of both lncRNAs and mRNAs with more than one exon (lnc-human2 and m-human2). Furthermore, some factors such as full-length transcript secondary structures played a contrary role in lnc-human1 and m-human2. Second, through the half-life comparisons of nucleus- and cytoplasm-specific and common lncRNAs and mRNAs, lncRNAs (mRNAs) in the nucleus were found to be less stable than those in the cytoplasm, which was derived from transcripts themselves rather than cellular location. Third, kmers-based protein−RNA or RNA−RNA interactions promoted lncRNA stability from lnc-human1 and decreased mRNA stability from m-human2 with high probability. Finally, through applying deep learning−based regression, a non-linear relationship was found to exist between the half-lives of lncRNAs (mRNAs) and related factors. The present study established lncRNA and mRNA half-life regulation networks in the A549 cell line and shed new light on the degradation behaviors of both lncRNAs and mRNAs.

Transcript stability is important for many biological processes. However, little is known about the features related to human lncRNA stability. Through quantitative analysis between the half-lives of lncRNAs (mRNAs) and various factors, we found a nonlinear relationship between the half-lives of lncRNAs (mRNAs) and the related factors and their combinations. Our research provided a comprehensive understanding of lncRNA stability. Further efforts are needed to develop an accurate quantitative prediction model for the half-lives of lncRNA (mRNA).

A Trans-Omics Comparison Reveals Common Gene Expression Strategies in Four Model Organisms and Exposes Similarities and Differences between Them

The ultimate goal of gene expression regulation is on the protein level. However, because the amounts of mRNAs and proteins are controlled by their synthesis and degradation rates, the cellular amount of a given protein can be attained by following different strategies. By studying omics data for six expression variables (mRNA and protein amounts, plus their synthesis and decay rates), we previously demonstrated the existence of common expression strategies (CESs) for functionally related genes in the yeast Saccharomyces cerevisiae. Here we extend that study to two other eukaryotes: the yeast Schizosaccharomyces pombe and cultured human HeLa cells. We also use genomic data from the model prokaryote Escherichia coli as an external reference. We show that six-variable profiles (6VPs) can be constructed for every gene and that these 6VPs are similar for genes with similar functions in all the studied organisms. The differences in 6VPs between organisms can be used to establish their phylogenetic relationships. The analysis of the correlations among the six variables supports the hypothesis that most gene expression control occurs in actively growing organisms at the transcription rate level, and that translation plays a minor role. We propose that living organisms use CESs for the genes acting on the same physiological pathways, especially for those belonging to stable macromolecular complexes, but CESs have been modeled by evolution to adapt to the specific life circumstances of each organism.

Loss of Cnot6l Impairs Inosine RNA Modifications in Mouse Oocytes

Mammalian oocytes must degrade maternal transcripts through a process called translational mRNA decay, in which maternal mRNA undergoes translational activation, followed by deadenylation and mRNA decay. Once a transcript is translationally activated, it becomes deadenylated by the CCR4-NOT complex. Knockout of CCR4-NOT Transcription Complex Subunit 6 Like (Cnot6l), a deadenylase within the CCR4-NOT complex, results in mRNA decay defects during metaphase I (MI) entry. Knockout of B-cell translocation gene-4 (Btg4), an adaptor protein of the CCR4-NOT complex, results in mRNA decay defects following fertilization. Therefore, mechanisms controlling mRNA turnover have significant impacts on oocyte competence and early embryonic development. Post-transcriptional inosine RNA modifications can impact mRNA stability, possibly through a translation mechanism. Here, we assessed inosine RNA modifications in oocytes, eggs, and embryos from Cnot6l^-/- and Btg4^-/- mice, which display stabilization of mRNA and over-translation of the stabilized transcripts. If inosine modifications have a role in modulating RNA stability, we hypothesize that in these mutant backgrounds, we would observe changes or a disruption in inosine mRNA modifications. To test this, we used a computational approach to identify inosine RNA modifications in total and polysomal RNA-seq data during meiotic maturation (GV, MI, and MII stages). We observed pronounced depletion of inosine mRNA modifications in samples from Cnot6l^-/-, but not in Btg4^-/- mice. Additionally, analysis of ribosome-associated RNA revealed clearance of inosine modified mRNA. These observations suggest a novel mechanism of mRNA clearance during oocyte maturation, in which inosine-containing transcripts decay in an independent, but parallel mechanism to CCR4-NOT deadenylation.

Crosstalk between codon optimality and cis-regulatory elements dictates mRNA stability

BACKGROUND

The regulation of messenger RNA (mRNA) stability has a profound impact on gene expression dynamics during embryogenesis. For example, in animals, maternally deposited mRNAs are degraded after fertilization to enable new developmental trajectories. Regulatory sequences in 3' untranslated regions (3'UTRs) have long been considered the central determinants of mRNA stability. However, recent work indicates that the coding sequence also possesses regulatory information. Specifically, translation in cis impacts mRNA stability in a codon-dependent manner. However, the strength of this mechanism during embryogenesis, as well as its relationship with other known regulatory elements, such as microRNA, remains unclear.

RESULTS

Here, we show that codon composition is a major predictor of mRNA stability in the early embryo. We show that this mechanism works in combination with other cis-regulatory elements to dictate mRNA stability in zebrafish and Xenopus embryos as well as in mouse and human cells. Furthermore, we show that microRNA targeting efficacy can be affected by substantial enrichment of optimal (stabilizing) or non-optimal (destabilizing) codons. Lastly, we find that one microRNA, miR-430, antagonizes the stabilizing effect of optimal codons during early embryogenesis in zebrafish.

CONCLUSIONS

By integrating the contributions of different regulatory mechanisms, our work provides a framework for understanding how combinatorial control of mRNA stability shapes the gene expression landscape.

Rqc1 and other yeast proteins containing highly positively charged sequences are not targets of the RQC complex

Previous work has suggested that highly positively charged protein segments coded by rare codons or poly (A) stretches induce ribosome stalling and translational arrest through electrostatic interactions with the negatively charged ribosome exit tunnel, leading to inefficient elongation. This arrest leads to the activation of the Ribosome Quality Control (RQC) pathway and results in low expression of these reporter proteins. However, the only endogenous yeast proteins known to activate the RQC are Rqc1, a protein essential for RQC function, and Sdd1, a protein with unknown function, both of which contain polybasic sequences. To explore the generality of this phenomenon, we investigated whether the RQC complex controls the expression of other proteins with polybasic sequences. We showed by ribosome profiling data analysis and western blot that proteins containing polybasic sequences similar to, or even more positively charged than those of Rqc1 and Sdd1, were not targeted by the RQC complex. We also observed that the previously reported Ltn1-dependent regulation of Rqc1 is posttranslational, independent of the RQC activity. Taken together, our results suggest that RQC should not be regarded as a general regulatory pathway for the expression of highly positively charged proteins in yeast.

Natural variation of the cardiac transcriptome in humans

We investigated the gene-expression variation among humans by analysing previously published mRNA-seq and ribosome footprint profiling of heart left-ventricles from healthy donors. We ranked the genes according to their coefficient of variation values and found that the top 5% most variable genes had special features compared to the rest of the genome, such as lower mRNA levels and shorter half-lives coupled to increased translation efficiency. We observed that these genes are mostly involved with immune response and have a pleiotropic effect on disease phenotypes, indicating that asymptomatic conditions contribute to the gene expression diversity of healthy individuals.

What determines eukaryotic translation elongation: recent molecular and quantitative analyses of protein synthesis

Protein synthesis from mRNA is an energy-intensive and tightly controlled cellular process. Translation elongation is a well-coordinated, multifactorial step in translation that undergoes dynamic regulation owing to cellular state and environmental determinants. Recent studies involving genome-wide approaches have uncovered some crucial aspects of translation elongation including the mRNA itself and the nascent polypeptide chain. Additionally, these studies have fuelled quantitative and mathematical modelling of translation elongation. In this review, we provide a comprehensive overview of the key determinants of translation elongation. We discuss consequences of ribosome stalling or collision, and how the cells regulate translation in case of such events. Next, we review theoretical approaches and widely used mathematical models that have become an essential ingredient to interpret complex molecular datasets and study translation dynamics quantitatively. Finally, we review recent advances in live-cell reporter and related analysis techniques, to monitor the translation dynamics of single cells and single-mRNA molecules in real time.

Natural Selection Shapes Codon Usage in the Human Genome

Synonymous codon usage has been identified as a determinant of translational efficiency and mRNA stability in model organisms and human cell lines. However, whether natural selection shapes human codon content to optimize translation efficiency is unclear. Furthermore, aside from those that affect splicing, synonymous mutations are typically ignored as potential contributors to disease. Using genetic sequencing data from nearly 200,000 individuals, we uncover clear evidence that natural selection optimizes codon content in the human genome. In deriving intolerance metrics to quantify gene-level constraint on synonymous variation, we discover that dosage-sensitive genes, DNA-damage-response genes, and cell-cycle-regulated genes are particularly intolerant to synonymous variation. Notably, we illustrate that reductions in codon optimality in BRCA1 can attenuate its function. Our results reveal that synonymous mutations most likely play an underappreciated role in human variation.

Machine learning using intrinsic genomic signatures for rapid classification of novel pathogens: COVID-19 case study

The 2019 novel coronavirus (renamed SARS-CoV-2, and generally referred to as the COVID-19 virus) has spread to 184 countries with over 1.5 million confirmed cases. Such major viral outbreaks demand early elucidation of taxonomic classification and origin of the virus genomic sequence, for strategic planning, containment, and treatment. This paper identifies an intrinsic COVID-19 virus genomic signature and uses it together with a machine learning-based alignment-free approach for an ultra-fast, scalable, and highly accurate classification of whole COVID-19 virus genomes. The proposed method combines supervised machine learning with digital signal processing (MLDSP) for genome analyses, augmented by a decision tree approach to the machine learning component, and a Spearman's rank correlation coefficient analysis for result validation. These tools are used to analyze a large dataset of over 5000 unique viral genomic sequences, totalling 61.8 million bp, including the 29 COVID-19 virus sequences available on January 27, 2020. Our results support a hypothesis of a bat origin and classify the COVID-19 virus as Sarbecovirus, within Betacoronavirus. Our method achieves 100% accurate classification of the COVID-19 virus sequences, and discovers the most relevant relationships among over 5000 viral genomes within a few minutes, ab initio, using raw DNA sequence data alone, and without any specialized biological knowledge, training, gene or genome annotations. This suggests that, for novel viral and pathogen genome sequences, this alignment-free whole-genome machine-learning approach can provide a reliable real-time option for taxonomic classification.

Pervasive Strong Selection at the Level of Codon Usage Bias in Drosophila melanogaster

Codon usage bias (CUB), where certain codons are used more frequently than expected by chance, is a ubiquitous phenomenon and occurs across the tree of life. The dominant paradigm is that the proportion of preferred codons is set by weak selection. While experimental changes in codon usage have at times shown large phenotypic effects in contrast to this paradigm, genome-wide population genetic estimates have supported the weak selection model. Here we use deep genomic population sequencing of two Drosophila melanogaster populations to measure selection on synonymous sites in a way that allowed us to estimate the prevalence of both weak and strong purifying selection. We find that selection in favor of preferred codons ranges from weak (|Nes| ∼ 1) to strong (|Nes| > 10), with strong selection acting on 10-20% of synonymous sites in preferred codons. While previous studies indicated that selection at synonymous sites could be strong, this is the first study to detect and quantify strong selection specifically at the level of CUB. Further, we find that CUB-associated polymorphism accounts for the majority of strong selection on synonymous sites, with secondary contributions of splicing (selection on alternatively spliced genes, splice junctions, and spliceosome-bound sites) and transcription factor binding. Our findings support a new model of CUB and indicate that the functional importance of CUB, as well as synonymous sites in general, have been underestimated.

From reporters to endogenous genes: the impact of the first five codons on translation efficiency in Escherichia coli

Translation initiation is a critical step in the regulation of protein synthesis, and it is subjected to different control mechanisms, such as 5' UTR secondary structure and initiation codon context, that can influence the rates at which initiation and consequentially translation occur. For some genes, translation elongation also affects the rate of protein synthesis. With a GFP library containing nearly all possible combinations of nucleotides from the 3rd to the 5th codon positions in the protein coding region of the mRNA, it was previously demonstrated that some nucleotide combinations increased GFP expression up to four orders of magnitude. While it is clear that the codon region from positions 3 to 5 can influence protein expression levels of artificial constructs, its impact on endogenous proteins is still unknown. Through bioinformatics analysis, we identified the nucleotide combinations of the GFP library in Escherichia coli genes and examined the correlation between the expected levels of translation according to the GFP data with the experimental measures of protein expression. We observed that E. coli genes were enriched with the nucleotide compositions that enhanced protein expression in the GFP library, but surprisingly, it seemed to affect the translation efficiency only marginally. Nevertheless, our data indicate that different enterobacteria present similar nucleotide composition enrichment as E. coli, suggesting an evolutionary pressure towards the conservation of short translational enhancer sequences.