Codon Usage and 3' UTR Length Determine Maternal mRNA Stability in Zebrafish

Deadenylation by the CCR4-NOT complex contributes to the turnover of hairy-related mRNAs in the zebrafish segmentation clock

In the zebrafish segmentation clock, hairy/enhancer of split-related genes her1, her7, and hes6 encodes components of core oscillators. Since the expression of cyclic genes proceeds rapidly in the presomitic mesoderm (PSM), these hairy-related mRNAs are subject to strict post-transcriptional regulation. In this study, we demonstrate that inhibition of the CCR4-NOT deadenylase complex lengthens poly(A) tails of hairy-related mRNAs and increases the amount of these mRNAs, which is accompanied by defective somite segmentation. In transgenic embryos, we show that EGFP mRNAs with 3'UTRs of hairy-related genes exhibit turnover similar to endogenous mRNAs. Our results suggest that turnover rates of her1, her7, and hes6 mRNAs are differently regulated by the CCR4-NOT deadenylase complex possibly through their 3'UTRs in the zebrafish PSM.

Codon usage regulates human KRAS expression at both transcriptional and translational levels

KRAS and HRAS are highly homologous oncogenic Ras GTPase family members that are mutated in a wide spectrum of human cancers. Despite having high amino acid identity, KRAS and HRAS have very different codon usage biases: the HRAS gene contains many common codons, and KRAS is enriched for rare codons. Rare codons in KRAS suppress its protein expression, which has been shown to affect both normal and cancer biology in mammals. Here, using HRAS or KRAS expression in different human cell lines and in vitro transcription and translation assays, we show that KRAS rare codons inhibit both translation efficiency and transcription and that the contribution of these two processes varies among different cell lines. We observed that codon usage regulates mRNA translation efficiency such that WT KRAS mRNA is poorly translated. On the other hand, common codons increased transcriptional rates by promoting activating histone modifications and recruitment of transcriptional coactivators. Finally, we found that codon usage also influences KRAS protein conformation, likely because of its effect on co-translational protein folding. Together, our results reveal that codon usage has multidimensional effects on protein expression, ranging from effects on transcription to protein folding in human cells.

Translation elongation and mRNA stability are coupled through the ribosomal A-site

Messenger RNA (mRNA) degradation plays a critical role in regulating transcript levels in eukaryotic cells. Previous work by us and others has shown that codon identity exerts a powerful influence on mRNA stability. In Saccharomyces cerevisiae, studies using a handful of reporter mRNAs show that optimal codons increase translation elongation rate, which in turn increases mRNA stability. However, a direct relationship between elongation rate and mRNA stability has not been established across the entire yeast transcriptome. In addition, there is evidence from work in higher eukaryotes that amino acid identity influences mRNA stability, raising the question as to whether the impact of translation elongation on mRNA decay is at the level of tRNA decoding, amino acid incorporation, or some combination of each. To address these questions, we performed ribosome profiling of wild-type yeast. In good agreement with other studies, our data showed faster codon-specific elongation over optimal codons and faster transcript-level elongation correlating with transcript optimality. At both the codon-level and transcript-level, faster elongation correlated with increased mRNA stability. These findings were reinforced by showing increased translation efficiency and kinetics for a panel of 11 HIS3 reporter mRNAs of increasing codon optimality. While we did observe that elongation measured by ribosome profiling is composed of both amino acid identity and synonymous codon effects, further analyses of these data establish that A-site tRNA decoding rather than other steps of translation elongation is driving mRNA decay in yeast.

Somatic and Germline MicroRNAs Form Distinct Silencing Complexes to Regulate Their Target mRNAs Differently

Animal germ cells possess a specific post-transcriptional regulatory context allowing the storage of maternal transcripts in the oocyte until their translation at a specific point in early development. As key regulators of gene expression, miRNAs repress translation mainly through mRNA destabilization. Thus, germline miRNAs likely use distinct ways to regulate their targets. Here, we use C. elegans to compare miRNA function within germline and somatic tissues. We show that the same miRNA displays tissue-specific gene regulatory mechanisms. While translational repression occurs in both tissues, targeted mRNAs are instead stabilized in the germline. Comparative analyses of miRNA silencing complexes (miRISC) demonstrate that their composition differs from germline to soma. We show that germline miRNA targets preferentially localize to perinuclear regions adjacent to P granules, and their repression is dependent on the core P granule component GLH-1. Together, our findings reveal the existence of different miRISC in animals that affect targeted mRNAs distinctively.

Poor codon optimality as a signal to degrade transcripts with frameshifts

Frameshifting errors are common and mRNA quality control pathways, such as nonsense-mediated decay (NMD), exist to degrade these aberrant transcripts. Recent work has shown the existence of a genetic link between NMD and codon-usage mediated mRNA decay. Here we present computational evidence that these pathways are synergic for removing frameshifts.

Repertoires of tRNAs: The Couplers of Genomics and Proteomics

The pool of transfer RNA (tRNA) molecules in cells allows the ribosome to decode genetic information. This repertoire of molecular decoders is positioned in the crossroad of the genome, the transcriptome, and the proteome. Omics and systems biology now allow scientists to explore the entire repertoire of tRNAs of many organisms, revealing basic exciting biology. The tRNA gene set of hundreds of species is now characterized, in addition to the tRNA genes of organelles and viruses. Genes encoding tRNAs for certain anticodon types appear in dozens of copies in a genome, while others are universally absent from any genome. Transcriptome measurement of tRNAs is challenging, but in recent years new technologies have allowed researchers to determine the dynamic expression patterns of tRNAs. These advances reveal that availability of ready-to-translate tRNA molecules is highly controlled by several transcriptional and posttranscriptional regulatory processes. This regulation shapes the proteome according to the cellular state. The tRNA pool profoundly impacts many aspects of cellular and organismal life, including protein expression level, translation accuracy, adequacy of folding, and even mRNA stability. As a result, the shape of the tRNA pool affects organismal health and may participate in causing conditions such as cancer and neurological conditions.

High resolution annotation of zebrafish transcriptome using long-read sequencing

With the emergence of zebrafish as an important model organism, a concerted effort has been made to study its transcriptome. This effort is limited, however, by gaps in zebrafish annotation, which are especially pronounced concerning transcripts dynamically expressed during zygotic genome activation (ZGA). To date, short-read sequencing has been the principal technology for zebrafish transcriptome annotation. In part because these sequence reads are too short for assembly methods to resolve the full complexity of the transcriptome, the current annotation is rudimentary. By providing direct observation of full-length transcripts, recently refined long-read sequencing platforms can dramatically improve annotation coverage and accuracy. Here, we leveraged the SMRT platform to study the transcriptome of zebrafish embryos before and after ZGA. Our analysis revealed additional novelty and complexity in the zebrafish transcriptome, identifying 2539 high-confidence novel transcripts that originated from previously unannotated loci and 1835 high-confidence new isoforms in previously annotated genes. We validated these findings using a suite of computational approaches including structural prediction, sequence homology, and functional conservation analyses, as well as by confirmatory transcript quantification with short-read sequencing data. Our analyses provided insight into new homologs and paralogs of functionally important proteins and noncoding RNAs, isoform switching occurrences, and different classes of novel splicing events. Several novel isoforms representing distinct splicing events were validated through PCR experiments, including the discovery and validation of a novel 8-kb transcript spanning multiple mir-430 elements, an important driver of early development. Our study provides a significantly improved zebrafish transcriptome annotation resource.

Decoupling the impact of microRNAs on translational repression versus RNA degradation in embryonic stem cells

Translation and mRNA degradation are intimately connected, yet the mechanisms that link them are not fully understood. Here, we studied these mechanisms in embryonic stem cells (ESCs). Transcripts showed a wide range of stabilities, which correlated with their relative translation levels and that did not change during early ESC differentiation. The protein DHH1 links translation to mRNA stability in yeast; however, loss of the mammalian homolog, DDX6, in ESCs did not disrupt the correlation across transcripts. Instead, the loss of DDX6 led to upregulated translation of microRNA targets, without concurrent changes in mRNA stability. The Ddx6 knockout cells were phenotypically and molecularly similar to cells lacking all microRNAs (Dgcr8 knockout ESCs). These data show that the loss of DDX6 can separate the two canonical functions of microRNAs: translational repression and transcript destabilization. Furthermore, these data uncover a central role for translational repression independent of transcript destabilization in defining the downstream consequences of microRNA loss.

The translational regulation of maternal mRNAs in time and space

Since their discovery, the study of maternal mRNAs has led to the identification of mechanisms underlying their spatiotemporal regulation within the context of oogenesis and early embryogenesis. Following synthesis in the oocyte, maternal mRNAs are translationally silenced and sequestered into storage in cytoplasmic granules. At the same time, their unique distribution patterns throughout the oocyte and embryo are tightly controlled and connected to their functions in downstream embryonic processes. At certain points in oogenesis and early embryogenesis, maternal mRNAs are translationally activated to perform their functions in a timely manner. The cytoplasmic polyadenylation machinery is responsible for the translational activation of maternal mRNAs, and its role in initiating the maternal to zygotic transition events has recently come to light. Here, we summarize the current knowledge on maternal mRNA regulation, with particular focus on cytoplasmic polyadenylation as a mechanism for translational regulation.

Dynamic transcriptomic m6A decoration: writers, erasers, readers and functions in RNA metabolism

N⁶-methyladenosine (m⁶A) is a chemical modification present in multiple RNA species, being most abundant in mRNAs. Studies on enzymes or factors that catalyze, recognize, and remove m⁶A have revealed its comprehensive roles in almost every aspect of mRNA metabolism, as well as in a variety of physiological processes. This review describes the current understanding of the m⁶A modification, particularly the functions of its writers, erasers, readers in RNA metabolism, with an emphasis on its role in regulating the isoform dosage of mRNAs.

Pervasive allele-specific regulation on RNA decay in hybrid mice

Cellular RNA abundance is determined by both RNA transcription and decay. Therefore, change in RNA abundance, which can drive phenotypic diversity between different species, could arise from genetic variants affecting either process. However, previous studies in the evolution of RNA expression have been largely focused on transcription. Here, to globally investigate the effects of cis-regulatory divergence on RNA decay in mammals for the first time, we quantified allele-specific differences in RNA decay rates (ASD) in an F1 hybrid mouse. Out of 8,815 genes with sufficient data, we identified 621 genes exhibiting significant cis-divergence. Systematic analysis of these genes revealed that the genetic variants affecting microRNA binding and RNA secondary structures contribute to the observed divergences. Finally, we demonstrated that although the divergences in RNA abundance were predominantly determined by allelic differences in RNA transcription, most genes with significant ASD did not exhibit significant difference in RNA abundance. For these genes, the apparently compensatory effect between the allelic differences in RNA transcription and ASD suggests that changes in RNA decay could serve as important means to stabilize RNA abundances during mammalian evolution.

The Interplay between the RNA Decay and Translation Machinery in Eukaryotes

RNA decay plays a major role in regulating gene expression and is tightly networked with other aspects of gene expression to effectively coordinate post-transcriptional regulation. The goal of this work is to provide an overview of the major factors and pathways of general messenger RNA (mRNA) decay in eukaryotic cells, and then discuss the effective interplay of this cytoplasmic process with the protein synthesis machinery. Given the transcript-specific and fluid nature of mRNA stability in response to changing cellular conditions, understanding the fundamental networking between RNA decay and translation will provide a foundation for a complete mechanistic understanding of this important aspect of cell biology.

Enhancing heterologous expression in Chlamydomonas reinhardtii by transcript sequence optimization

Various species of microalgae have recently emerged as promising host-organisms for use in biotechnology industries due to their unique properties. These include efficient conversion of sunlight into organic compounds, the ability to grow in extreme conditions and the occurrence of numerous post-translational modification pathways. However, the inability to obtain high levels of nuclear heterologous gene expression in microalgae hinders the development of the entire field. To overcome this limitation, we analyzed different sequence optimization algorithms while studying the effect of transcript sequence features on heterologous expression in the model microalga Chlamydomonas reinhardtii, whose genome consists of rare features such as a high GC content. Based on the analysis of genomic data, we created eight unique sequences coding for a synthetic ferredoxin-hydrogenase enzyme, used here as a reporter gene. Following in silico design, these synthetic genes were transformed into the C. reinhardtii nucleus, after which gene expression levels were measured. The empirical data, measured in vivo show a discrepancy of up to 65-fold between the different constructs. In this work we demonstrate how the combination of computational methods and our empirical results enable us to learn about the way gene expression is encoded in the C. reinhardtii transcripts. We describe the deleterious effect on overall expression of codons encoding for splicing signals. Subsequently, our analysis shows that utilization of a frequent subset of preferred codons results in elevated transcript levels, and that mRNA folding energy in the vicinity of translation initiation significantly affects gene expression.

RNA sequencing of early round goby embryos reveals that maternal experiences can shape the maternal RNA contribution in a wild vertebrate

Background

It has been proposed that non-genetic inheritance could promote species fitness. Non-genetic inheritance could allow offspring to benefit from the experience of their parents, and could advocate pre-adaptation to prevailing and potentially selective conditions. Indeed, adaptive parental effects have been modeled and observed, but the molecular mechanisms behind them are far from understood.

Results

In the present study, we investigated whether maternal RNA can carry information about environmental conditions experienced by the mother in a wild vertebrate. Maternal RNA directs the development of the early embryo in many non-mammalian vertebrates and invertebrates. However, it is not known whether vertebrate maternal RNA integrates information about the parental environment. We sequenced the maternal RNA contribution from a model that we expected to rely on parental effects: the invasive benthic fish species Neogobius melanostomus (Round Goby). We found that maternal RNA expression levels correlated with the water temperature experienced by the mother before oviposition, and identified temperature-responsive gene groups such as core nucleosome components or the microtubule cytoskeleton.

Conclusions

Our findings suggest that the maternal RNA contribution may incorporate environmental information. Maternal RNA should therefore be considered a potentially relevant pathway for non-genetic inheritance. Also, the ability of a species to integrate environmental information in the maternal RNA contribution could potentially contribute to species fitness and may also play a role in extraordinary adaptive success stories of invasive species such as the round goby.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1132-2) contains supplementary material, which is available to authorized users.

Codon-Resolution Analysis Reveals a Direct and Context-Dependent Impact of Individual Synonymous Mutations on mRNA Level

Codon usage bias (CUB) refers to the observation that synonymous codons are not used equally frequently in a genome. CUB is stronger in more highly expressed genes, a phenomenon commonly explained by stronger natural selection on translational accuracy and/or efficiency among these genes. Nevertheless, this phenomenon could also occur if CUB regulates gene expression at the mRNA level, a hypothesis that has not been tested until recently. Here, we attempt to quantify the impact of synonymous mutations on mRNA level in yeast using 3,556 synonymous variants of a heterologous gene encoding green fluorescent protein (GFP) and 523 synonymous variants of an endogenous gene TDH3. We found that mRNA level was positively correlated with CUB among these synonymous variants, demonstrating a direct role of CUB in regulating transcript concentration, likely via regulating mRNA degradation rate, as our additional experiments suggested. More importantly, we quantified the effects of individual synonymous mutations on mRNA level and found them dependent on 1) CUB and 2) mRNA secondary structure, both in proximal sequence contexts. Our study reveals the pleiotropic effects of synonymous codon usage and provides an additional explanation for the well-known correlation between CUB and gene expression level.

Codon usage biases co-evolve with transcription termination machinery to suppress premature cleavage and polyadenylation

Codon usage biases are found in all genomes and influence protein expression levels. The codon usage effect on protein expression was thought to be mainly due to its impact on translation. Here, we show that transcription termination is an important driving force for codon usage bias in eukaryotes. Using Neurospora crassa as a model organism, we demonstrated that introduction of rare codons results in premature transcription termination (PTT) within open reading frames and abolishment of full-length mRNA. PTT is a wide-spread phenomenon in Neurospora, and there is a strong negative correlation between codon usage bias and PTT events. Rare codons lead to the formation of putative poly(A) signals and PTT. A similar role for codon usage bias was also observed in mouse cells. Together, these results suggest that codon usage biases co-evolve with the transcription termination machinery to suppress premature termination of transcription and thus allow for optimal gene expression.

Codon choice directs constitutive mRNA levels in trypanosomes

Selective transcription of individual protein coding genes does not occur in trypanosomes and the cellular copy number of each mRNA must be determined post-transcriptionally. Here, we provide evidence that codon choice directs the levels of constitutively expressed mRNAs. First, a novel codon usage metric, the gene expression codon adaptation index (geCAI), was developed that maximised the relationship between codon choice and the measured abundance for a transcriptome. Second, geCAI predictions of mRNA levels were tested using differently coded GFP transgenes and were successful over a 25-fold range, similar to the variation in endogenous mRNAs. Third, translation was necessary for the accelerated mRNA turnover resulting from codon choice. Thus, in trypanosomes, the information determining the levels of most mRNAs resides in the open reading frame and translation is required to access this information.

Acute multi-sgRNA knockdown of KEOPS complex genes reproduces the microcephaly phenotype of the stable knockout zebrafish model

Until recently, morpholino oligonucleotides have been widely employed in zebrafish as an acute and efficient loss-of-function assay. However, off-target effects and reproducibility issues when compared to stable knockout lines have compromised their further use. Here we employed an acute CRISPR/Cas approach using multiple single guide RNAs targeting simultaneously different positions in two exemplar genes (osgep or tprkb) to increase the likelihood of generating mutations on both alleles in the injected F0 generation and to achieve a similar effect as morpholinos but with the reproducibility of stable lines. This multi single guide RNA approach resulted in median likelihoods for at least one mutation on each allele of >99% and sgRNA specific insertion/deletion profiles as revealed by deep-sequencing. Immunoblot showed a significant reduction for Osgep and Tprkb proteins. For both genes, the acute multi-sgRNA knockout recapitulated the microcephaly phenotype and reduction in survival that we observed previously in stable knockout lines, though milder in the acute multi-sgRNA knockout. Finally, we quantify the degree of mutagenesis by deep sequencing, and provide a mathematical model to quantitate the chance for a biallelic loss-of-function mutation. Our findings can be generalized to acute and stable CRISPR/Cas targeting for any zebrafish gene of interest.

RNA tales – how embryos read and discard messages from mom

Following fertilization, embryos develop for a substantial amount of time with a transcriptionally silent genome. Thus, early development is maternally programmed, as it solely relies on RNAs and proteins that are provided by the female gamete. However, these maternal instructions are not sufficient to support later steps of embryogenesis and are therefore gradually replaced by novel products synthesized from the zygotic genome. This switch in the origin of molecular players that drive early development is known as the maternal-to-zygotic transition (MZT). MZT is a universal phenomenon among all metazoans and comprises two interconnected processes: maternal mRNA degradation and the transcriptional awakening of the zygotic genome. The recent adaptation of high-throughput methods for use in embryos has deepened our knowledge of the molecular principles underlying MZT. These mechanisms comprise conserved strategies for RNA regulation that operate in many well-studied cellular contexts but that have adapted differently to early development. In this Review, we will discuss advances in our understanding of post-transcriptional regulatory pathways that drive maternal mRNA clearance during MZT, with an emphasis on recent data in zebrafish embryos on codon-mediated mRNA decay, the contributions of microRNAs (miRNAs) and RNA-binding proteins to this process, and the roles of RNA modifications in the stability control of maternal mRNAs.

Quantitating translational control: mRNA abundance-dependent and independent contributions and the mRNA sequences that specify them

Translation rate per mRNA molecule correlates positively with mRNA abundance. As a result, protein levels do not scale linearly with mRNA levels, but instead scale with the abundance of mRNA raised to the power of an ‘amplification exponent’. Here we show that to quantitate translational control, the translation rate must be decomposed into two components. One, TR_mD, depends on the mRNA level and defines the amplification exponent. The other, TR_mIND, is independent of mRNA amount and impacts the correlation coefficient between protein and mRNA levels. We show that in Saccharomyces cerevisiae TR_mD represents ∼20% of the variance in translation and directs an amplification exponent of 1.20 with a 95% confidence interval [1.14, 1.26]. TR_mIND constitutes the remaining ∼80% of the variance in translation and explains ∼5% of the variance in protein expression. We also find that TR_mD and TR_mIND are preferentially determined by different mRNA sequence features: TR_mIND by the length of the open reading frame and TR_mD both by a ∼60 nucleotide element that spans the initiating AUG and by codon and amino acid frequency. Our work provides more appropriate estimates of translational control and implies that TR_mIND is under different evolutionary selective pressures than TR_mD.

A Massively Parallel Reporter Assay of 3' UTR Sequences Identifies In Vivo Rules for mRNA Degradation

The stability of mRNAs is regulated by signals within their sequences, but a systematic and predictive understanding of the underlying sequence rules remains elusive. Here we introduce UTR-seq, a combination of massively parallel reporter assays and regression models, to survey the dynamics of tens of thousands of 3' UTR sequences during early zebrafish embryogenesis. UTR-seq revealed two temporal degradation programs: a maternally encoded early-onset program and a late-onset program that accelerated degradation after zygotic genome activation. Three signals regulated early-onset rates: stabilizing poly-U and UUAG sequences and destabilizing GC-rich signals. Three signals explained late-onset degradation: miR-430 seeds, AU-rich sequences, and Pumilio recognition sites. Sequence-based regression models translated 3' UTRs into their unique decay patterns and predicted the in vivo effect of sequence signals on mRNA stability. Their application led to the successful design of artificial 3' UTRs that conferred specific mRNA dynamics. UTR-seq provides a general strategy to uncover the rules of RNA cis regulation.

Short poly(A) tails are a conserved feature of highly expressed genes

Poly(A) tails are important elements in mRNA translation and stability. However, recent genome-wide studies concluded that poly(A) tail length was generally not associated with translational efficiency in non-embryonic cells. To investigate if poly(A) tail size might be coupled to gene expression in an intact organism, we used an adapted TAIL-seq protocol to measure poly(A) tails in Caenorhabditis elegans. Surprisingly, we found that well-expressed transcripts contain relatively short, well-defined tails. This attribute appears dependent on translational efficiency, as transcripts enriched for optimal codons and ribosome association had the shortest tail sizes, while non-coding RNAs retained long tails. Across eukaryotes, short tails were a feature of abundant and well-translated mRNAs. Although this seems to contradict the dogma that deadenylation induces translational inhibition and mRNA decay, it instead suggests that well-expressed mRNAs accumulate with pruned tails that accommodate a minimal number of poly(A) binding proteins, which may be ideal for protective and translational functions.

Translation efficiency is a determinant of the magnitude of miRNA-mediated repression

MicroRNAs are well known regulators of mRNA stability and translation. However, the magnitude of both translational repression and mRNA decay induced by miRNA binding varies greatly between miRNA targets. This can be the result of cis and trans factors that affect miRNA binding or action. We set out to address this issue by studying how various mRNA characteristics affect miRNA-mediated repression. Using a dual luciferase reporter system, we systematically analyzed the ability of selected mRNA elements to modulate miRNA-mediated repression. We found that changing the 3'UTR of a miRNA-targeted reporter modulates translational repression by affecting the translation efficiency. This 3'UTR dependent modulation can be further altered by changing the codon-optimality or 5'UTR of the luciferase reporter. We observed maximal repression with intermediate codon optimality and weak repression with very high or low codon optimality. Analysis of ribosome profiling and RNA-seq data for endogenous miRNA targets revealed translation efficiency as a key determinant of the magnitude of miRNA-mediated translational repression. Messages with high translation efficiency were more robustly repressed. Together our results reveal modulation of miRNA-mediated repression by characteristics and features of the 5'UTR, CDS and 3'UTR.

Cis-regulatory elements explain most of the mRNA stability variation across genes in yeast

The stability of mRNA is one of the major determinants of gene expression. Although a wealth of sequence elements regulating mRNA stability has been described, their quantitative contributions to half-life are unknown. Here, we built a quantitative model for Saccharomyces cerevisiae based on functional mRNA sequence features that explains 59% of the half-life variation between genes and predicts half-life at a median relative error of 30%. The model revealed a new destabilizing 3' UTR motif, ATATTC, which we functionally validated. Codon usage proves to be the major determinant of mRNA stability. Nonetheless, single-nucleotide variations have the largest effect when occurring on 3' UTR motifs or upstream AUGs. Analyzing mRNA half-life data of 34 knockout strains showed that the effect of codon usage not only requires functional decapping and deadenylation, but also the 5'-to-3' exonuclease Xrn1, the nonsense-mediated decay genes, but not no-go decay. Altogether, this study quantitatively delineates the contributions of mRNA sequence features on stability in yeast, reveals their functional dependencies on degradation pathways, and allows accurate prediction of half-life from mRNA sequence.

Codon usage regulates protein structure and function by affecting translation elongation speed in Drosophila cells

Codon usage biases are found in all eukaryotic and prokaryotic genomes and have been proposed to regulate different aspects of translation process. Codon optimality has been shown to regulate translation elongation speed in fungal systems, but its effect on translation elongation speed in animal systems is not clear. In this study, we used a Drosophila cell-free translation system to directly compare the velocity of mRNA translation elongation. Our results demonstrate that optimal synonymous codons speed up translation elongation while non-optimal codons slow down translation. In addition, codon usage regulates ribosome movement and stalling on mRNA during translation. Finally, we show that codon usage affects protein structure and function in vitro and in Drosophila cells. Together, these results suggest that the effect of codon usage on translation elongation speed is a conserved mechanism from fungi to animals that can affect protein folding in eukaryotic organisms.

Codon optimality, bias and usage in translation and mRNA decay

The advent of ribosome profiling and other tools to probe mRNA translation has revealed that codon bias - the uneven use of synonymous codons in the transcriptome - serves as a secondary genetic code: a code that guides the efficiency of protein production, the fidelity of translation and the metabolism of mRNAs. Recent advancements in our understanding of mRNA decay have revealed a tight coupling between ribosome dynamics and the stability of mRNA transcripts; this coupling integrates codon bias into the concept of codon optimality, or the effects that specific codons and tRNA concentrations have on the efficiency and fidelity of the translation machinery. In this Review, we first discuss the evidence for codon-dependent effects on translation, beginning with the basic mechanisms through which translation perturbation can affect translation efficiency, protein folding and transcript stability. We then discuss how codon effects are leveraged by the cell to tailor the proteome to maintain homeostasis, execute specific gene expression programmes of growth or differentiation and optimize the efficiency of protein production.

Determinants and Regulation of Protein Turnover in Yeast

Protein turnover maintains the recycling needs of the proteome, and its malfunction has been linked to aging and age-related diseases. However, not all proteins turnover equally, and the factors that contribute to accelerate or slow down turnover are mostly unknown. We measured turnover rates for 3,160 proteins in exponentially growing yeast and analyzed their dependence on physical, functional, and genetic properties. We found that functional characteristics, including protein localization, complex membership, and connectivity, have greater effect on turnover than sequence elements. We also found that protein turnover and mRNA turnover are correlated. Analysis under nutrient perturbation and osmotic stress revealed that protein turnover highly depends on cellular state and is faster when proteins are being actively used. Finally, stress-induced changes in protein and transcript abundance correlated with changes in protein turnover. This study provides a resource of protein turnover rates and principles to understand the recycling needs of the proteome under basal conditions and perturbation.

When LIN41 Comes to a Fork in the Road, It Takes BOTH Paths: Translational Repression OR mRNA Decay, Depending on the Target Site Position

In this issue, Aeschimann et al. (2017) demonstrate that, depending on the target location site (5'UTR or 3'UTR), LIN41 triggers repression of translation or mRNA decay, suggesting that one factor may use two independent pathways of post-transcriptional gene regulation.

RNA Methylation Clears the Way

During the maternal-to-zygotic transition, maternal mRNAs are cleared by multiple distinct but interrelated pathways. A recent study in Nature by Zhao et al. (2017) finds that YTHDF2, a reader of N⁶- methylation, facilitates maternal mRNA decay, introducing an additional facet of control over transcript fate and developmental reprogramming.

Pnrc2 regulates 3'UTR-mediated decay of segmentation clock-associated transcripts during zebrafish segmentation

Vertebrate segmentation is controlled by the segmentation clock, a molecular oscillator that regulates gene expression and cycles rapidly. The expression of many genes oscillates during segmentation, including hairy/Enhancer of split-related (her or Hes) genes, which encode transcriptional repressors that auto-inhibit their own expression, and deltaC (dlc), which encodes a Notch ligand. We previously identified the tortuga (tor) locus in a zebrafish forward genetic screen for genes involved in cyclic transcript regulation and showed that cyclic transcripts accumulate post-splicing in tor mutants. Here we show that cyclic mRNA accumulation in tor mutants is due to loss of pnrc2, which encodes a proline-rich nuclear receptor co-activator implicated in mRNA decay. Using an inducible in vivo reporter system to analyze transcript stability, we find that the her1 3'UTR confers Pnrc2-dependent instability to a heterologous transcript. her1 mRNA decay is Dicer-independent and likely employs a Pnrc2-Upf1-containing mRNA decay complex. Surprisingly, despite accumulation of cyclic transcripts in pnrc2-deficient embryos, we find that cyclic protein is expressed normally. Overall, we show that Pnrc2 promotes 3'UTR-mediated decay of developmentally-regulated segmentation clock transcripts and we uncover an additional post-transcriptional regulatory layer that ensures oscillatory protein expression in the absence of cyclic mRNA decay.

Pervasive yet nonuniform contributions of Dcp2 and Cnot7 to maternal mRNA clearance in zebrafish

mRNA degradation is a fundamental biological process that erases transcribed genetic information from cells. During maternal-to-zygotic transition of animal development, thousands of maternal mRNAs are degraded by multiple mechanisms including microRNAs and codon-mediated decay. Enzymatic requirements for maternal mRNA clearance, however, are not fully understood. Here, we analyzed a contribution of the decapping enzyme Dcp2 to maternal mRNA clearance in zebrafish by over-expressing catalytically inactive Dcp2 and performing RNA-seq analysis. As expected, Dcp2 had a widespread role in maternal mRNA clearance. Interestingly, each mRNA showed differential dependency on Dcp2-mediated decapping and Cnot7-mediated deadenylation for degradation. Correlation analysis identified several mRNA features that were associated with the observed differential dependency. Our results show pervasive yet nonuniform contributions of the decapping enzyme Dcp2 and the deadenylase Cnot7 to maternal mRNA clearance.

The link between adjacent codon pairs and mRNA stability

Background

Evidence in diverse organisms suggests that codon optimality is a major determinant of mRNA translation and degradation. Codon optimality is thought to act by modulating the efficiency of ribosome elongation. In Saccharomyces cerevisiae, a recent study has identified 17 adjacent codon pairs that mediate strong inhibition of translation elongation. However, relationships between the inhibitory codon pairs and other aspects of gene expression are unknown.

Results

To gain insights into how the inhibitory codon pairs may affect aspects of gene expression, we utilized existing datasets to conduct genome-scale analyses in S. cerevisiae. Our analysis revealed the following points. First, the inhibitory codon pairs are significantly associated with faster mRNA decay. The association is not solely due to the content of nucleotides, individual codons, or dipeptides encoded by the inhibitory codon pairs. Second, the inhibitory codon pairs cannot fully explain the previously known relationship of codon optimality with mRNA stability, suggesting that optimality of individual codons and properties of adjacent codon pairs both contribute to gene regulation. Finally, although the inhibitory codon pairs are associated with slower mRNA synthesis and protein instability, the associations can be attributed to usage bias in individual codons.

Conclusions

This study suggests an association of inhibitory codon pairs with mRNA stability and thus another layer of complexity in the codon-mediated gene regulation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3749-8) contains supplementary material, which is available to authorized users.

When mRNA translation meets decay

Messenger RNA (mRNA) translation and mRNA degradation are important determinants of protein output, and they are interconnected. Previously, it was thought that translation of an mRNA, as a rule, prevents its degradation. mRNA surveillance mechanisms, which degrade mRNAs as a consequence of their translation, were considered to be exceptions to this rule. Recently, however, it has become clear that many mRNAs are degraded co-translationally, and it has emerged that codon choice, by influencing the rate of ribosome elongation, affects the rate of mRNA decay. In this review, we discuss the links between translation and mRNA stability, with an emphasis on emerging data suggesting that codon optimality may regulate mRNA degradation.

Dynamic RNA-protein interactions underlie the zebrafish maternal-to-zygotic transition

During the maternal-to-zygotic transition (MZT), transcriptionally silent embryos rely on post-transcriptional regulation of maternal mRNAs until zygotic genome activation (ZGA). RNA-binding proteins (RBPs) are important regulators of post-transcriptional RNA processing events, yet their identities and functions during developmental transitions in vertebrates remain largely unexplored. Using mRNA interactome capture, we identified 227 RBPs in zebrafish embryos before and during ZGA, hereby named the zebrafish MZT mRNA-bound proteome. This protein constellation consists of many conserved RBPs, some of which are potential stage-specific mRNA interactors that likely reflect the dynamics of RNA-protein interactions during MZT. The enrichment of numerous splicing factors like hnRNP proteins before ZGA was surprising, because maternal mRNAs were found to be fully spliced. To address potentially unique roles of these RBPs in embryogenesis, we focused on Hnrnpa1. iCLIP and subsequent mRNA reporter assays revealed a function for Hnrnpa1 in the regulation of poly(A) tail length and translation of maternal mRNAs through sequence-specific association with 3' UTRs before ZGA. Comparison of iCLIP data from two developmental stages revealed that Hnrnpa1 dissociates from maternal mRNAs at ZGA and instead regulates the nuclear processing of pri-mir-430 transcripts, which we validated experimentally. The shift from cytoplasmic to nuclear RNA targets was accompanied by a dramatic translocation of Hnrnpa1 and other pre-mRNA splicing factors to the nucleus in a transcription-dependent manner. Thus, our study identifies global changes in RNA-protein interactions during vertebrate MZT and shows that Hnrnpa1 RNA-binding activities are spatially and temporally coordinated to regulate RNA metabolism during early development.

Tubulins and brain development - The origins of functional specification

The development of the vertebrate central nervous system is reliant on a complex cascade of biological processes that include mitotic division, relocation of migrating neurons, and the extension of dendritic and axonal processes. Each of these cellular events requires the diverse functional repertoire of the microtubule cytoskeleton for the generation of forces, assembly of macromolecular complexes and transport of molecules and organelles. The tubulins are a multi-gene family that encode for the constituents of microtubules, and have been implicated in a spectrum of neurological disorders. Evidence is building that different tubulins tune the functional properties of the microtubule cytoskeleton dependent on the cell type, developmental profile and subcellular localisation. Here we review of the origins of the functional specification of the tubulin gene family in the developing brain at a transcriptional, translational, and post-transcriptional level. We remind the reader that tubulins are not just loading controls for your average Western blot.

Synonymous Codons: Choose Wisely for Expression

The genetic code, which defines the amino acid sequence of a protein, also contains information that influences the rate and efficiency of translation. Neither the mechanisms nor functions of codon-mediated regulation were well understood. The prevailing model was that the slow translation of codons decoded by rare tRNAs reduces efficiency. Recent genome-wide analyses have clarified several issues. Specific codons and codon combinations modulate ribosome speed and facilitate protein folding. However, tRNA availability is not the sole determinant of rate; rather, interactions between adjacent codons and wobble base pairing are key. One mechanism linking translation efficiency and codon use is that slower decoding is coupled to reduced mRNA stability. Changes in tRNA supply mediate biological regulationfor instance,, changes in tRNA amounts facilitate cancer metastasis.

Dynamics of miRNA transcriptome during gonadal development of zebrafish

Studies in non-teleost vertebrates have found microRNAs (miRNAs) to be essential for proper gonadal development. However, comparatively little is known about their role during gonadal development in teleost fishes. So far in zebrafish, a model teleost, transcript profiling throughout gonadal development has not been established because of a tiny size of an organ in juvenile stages and its poor distinguishability from surrounding tissues. We performed small RNA sequencing on isolated gonads of See-Thru-Gonad line, from the undifferentiated state at 3 weeks post fertilization (wpf) to fully mature adults at 24 wpf. We identified 520 gonadal mature miRNAs; 111 of them had significant changes in abundance over time, while 50 miRNAs were either testis- or ovary-enriched significantly in at least one developmental stage. We characterized patterns of miRNA abundance over time including isomiR variants. We identified putative germline versus gonadal somatic miRNAs through differential small RNA sequencing of isolated gametes versus the whole gonads. This report is the most comprehensive analysis of the miRNA repertoire in zebrafish gonads during the sexual development to date and provides an important database from which functional studies can be performed.

m6A-dependent maternal mRNA clearance facilitates zebrafish maternal-to-zygotic transition

The maternal-to-zygotic transition (MZT) is one of the most profound and orchestrated processes during the early life of embryos, yet factors that shape the temporal pattern of vertebrate MZT are largely unknown. Here we show that over one-third of zebrafish maternal messenger RNAs (mRNAs) can be N⁶-methyladenosine (m⁶A) modified, and the clearance of these maternal mRNAs is facilitated by an m⁶A-binding protein, Ythdf2. Removal of Ythdf2 in zebrafish embryos decelerates the decay of m⁶A-modified maternal mRNAs and impedes zygotic genome activation. These embryos fail to initiate timely MZT, undergo cell cycle pause, and remain developmentally delayed throughout larval life. Our study reveals m⁶A-dependent RNA decay as a previously unidentified maternally driven mechanism that regulates maternal mRNA clearance during zebrafish MZT, highlighting the critical role of m⁶A mRNA methylation in transcriptome switching and animal development.

Expression of the novel maternal centrosome assembly factor Wdr8 is required for vertebrate embryonic mitoses

The assembly of the first centrosome occurs upon fertilisation when male centrioles recruit pericentriolar material (PCM) from the egg cytoplasm. The mechanisms underlying the proper assembly of centrosomes during early embryogenesis remain obscure. We identify Wdr8 as a novel maternally essential protein that is required for centrosome assembly during embryonic mitoses of medaka (Oryzias latipes). By CRISPR-Cas9-mediated knockout, maternal/zygotic Wdr8-null (m/zWdr8^-/-) blastomeres exhibit severe defects in centrosome structure that lead to asymmetric division, multipolar mitotic spindles and chromosome alignment errors. Via its WD40 domains, Wdr8 interacts with the centriolar satellite protein SSX2IP. Combining targeted gene knockout and in vivo reconstitution of the maternally essential Wdr8-SSX2IP complex reveals an essential link between maternal centrosome proteins and the stability of the zygotic genome for accurate vertebrate embryogenesis. Our approach provides a way of distinguishing maternal from paternal effects in early embryos and should contribute to understanding molecular defects in human infertility.

RESA identifies mRNA-regulatory sequences at high resolution

Gene expression is regulated extensively at the level of mRNA stability, localization, and translation. However, decoding functional RNA regulatory features remains a limitation to understanding post-transcriptional regulation in vivo. Here, we developed RNA Element Selection Assay (RESA), a method that selects RNA elements based on their activity in vivo and uses high-throughput sequencing to provide quantitative measurement of their regulatory function with near nucleotide resolution. We implemented RESA to identify sequence elements modulating mRNA stability during zebrafish embryogenesis. RESA provides a sensitive and quantitative measure of microRNA activity in vivo and also identifies novel regulatory sequences. To uncover specific sequence requirements within regulatory elements, we developed a bisulfite-mediated nucleotide conversion strategy for large-scale mutational analysis (RESA-bisulfite). Finally, we used the versatile RESA platform to map candidate protein-RNA interactions in vivo (RESA-CLIP). The RESA platform can be broadly applicable to uncover the regulatory features shaping gene expression and cellular function.

Translation repression by maternal RNA binding protein Zar1 is essential for early oogenesis in zebrafish

A large amount of maternal RNA is deposited in oocytes and is reserved for later development. Control of maternal RNA translation during oocyte maturation has been extensively investigated and its regulatory mechanisms are well documented. However, translational regulation of maternal RNA in early oogenesis is largely unexplored. In this study, we generated zebrafish zar1 mutants that result in early oocyte apoptosis and fully penetrant male development. Loss of p53 suppresses the apoptosis in zar1 mutants and restores oocyte development. zar1 immature ovaries show upregulation of proteins implicated in endoplasmic reticulum (ER) stress and the unfolded protein response (UPR). More importantly, loss of Zar1 causes marked upregulation of zona pellucida (ZP) family proteins, while overexpression of ZP proteins in oocytes causes upregulation of stress-related activating transcription factor 3 (atf3), arguing that tightly controlled translation of ZP proteins is essential for ER homeostasis during early oogenesis. Furthermore, Zar1 binds to ZP gene mRNAs and represses their translation. Together, our results indicate that regulation of translational repression and de-repression are essential for precisely controlling protein expression during early oogenesis.

Activation of transcription enforces the formation of distinct nuclear bodies in zebrafish embryos

Nuclear bodies are cellular compartments that lack lipid bilayers and harbor specific RNAs and proteins. Recent proposals that nuclear bodies form through liquid-liquid phase separation leave the question of how different nuclear bodies maintain their distinct identities unanswered. Here we investigate Cajal bodies (CBs), histone locus bodies (HLBs) and nucleoli - involved in assembly of the splicing machinery, histone mRNA 3' end processing, and rRNA processing, respectively - in the embryos of the zebrafish, Danio rerio. We take advantage of the transcriptional silence of the 1-cell embryo and follow nuclear body appearance as zygotic transcription becomes activated. CBs are present from fertilization onwards, while HLB and nucleolar components formed foci several hours later when histone genes and rDNA became active. HLB formation was blocked by transcription inhibition, suggesting nascent histone transcripts recruit HLB components like U7 snRNP. Surprisingly, we found that U7 base-pairing with nascent histone transcripts was not required for localization to HLBs. Rather, the type of Sm ring assembled on U7 determined its targeting to HLBs or CBs; the spliceosomal Sm ring targeted snRNAs to CBs while the specialized U7 Sm-ring localized to HLBs, demonstrating the contribution of protein constituents to the distinction among nuclear bodies. Thus, nucleolar, HLB, and CB components can mix in early embryogenesis when transcription is naturally or artificially silenced. These data support a model in which transcription of specific gene loci nucleates nuclear body components with high specificity and fidelity to perform distinct regulatory functions.

Analysis of the association between codon optimality and mRNA stability in Schizosaccharomyces pombe

Background

Recent experiments have shown that codon optimality is a major determinant of mRNA stability in Saccharomyces cerevisiae and that this phenomenon may be conserved in Escherichia coli and some metazoans, although work in Neurospora crassa is not consistent with this model.

Results

We examined the association between codon optimality and mRNA stability in the fission yeast Schizosaccharomyces pombe. Our analysis revealed the following points. First, we observe a genome-wide association between codon optimality and mRNA stability also in S. pombe, suggesting evolutionary conservation of the phenomenon. Second, in both S. pombe and S. cerevisiae, mRNA synthesis rates are also correlated at the genome-wide analysis with codon optimality, suggesting that the long-appreciated association between codon optimality and mRNA abundance is due to regulation of both mRNA synthesis and degradation. However, when we examined correlation of codon optimality and either mRNA half-lives or synthesis rates controlling for mRNA abundance, codon optimality was still positively correlated with mRNA half-lives in S. cerevisiae, but the association was no longer significant for mRNA half-lives in S. pombe or for synthesis rates in either organism. This illustrates how only the pairwise analysis of multiple correlating variables may limit these types of analyses. Finally, in S. pombe, codon optimality is associated with known DNA/RNA sequence motifs that are associated with mRNA production/stability, suggesting these two features have been under similar selective pressures for optimal gene expression.

Conclusions

Consistent with the emerging body of studies, this study suggests that the association between codon optimality and mRNA stability may be a broadly conserved phenomenon. It also suggests that the association can be explained at least in part by independent adaptations of codon optimality and other transcript features for elevated expression during evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3237-6) contains supplementary material, which is available to authorized users.

Codon optimality and mRNA decay

Recent evidence indicates that codon optimality is a broad determinant of mRNA stability. A study by Radhakrishnan et al. in Cell raises the possibility that the conserved DEAD-box protein Dhh1 underlies the phenomenon.

A Universal Code for mRNA Stability?

Precise elimination of maternal mRNAs plays a critical role during the maternal-to-zygotic transition (MZT) to promote developmental processing. Two new studies demonstrate that, in eukaryotes, codon-mediated decay is a conserved mechanism to shape maternal mRNA stability by affecting deadenylation rate in a translation-dependent manner. These studies add to a growing body of literature suggesting that translational elongation rates are a major determinant of mRNA stability.

Codon usage is an important determinant of gene expression levels largely through its effects on transcription

Codon usage biases are found in all eukaryotic and prokaryotic genomes, and preferred codons are more frequently used in highly expressed genes. The effects of codon usage on gene expression were previously thought to be mainly mediated by its impacts on translation. Here, we show that codon usage strongly correlates with both protein and mRNA levels genome-wide in the filamentous fungus Neurospora Gene codon optimization also results in strong up-regulation of protein and RNA levels, suggesting that codon usage is an important determinant of gene expression. Surprisingly, we found that the impact of codon usage on gene expression results mainly from effects on transcription and is largely independent of mRNA translation and mRNA stability. Furthermore, we show that histone H3 lysine 9 trimethylation is one of the mechanisms responsible for the codon usage-mediated transcriptional silencing of some genes with nonoptimal codons. Together, these results uncovered an unexpected important role of codon usage in ORF sequences in determining transcription levels and suggest that codon biases are an adaptation of protein coding sequences to both transcription and translation machineries. Therefore, synonymous codons not only specify protein sequences and translation dynamics, but also help determine gene expression levels.

The DEAD-Box Protein Dhh1p Couples mRNA Decay and Translation by Monitoring Codon Optimality

A major determinant of mRNA half-life is the codon-dependent rate of translational elongation. How the processes of translational elongation and mRNA decay communicate is unclear. Here, we establish that the DEAD-box protein Dhh1p is a sensor of codon optimality that targets an mRNA for decay. First, we find mRNAs whose translation elongation rate is slowed by inclusion of non-optimal codons are specifically degraded in a Dhh1p-dependent manner. Biochemical experiments show Dhh1p is preferentially associated with mRNAs with suboptimal codon choice. We find these effects on mRNA decay are sensitive to the number of slow-moving ribosomes on an mRNA. Moreover, we find Dhh1p overexpression leads to the accumulation of ribosomes specifically on mRNAs (and even codons) of low codon optimality. Lastly, Dhh1p physically interacts with ribosomes in vivo. Together, these data argue that Dhh1p is a sensor for ribosome speed, targeting an mRNA for repression and subsequent decay.

Codon identity regulates mRNA stability and translation efficiency during the maternal-to-zygotic transition

Cellular transitions require dramatic changes in gene expression that are supported by regulated mRNA decay and new transcription. The maternal-to-zygotic transition is a conserved developmental progression during which thousands of maternal mRNAs are cleared by post-transcriptional mechanisms. Although some maternal mRNAs are targeted for degradation by microRNAs, this pathway does not fully explain mRNA clearance. We investigated how codon identity and translation affect mRNA stability during development and homeostasis. We show that the codon triplet contains translation-dependent regulatory information that influences transcript decay. Codon composition shapes maternal mRNA clearance during the maternal-to-zygotic transition in zebrafish, Xenopus, mouse, and Drosophila, and gene expression during homeostasis across human tissues. Some synonymous codons show consistent stabilizing or destabilizing effects, suggesting that amino acid composition influences mRNA stability. Codon composition affects both polyadenylation status and translation efficiency. Thus, the ribosome interprets two codes within the mRNA: the genetic code which specifies the amino acid sequence and a conserved "codon optimality code" that shapes mRNA stability and translation efficiency across vertebrates.