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

Global regulation via modulation of ribosome pausing by the ABC-F protein EttA

Having multiple rounds of translation of the same mRNA creates dynamic complexities along with opportunities for regulation related to ribosome pausing and stalling at specific sequences. Yet, mechanisms controlling these critical processes and the principles guiding their evolution remain poorly understood. Through genetic, genomic, physiological, and biochemical approaches, we demonstrate that regulating ribosome pausing at specific amino acid sequences can produce ~2-fold changes in protein expression levels which strongly influence cell growth and therefore evolutionary fitness. We demonstrate, both in vivo and in vitro, that the ABC-F protein EttA directly controls the translation of mRNAs coding for a subset of enzymes in the tricarboxylic acid (TCA) cycle and its glyoxylate shunt, which modulates growth in some chemical environments. EttA also modulates expression of specific proteins involved in metabolically related physiological and stress-response pathways. These regulatory activities are mediated by EttA rescuing ribosomes paused at specific patterns of negatively charged residues within the first 30 amino acids of nascent proteins. We thus establish a unique global regulatory paradigm based on sequence-specific modulation of translational pausing.

Hypoxia and inflammation induce synergistic transcriptome turnover in macrophages

Macrophages are effector immune cells that experience substantial changes to oxygenation when transiting through tissues, especially when entering tumors or infected wounds. How hypoxia alters gene expression and macrophage effector function at the post-transcriptional level remains poorly understood. Here, we use TimeLapse-seq to measure how inflammatory activation modifies the hypoxic response in primary macrophages. Nucleoside recoding sequencing allows the derivation of steady-state transcript levels, degradation rates, and transcriptional synthesis rates from the same dataset. We find that hypoxia produces distinct responses from resting and inflammatory macrophages. Hypoxia induces destabilization of mRNA transcripts, though inflammatory macrophages substantially increase mRNA degradation compared to resting macrophages. Increased RNA turnover results in the upregulation of ribosomal protein genes and downregulation of extracellular matrix components in inflammatory macrophages. Pathways regulated by mRNA decay in vitro are differentially regulated in tumor-associated macrophages implying that mixed stimuli could induce post-transcriptional regulation of macrophage function in solid tumors.

Dengue virus preferentially uses human and mosquito non-optimal codons

Codon optimality refers to the effect that codon composition has on messenger RNA (mRNA) stability and translation level and implies that synonymous codons are not silent from a regulatory point of view. Here, we investigated the adaptation of virus genomes to the host optimality code using mosquito-borne dengue virus (DENV) as a model. We demonstrated that codon optimality exists in mosquito cells and showed that DENV preferentially uses nonoptimal (destabilizing) codons and avoids codons that are defined as optimal (stabilizing) in either human or mosquito cells. Human genes enriched in the codons preferentially and frequently used by DENV are upregulated during infection, and so is the tRNA decoding the nonoptimal and DENV preferentially used codon for arginine. We found that adaptation during single-host passaging in human or mosquito cells results in the selection of synonymous mutations towards DENV's preferred nonoptimal codons that increase virus fitness. Finally, our analyses revealed that hundreds of viruses preferentially use nonoptimal codons, with those infecting a single host displaying an even stronger bias, suggesting that host-pathogen interaction shapes virus-synonymous codon choice.

Relationship between circulating thrombospondin-1 messenger ribonucleic acid and microribonucleic acid-194 levels in Chinese patients with type 2 diabetic kidney disease: The outcomes of a case-control study

AIMS/INTRODUCTION

We investigated the relationship of circulating TSP-1 mRNA and miR-194 with diabetic kidney disease's degree.

MATERIALS AND METHODS

We enrolled 167 hospitalized type 2 diabetes patients in the endocrinology department. Patients were split into three groups according to urinary microalbumin: A, B and C. The control group comprised healthy outpatients (n = 163). The quantities of microribonucleic acid (miR)-194 and thrombospondin-1 (TSP-1) messenger ribonucleic acid (mRNA) in the participants' circulation were measured using a quantitative real-time polymerase chain reaction.

RESULTS

Circulating TSP-1 mRNA (P = 0.024) and miR-194 (P = 0.029) expressions significantly increased in type 2 diabetes patients. Circulating TSP-1 mRNA (P = 0.040) and miR-194 (P = 0.007) expression levels differed significantly among the three groups; circulating TSP-1 mRNA expression increased with urinary microalbumin. However, miR-194 declined in group B and increased in group C. Circulating TSP-1 mRNA was positively correlated with cystatin-c (r = 0.281; P = 0.021) and microalbumin/creatinine ratio (UmALB/Cr; r = 0.317; P = 0.009); miR-194 was positively correlated with UmALB/Cr (r = 0.405; P = 0.003). Stepwise multivariate linear regression analysis showed cystatin-c (β = 0.578; P = 0.021) and UmALB/Cr (β = 0.001; P = 0.009) as independent factors for TSP-1 mRNA; UmALB/Cr (β = 0.005; P = 0.028) as an independent factor for miR194. Areas under the curve for circulating TSP-1 mRNA and miR194 were 0.756 (95% confidence interval 0.620-0.893; sensitivity 0.69 and specificity 0.71, P < 0.01) and 0.584 (95% confidence interval 0.421-0.748; sensitivity 0.54 and specificity 0.52, P < 0.01), respectively.

CONCLUSIONS

Circulating TSP-1 mRNA and miR-194 expressions significantly increased in type 2 diabetes patients. The microalbumin group had lower levels of miR-194 (a risk factor that is valuable for type 2 diabetes kidney disease evaluation).

2'-O-methylation at internal sites on mRNA promotes mRNA stability

2'-O-methylation (Nm) is a prominent RNA modification well known in noncoding RNAs and more recently also found at many mRNA internal sites. However, their function and base-resolution stoichiometry remain underexplored. Here, we investigate the transcriptome-wide effect of internal site Nm on mRNA stability. Combining nanopore sequencing with our developed machine learning method, NanoNm, we identify thousands of Nm sites on mRNAs with a single-base resolution. We observe a positive effect of FBL-mediated Nm modification on mRNA stability and expression level. Elevated FBL expression in cancer cells is associated with increased expression levels for 2'-O-methylated mRNAs of cancer pathways, implying the role of FBL in post-transcriptional regulation. Lastly, we find that FBL-mediated 2'-O-methylation connects to widespread 3' UTR shortening, a mechanism that globally increases RNA stability. Collectively, we demonstrate that FBL-mediated Nm modifications at mRNA internal sites regulate gene expression by enhancing mRNA stability.

Blood flow regulates acvrl1 transcription via ligand-dependent Alk1 activity

Hereditary hemorrhagic telangiectasia (HHT) is an autosomal dominant disease characterized by the development of arteriovenous malformations (AVMs) that can result in significant morbidity and mortality. HHT is caused primarily by mutations in bone morphogenetic protein receptors ACVRL1/ALK1, a signaling receptor, or endoglin (ENG), an accessory receptor. Because overexpression of Acvrl1 prevents AVM development in both Acvrl1 and Eng null mice, enhancing ACVRL1 expression may be a promising approach to development of targeted therapies for HHT. Therefore, we sought to understand the molecular mechanism of ACVRL1 regulation. We previously demonstrated in zebrafish embryos that acvrl1 is predominantly expressed in arterial endothelial cells and that expression requires blood flow. Here, we document that flow dependence exhibits regional heterogeneity and that acvrl1 expression is rapidly restored after reinitiation of flow. Furthermore, we find that acvrl1 expression is significantly decreased in mutants that lack the circulating Alk1 ligand, Bmp10, and that, in the absence of flow, intravascular injection of BMP10 or the related ligand, BMP9, restores acvrl1 expression in an Alk1-dependent manner. Using a transgenic acvrl1:egfp reporter line, we find that flow and Bmp10 regulate acvrl1 at the level of transcription. Finally, we observe similar ALK1 ligand-dependent increases in ACVRL1 in human endothelial cells subjected to shear stress. These data suggest that ligand-dependent Alk1 activity acts downstream of blood flow to maintain or enhance acvrl1 expression via a positive feedback mechanism, and that ALK1 activating therapeutics may have dual functionality by increasing both ALK1 signaling flux and ACVRL1 expression.

Comparative Analysis of Maternal Gene Expression Patterns Unravels Evolutionary Signatures Across Reproductive Modes

Maternal genes have a pivotal role in regulating metazoan early development. As such their functions have been extensively studied since the dawn of developmental biology. The temporal and spatial dynamics of their transcripts have been thoroughly described in model organisms and their functions have been undergoing heavy investigations. Yet, less is known about the evolutionary changes shaping their presence within diverse oocytes. Due to their unique maternal inheritance pattern, a high degree is predicted to be present when it comes to their expression. Insofar only limited and conflicting results have emerged around it. Here, we set out to elucidate which evolutionary changes could be detected in the maternal gene expression patterns using phylogenetic comparative methods on RNAseq data from 43 species. Using normalized gene expression values and fold change information throughout early development we set out to find the best-fitting evolutionary model. Through modeling, we find evidence supporting both the high degree of divergence and constraint on gene expression values, together with their temporal dynamics. Furthermore, we find that maternal gene expression alone can be used to explain the reproductive modes of different species. Together, these results suggest a highly dynamic evolutionary landscape of maternal gene expression. We also propose a possible functional dichotomy of maternal genes which is influenced by the reproductive strategy undertaken by examined species.

Alter codon bias of the P. pastoris genome to overcome a bottleneck in codon optimization strategy development and improve protein expression

Apart from its role in translation, codon bias is also an important mechanism to regulate mRNA levels. The traditional frequency-based codon optimization strategy is rather efficient in organisms such as N. crassa, but much less in yeast P. pastoris which is a popular host for heterologous protein expression. This is because that unlike N. crassa, the preferred codons of P. pastoris are actually AU-rich and hence codon optimization for extremely low GC content comes with issues of pre-mature transcriptional termination or low RNA stability in spite of translational advantages. To overcome this bottleneck, we focused on three reporter genes in P. pastoris first and confirmed the great advantage of GC-prone codon optimization on mRNA levels. Then we altered the codon bias profile of P. pastoris by introducing additional rare tRNA gene copies. Prior to that we constructed IPTG-regulated tRNA species to enable chassis cells to switch between different codon bias status. As demonstrated again with reporter genes, protein yield of luc and 0788 was successfully increased by 4-5 folds in chassis cells. In summary, here we provide an alternative codon optimization strategy for genes with unsatisfactory performance under traditional codon frequency-based optimization.

UPF1 regulates mRNA stability by sensing poorly translated coding sequences

Post-transcriptional mRNA regulation shapes gene expression, yet how cis-elements and mRNA translation interface to regulate mRNA stability is poorly understood. We find that the strength of translation initiation, upstream open reading frame (uORF) content, codon optimality, AU-rich elements, microRNA binding sites, and open reading frame (ORF) length function combinatorially to regulate mRNA stability. Machine-learning analysis identifies ORF length as the most important conserved feature regulating mRNA decay. We find that Upf1 binds poorly translated and untranslated ORFs, which are associated with a higher decay rate, including mRNAs with uORFs and those with exposed ORFs after stop codons. Our study emphasizes Upf1's converging role in surveilling mRNAs with exposed ORFs that are poorly translated, such as mRNAs with long ORFs, ORF-like 3' UTRs, and mRNAs containing uORFs. We propose that Upf1 regulation of poorly/untranslated ORFs provides a unifying mechanism of surveillance in regulating mRNA stability and homeostasis in an exon-junction complex (EJC)-independent nonsense-mediated decay (NMD) pathway that we term ORF-mediated decay (OMD).

miR-430 regulates zygotic mRNA during zebrafish embryogenesis

BACKGROUND

Early embryonic developmental programs are guided by the coordinated interplay between maternally inherited and zygotically manufactured RNAs and proteins. Although these processes happen concomitantly and affecting gene function during this period is bound to affect both pools of mRNAs, it has been challenging to study their expression dynamics separately.

RESULTS

By employing SLAM-seq, a nascent mRNA labeling transcriptomic approach, in a developmental time series we observe that over half of the early zebrafish embryo transcriptome consists of maternal-zygotic genes, emphasizing their pivotal role in early embryogenesis. We provide an hourly resolution of de novo transcriptional activation events and follow nascent mRNA trajectories, finding that most de novo transcriptional events are stable throughout this period. Additionally, by blocking microRNA-430 function, a key post transcriptional regulator during zebrafish embryogenesis, we directly show that it destabilizes hundreds of de novo transcribed mRNAs from pure zygotic as well as maternal-zygotic genes. This unveils a novel miR-430 function during embryogenesis, fine-tuning zygotic gene expression.

CONCLUSION

These insights into zebrafish early embryo transcriptome dynamics emphasize the significance of post-transcriptional regulators in zygotic genome activation. The findings pave the way for future investigations into the coordinated interplay between transcriptional and post-transcriptional landscapes required for the establishment of animal cell identities and functions.

Exploiting Systematic Engineering of the Expression Cassette as a Powerful Tool to Enhance Heterologous Gene Expression in Trichoderma reesei

Many endeavors in expressing a heterologous gene in microbial hosts rely on simply placing the gene of interest between a selected pair of promoters and terminator. However, although the expression efficiency could be improved by engineering the host cell, how modifying the expression cassette itself systematically would affect heterologous gene expression remains largely unknown. As the promoter and terminator bear plentiful cis-elements, herein using the Aspergillus niger mannanase with high application value in animal feeds and the eukaryotic filamentous fungus workhorse Trichoderma reesei as a model gene/host, systematic engineering of an expression cassette was investigated to decipher the effect of its mutagenesis on heterologous gene expression. Modifying the promoter, signal peptide, the eukaryotic-specific Kozak sequence, and the 3'-UTR could stepwise improve extracellular mannanase production from 17 U/mL to an ultimate 471 U/mL, representing a 27.7-fold increase in expression. The strategies can be generally applied in improving the production of heterologous proteins in eukaryotic microbial hosts.

Computational design of mRNA vaccines

mRNA technology has emerged as a successful vaccine platform that offered a swift response to the COVID-19 pandemic. Accumulating evidence shows that vaccine efficacy, thermostability, and other important properties, are largely impacted by intrinsic properties of the mRNA molecule, such as RNA sequence and structure, both of which can be optimized. Designing mRNA sequence for vaccines presents a combinatorial problem due to an extremely large selection space. For instance, due to the degeneracy of the genetic code, there are over 10632 possible mRNA sequences that could encode the spike protein, the COVID-19 vaccines' target. Moreover, designing different elements of the mRNA sequence simultaneously against multiple objectives such as translational efficiency, reduced reactogenicity, and improved stability requires an efficient and sophisticated optimization strategy. Recently, there has been a growing interest in utilizing computational tools to redesign mRNA sequences to improve vaccine characteristics and expedite discovery timelines. In this review, we explore important biophysical features of mRNA to be considered for vaccine design and discuss how computational approaches can be applied to rapidly design mRNA sequences with desirable characteristics.

Transcriptomic Profiling at the Maternal-to-Zygotic Transition in Leech, Helobdella austinensis

The glossiphoniid leech, Helobdella austinensis, is an experimentally tractable member of the superphylum, Lophotrochozoa. Its large embryonic cells, stereotyped asymmetric cell divisions and ex vivo development capabilities makes it a favorable model for studying the molecular and cellular events of a representative spiralian. In this study, we focused on a narrow developmental time window of ~6-8 h, comprising stages just prior to and immediately following zygote deposition. Employing RNA-Seq methodology, we identified differentially expressed transcripts at this fundamental ontogenic boundary, known as the maternal-to-zygotic transition (MZT). Gene expression changes were characterized by the massive degradation of maternal RNAs (~45%) coupled with the rapid transcription of ~5000 zygotic genes (~20% of the genome) in the first mitotic cell cycle. The latter transcripts encoded a mixture of cell maintenance and regulatory proteins that predictably influence downstream developmental events.

Decoding the Heterogeneity and Specialized Function of Translation Machinery Through Ribosome Profiling in Yeast Mutants of Initiation Factors

The nuanced heterogeneity and specialized functions of translation machinery are increasingly recognized as crucial for precise translational regulation. Here, high-throughput ribosomal profiling (ribo-seq) is used to analyze the specialized roles of eukaryotic initiation factors (eIFs) in the budding yeast. By examining changes in ribosomal distribution across the genome resulting from knockouts of eIF4A, eIF4B, eIF4G1, CAF20, or EAP1, or knockdowns of eIF1, eIF1A, eIF4E, or PAB1, two distinct initiation-factor groups, the "looping" and "scanning" groups are discerned, based on similarities in the ribosomal landscapes their perturbation induced. The study delves into the cis-regulatory sequence features of genes influenced predominantly by each group, revealing that genes more dependent on the looping-group factors generally have shorter transcripts and poly(A) tails. In contrast, genes more dependent on the scanning-group factors often possess upstream open reading frames and exhibit a higher GC content in their 5' untranslated regions. From the ribosomal RNA fragments identified in the ribo-seq data, ribosomal heterogeneity associated with perturbation of specific initiation factors is further identified, suggesting their potential roles in regulating ribosomal components. Collectively, the study illuminates the complexity of translational regulation driven by heterogeneity and specialized functions of translation machinery, presenting potential approaches for targeted gene translation manipulation.

Blood flow regulates acvrl1 transcription via ligand-dependent Alk1 activity

Hereditary hemorrhagic telangiectasia (HHT) is an autosomal dominant disease characterized by the development of arteriovenous malformations (AVMs) that can result in significant morbidity and mortality. HHT is caused primarily by mutations in bone morphogenetic protein receptors ACVRL1/ALK1, a signaling receptor, or endoglin (ENG), an accessory receptor. Because overexpression of Acvrl1 prevents AVM development in both Acvrl1 and Eng null mice, enhancing ACVRL1 expression may be a promising approach to development of targeted therapies for HHT. Therefore, we sought to understand the molecular mechanism of ACVRL1 regulation. We previously demonstrated in zebrafish embryos that acvrl1 is predominantly expressed in arterial endothelial cells and that expression requires blood flow. Here, we document that flow dependence exhibits regional heterogeneity and that acvrl1 expression is rapidly restored after reinitiation of flow. Furthermore, we find that acvrl1 expression is significantly decreased in mutants that lack the circulating Alk1 ligand, Bmp10, and that BMP10 microinjection into the vasculature in the absence of flow enhances acvrl1 expression in an Alk1-dependent manner. Using a transgenic acvrl1:egfp reporter line, we find that flow and Bmp10 regulate acvrl1 at the level of transcription. Finally, we observe similar ALK1 ligand-dependent increases in ACVRL1 in human endothelial cells subjected to shear stress. These data suggest that Bmp10 acts downstream of blood flow to maintain or enhance acvrl1 expression via a positive feedback mechanism, and that ALK1 activating therapeutics may have dual functionality by increasing both ALK1 signaling flux and ACVRL1 expression.

A spatiotemporally resolved atlas of mRNA decay in the C. elegans embryo reveals differential regulation of mRNA stability across stages and cell types

During embryonic development, cells undergo dynamic changes in gene expression that are required for appropriate cell fate specification. Although both transcription and mRNA degradation contribute to gene expression dynamics, patterns of mRNA decay are less well-understood. Here we directly measured spatiotemporally resolved mRNA decay rates transcriptome-wide throughout C. elegans embryogenesis by transcription inhibition followed by bulk and single-cell RNA-sequencing. This allowed us to calculate mRNA half-lives within specific cell types and developmental stages and identify differentially regulated mRNA decay throughout embryonic development. We identified transcript features that are correlated with mRNA stability and found that mRNA decay rates are associated with distinct peaks in gene expression over time. Moreover, we provide evidence that, on average, mRNA is more stable in the germline compared to in the soma and in later embryonic stages compared to in earlier stages. This work suggests that differential mRNA decay across cell states and time helps to shape developmental gene expression, and it provides a valuable resource for studies of mRNA turnover regulatory mechanisms.

Predictive modeling of oocyte maternal mRNA features for five mammalian species reveals potential shared and species-restricted regulators during maturation

Oocyte maturation is accompanied by changes in abundances of thousands of mRNAs, many degraded and many preferentially stabilized. mRNA stability can be regulated by diverse features including GC content, codon bias, and motifs within the 3'-untranslated region (UTR) interacting with RNA binding proteins (RBPs) and miRNAs. Many studies have identified factors participating in mRNA splicing, bulk mRNA storage, and translational recruitment in mammalian oocytes, but the roles of potentially hundreds of expressed factors, how they regulate cohorts of thousands of mRNAs, and to what extent their functions are conserved across species has not been determined. We performed an extensive in silico cross-species analysis of features associated with mRNAs of different stability classes during oocyte maturation (stable, moderately degraded, and highly degraded) for five mammalian species. Using publicly available RNA sequencing data for germinal vesicle (GV) and MII oocyte transcriptomes, we determined that 3'-UTR length and synonymous codon usage are positively associated with stability, while greater GC content is negatively associated with stability. By applying machine learning and feature selection strategies, we identified RBPs and miRNAs that are predictive of mRNA stability, including some across multiple species and others more species-restricted. The results provide new insight into the mechanisms regulating maternal mRNA stabilization or degradation.NEW & NOTEWORTHY Conservation across species of mRNA features regulating maternal mRNA stability during mammalian oocyte maturation was analyzed. 3'-Untranslated region length and synonymous codon usage are positively associated with stability, while GC content is negatively associated. Just three RNA binding protein motifs were predicted to regulate mRNA stability across all five species examined, but associated pathways and functions are shared, indicating oocytes of different species arrive at comparable physiological destinations via different routes.

Synonymous codon usage regulates translation initiation

Nonoptimal synonymous codons repress gene expression, but the underlying mechanisms are poorly understood. We and others have previously shown that nonoptimal codons slow translation elongation speeds and thereby trigger messenger RNA (mRNA) degradation. Nevertheless, transcript levels are often insufficient to explain protein levels, suggesting additional mechanisms by which codon usage regulates gene expression. Using reporters in human and Drosophila cells, we find that transcript levels account for less than half of the variation in protein abundance due to codon usage. This discrepancy is explained by translational differences whereby nonoptimal codons repress translation initiation. Nonoptimal transcripts are also less bound by the translation initiation factors eIF4E and eIF4G1, providing a mechanistic explanation for their reduced initiation rates. Importantly, translational repression can occur without mRNA decay and deadenylation, and it does not depend on the known nonoptimality sensor, CNOT3. Our results reveal a potent mechanism of regulation by codon usage where nonoptimal codons repress further rounds of translation.

Glycosylated queuosines in tRNAs optimize translational rate and post-embryonic growth

Transfer RNA (tRNA) modifications are critical for protein synthesis. Queuosine (Q), a 7-deaza-guanosine derivative, is present in tRNA anticodons. In vertebrate tRNAs for Tyr and Asp, Q is further glycosylated with galactose and mannose to generate galQ and manQ, respectively. However, biogenesis and physiological relevance of Q-glycosylation remain poorly understood. Here, we biochemically identified two RNA glycosylases, QTGAL and QTMAN, and successfully reconstituted Q-glycosylation of tRNAs using nucleotide diphosphate sugars. Ribosome profiling of knockout cells revealed that Q-glycosylation slowed down elongation at cognate codons, UAC and GAC (GAU), respectively. We also found that galactosylation of Q suppresses stop codon readthrough. Moreover, protein aggregates increased in cells lacking Q-glycosylation, indicating that Q-glycosylation contributes to proteostasis. Cryo-EM of human ribosome-tRNA complex revealed the molecular basis of codon recognition regulated by Q-glycosylations. Furthermore, zebrafish qtgal and qtman knockout lines displayed shortened body length, implying that Q-glycosylation is required for post-embryonic growth in vertebrates.

Znf598-mediated Rps10/eS10 ubiquitination contributes to the ribosome ubiquitination dynamics during zebrafish development

The ribosome is a translational apparatus that comprises about 80 ribosomal proteins and four rRNAs. Recent studies reported that ribosome ubiquitination is crucial for translational regulation and ribosome-associated quality control (RQC). However, little is known about the dynamics of ribosome ubiquitination under complex biological processes of multicellular organisms. To explore ribosome ubiquitination during animal development, we generated a zebrafish strain that expresses a FLAG-tagged ribosomal protein Rpl36/eL36 from its endogenous locus. We examined ribosome ubiquitination during zebrafish development by combining affinity purification of ribosomes from rpl36-FLAG zebrafish embryos with immunoblotting analysis. Our findings showed that the ubiquitination of ribosomal proteins dynamically changed as development proceeded. We also showed that during zebrafish development, the ribosome was ubiquitinated by Znf598, an E3 ubiquitin ligase that activates RQC. Ribosomal protein Rps10/eS10 was found to be a key ubiquitinated protein during development. Furthermore, we showed that Rps10/eS10 ubiquitination-site mutations reduced the overall ubiquitination pattern of the ribosome. These results demonstrate the complexity and dynamics of ribosome ubiquitination during zebrafish development.

Translational control of furina by an RNA regulon is important for left-right patterning, heart morphogenesis and cardiac valve function

Heart development is a complex process that requires asymmetric positioning of the heart, cardiac growth and valve morphogenesis. The mechanisms controlling heart morphogenesis and valve formation are not fully understood. The pro-convertase FurinA functions in heart development across vertebrates. How FurinA activity is regulated during heart development is unknown. Through computational analysis of the zebrafish transcriptome, we identified an RNA motif in a variant FurinA transcript harbouring a long 3' untranslated region (3'UTR). The alternative 3'UTR furina isoform is expressed prior to organ positioning. Somatic deletions in the furina 3'UTR lead to embryonic left-right patterning defects. Reporter localisation and RNA-binding assays show that the furina 3'UTR forms complexes with the conserved RNA-binding translational repressor, Ybx1. Conditional ybx1 mutant embryos show premature and increased Furin reporter expression, abnormal cardiac morphogenesis and looping defects. Mutant ybx1 hearts have an expanded atrioventricular canal, abnormal sino-atrial valves and retrograde blood flow from the ventricle to the atrium. This is similar to observations in humans with heart valve regurgitation. Thus, the furina 3'UTR element/Ybx1 regulon is important for translational repression of FurinA and regulation of heart development.

Mature mRNA processing that deletes 3' end sequences directs translational activation and embryonic development

Eggs accumulate thousands of translationally repressed mRNAs that are translated into proteins after fertilization to direct diverse developmental processes. However, molecular mechanisms underlying the translation of stored mRNAs after fertilization remain unclear. Here, we report a previously unknown RNA processing of 3' end sequences of mature mRNAs that activates the translation of stored mRNAs. Specifically, 9 to 72 nucleotides at the 3' ends of zebrafish pou5f3 and mouse Pou5f1 mRNAs were deleted in the early stages of development. Reporter assays illustrated the effective translation of the truncated forms of mRNAs. Moreover, promotion and inhibition of the shortening of 3' ends accelerated and attenuated Pou5f3 accumulation, respectively, resulting in defective development. Identification of proteins binding to unprocessed and/or processed mRNAs revealed that mRNA shortening acts as molecular switches. Comprehensive analysis revealed that >250 mRNAs underwent this processing. Therefore, our results provide a molecular principle that triggers the translational activation and directs development.

Eukaryotic mRNA decapping factors: molecular mechanisms and activity

Decapping is the enzymatic removal of 5' cap structures from mRNAs in eukaryotic cells. Cap structures normally enhance mRNA translation and stability, and their excision commits an mRNA to complete 5'-3' exoribonucleolytic digestion and generally ends the physical and functional cellular presence of the mRNA. Decapping plays a pivotal role in eukaryotic cytoplasmic mRNA turnover and is a critical and highly regulated event in multiple 5'-3' mRNA decay pathways, including general 5'-3' decay, nonsense-mediated mRNA decay (NMD), AU-rich element-mediated mRNA decay, microRNA-mediated gene silencing, and targeted transcript-specific mRNA decay. In the yeast Saccharomyces cerevisiae, mRNA decapping is carried out by a single Dcp1-Dcp2 decapping enzyme in concert with the accessory activities of specific regulators commonly known as decapping activators or enhancers. These regulatory proteins include the general decapping activators Edc1, 2, and 3, Dhh1, Scd6, Pat1, and the Lsm1-7 complex, as well as the NMD-specific factors, Upf1, 2, and 3. Here, we focus on in vivo mRNA decapping regulation in yeast. We summarize recently uncovered molecular mechanisms that control selective targeting of the yeast decapping enzyme and discuss new roles for specific decapping activators in controlling decapping enzyme targeting, assembly of target-specific decapping complexes, and the monitoring of mRNA translation. Further, we discuss the kinetic contribution of mRNA decapping for overall decay of different substrate mRNAs and highlight experimental evidence pointing to the functional coordination and physical coupling between events in mRNA deadenylation, decapping, and 5'-3' exoribonucleolytic decay.

microRNAs slow translating ribosomes to prevent protein misfolding in eukaryotes

Slower translation rates reduce protein misfolding. Such reductions in speed can be mediated by the presence of non-optimal codons, which allow time for proper folding to occur. Although this phenomenon is conserved from bacteria to humans, it is not known whether there are additional eukaryote-specific mechanisms which act in the same way. MicroRNAs (miRNAs), not present in prokaryotes, target both coding sequences (CDS) and 3' untranslated regions (UTR). Given their low suppressive efficiency, it has been unclear why miRNAs are equally likely to bind to a CDS. Here, we show that miRNAs transiently stall translating ribosomes, preventing protein misfolding with little negative effect on protein abundance. We first analyzed ribosome profiles and miRNA binding sites to examine whether miRNAs stall ribosomes. Furthermore, either global or specific miRNA deficiency accelerated ribosomes and induced aggregation of a misfolding-prone polypeptide reporter. These defects were rescued by slowing ribosomes using non-cleaving shRNAs as miRNA mimics. We finally show that proinsulin misfolding, associated with type II diabetes, was resolved by non-cleaving shRNAs. Our findings provide a eukaryote-specific mechanism of co-translational protein folding and a previously unknown mechanism of action to target protein misfolding diseases.

Specific recognition and ubiquitination of translating ribosomes by mammalian CCR4-NOT

Translation affects messenger RNA stability and, in yeast, this is mediated by the Ccr4-Not deadenylation complex. The details of this process in mammals remain unclear. Here, we use cryogenic electron microscopy (cryo-EM) and crosslinking mass spectrometry to show that mammalian CCR4-NOT specifically recognizes ribosomes that are stalled during translation elongation in an in vitro reconstituted system with rabbit and human components. Similar to yeast, mammalian CCR4-NOT inserts a helical bundle of its CNOT3 subunit into the empty E site of the ribosome. Our cryo-EM structure shows that CNOT3 also locks the L1 stalk in an open conformation to inhibit further translation. CCR4-NOT is required for stable association of the nonconstitutive subunit CNOT4, which ubiquitinates the ribosome, likely to signal stalled translation elongation. Overall, our work shows that human CCR4-NOT not only detects but also enforces ribosomal stalling to couple translation and mRNA decay.

Parallels and contrasts between the cnidarian and bilaterian maternal-to-zygotic transition are revealed in Hydractinia embryos

Embryogenesis requires coordinated gene regulatory activities early on that establish the trajectory of subsequent development, during a period called the maternal-to-zygotic transition (MZT). The MZT comprises transcriptional activation of the embryonic genome and post-transcriptional regulation of egg-inherited maternal mRNA. Investigation into the MZT in animals has focused almost exclusively on bilaterians, which include all classical models such as flies, worms, sea urchin, and vertebrates, thus limiting our capacity to understand the gene regulatory paradigms uniting the MZT across all animals. Here, we elucidate the MZT of a non-bilaterian, the cnidarian Hydractinia symbiolongicarpus. Using parallel poly(A)-selected and non poly(A)-dependent RNA-seq approaches, we find that the Hydractinia MZT is composed of regulatory activities similar to many bilaterians, including cytoplasmic readenylation of maternally contributed mRNA, delayed genome activation, and separate phases of maternal mRNA deadenylation and degradation that likely depend on both maternally and zygotically encoded clearance factors, including microRNAs. But we also observe massive upregulation of histone genes and an expanded repertoire of predicted H4K20 methyltransferases, aspects thus far particular to the Hydractinia MZT and potentially underlying a novel mode of early embryonic chromatin regulation. Thus, similar regulatory strategies with taxon-specific elaboration underlie the MZT in both bilaterian and non-bilaterian embryos, providing insight into how an essential developmental transition may have arisen in ancestral animals.

Translation and mRNA Stability Control

Messenger RNA (mRNA) stability and translational efficiency are two crucial aspects of the post-transcriptional process that profoundly impact protein production in a cell. While it is widely known that ribosomes produce proteins, studies during the past decade have surprisingly revealed that ribosomes also control mRNA stability in a codon-dependent manner, a process referred to as codon optimality. Therefore, codons, the three-nucleotide words read by the ribosome, have a potent effect on mRNA stability and provide cis-regulatory information that extends beyond the amino acids they encode. While the codon optimality molecular mechanism is still unclear, the translation elongation rate appears to trigger mRNA decay. Thus, transfer RNAs emerge as potential master gene regulators affecting mRNA stability. Furthermore, while few factors related to codon optimality have been identified in yeast, the orthologous genes in vertebrates do not necessary share the same functions. Here, we discuss codon optimality findings and gene regulation layers related to codon composition in different eukaryotic species.

Quality control of protein synthesis in the early elongation stage

In the early stage of bacterial translation, peptidyl-tRNAs frequently dissociate from the ribosome (pep-tRNA drop-off) and are recycled by peptidyl-tRNA hydrolase. Here, we establish a highly sensitive method for profiling of pep-tRNAs using mass spectrometry, and successfully detect a large number of nascent peptides from pep-tRNAs accumulated in Escherichia coli pthts strain. Based on molecular mass analysis, we found about 20% of the peptides bear single amino-acid substitutions of the N-terminal sequences of E. coli ORFs. Detailed analysis of individual pep-tRNAs and reporter assay revealed that most of the substitutions take place at the C-terminal drop-off site and that the miscoded pep-tRNAs rarely participate in the next round of elongation but dissociate from the ribosome. These findings suggest that pep-tRNA drop-off is an active mechanism by which the ribosome rejects miscoded pep-tRNAs in the early elongation, thereby contributing to quality control of protein synthesis after peptide bond formation.

Parallels and contrasts between the cnidarian and bilaterian maternal-to-zygotic transition are revealed in Hydractinia embryos

Embryogenesis requires coordinated gene regulatory activities early on that establish the trajectory of subsequent development, during a period called the maternal-to-zygotic transition (MZT). The MZT comprises transcriptional activation of the embryonic genome and post-transcriptional regulation of egg-inherited maternal mRNA. Investigation into the MZT in animals has focused almost exclusively on bilaterians, which include all classical models such as flies, worms, sea urchin, and vertebrates, thus limiting our capacity to understand the gene regulatory paradigms uniting the MZT across all animals. Here, we elucidate the MZT of a non-bilaterian, the cnidarian Hydractinia symbiolongicarpus . Using parallel poly(A)-selected and non poly(A)-dependent RNA-seq approaches, we find that the Hydractinia MZT is composed of regulatory activities analogous to many bilaterians, including cytoplasmic readenylation of maternally contributed mRNA, delayed genome activation, and separate phases of maternal mRNA deadenylation and degradation that likely depend on both maternally and zygotically encoded clearance factors, including microRNAs. But we also observe massive upregulation of histone genes and an expanded repertoire of predicted H4K20 methyltransferases, aspects thus far unique to the Hydractinia MZT and potentially underlying a novel mode of early embryonic chromatin regulation. Thus, similar regulatory strategies with taxon-specific elaboration underlie the MZT in both bilaterian and non-bilaterian embryos, providing insight into how an essential developmental transition may have arisen in ancestral animals.

mRNA transport, translation, and decay in adult mammalian central nervous system axons

Localized mRNA translation regulates synapse function and axon maintenance, but how compartment-specific mRNA repertoires are regulated is largely unknown. We developed an axonal transcriptome capture method that allows deep sequencing of metabolically labeled mRNAs from retinal ganglion cell axon terminals in mouse. Comparing axonal-to-somal transcriptomes and axonal translatome-to-transcriptome enables genome-wide visualization of mRNA transport and translation and unveils potential regulators tuned to each process. FMRP and TDP-43 stand out as key regulators of transport, and experiments in Fmr1 knockout mice validate FMRP's role in the axonal transportation of synapse-related mRNAs. Pulse-and-chase experiments enable genome-wide assessment of mRNA stability in axons and reveal a strong coupling between mRNA translation and decay. Measuring the absolute mRNA abundance per axon terminal shows that the adult axonal transcriptome is stably maintained by persistent transport. Our datasets provide a rich resource for unique insights into RNA-based mechanisms in maintaining presynaptic structure and function in vivo.

Transcriptomic, proteomic, and functional consequences of codon usage bias in human cells during heterologous gene expression

Differences in codon frequency between genomes, genes, or positions along a gene, modulate transcription and translation efficiency, leading to phenotypic and functional differences. Here, we present a multiscale analysis of the effects of synonymous codon recoding during heterologous gene expression in human cells, quantifying the phenotypic consequences of codon usage bias at different molecular and cellular levels, with an emphasis on translation elongation. Six synonymous versions of an antibiotic resistance gene were generated, fused to a fluorescent reporter, and independently expressed in HEK293 cells. Multiscale phenotype was analyzed by means of quantitative transcriptome and proteome assessment, as proxies for gene expression; cellular fluorescence, as a proxy for single-cell level expression; and real-time cell proliferation in absence or presence of antibiotic, as a proxy for the cell fitness. We show that differences in codon usage bias strongly impact the molecular and cellular phenotype: (i) they result in large differences in mRNA levels and protein levels, leading to differences of over 15 times in translation efficiency; (ii) they introduce unpredicted splicing events; (iii) they lead to reproducible phenotypic heterogeneity; and (iv) they lead to a trade-off between the benefit of antibiotic resistance and the burden of heterologous expression. In human cells in culture, codon usage bias modulates gene expression by modifying mRNA availability and suitability for translation, leading to differences in protein levels and eventually eliciting functional phenotypic changes.

SLAMseq resolves the kinetics of maternal and zygotic gene expression during early zebrafish embryogenesis

The maternal-to-zygotic transition (MZT) is a key developmental process in metazoan embryos that involves the activation of zygotic transcription (ZGA) and degradation of maternal transcripts. We employed metabolic mRNA sequencing (SLAMseq) to deconvolute the compound embryonic transcriptome in zebrafish. While mitochondrial zygotic transcripts prevail prior to MZT, we uncover the spurious transcription of hundreds of short and intron-poor genes as early as the 2-cell stage. Upon ZGA, most zygotic transcripts originate from thousands of maternal-zygotic (MZ) genes that are transcribed at rates comparable to those of hundreds of purely zygotic genes and replenish maternal mRNAs at distinct timescales. Rapid replacement of MZ transcripts involves transcript decay features unrelated to major maternal degradation pathways and promotes de novo synthesis of the core gene expression machinery by increasing poly(A)-tail length and translation efficiency. SLAMseq hence provides insights into the timescales, molecular features, and regulation of MZT during zebrafish embryogenesis.

Maternal mRNA deadenylation and allocation via Rbm14 condensates facilitate vertebrate blastula development

Early embryonic development depends on proper utilization and clearance of maternal transcriptomes. How these processes are spatiotemporally regulated remains unclear. Here we show that nuclear RNA-binding protein Rbm14 and maternal mRNAs co-phase separate into cytoplasmic condensates to facilitate vertebrate blastula-to-gastrula development. In zebrafish, Rbm14 condensates were highly abundant in blastomeres and markedly reduced after prominent activation of zygotic transcription. They concentrated at spindle poles by associating with centrosomal γ-tubulin puncta and displayed mainly asymmetric divisions with a global symmetry across embryonic midline in 8- and 16-cell embryos. Their formation was dose-dependently stimulated by m6 A, but repressed by m5 C modification of the maternal mRNA. Furthermore, deadenylase Parn co-phase separated with these condensates, and this was required for deadenylation of the mRNAs in early blastomeres. Depletion of Rbm14 impaired embryonic cell differentiations and full activations of the zygotic genome in both zebrafish and mouse and resulted in developmental arrest at the blastula stage. Our results suggest that cytoplasmic Rbm14 condensate formation regulates early embryogenesis by facilitating deadenylation, protection, and mitotic allocation of m6 A-modified maternal mRNAs, and by releasing the poly(A)-less transcripts upon regulated disassembly to allow their re-polyadenylation and translation or clearance.

Codon modification of Tuba1a alters mRNA levels and causes a severe neurodevelopmental phenotype in mice

The tubulinopathies are an umbrella of rare diseases that result from mutations in tubulin genes and are frequently characterised by severe brain malformations. The characteristics of a given disease reflect the expression pattern of the transcript, the function of a given tubulin gene, and the role microtubules play in a particular cell type. Mouse models have proved to be valuable tools that have provided insight into the molecular and cellular mechanisms that underlie the disease state. In this manuscript we compare two Tuba1a mouse models, both of which express wild-type TUBA1A protein but employ different codon usage. We show that modification of the Tuba1a mRNA sequence results in homozygous lethality and a severe neurodevelopmental phenotype. This is associated with a decrease in the number of post-mitotic neurons, PAX6 positive progenitors, and an increase in the number of apoptotic cells. We attribute this to a decrease in the stability of the modified Tuba1a transcript, and the absence of compensation by the other neurogenic tubulins. Our findings highlight the importance of maintaining the wild-type coding sequence when engineering mouse lines and the impact of synonymous genetic variation.

The role of post-transcriptional modifications during development

While the existence of post-transcriptional modifications of RNA nucleotides has been known for decades, in most RNA species the exact positions of these modifications and their physiological function have been elusive until recently. Technological advances, such as high-throughput next-generation sequencing (NGS) methods and nanopore-based mapping technologies, have made it possible to map the position of these modifications with single nucleotide accuracy, and genetic screens have uncovered the “writer”, “reader” and “eraser” proteins that help to install, interpret and remove such modifications, respectively. These discoveries led to intensive research programmes with the aim of uncovering the roles of these modifications during diverse biological processes. In this review, we assess novel discoveries related to the role of post-transcriptional modifications during animal development, highlighting how these discoveries can affect multiple aspects of development from fertilization to differentiation in many species.

A nascent peptide code for translational control of mRNA stability in human cells

Stability of eukaryotic mRNAs is associated with their codon, amino acid, and GC content. Yet, coding sequence motifs that predictably alter mRNA stability in human cells remain poorly defined. Here, we develop a massively parallel assay to measure mRNA effects of thousands of synthetic and endogenous coding sequence motifs in human cells. We identify several families of simple dipeptide repeats whose translation triggers mRNA destabilization. Rather than individual amino acids, specific combinations of bulky and positively charged amino acids are critical for the destabilizing effects of dipeptide repeats. Remarkably, dipeptide sequences that form extended β strands in silico and in vitro slowdown ribosomes and reduce mRNA levels in vivo. The resulting nascent peptide code underlies the mRNA effects of hundreds of endogenous peptide sequences in the human proteome. Our work suggests an intrinsic role for the ribosome as a selectivity filter against the synthesis of bulky and aggregation-prone peptides.

Codon Usage and mRNA Stability are Translational Determinants of Cellular Response to Canonical Ferroptosis Inducers

Ferroptosis is a non-apoptotic cell death mechanism characterized by the generation of lipid peroxides. While many effectors in the ferroptosis pathway have been mapped, its epitranscriptional regulation is not yet fully understood. Ferroptosis can be induced via system xCT inhibition (Class I) or GPX4 inhibition (Class II). Previous works have revealed important differences in cellular response to different ferroptosis inducers. Importantly, blocking mRNA transcription or translation appears to protect cells against Class I ferroptosis inducing agents but not Class II. In this work, we examined the impact of blocking transcription (via Actinomycin D) or translation (via Cycloheximide) on Erastin (Class I) or RSL3 (Class II) induced ferroptosis. Blocking transcription or translation protected cells against Erastin but was detrimental against RSL3. Cycloheximide led to increased levels of GSH alone or when co-treated with Erastin via the activation of the reverse transsulfuration pathway. RNA sequencing analysis revealed early activation of a strong alternative splice program before observed changes in transcription. mRNA stability analysis revealed divergent mRNA stability changes in cellular response to Erastin or RSL3. Importantly, codon optimality biases were drastically different in either condition. Our data also implicated translation repression and rate as an important determinant of the cellular response to ferroptosis inducers. Given that mRNA stability and codon usage can be influenced via the tRNA epitranscriptome, we evaluated the role of a tRNA modifying enzyme in ferroptosis stress response. Alkbh1, a tRNA demethylase, led to translation repression and increased the resistance to Erastin but made cells more sensitive to RSL3.

Translatome and transcriptome co-profiling reveals a role of TPRXs in human zygotic genome activation

Translational regulation plays a critical role during the oocyte-to-embryo transition (OET) and zygotic genome activation (ZGA). Here, we integrated ultra-low-input Ribo-seq with mRNA-seq to co-profile the translatome and transcriptome in human oocytes and early embryos. Comparison with mouse counterparts identified widespread differentially translated genes functioning in epigenetic reprogramming, transposon defense, and small RNA biogenesis, in part driven by species-specific regulatory elements in 3' untranslated regions. Moreover, PRD-like homeobox transcription factors, including TPRXL, TPRX1, and TPRX2, are highly translated around ZGA. TPRX1/2/L knockdown leads to defective ZGA and preimplantation development. Ectopically expressed TPRXs bind and activate key ZGA genes in human embryonic stem cells. These data reveal the conservation and divergence of translation landscapes during OET and identify critical regulators of human ZGA.

Geneticin reduces mRNA stability

Messenger RNA (mRNA) translation can lead to higher rates of mRNA decay, suggesting the ribosome plays a role in mRNA destruction. Furthermore, mRNA features, such as codon identities, which are directly probed by the ribosome, correlate with mRNA decay rates. Many amino acids are encoded by synonymous codons, some of which are decoded by more abundant tRNAs leading to more optimal translation and increased mRNA stability. Variable translation rates for synonymous codons can lead to ribosomal collisions as ribosomes transit regions with suboptimal codons, and ribosomal collisions can promote mRNA decay. In addition to different translation rates, the presence of certain codons can also lead to higher or lower rates of amino acid misincorporation which could potentially lead to protein misfolding if a substituted amino acid fails to make critical contacts in a structure. Here, we test whether Geneticin—G418, an aminoglycoside antibiotic known to promote amino acid misincorporation—affects mRNA stability. We observe that G418 decreases firefly luciferase mRNA stability in an in vitro translation system and also reduces mRNA stability in mouse embryonic stem cells (mESCs). G418-sensitive mRNAs are enriched for certain optimal codons that contain G or C in the wobble position, arguing that G418 blunts the stabilizing effects of codon optimality.

Distinct elongation stalls during translation are linked with distinct pathways for mRNA degradation

Key protein adapters couple translation to mRNA decay on specific classes of problematic mRNAs in eukaryotes. Slow decoding on non-optimal codons leads to codon-optimality-mediated decay (COMD) and prolonged arrest at stall sites leads to no-go decay (NGD). The identities of the decay factors underlying these processes and the mechanisms by which they respond to translational distress remain open areas of investigation. We use carefully designed reporter mRNAs to perform genetic screens and functional assays in Saccharomyces cerevisiae. We characterize the roles of Hel2, Syh1, and Smy2 in coordinating translational repression and mRNA decay on NGD reporter mRNAs, finding that Syh1 and, to a lesser extent its paralog Smy2, act in a distinct pathway from Hel2. This Syh1/Smy2-mediated pathway acts as a redundant, compensatory pathway to elicit NGD when Hel2-dependent NGD is impaired. Importantly, we observe that these NGD factors are not involved in the degradation of mRNAs enriched in non-optimal codons. Further, we establish that a key factor previously implicated in COMD, Not5, contributes modestly to the degradation of an NGD-targeted mRNA. Finally, we use ribosome profiling to reveal distinct ribosomal states associated with each reporter mRNA that readily rationalize the contributions of NGD and COMD factors to degradation of these reporters. Taken together, these results provide new insight into the role of Syh1 and Smy2 in NGD and into the ribosomal states that correlate with the activation of distinct pathways targeting mRNAs for degradation in yeast.

iCodon customizes gene expression based on the codon composition

Messenger RNA (mRNA) stability substantially impacts steady-state gene expression levels in a cell. mRNA stability is strongly affected by codon composition in a translation-dependent manner across species, through a mechanism termed codon optimality. We have developed iCodon (www.iCodon.org), an algorithm for customizing mRNA expression through the introduction of synonymous codon substitutions into the coding sequence. iCodon is optimized for four vertebrate transcriptomes: mouse, human, frog, and fish. Users can predict the mRNA stability of any coding sequence based on its codon composition and subsequently generate more stable (optimized) or unstable (deoptimized) variants encoding for the same protein. Further, we show that codon optimality predictions correlate with both mRNA stability using a massive reporter library and expression levels using fluorescent reporters and analysis of endogenous gene expression in zebrafish embryos and/or human cells. Therefore, iCodon will benefit basic biological research, as well as a wide range of applications for biotechnology and biomedicine.

METTL3-Mediated m6A RNA Methylation of ZBTB4 Interferes With Trophoblast Invasion and Maybe Involved in RSA

N⁶-methyladenosine (m⁶A) was the most abundant modification of mRNA and lncRNA in mammalian cells and played an important role in many biological processes. However, whether m⁶A modification was associated with recurrent spontaneous abortion (RSA) and its roles were still unclear.

Methods: Methylated RNA immunoprecipitation sequencing (MeRIP-Seq) was used to study the global m⁶A modification pattern in RSAs and controls. RNA sequencing (RNA-Seq) was used to study the level of global mRNA in two groups. Real-time quantitative PCR (RT-qPCR) was used to verify the level of mRNA of METTL3 and ZBTB4. MeRIP–qPCR was conducted to test the level of ZBTB4 m⁶A modification in two groups. In order to further explore whether ZBTB4 was the substrate of METTL3, the HTR-8/SVneo (HTR-8) cell line was selected for the knockdown and overexpression of METTL3. To study whether METTL3 regulated the ZBTB4 expression by recognizing ZBTB4 mRNA m⁶A motifs in coding sequences (CDS), dual-luciferase reporter assay was conducted. RNA stability assays using actinomycin D were conducted to study the RNA stability of the HTR-8 cell line with METTL3 overexpression and knockdown. To illustrate the role of METTL3 in the invasion of trophoblast, matrigel invasion assays and transwell migration assays were conducted using the HTR-8 cell line with METTL3 overexpression and knockdown.

Results: A total of 65 genes were found with significant differences both in m⁶A modification and mRNA expression. We found m⁶A methyltransferase METTL3 was significantly down-regulated in the RSA group. Through gene function analysis, RT-qPCR, MeRIP–qPCR validation experiment, knockdown, and overexpression of METTL3 in the HTR-8 cell line, ZBTB4 was selected as one target of METTL3. Furthermore, we clarified that METTL3 regulated the expression of ZBTB4 by recognizing ZBTB4 mRNA m⁶A motifs in the CDS using the dual-luciferase reporter assay and METTL3 regulated the invasion of trophoblast by altering the stability and expression of ZBTB4 by RNA stability, matrigel invasion, and transwell migration assays.

Conclusion: Our study revealed the mechanism by which METTL3 regulated the stability and expression of ZBTB4 and the trophoblast migration ability of RSA. A new perspective was provided for exploring the mechanism of embryonic development in RSA patients.

Phase separation of Ddx3xb helicase regulates maternal-to-zygotic transition in zebrafish

Vertebrate embryogenesis involves a conserved and fundamental process, called the maternal-to-zygotic transition (MZT), which marks the switch from a maternal factors-dominated state to a zygotic factors-driven state. Yet the precise mechanism underlying MZT remains largely unknown. Here we report that the RNA helicase Ddx3xb in zebrafish undergoes liquid-liquid phase separation (LLPS) via its N-terminal intrinsically disordered region (IDR), and an increase in ATP content promotes the condensation of Ddx3xb during MZT. Mutant form of Ddx3xb losing LLPS ability fails to rescue the developmental defect of Ddx3xb-deficient embryos. Interestingly, the IDR of either FUS or hnRNPA1 can functionally replace the N-terminal IDR in Ddx3xb. Phase separation of Ddx3xb facilitates the unwinding of 5' UTR structures of maternal mRNAs to enhance their translation. Our study reveals an unprecedent mechanism whereby the Ddx3xb phase separation regulates MZT by promoting maternal mRNA translation.

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.

Sorting and packaging of RNA into extracellular vesicles shape intracellular transcript levels

BACKGROUND

Extracellular vesicles (EVs) are released by nearly every cell type and have attracted much attention for their ability to transfer protein and diverse RNA species from donor to recipient cells. Much attention has been given so far to the features of EV short RNAs such as miRNAs. However, while the presence of mRNA and long noncoding RNA (lncRNA) transcripts in EVs has also been reported by multiple different groups, the properties and function of these longer transcripts have been less thoroughly explored than EV miRNA. Additionally, the impact of EV export on the transcriptome of exporting cells has remained almost completely unexamined. Here, we globally investigate mRNA and lncRNA transcripts in endothelial EVs in multiple different conditions.

RESULTS

In basal conditions, long RNA transcripts enriched in EVs have longer than average half-lives and distinctive stability-related sequence and structure characteristics including shorter transcript length, higher exon density, and fewer 3' UTR A/U-rich elements. EV-enriched long RNA transcripts are also enriched in HNRNPA2B1 binding motifs and are impacted by HNRNPA2B1 depletion, implicating this RNA-binding protein in the sorting of long RNA to EVs. After signaling-dependent modification of the cellular transcriptome, we observed that, unexpectedly, the rate of EV enrichment relative to cells was altered for many mRNA and lncRNA transcripts. This change in EV enrichment was negatively correlated with intracellular abundance, with transcripts whose export to EVs increased showing decreased abundance in cells and vice versa. Correspondingly, after treatment with inhibitors of EV secretion, levels of mRNA and lncRNA transcripts that are normally highly exported to EVs increased in cells, indicating a measurable impact of EV export on the long RNA transcriptome of the exporting cells. Compounds with different mechanisms of inhibition of EV secretion affected the cellular transcriptome differently, suggesting the existence of multiple EV subtypes with different long RNA profiles.

CONCLUSIONS

We present evidence for an impact of EV physiology on the characteristics of EV-producing cell transcriptomes. Our work suggests a new paradigm in which the sorting and packaging of transcripts into EVs participate, together with transcription and RNA decay, in controlling RNA homeostasis and shape the cellular long RNA abundance profile.

Genetic variants associated mRNA stability in lung

Background

Expression quantitative trait loci (eQTLs) analyses have been widely used to identify genetic variants associated with gene expression levels to understand what molecular mechanisms underlie genetic traits. The resultant eQTLs might affect the expression of associated genes through transcriptional or post-transcriptional regulation. In this study, we attempt to distinguish these two types of regulation by identifying genetic variants associated with mRNA stability of genes (stQTLs).

Results

Here, we presented a computational framework that takes advantage of recently developed methods to infer the mRNA stability of genes based on RNA-seq data and performed association analysis to identify stQTLs. Using the Genotype-Tissue Expression (GTEx) lung RNA-Seq data, we identified a total of 142,801 stQTLs for 3942 genes and 186,132 eQTLs for 4751 genes from 15,122,700 genetic variants for 13,476 genes on the autosomes, respectively. Interestingly, our results indicated that stQTLs were enriched in the CDS and 3’UTR regions, while eQTLs are enriched in the CDS, 3’UTR, 5’UTR, and upstream regions. We also found that stQTLs are more likely than eQTLs to overlap with RNA binding protein (RBP) and microRNA (miRNA) binding sites. Our analyses demonstrate that simultaneous identification of stQTLs and eQTLs can provide more mechanistic insight on the association between genetic variants and gene expression levels.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08405-y.

Ribosome slowdown triggers codon-mediated mRNA decay independently of ribosome quality control

The control of mRNA stability plays a central role in regulating gene expression patterns. Recent studies have revealed that codon composition in the open reading frame determines mRNA stability in multiple organisms. Based on genome-wide correlation approaches, this previously unrecognized role for the genetic code is attributable to the kinetics of the codon-decoding process by the ribosome. However, complementary experimental analyses are required to clarify the codon effects on mRNA stability and the related cotranslational mRNA decay pathways, for example, those triggered by aberrant ribosome stalling. In the current study, we performed a set of reporter-based analyses to define codon-mediated mRNA decay and ribosome stall-dependent mRNA decay in zebrafish embryos. Our analysis showed that the effect of codons on mRNA stability stems from the decoding process, independent of the ribosome quality control factor Znf598 and stalling-dependent mRNA decay. We propose that codon-mediated mRNA decay is rather triggered by transiently slowed ribosomes engaging in a productive translation cycle in zebrafish embryos.

The Reversible Methylation of m6A Is Involved in Plant Virus Infection

Simple Summary

N6-methyladenosine (m6A) is the most prevalent modification in the mRNAs of many eukaryotic species. The abundance and effects of m6A are determined by dynamic interactions between its methyltransferases (“writers”), demethylases (“erasers”), and binding proteins (“readers”). It has been indicated that there is a strong correlation between m6A and virus infection in mammals. In the case of plant virus infection, it appears that m6A plays a dual role. On the one hand, m6A acts as a plant immune response induced by virus infection, inhibiting viral replication or translation through methylation of viral genome RNAs. On the other hand, m6A acts as part of an infection strategy employed by plant viruses to overcome the host immune system by interacting with m6A-related proteins. We proposed that antagonists of m6A-related proteins might be used to design new strategies for plant virus control in the future.

Abstract

In recent years, m6A RNA methylation has attracted broad interest and is becoming a hot research topic. It has been demonstrated that there is a strong association between m6A and viral infection in the human system. The life cycles of plant RNA viruses are often coordinated with the mechanisms of their RNA modification. Here, we reviewed recent advances in m6A methylation in plant viruses. It appears that m6A methylation plays a dual role during viral infection in plants. On the one hand, m6A methylation acts as an antiviral immune response induced by virus infection, which inhibits viral replication or translation through the methylation of viral genome RNAs. On the other hand, plant viruses could disrupt the m6A methylation through interacting with the key proteins of the m6A pathway to avoid modification. Those plant viruses containing ALKB domain are discussed as well. Based on this mechanism, we propose that new strategies for plant virus control could be designed with competitive antagonists of m6A-associated proteins.

Roles of mRNA poly(A) tails in regulation of eukaryotic gene expression

In eukaryotes, poly(A) tails are present on almost every mRNA. Early experiments led to the hypothesis that poly(A) tails and the cytoplasmic polyadenylate-binding protein (PABPC) promote translation and prevent mRNA degradation, but the details remained unclear. More recent data suggest that the role of poly(A) tails is much more complex: poly(A)-binding protein can stimulate poly(A) tail removal (deadenylation) and the poly(A) tails of stable, highly translated mRNAs at steady state are much shorter than expected. Furthermore, the rate of translation elongation affects deadenylation. Consequently, the interplay between poly(A) tails, PABPC, translation and mRNA decay has a major role in gene regulation. In this Review, we discuss recent work that is revolutionizing our understanding of the roles of poly(A) tails in the cytoplasm. Specifically, we discuss the roles of poly(A) tails in translation and control of mRNA stability and how poly(A) tails are removed by exonucleases (deadenylases), including CCR4-NOT and PAN2-PAN3. We also discuss how deadenylation rate is determined, the integration of deadenylation with other cellular processes and the function of PABPC. We conclude with an outlook for the future of research in this field.

Computational modeling of mRNA degradation dynamics using deep neural networks

MOTIVATION

messenger RNA (mRNA) degradation plays critical roles in post-transcriptional gene regulation. A major component of mRNA degradation is determined by 3'-UTR elements. Hence, researchers are interested in studying mRNA dynamics as a function of 3'-UTR elements. A recent study measured the mRNA degradation dynamics of tens of thousands of 3'-UTR sequences using a massively parallel reporter assay. However, the computational approach used to model mRNA degradation was based on a simplifying assumption of a linear degradation rate. Consequently, the underlying mechanism of 3'-UTR elements is still not fully understood.

RESULTS

Here, we developed deep neural networks to predict mRNA degradation dynamics and interpreted the networks to identify regulatory elements in the 3'-UTR and their positional effect. Given an input of a 110 nt-long 3'-UTR sequence and an initial mRNA level, the model predicts mRNA levels of eight consecutive time points. Our deep neural networks significantly improved prediction performance of mRNA degradation dynamics compared with extant methods for the task. Moreover, we demonstrated that models predicting the dynamics of two identical 3'-UTR sequences, differing by their poly(A) tail, performed better than single-task models. On the interpretability front, by using Integrated Gradients, our convolutional neural networks (CNNs) models identified known and novel cis-regulatory sequence elements of mRNA degradation. By applying a novel systematic evaluation of model interpretability, we demonstrated that the recurrent neural network models are inferior to the CNN models in terms of interpretability and that random initialization ensemble improves both prediction and interoperability performance. Moreover, using a mutagenesis analysis, we newly discovered the positional effect of various 3'-UTR elements.

AVAILABILITY AND IMPLEMENTATION

All the code developed through this study is available at github.com/OrensteinLab/DeepUTR/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.