Click here to close
Hello! We notice that you are using Internet Explorer, which is not supported by Xenbase and may cause the site to display incorrectly.
We suggest using a current version of Chrome,
FireFox, or Safari.
Genome Biol
2021 Jan 05;221:14. doi: 10.1186/s13059-020-02251-5.
Show Gene links
Show Anatomy links
Crosstalk between codon optimality and cis-regulatory elements dictates mRNA stability.
Medina-Muñoz SG
,
Kushawah G
,
Castellano LA
,
Diez M
,
DeVore ML
,
Salazar MJB
,
Bazzini AA
.
???displayArticle.abstract???
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.
Fig. 1. Codon composition predicts mRNA stability in vertebrates. a Scheme of the procedure to train a predictive model of mRNA stability. For each endogenous mRNA, the codon frequencies and the 3′UTR length are used as predictors to train a lasso regression model [28]. The scatter plot shows the point density of predicted and observed mRNA stability (test set genes n = 7576, Pearson correlation test). b Scheme of the 1 nucleotide out of frame reporters (optimal and non-optimal): two mRNAs that differ in the codon composition due to a single nucleotide deletion (G in red, highlight with green) which creates a frameshift. The encoding mCherry fluorescent protein was followed by a cis-acting hydrolase element (P2A) and then by a coding region enriched in optimal or non-optimal codons due to the frameshift. P2A causes ribosome skipping; therefore, the mCherry is not fused to the optimal or non-optimal encoded proteins. The mRNA reporter pairs were co-injected with mRNA encoding for GFP as an internal control [24]. c Fluorescence microscopy images of representative embryos at 8 h post-injection (hpi) with the indicated 1 nt out of frame reporter and GFP. Box plot displays fluorescence quantification at 8 hpi with each reporter. The mCherry fluorescence intensity was normalized to GFP intensity in each embryo (p = 0.041, paired t test). d mRNA stability predictions for 1 nucleotide out of frame reporters in fish and human cells. In all cases, the prediction for the optimal reporter is higher than that for the non-optimal (p = 0.007, binomial test)
Fig. 2. Codon optimality is the major determinant of mRNA stability. a Diagram depicting the maternal to zygotic transition in zebrafish. b Scatter plots of predicted and observed mRNA stability, during MZT, for maternal mRNAs in zebrafish and Xenopus. The gradient of color represents the content of optimal codons [22]. The predicted mRNA stability correlates with the proportion of optimal codons (p < 2 × 10−16, Pearson correlation). The mRNA stability median of mRNAs enriched in optimal (red) or non-optimal (blue) codons, as well as mRNAs with miR-430/-427 target sites in the 3′UTR (green) are shown. c Codon content explains most of the mRNA decay during MZT. The x-axis shows the fraction of coding genes that can be regulated by different mRNA stability pathways. For microRNAs, this fraction corresponds to the number of seed sites (GCACTT) in the 3′UTR, and for m6A, the fraction is the number of target genes reported [10]. The y-axis shows the Bayesian model comparison weights [32]. These weights represent which model is more likely to predict unobserved data better, higher values indicate stronger regulatory effects
Fig. 3. Dissecting cis-regulatory elements after accounting for codon-mediated regulation. a Diagram describing the residual score (observed − predicted). mRNAs that decay more than expected by the model show negative residuals and potentially contain destabilizing cis-regulatory elements. The mRNAs with positive residual scores might have stabilizing cis-regulatory elements. b, c The model overestimates the mRNA stability of miR-430 targets. Sinaplot showing the distribution of the residual scores for targets and not targets of miR-430. In this type of plot, each dot represents an individual mRNA. The targets of miR-430 are grouped by the type of miR-430 seed (b) (p < 2 × 10−16, ANOVA test) or by the number of target seeds (c) (p < 2 × 10−16, ANOVA test) present in the 3′UTR during the MZT in zebrafish (Additional file 1: Fig. S1a-b) and Xenopus [33]. d The model also overestimates the stability of mRNAs that contain m6A methylation mark. Sinaplot showing that m6A targets [10] display lower residual score distribution than not targets during the MZT in zebrafish (p < 2 × 10−16, t test). e Sylamer landscape plot that tracks occurrence biases of 6-nucleotide words in the 3′UTRs using hypergeometric p values for all words across the mRNAs ranking based on residual values [35]. The highlighted 6-mers are significantly associated with mRNA stabilization (TATCTA, CTATCT, and TCTATC) and destabilization (GCACTT, TAGGAC, and GGACTT). The color shows the putative regulatory pathway that recognizes these 6-mers. The dotted line shows a hypergeometric p = 0.01
Fig. 4. Codon optimality affects mRNA stability and gene expression of microRNA targets and m6A targets in vertebrates. mRNAs were divided into targets (microRNA or m6A) and nontargets. Each group was divided into four equal groups with decreasing levels of optimal codons. All p values were computed with a linear model. a Sinaplot showing the distribution of mRNA stability during MZT (log2-fold change 6 h post-fertilization (hpf)/2 hpf) for targets and not target genes of miR-430/-427 in zebrafish (Additional file 1: Fig. S1a-b) and Xenopus [33] embryos. b Sinaplot showing the distribution of mRNA stability during MZT in zebrafish for methylated (m6A) and non-methylated mRNAs [10]. c Sinaplot showing the mRNA stability distribution of m6A and miR-291a targets [29] and not targets in mouse embryonic stem cells. d Sinaplot showing the RNA-level distribution of miR-1 and miR-155 targets [37] after microRNA transfection in human cells. e Diagram depicting reporter genes containing almost identical nucleotide sequence but different codon composition (enriched in optimal or non-optimal codons) due to a single insertion changing the frame. For each coding region, a 3′UTR sequence containing miR-17 seed sites (7-mer AGCACTT) or a mutant version with two nucleotides mutated disrupting the miR-17 seed sites were cloned. The boxplot shows the distribution of scaled mCherry/GFP intensity for reporters with and without miR-17 seed site in transfected 293T human cells. Both coding sequence and miR-17 affect the reporter expression (p values computed with paired t test)
Fig. 5. Targeting efficacy of miR-430/-427 can be affected by the coding sequence. a Scatter plot comparing codon optimality level and change in expression between wildtype zebrafish embryos vs maternal/zygotic Dicer mutant embryos at 6 h post-fertilization [39]. Only mRNAs with miR-430 seed in the 3′UTR (GCACTT) are shown. The line represents the average change in expression (log2-fold WT/Dicer at 6 hpf) as a function of the codon optimality level (Additional file 1: Table S5). The effect of miR-430 is not constant across different levels of codon optimality (p value was obtained using an F test comparing a model with a non-linear effect on codon optimality vs constant effect). The confidence interval was determined with bootstrap replicates (n = 100) [40]. b Line plot showing the expected decrease in gene expression due to miR-430 during zebrafish MZT (log2-fold change 6 vs 2 hpf) (Additional file 1: Fig. S1a-b). The y-axis represents a measure, estimated from the data, of the miR-430 repressive strength with respect to codon optimality (x-axis). For example, for an mRNA that is very repressed by miR-430, the miR-430 component will be larger (higher negative value in the y-axis). However, for another gene, which has a weak miR-430, the miR-430 component is smaller (closer to 0 in the y-axis). The p value denotes the statistical significance of the non-linear interaction between codon optimality and miR-430 presence (F test) obtained with a generalized additive model [41]. The confidence interval was determined with bootstrap replicates (n = 100) [40]. c Scheme of the reporter library which includes random fragments of the zebrafish transcriptome [22]. Transcripts share the same 5′ and 3′UTR but some sequences contain a stop codon in the coding region. These stop codons create a random and longer 3′UTR sequence. mRNAs were injected at the one‐cell stage in zebrafish [22], and the reporter library is analyzed at 2 and 8 hpi using high-throughput sequencing. To analyze the depletion of miR-430 with respect to the codon content of the transcripts, we filtered those sequences that contain a coding region of at least 350 nucleotides and a random 3′UTR length of at least 75 nucleotides. d, e Analysis of miR-430/427 depletion in the reporter library. The reporters were grouped into 7-tiles, equal size, with increasing levels of codon optimality (Additional file 1: Fig. S5c). For each tile, we computed the depletion of the miR-430/427 seed (GCACTT) in the 3′UTR (log2-fold 8 h/2 h for zebrafish and 9h/1 h for Xenopus). Each boxplot is formed by bootstrap replicates (n = 100) [40]
Fig. 6. MicroRNAs antagonize codon optimality effect on mRNA stability during the MZT. a Sinaplot showing the codon optimality distribution (Table S5) for the top 1000 most unstable maternal genes during zebrafish and Xenopus MZT. The stability was defined based on the log2-fold change of early (2 h) vs late time points (fish = 6 h, Xenopus = 9 h) in fish (Additional file 1: Fig. S1a-b) and Xenopus [33]. The genes were divided into groups based on the numbers or seed type of miR-430/427. The content of optimal codons increases with the miR-430 regulation strength (p = 2 × 10−4 zebrafish, p = 2 × 10−3
Xenopus, linear regression). The p value was obtained by comparing the difference in the mean level of optimal codons (PLS1) between genes with and without miR-430/427 sites using a linear model. b Heatmap of miR-430 enrichment in the 3′UTR as a function of mRNA stability and codon optimality level. The miR-430 enrichment was estimated with a generalized linear model. Unstable mRNAs enriched in optimal codons temp to contain miR-430 sites (e.g., smarca2). c Scatter plot comparing the RNA stability, during MZT, of ortholog genes in zebrafish (Additional file 1: Fig. S1a-b) and Xenopus [33]. d Scatter plot comparing the content of optimal codons in zebrafish and Xenopus for ortholog genes [22]. e Sinaplot showing the codon optimality distribution for unstable mRNAs in panel a that are orthologs (n = 280). Messenger RNAs were divided into four categories according to the presence or absence of miR-430/-427 seeds in both species. The orthologous mRNA with miR-430/-427 seeds in both species are the most enriched in optimal codons (p = 0.0868, one-way ANOVA test)
Fig. 7. Model showing that the mRNA stability depends on the regulatory elements of the coding and the 3′UTR, suggesting that to fully understand mRNA stability, the regulatory information across the entire mRNA sequence needs to be integrated, rather than focusing solely on the 3′UTR or in the coding sequence
Aanes,
Zebrafish mRNA sequencing deciphers novelties in transcriptome dynamics during maternal to zygotic transition.
2011, Pubmed
Aanes,
Zebrafish mRNA sequencing deciphers novelties in transcriptome dynamics during maternal to zygotic transition.
2011,
Pubmed
Ambros,
MicroRNAs and developmental timing.
2011,
Pubmed
Bartel,
Metazoan MicroRNAs.
2018,
Pubmed
Bartel,
MicroRNAs: genomics, biogenesis, mechanism, and function.
2004,
Pubmed
Batista,
m(6)A RNA modification controls cell fate transition in mammalian embryonic stem cells.
2014,
Pubmed
Bazzini,
Identification of small ORFs in vertebrates using ribosome footprinting and evolutionary conservation.
2014,
Pubmed
Bazzini,
Ribosome profiling shows that miR-430 reduces translation before causing mRNA decay in zebrafish.
2012,
Pubmed
Bazzini,
Codon identity regulates mRNA stability and translation efficiency during the maternal-to-zygotic transition.
2016,
Pubmed
,
Xenbase
Boël,
Codon influence on protein expression in E. coli correlates with mRNA levels.
2016,
Pubmed
Bushati,
Temporal reciprocity of miRNAs and their targets during the maternal-to-zygotic transition in Drosophila.
2008,
Pubmed
Carneiro,
Codon stabilization coefficient as a metric to gain insights into mRNA stability and codon bias and their relationships with translation.
2019,
Pubmed
Chen,
Global regulation of mRNA translation and stability in the early Drosophila embryo by the Smaug RNA-binding protein.
2014,
Pubmed
Cheng,
Cis-regulatory elements explain most of the mRNA stability variation across genes in yeast.
2017,
Pubmed
Cifuentes,
A novel miRNA processing pathway independent of Dicer requires Argonaute2 catalytic activity.
2010,
Pubmed
Despic,
RNA tales - how embryos read and discard messages from mom.
2018,
Pubmed
Flicek,
Ensembl 2014.
2014,
Pubmed
Giraldez,
Zebrafish MiR-430 promotes deadenylation and clearance of maternal mRNAs.
2006,
Pubmed
Guo,
Mammalian microRNAs predominantly act to decrease target mRNA levels.
2010,
Pubmed
Hanson,
Codon optimality, bias and usage in translation and mRNA decay.
2018,
Pubmed
Herzog,
Thiol-linked alkylation of RNA to assess expression dynamics.
2017,
Pubmed
Ivanova,
The RNA m6A Reader YTHDF2 Is Essential for the Post-transcriptional Regulation of the Maternal Transcriptome and Oocyte Competence.
2017,
Pubmed
Johnstone,
Upstream ORFs are prevalent translational repressors in vertebrates.
2016,
Pubmed
Kushawah,
CRISPR-Cas13d Induces Efficient mRNA Knockdown in Animal Embryos.
2020,
Pubmed
Langmead,
Fast gapped-read alignment with Bowtie 2.
2012,
Pubmed
Laue,
The maternal to zygotic transition regulates genome-wide heterochromatin establishment in the zebrafish embryo.
2019,
Pubmed
Laver,
Regulation and Function of Maternal Gene Products During the Maternal-to-Zygotic Transition in Drosophila.
2015,
Pubmed
Laver,
Brain tumor is a sequence-specific RNA-binding protein that directs maternal mRNA clearance during the Drosophila maternal-to-zygotic transition.
2015,
Pubmed
Laver,
The RNA-Binding Protein Rasputin/G3BP Enhances the Stability and Translation of Its Target mRNAs.
2020,
Pubmed
Lee,
Nanog, Pou5f1 and SoxB1 activate zygotic gene expression during the maternal-to-zygotic transition.
2013,
Pubmed
Li,
The Sequence Alignment/Map format and SAMtools.
2009,
Pubmed
Lichner,
The miR-290-295 cluster promotes pluripotency maintenance by regulating cell cycle phase distribution in mouse embryonic stem cells.
2011,
Pubmed
Lund,
Deadenylation of maternal mRNAs mediated by miR-427 in Xenopus laevis embryos.
2009,
Pubmed
,
Xenbase
Mayr,
Regulation by 3'-Untranslated Regions.
2017,
Pubmed
McGeary,
The biochemical basis of microRNA targeting efficacy.
2019,
Pubmed
Mishima,
Codon Usage and 3' UTR Length Determine Maternal mRNA Stability in Zebrafish.
2016,
Pubmed
Owens,
Measuring Absolute RNA Copy Numbers at High Temporal Resolution Reveals Transcriptome Kinetics in Development.
2016,
Pubmed
,
Xenbase
Pinder,
Smaug: an unexpected journey into the mechanisms of post-transcriptional regulation.
2013,
Pubmed
Presnyak,
Codon optimality is a major determinant of mRNA stability.
2015,
Pubmed
Quinlan,
BEDTools: a flexible suite of utilities for comparing genomic features.
2010,
Pubmed
Rabani,
A Massively Parallel Reporter Assay of 3' UTR Sequences Identifies In Vivo Rules for mRNA Degradation.
2017,
Pubmed
Radhakrishnan,
The DEAD-Box Protein Dhh1p Couples mRNA Decay and Translation by Monitoring Codon Optimality.
2016,
Pubmed
Ross,
mRNA stability in mammalian cells.
1995,
Pubmed
Rutkowski,
Widespread disruption of host transcription termination in HSV-1 infection.
2015,
Pubmed
Schoenberg,
Regulation of cytoplasmic mRNA decay.
2012,
Pubmed
Seok,
MicroRNA Target Recognition: Insights from Transcriptome-Wide Non-Canonical Interactions.
2016,
Pubmed
Sheu-Gruttadauria,
Beyond the seed: structural basis for supplementary microRNA targeting by human Argonaute2.
2019,
Pubmed
Tadros,
The maternal-to-zygotic transition: a play in two acts.
2009,
Pubmed
Vastenhouw,
The maternal-to-zygotic transition revisited.
2019,
Pubmed
Vejnar,
Genome wide analysis of 3' UTR sequence elements and proteins regulating mRNA stability during maternal-to-zygotic transition in zebrafish.
2019,
Pubmed
Wu,
Translation affects mRNA stability in a codon-dependent manner in human cells.
2019,
Pubmed
Wu,
Translation of small downstream ORFs enhances translation of canonical main open reading frames.
2020,
Pubmed
Wu,
Systems to study codon effect on post-transcriptional regulation of gene expression.
2018,
Pubmed
Yang,
RNA 5-Methylcytosine Facilitates the Maternal-to-Zygotic Transition by Preventing Maternal mRNA Decay.
2019,
Pubmed
Yartseva,
RESA identifies mRNA-regulatory sequences at high resolution.
2017,
Pubmed
Zhao,
m6A-dependent maternal mRNA clearance facilitates zebrafish maternal-to-zygotic transition.
2017,
Pubmed
van Dongen,
Detecting microRNA binding and siRNA off-target effects from expression data.
2008,
Pubmed