Global analysis of mRNA isoform half-lives reveals stabilizing and destabilizing elements in yeast

Advances and perspectives in genetic expression and operation for the oleaginous yeast Yarrowia lipolytica

The utilization of industrial biomanufacturing has emerged as a viable and sustainable alternative to fossil-based resources for producing functional chemicals. Moreover, advancements in synthetic biology have created new opportunities for the development of innovative cell factories. Notably, Yarrowia lipolytica, an oleaginous yeast that is generally regarded as safe, possesses several advantageous characteristics, including the ability to utilize inexpensive renewable carbon sources, well-established genetic backgrounds, and mature genetic manipulation methods. Consequently, there is increasing interest in manipulating the metabolism of this yeast to enhance its potential as a biomanufacturing platform. Here, we reviewed the latest developments in genetic expression strategies and manipulation tools related to Y. lipolytica, particularly focusing on gene expression, chromosomal operation, CRISPR-based tool, and dynamic biosensors. The purpose of this review is to serve as a valuable reference for those interested in the development of a Y. lipolytica microbial factory.

Chromatin regulates alternative polyadenylation via the RNA polymerase II elongation rate

The RNA polymerase II (Pol II) elongation rate influences poly(A) site selection, with slow and fast Pol II derivatives causing upstream and downstream shifts, respectively, in poly(A) site utilization. In yeast, depletion of either of the histone chaperones FACT or Spt6 causes an upstream shift of poly(A) site use that strongly resembles the poly(A) profiles of slow Pol II mutant strains. Like slow Pol II mutant strains, FACT- and Spt6-depleted cells exhibit Pol II processivity defects, indicating that both Spt6 and FACT stimulate the Pol II elongation rate. Poly(A) profiles of some genes show atypical downstream shifts; this subset of genes overlaps well for FACT- or Spt6-depleted strains but is different from the atypical genes in Pol II speed mutant strains. In contrast, depletion of histone H3 or H4 causes a downstream shift of poly(A) sites for most genes, indicating that nucleosomes inhibit the Pol II elongation rate in vivo. Thus, chromatin-based control of the Pol II elongation rate is a potential mechanism, distinct from direct effects on the cleavage/polyadenylation machinery, to regulate alternative polyadenylation in response to genetic or environmental changes.

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).

Nanocage-incorporated engineered destabilized 3'UTR ARE of ERBB2 inhibits tumor growth and liver and lung metastasis in EGFR T790M osimertinib- and trastuzumab-resistant and ERBB2-expressing NSCLC via the reduction of ERBB2

Non-small cell lung cancer (NSCLC) caused more deaths in 2017 than breast cancer, prostate, and brain cancers combined. This is primarily due to their aggressive metastatic nature, leading to more fatal rates of cancer patients. Despite this condition, there are no clinically approved drugs that can target metastasis. The NSCLC with EGFR T790M-overexpressing HER2 shows the resistance to osimertinib and trastuzumab starting 10-18 months after the therapy, and thus prospects are grim to these patients. To target the recalcitrant ERBB2 driver oncogene, we developed two engineered destabilizing 3'UTR ERBB2 constructs that degrade the endogenous ERBB2 transcript and proteins by overwriting the encoded endogenous ERBB2 mRNA with the destabilizing message. When iron oxide nanocages (IO nanocages) were used as vehicles to deliver them to tumors and whole tissues in mice bearing tumors, it was well tolerated and safe and caused no genome rearrangement whereas they were integrated into genome deserts (non-coding regions). We achieved significant reduction of the primary tumor volume with desARE3'UTRERBB2-30, achieving 50% complete tumor lysis and inhibiting 60%-80% of liver metastasis, hepatomegaly, and 90% of lung metastasis, through ERBB2 downregulation. These constructs were distributed robustly into tumors, livers, lungs, kidneys, and spleen and mildly in the brain and not in the heart. They caused no abnormality in both short- and long-term administrations as well as in healthy mice. In summary, we accomplished significant breakthrough for the therapeutics of intractable lung cancer patients whose cancers become resistant and metastasize.

Metabolic Pathway Coupled with Fermentation Process Optimization for High-Level Production of Retinol in Yarrowia lipolytica

Retinol is a lipid-soluble form of vitamin A that is crucial for human visual and immune functions. The production of retinol through microbial fermentation has been the focus of recent exploration. However, the obtained titer remains limited and the product is often a mixture of retinal, retinol, and retinoic acid, necessitating purification. To achieve efficient biosynthesis of retinol in Yarrowia lipolytica, we improved the metabolic flux of β-carotene to provide sufficient precursors for retinol in this study. Coupled with the optimization of the expression level of β-carotene 15,15'-dioxygenase, de novo production of retinol was achieved. Furthermore, Tween 80 was used as an extractant and butylated hydroxytoluene as an antioxidant to extract intracellular retinol and prevent retinol oxidation, respectively. This strategy significantly increased the level of retinol production. By optimizing the enzymes converting retinal to retinol, the proportion of extracellular retinol in the produced retinoids reached 100%, totaling 1042.3 mg/L. Finally, total retinol production reached 5.4 g/L through fed-batch fermentation in a 5 L bioreactor, comprising 4.2 g/L extracellular retinol and 1.2 g/L intracellular retinol. This achievement represents the highest reported titer so far and advances the industrial production of retinol.

Comprehensive chromatographic assessment of forced degraded in vitro transcribed mRNA

Heightened interest in messenger RNA (mRNA) therapeutics has accelerated the need for analytical methodologies that facilitate the production of supplies for clinical trials. Forced degradation studies are routinely conducted to provide an understanding of potential weak spots in the molecule that are exploited by stresses encountered during bulk purification, production, shipment, and storage. Consequently, temperature fluctuations and excursions are often experienced during these unit operations and may accelerate mRNA degradation. Here, we present a concise panel of chromatography-based stability-indicating assays for evaluating thermally stressed in vitro transcribed (IVT) mRNA as part of a forced degradation study. We found that addition of EDTA to the mRNAs prior to heat exposure reduced the extent of degradation, suggesting that transcripts may be fragmenting via a divalent metal-ion mediated pathway. Trace divalent metal contamination that can accelerate RNA instability is likely carried over from upstream steps. We demonstrate the application of these methods to evaluate the critical quality attributes (CQAs) of mRNAs as well as to detect intrinsic process- and product-related impurities.

Mathematical model for the role of multiple pericentromeric repeats on heterochromatin assembly

Although the length and constituting sequences for pericentromeric repeats are highly variable across eukaryotes, the presence of multiple pericentromeric repeats is one of the conserved features of the eukaryotic chromosomes. Pericentromeric heterochromatin is often misregulated in human diseases, with the expansion of pericentromeric repeats in human solid cancers. In this article, we have developed a mathematical model of the RNAi-dependent methylation of H3K9 in the pericentromeric region of fission yeast. Our model, which takes copy number as an explicit parameter, predicts that the pericentromere is silenced only if there are many copies of repeats. It becomes bistable or desilenced if the copy number of repeats is reduced. This suggests that the copy number of pericentromeric repeats alone can determine the fate of heterochromatin silencing in fission yeast. Through sensitivity analysis, we identified parameters that favor bistability and desilencing. Stochastic simulation shows that faster cell division and noise favor the desilenced state. These results show the unexpected role of pericentromeric repeat copy number in gene silencing and provide a quantitative basis for how the copy number allows or protects repetitive and unique parts of the genome from heterochromatin silencing, respectively.

How is polyadenylation restricted to 3'-untranslated regions?

Polyadenylation occurs at numerous sites within 3'-untranslated regions (3'-UTRs) but rarely within coding regions. How does Pol II travel through long coding regions without generating poly(A) sites, yet then permits promiscuous polyadenylation once it reaches the 3'-UTR? The cleavage/polyadenylation (CpA) machinery preferentially associates with 3'-UTRs, but it is unknown how its recruitment is restricted to 3'-UTRs during Pol II elongation. Unlike coding regions, 3'-UTRs have long AT-rich stretches of DNA that may be important for restricting polyadenylation to 3'-UTRs. Recognition of the 3'-UTR could occur at the DNA (AT-rich), RNA (AU-rich), or RNA:DNA hybrid (rU:dA- and/or rA:dT-rich) level. Based on the nucleic acid critical for 3'-UTR recognition, there are three classes of models, not mutually exclusive, for how the CpA machinery is selectively recruited to 3'-UTRs, thereby restricting where polyadenylation occurs: (1) RNA-based models suggest that the CpA complex directly (or indirectly through one or more intermediary proteins) binds long AU-rich stretches that are exposed after Pol II passes through these regions. (2) DNA-based models suggest that the AT-rich sequence affects nucleosome depletion or the elongating Pol II machinery, resulting in dissociation of some elongation factors and subsequent recruitment of the CpA machinery. (3) RNA:DNA hybrid models suggest that preferential destabilization of the Pol II elongation complex at rU:dA- and/or rA:dT-rich duplexes bridging the nucleotide addition and RNA exit sites permits preferential association of the CpA machinery with 3'-UTRs. Experiments to provide evidence for one or more of these models are suggested.

The fitness trade-off between growth and stress resistance determines the phenotypic landscape

BACKGROUND

A central challenge in biology is to discover a principle that determines individual phenotypic differences within a species. The growth rate is particularly important for a unicellular organism, and the growth rate under a certain condition is negatively associated with that of another condition, termed fitness trade-off. Therefore, there should exist a common molecular mechanism that regulates multiple growth rates under various conditions, but most studies so far have focused on discovering those genes associated with growth rates under a specific condition.

RESULTS

In this study, we found that there exists a recurrent gene expression signature whose expression levels are related to the fitness trade-off between growth preference and stress resistance across various yeast strains and multiple conditions. We further found that the genomic variation of stress-response, ribosomal, and cell cycle regulators are potential causal genes that determine the sensitivity between growth and survival. Intriguingly, we further observed that the same principle holds for human cells using anticancer drug sensitivities across multiple cancer cell lines.

CONCLUSIONS

Together, we suggest that the fitness trade-off is an evolutionary trait that determines individual growth phenotype within a species. By using this trait, we can possibly overcome anticancer drug resistance in cancer cells.

A CUG codon-adapted anchor-away toolkit for functional analysis of genes in Candida albicans

Promoter shutoff of essential genes in the diploid Candida albicans has often been insufficient to create tight, conditional null alleles due to leaky expression and has been a stumbling block in pathogenesis research. Moreover, homozygous deletion of non-essential genes has often been problematic due to the frequent aneuploidy in the mutant strains. Rapid, conditional depletion of essential genes by the anchor-away strategy has been successfully employed in Saccharomyces cerevisiae and other model organisms. Here, rapamycin mediates the dimerization of human FK506-binding protein (FKBP12) and FKBP12-rapamycin-binding (FRB) domain-containing target protein, resulting in relocalization to altered sub-cellular locations. In this work, we used the ribosomal protein Rpl13 as the anchor and took two nuclear proteins as targets to construct a set of mutants in a proof-of-principle approach. We first constructed a rapamycin-resistant C. albicans strain by introducing a dominant mutation in the CaTOR1 gene and a homozygous deletion of RBP1, the ortholog of FKBP12, a primary target of rapamycin. The FKBP12 and the FRB coding sequences were then CUG codon-adapted for C. albicans by site-directed mutagenesis. Anchor-away strains expressing the essential TBP1 gene or the non-essential SPT8 gene as FRB fusions were constructed. We found that rapamycin caused rapid cessation of growth of the TBP-AA strain within 15 minutes and the SPT8-AA strain phenocopied the constitutive filamentous phenotype of the spt8Δ/spt8Δ mutant. Thus, the anchor-away toolbox for C. albicans developed here can be employed for genome-wide analysis to identify gene function in a rapid and reliable manner, further accelerating anti-fungal drug development in C. albicans.

IMPORTANCE

Molecular genetic studies thus far have identified ~27% open-reading frames as being essential for the vegetative growth of Candida albicans in rich medium out of a total 6,198 haploid set of open reading frames. However, a major limitation has been to construct rapid conditional alleles of essential C. albicans genes with near quantitative depletion of encoded proteins. Here, we have developed a toolbox for rapid and conditional depletion of genes that would aid studies of gene function of both essential and non-essential genes.

Structured 3' UTRs destabilize mRNAs in plants

BACKGROUND

RNA secondary structure (RSS) can influence the regulation of transcription, RNA processing, and protein synthesis, among other processes. 3' untranslated regions (3' UTRs) of mRNA also hold the key for many aspects of gene regulation. However, there are often contradictory results regarding the roles of RSS in 3' UTRs in gene expression in different organisms and/or contexts.

RESULTS

Here, we incidentally observe that the primary substrate of miR159a (pri-miR159a), when embedded in a 3' UTR, could promote mRNA accumulation. The enhanced expression is attributed to the earlier polyadenylation of the transcript within the hybrid pri-miR159a-3' UTR and, resultantly, a poorly structured 3' UTR. RNA decay assays indicate that poorly structured 3' UTRs could promote mRNA stability, whereas highly structured 3' UTRs destabilize mRNA in vivo. Genome-wide DMS-MaPseq also reveals the prevailing inverse relationship between 3' UTRs' RSS and transcript accumulation in the transcriptomes of Arabidopsis, rice, and even human. Mechanistically, transcripts with highly structured 3' UTRs are preferentially degraded by 3'-5' exoribonuclease SOV and 5'-3' exoribonuclease XRN4, leading to decreased expression in Arabidopsis. Finally, we engineer different structured 3' UTRs to an endogenous FT gene and alter the FT-regulated flowering time in Arabidopsis.

CONCLUSIONS

We conclude that highly structured 3' UTRs typically cause reduced accumulation of the harbored transcripts in Arabidopsis. This pattern extends to rice and even mammals. Furthermore, our study provides a new strategy of engineering the 3' UTRs' RSS to modify plant traits in agricultural production and mRNA stability in biotechnology.

Synthetic Biology Tools for Engineering Aspergillus oryzae

For more than a thousand years, Aspergillus oryzae has been used in traditional culinary industries, including for food fermentation, brewing, and flavoring. In recent years, A. oryzae has been extensively used in deciphering the pathways of natural product synthesis and value-added compound bioproduction. Moreover, it is increasingly being used in modern biotechnology industries, such as for the production of enzymes and recombinant proteins. The investigation of A. oryzae has been significantly accelerated through the successive application of a diverse array of synthetic biology techniques and methodologies. In this review, the advancements in biological tools for the synthesis of A. oryzae, including DNA assembly technologies, gene expression regulatory elements, and genome editing systems, are discussed. Additionally, the challenges associated with the heterologous expression of A. oryzae are addressed.

Alternative Transcripts Diversify Genome Function for Phenome Relevance to Health and Diseases

Manipulation using alternative exon splicing (AES), alternative transcription start (ATS), and alternative polyadenylation (APA) sites are key to transcript diversity underlying health and disease. All three are pervasive in organisms, present in at least 50% of human protein-coding genes. In fact, ATS and APA site use has the highest impact on protein identity, with their ability to alter which first and last exons are utilized as well as impacting stability and translation efficiency. These RNA variants have been shown to be highly specific, both in tissue type and stage, with demonstrated importance to cell proliferation, differentiation and the transition from fetal to adult cells. While alternative exon splicing has a limited effect on protein identity, its ubiquity highlights the importance of these minor alterations, which can alter other features such as localization. The three processes are also highly interwoven, with overlapping, complementary, and competing factors, RNA polymerase II and its CTD (C-terminal domain) chief among them. Their role in development means dysregulation leads to a wide variety of disorders and cancers, with some forms of disease disproportionately affected by specific mechanisms (AES, ATS, or APA). Challenges associated with the genome-wide profiling of RNA variants and their potential solutions are also discussed in this review.

Ceg1 depletion reveals mechanisms governing degradation of non-capped RNAs in Saccharomyces cerevisiae

Most functional eukaryotic mRNAs contain a 5' 7-methylguanosine (m7G) cap. Although capping is essential for many biological processes including mRNA processing, export and translation, the fate of uncapped transcripts has not been studied extensively. Here, we employed fast nuclear depletion of the capping enzymes in Saccharomyces cerevisiae to uncover the turnover of the transcripts that failed to be capped. We show that although the degradation of cap-deficient mRNA is dominant, the levels of hundreds of non-capped mRNAs increase upon depletion of the capping enzymes. Overall, the abundance of non-capped mRNAs is inversely correlated to the expression levels, altogether resembling the effects observed in cells lacking the cytoplasmic 5'-3' exonuclease Xrn1 and indicating differential degradation fates of non-capped mRNAs. The inactivation of the nuclear 5'-3' exonuclease Rat1 does not rescue the non-capped mRNA levels indicating that Rat1 is not involved in their degradation and consequently, the lack of the capping does not affect the distribution of RNA Polymerase II on the chromatin. Our data indicate that the cap presence is essential to initiate the Xrn1-dependent degradation of mRNAs underpinning the role of 5' cap in the Xrn1-dependent buffering of the cellular mRNA levels.

Elongation rate of RNA polymerase II affects pausing patterns across 3' UTRs

Yeast mRNAs are polyadenylated at multiple sites in their 3' untranslated regions (3' UTRs), and poly(A) site usage is regulated by the rate of transcriptional elongation by RNA polymerase II (Pol II). Slow Pol II derivatives favor upstream poly(A) sites, and fast Pol II derivatives favor downstream poly(A) sites. Transcriptional elongation and polyadenylation are linked at the nucleotide level, presumably reflecting Pol II dwell time at each residue that influences the level of polyadenylation. Here, we investigate the effect of Pol II elongation rate on pausing patterns and the relationship between Pol II pause sites and poly(A) sites within 3' UTRs. Mutations that affect Pol II elongation rate alter sequence preferences at pause sites within 3' UTRs, and pausing preferences differ between 3' UTRs and coding regions. In addition, sequences immediately flanking the pause sites show preferences that are largely independent of Pol II speed. In wild-type cells, poly(A) sites are preferentially located < 50 nucleotides upstream from Pol II pause sites, but this spatial relationship is diminished in cells harboring Pol II speed mutants. Based on a random forest classifier, Pol II pause sites are modestly predicted by the distance to poly(A) sites but are better predicted by the chromatin landscape in Pol II speed derivatives. Transcriptional regulatory proteins can influence the relationship between Pol II pausing and polyadenylation but in a manner distinct from Pol II elongation rate derivatives. These results indicate a complex relationship between Pol II pausing and polyadenylation.

The LARP1 homolog Slr1p controls the stability and expression of proto-5'TOP mRNAs in fission yeast

Messenger RNAs (mRNAs) in higher eukaryotes that encode proteins important for the assembly of the translational apparatus (e.g., ribosomal proteins) often harbor a pyrimidine-rich motif at the extreme 5' end known as a 5' terminal oligopyrimidine (5'TOP) sequence. Members of the La-related protein 1 (LARP1) family control 5'TOP expression through a conserved DM15 motif, but the mechanism is not well understood. 5'TOP motifs have not been described in many lower organisms, and fission yeast harbors a LARP1 homolog that also lacks a DM15 motif. In this work, we show that the fission yeast LARP1 homolog, Slr1p, controls the translation and stability of mRNAs encoding proteins analogous to 5'TOP mRNAs in higher eukaryotes, which we thus refer to as proto-5'TOPs. Our data suggest that the LARP1 DM15 motif and the mRNA 5'TOP motif may be features that were scaffolded over a more fundamental mechanism of LARP1-associated control of gene expression.

Delineating yeast cleavage and polyadenylation signals using deep learning

3'-end cleavage and polyadenylation is an essential process for eukaryotic mRNA maturation. In yeast species, the polyadenylation signals that recruit the processing machinery are degenerate and remain poorly characterized compared to well-defined regulatory elements in mammals. Especially, recent deep sequencing experiments showed extensive cleavage heterogeneity for some mRNAs in Saccharomyces cerevisiae and uncovered the polyA motif differences between S. cerevisiae vs. Schizosaccharomyces pombe . The findings raised the fundamental question of how polyadenylation signals are formed in yeast. Here we addressed this question by developing deep learning models to deconvolute degenerate cis -regulatory elements and quantify their positional importance in mediating yeast polyA site formation, cleavage heterogeneity, and strength. In S. cerevisiae , cleavage heterogeneity is promoted by the depletion of U-rich elements around polyA sites as well as multiple occurrences of upstream UA-rich elements. Sites with high cleavage heterogeneity show overall lower strength. The site strength and tandem site distances modulate alternative polyadenylation (APA) under the diauxic stress. Finally, we developed a deep learning model to reveal the distinct motif configuration of S. pombe polyA sites which show more precise cleavage than S. cerevisiae . Altogether, our deep learning models provide unprecedented insights into polyA site formation across yeast species.

Simultaneous estimation of gene regulatory network structure and RNA kinetics from single cell gene expression

Cells respond to environmental and developmental stimuli by remodeling their transcriptomes through regulation of both mRNA transcription and mRNA decay. A central goal of biology is identifying the global set of regulatory relationships between factors that control mRNA production and degradation and their target transcripts and construct a predictive model of gene expression. Regulatory relationships are typically identified using transcriptome measurements and causal inference algorithms. RNA kinetic parameters are determined experimentally by employing run-on or metabolic labeling (e.g. 4-thiouracil) methods that allow transcription and decay rates to be separately measured. Here, we develop a deep learning model, trained with single-cell RNA-seq data, that both infers causal regulatory relationships and estimates RNA kinetic parameters. The resulting in silico model predicts future gene expression states and can be perturbed to simulate the effect of transcription factor changes. We acquired model training data by sequencing the transcriptomes of 175,000 individual Saccharomyces cerevisiae cells that were subject to an external perturbation and continuously sampled over a one hour period. The rate of change for each transcript was calculated on a per-cell basis to estimate RNA velocity. We then trained a deep learning model with transcriptome and RNA velocity data to calculate time-dependent estimates of mRNA production and decay rates. By separating RNA velocity into transcription and decay rates, we show that rapamycin treatment causes existing ribosomal protein transcripts to be rapidly destabilized, while production of new transcripts gradually slows over the course of an hour. The neural network framework we present is designed to explicitly model causal regulatory relationships between transcription factors and their genes, and shows superior performance to existing models on the basis of recovery of known regulatory relationships. We validated the predictive power of the model by perturbing transcription factors in silico and comparing transcriptome-wide effects with experimental data. Our study represents the first step in constructing a complete, predictive, biophysical model of gene expression regulation.

Destabilized 3'UTR elements therapeutically degrade ERBB2 mRNA in drug-resistant ERBB2+ cancer models

Breast, lung, and colorectal cancer resistance to molecular targeted therapy is a major challenge that unfavorably impacts clinical outcomes leading to hundreds of thousands of deaths annually. In ERBB2+ cancers regardless of the tissue of origin, many ERBB2+ cancers are resistant to ERBB2-targeted therapy. We discovered that ERBB2+ cancer cells are enriched with poly U sequences on their 3'UTR which are mRNA-stabilizing sequences. We developed a novel technology, in which we engineered these ERBB2 mRNA-stabilizing sequences to unstable forms that successfully overwrote and outcompeted the endogenous ERBB2 mRNA-encoded message and degraded ERBB2 transcripts which led to the loss of the protein across multiple cancer cell types both in the wildtype and drug-resistance settings in vitro and in vivo, offering a unique safe novel modality to control ERBB2 mRNA and other pervasive oncogenic signals where current targeted therapies fail.

Oligonucleotide mapping via mass spectrometry to enable comprehensive primary structure characterization of an mRNA vaccine against SARS-CoV-2

Oligonucleotide mapping via liquid chromatography with UV detection coupled to tandem mass spectrometry (LC-UV-MS/MS) was recently developed to support development of Comirnaty, the world's first commercial mRNA vaccine which immunizes against the SARS-CoV-2 virus. Analogous to peptide mapping of therapeutic protein modalities, oligonucleotide mapping described here provides direct primary structure characterization of mRNA, through enzymatic digestion, accurate mass determinations, and optimized collisionally-induced fragmentation. Sample preparation for oligonucleotide mapping is a rapid, one-pot, one-enzyme digestion. The digest is analyzed via LC-MS/MS with an extended gradient and resulting data analysis employs semi-automated software. In a single method, oligonucleotide mapping readouts include a highly reproducible and completely annotated UV chromatogram with 100% maximum sequence coverage, and a microheterogeneity assessment of 5' terminus capping and 3' terminus poly(A)-tail length. Oligonucleotide mapping was pivotal to ensure the quality, safety, and efficacy of mRNA vaccines by providing: confirmation of construct identity and primary structure and assessment of product comparability following manufacturing process changes. More broadly, this technique may be used to directly interrogate the primary structure of RNA molecules in general.

Evolutionary trade-off and mutational bias could favor transcriptional over translational divergence within paralog pairs

How changes in the different steps of protein synthesis-transcription, translation and degradation-contribute to differences of protein abundance among genes is not fully understood. There is however accumulating evidence that transcriptional divergence might have a prominent role. Here, we show that yeast paralogous genes are more divergent in transcription than in translation. We explore two causal mechanisms for this predominance of transcriptional divergence: an evolutionary trade-off between the precision and economy of gene expression and a larger mutational target size for transcription. Performing simulations within a minimal model of post-duplication evolution, we find that both mechanisms are consistent with the observed divergence patterns. We also investigate how additional properties of the effects of mutations on gene expression, such as their asymmetry and correlation across levels of regulation, can shape the evolution of paralogs. Our results highlight the importance of fully characterizing the distributions of mutational effects on transcription and translation. They also show how general trade-offs in cellular processes and mutation bias can have far-reaching evolutionary impacts.

Maximizing protein production by keeping cells at optimal secretory stress levels using real-time control approaches

Optimizing the production of recombinant proteins is a problem of major industrial and pharmaceutical importance. Secretion of the protein by the host cell considerably simplifies downstream purification processes. However, for many proteins, this is also the limiting production step. Current solutions involve extensive engineering of the chassis cell to facilitate protein trafficking and limit protein degradation triggered by excessive secretion-associated stress. Here, we propose instead a regulation-based strategy in which induction is dynamically adjusted to an optimal strength based on the current stress level of the cells. Using a small collection of hard-to-secrete proteins, a bioreactor-based platform with automated cytometry measurements, and a systematic assay to quantify secreted protein levels, we demonstrate that the secretion sweet spot is indicated by the appearance of a subpopulation of cells that accumulate high amounts of proteins, decrease growth, and face significant stress, that is, experience a secretion burnout. In these cells, adaptations capabilities are overwhelmed by a too strong production. Using these notions, we show for a single-chain antibody variable fragment that secretion levels can be improved by 70% by dynamically keeping the cell population at optimal stress levels using real-time closed-loop control.

Condition-specific 3' mRNA isoform half-lives and stability elements in yeast

Alternative polyadenylation generates numerous 3' mRNA isoforms that can differ in their stability, structure, and function. These isoforms can be used to map mRNA stabilizing and destabilizing elements within 3' untranslated regions (3'UTRs). Here, we examine how environmental conditions affect 3' mRNA isoform turnover and structure in yeast cells on a transcriptome scale. Isoform stability broadly increases when cells grow more slowly, with relative half-lives of most isoforms being well correlated across multiple conditions. Surprisingly, dimethyl sulfate probing reveals that individual 3' isoforms have similar structures across different conditions, in contrast to the extensive structural differences that can exist between closely related isoforms in an individual condition. Unexpectedly, most mRNA stabilizing and destabilizing elements function only in a single growth condition. The genes associated with some classes of condition-specific stability elements are enriched for different functional categories, suggesting that regulated mRNA stability might contribute to adaptation to different growth environments. Condition-specific stability elements do not result in corresponding condition-specific changes in steady-state mRNA isoform levels. This observation is consistent with a compensatory mechanism between polyadenylation and stability, and it suggests that condition-specific mRNA stability elements might largely reflect condition-specific regulation of mRNA 3' end formation.

The emerging roles of CFIm25 (NUDT21/CPSF5) in human biology and disease

The mammalian cleavage factor I subunit CFIm25 (NUDT21) binds to the UGUA sequences of precursor RNAs. Traditionally, CFIm25 is known to facilitate 3' end formation of pre-mRNAs resulting in the formation of polyadenylated transcripts. Recent studies suggest that CFIm25 may be involved in the cyclization and hence generation of circular RNAs (circRNAs) that contain UGUA motifs. These circRNAs act as competing endogenous RNAs (ceRNAs) that disrupt the ceRNA-miRNA-mRNA axis. Other emerging roles of CFIm25 include regulating both alternative splicing and alternative polyadenylation (APA). APA generates different sized transcripts that may code for different proteins, or more commonly transcripts that code for the same protein but differ in the length and sequence content of their 3' UTRs (3' UTR-APA). CFIm25 mediated global changes in 3' UTR-APA affect human physiology including spermatogenesis and the determination of cell fate. Deregulation of CFIm25 and changes in 3' UTR-APA have been implicated in several human diseases including cancer. In many cancers, CFIm25 acts as a tumor suppressor. However, there are some cancers where CFIm25 has the opposite effect. Alterations in CFIm25-driven 3' UTR-APA may also play a role in neural dysfunction and fibrosis. CFIm25 mediated 3' UTR-APA changes can be used to generate specific signatures that can be used as potential biomarkers in development and disease. Due to the emerging role of CFIm25 as a regulator of the aforementioned RNA processing events, modulation of CFIm25 levels may be a novel viable therapeutic approach. This article is categorized under: RNA Processing > 3' End Processing RNA in Disease and Development > RNA in Disease.

Alternative 3'-untranslated regions regulate high-salt tolerance of Spartina alterniflora

High-salt stress continues to challenge the growth and survival of many plants. Alternative polyadenylation (APA) produces mRNAs with different 3'-untranslated regions (3' UTRs) to regulate gene expression at the post-transcriptional level. However, the roles of alternative 3' UTRs in response to salt stress remain elusive. Here, we report the function of alternative 3' UTRs in response to high-salt stress in S. alterniflora (Spartina alterniflora), a monocotyledonous halophyte tolerant of high-salt environments. We found that high-salt stress induced global APA dynamics, and ∼42% of APA genes responded to salt stress. High-salt stress led to 3' UTR lengthening of 207 transcripts through increasing the usage of distal poly(A) sites. Transcripts with alternative 3' UTRs were mainly enriched in salt stress-related ion transporters. Alternative 3' UTRs of HIGH-AFFINITY K+ TRANSPORTER 1 (SaHKT1) increased RNA stability and protein synthesis in vivo. Regulatory AU-rich elements were identified in alternative 3' UTRs, boosting the protein level of SaHKT1. RNAi-knock-down experiments revealed that the biogenesis of 3' UTR lengthening in SaHKT1 was controlled by the poly(A) factor CLEAVAGE AND POLYADENYLATION SPECIFICITY FACTOR 30 (SaCPSF30). Over-expression of SaHKT1 with an alternative 3' UTR in rice (Oryza sativa) protoplasts increased mRNA accumulation of salt-tolerance genes in an AU-rich element-dependent manner. These results suggest that mRNA 3' UTR lengthening is a potential mechanism in response to high-salt stress. These results also reveal complex regulatory roles of alternative 3' UTRs coupling APA and regulatory elements at the post-transcriptional level in plants.

Alternative Polyadenylation Is a Novel Strategy for the Regulation of Gene Expression in Response to Stresses in Plants

Precursor message RNA requires processing to generate mature RNA. Cleavage and polyadenylation at the 3'-end in the maturation of mRNA is one of key processing steps in eukaryotes. The polyadenylation (poly(A)) tail of mRNA is an essential feature that is required to mediate its nuclear export, stability, translation efficiency, and subcellular localization. Most genes have at least two mRNA isoforms via alternative splicing (AS) or alternative polyadenylation (APA), which increases the diversity of transcriptome and proteome. However, most previous studies have focused on the role of alternative splicing on the regulation of gene expression. In this review, we summarize the recent advances concerning APA in the regulation of gene expression and in response to stresses in plants. We also discuss the mechanisms for the regulation of APA for plants in the adaptation to stress responses, and suggest that APA is a novel strategy for the adaptation to environmental changes and response to stresses in plants.

Translation of in vitro-transcribed RNA therapeutics

In vitro transcribed, modified messenger RNAs (IVTmRNAs) have been used to vaccinate billions of individuals against the SARS-CoV-2 virus, and are currently being developed for many additional therapeutic applications. IVTmRNAs must be translated into proteins with therapeutic activity by the same cellular machinery that also translates native endogenous transcripts. However, different genesis pathways and routes of entry into target cells as well as the presence of modified nucleotides mean that the way in which IVTmRNAs engage with the translational machinery, and the efficiency with which they are being translated, differs from native mRNAs. This review summarises our current knowledge of commonalities and differences in translation between IVTmRNAs and cellular mRNAs, which is key for the development of future design strategies that can generate IVTmRNAs with improved activity in therapeutic applications.

Advances of mRNA vaccine in tumor: a maze of opportunities and challenges

High-frequency mutations in tumor genomes could be exploited as an asset for developing tumor vaccines. In recent years, with the tremendous breakthrough in genomics, intelligence algorithm, and in-depth insight of tumor immunology, it has become possible to rapidly target genomic alterations in tumor cell and rationally select vaccine targets. Among a variety of candidate vaccine platforms, the early application of mRNA was limited by instability low efficiency and excessive immunogenicity until the successful development of mRNA vaccines against SARS-COV-2 broken of technical bottleneck in vaccine preparation, allowing tumor mRNA vaccines to be prepared rapidly in an economical way with good performance of stability and efficiency. In this review, we systematically summarized the classification and characteristics of tumor antigens, the general process and methods for screening neoantigens, the strategies of vaccine preparations and advances in clinical trials, as well as presented the main challenges in the current mRNA tumor vaccine development.

Decapping factor Dcp2 controls mRNA abundance and translation to adjust metabolism and filamentation to nutrient availability

Degradation of most yeast mRNAs involves decapping by Dcp1/Dcp2. DEAD-box protein Dhh1 has been implicated as an activator of decapping, in coupling codon non-optimality to enhanced degradation, and as a translational repressor, but its functions in cells are incompletely understood. RNA-Seq analyses coupled with CAGE sequencing of all capped mRNAs revealed increased abundance of hundreds of mRNAs in dcp2 Δ cells that appears to result directly from impaired decapping rather than elevated transcription, which was confirmed by ChIP-Seq analysis of RNA Polymerase II occupancies genome-wide. Interestingly, only a subset of mRNAs requires Dhh1 for targeting by Dcp2, and also generally requires the other decapping activators Pat1, Lsm2, Edc3 or Scd6; whereas most of the remaining transcripts utilize NMD factors for Dcp2-mediated turnover. Neither inefficient translation initiation nor stalled elongation appears to be a major driver of Dhh1-enhanced mRNA degradation. Surprisingly, ribosome profiling revealed that dcp2 Δ confers widespread changes in relative TEs that generally favor well-translated mRNAs. Because ribosome biogenesis is reduced while capped mRNA abundance is increased by dcp2 Δ, we propose that an increased ratio of mRNA to ribosomes increases competition among mRNAs for limiting ribosomes to favor efficiently translated mRNAs in dcp2 Δ cells. Interestingly, genes involved in respiration or utilization of alternative carbon or nitrogen sources are derepressed, and both mitochondrial function and cell filamentation (a strategy for nutrient foraging) are elevated by dcp2 Δ, suggesting that mRNA decapping sculpts gene expression post-transcriptionally to fine-tune metabolic pathways and morphological transitions according to nutrient availability.

Identifying transcript-level differential expression in primary human immune cells

BACKGROUND

Multipotential hematopoietic stem cells differentiate into a wide variety of immune cells with a diversity of functions, including the ability to respond to a variety of stimuli. Importantly, numerous studies have demonstrated the importance of gene transcription in defining cell identity and functions. While these studies have primarily been performed at the level of the gene, it is known that key immune genes such as CD44 and CD45 generate multiple different transcripts that are differentially expressed across different immune cells, and that encode proteins with different sequences and functions. Prior genomic surveys have shown that the mechanisms for generating diversity in expressed transcripts (alternate splicing, alternate transcription start sites, etc.) are very active in immune cells, but have been lacking in terms of identifying genes with multiple transcripts, that are differentially expressed, and likely to affect cell functions.

METHODS

We first identified the set of genes that had at least two transcripts expressed in our RNA sequencing dataset generated from purified populations of neutrophils, monocytes and five lymphocyte populations (B, NK, γδ T, CD4 + T and CD8 + T) from twelve healthy donors. Next, we developed a heuristic approach to identify genes where two or more transcripts have distinct expression patterns across lymphoid and/or myeloid populations. We then focused our annotation and interpretation on differentially expressed transcripts that affect the coding sequence. This process was repeated to identify transcripts that were differentially expressed between monocytes and populations of macrophages and LPS-stimulated macrophages derived from these monocytes in vitro.

RESULTS

We found that over 55 % of genes had two or more expressed transcripts, with an average ∼3 transcripts per gene, and that 70 % of these had at least two of the transcripts that encoded proteins with different sequences. As expected, we identified a complex pattern of differential expression for multiple transcripts encoding the CD45 transmembrane protein, but we also found similar evidence for ten other genes (CD300A, FYB1, GPI, LITAF, PSMA1, PTMA, RPL32, SEPTIN9, SH3BP2, SH3KBP1) when comparing the expression patterns of transcripts within myeloid and lymphoid cells. We also identified five genes with differentially expressed transcripts associated with the transition from monocytes to macrophages (FNBP1, KLF6, and SEPTIN9) or between macrophages and LPS-stimulated macrophages (CD44, OAZ2, and SEPTIN9). For the most part, we found that the different transcripts of these genes are expected to impact specific biological functions, for example the different transcripts of SEPTIN9 likely regulate the cytoskeleton in immune cells via their interactions with actins filaments and microtubules.

CONCLUSIONS

This analytic approach successfully identified multi-transcript genes that are differentially expressed across immune cells and could be applied to other transcriptomic data.

DATA AVAILABILITY STATEMENT

Researchers can request access to the individual-level data from the current study by contacting the Montreal Heart Institute ethics committee at the following institutional email address: cer.icm@icm-mhi.org.

Transcription factor binding process is the primary driver of noise in gene expression

Noise in expression of individual genes gives rise to variations in activity of cellular pathways and generates heterogeneity in cellular phenotypes. Phenotypic heterogeneity has important implications for antibiotic persistence, mutation penetrance, cancer growth and therapy resistance. Specific molecular features such as the presence of the TATA box sequence and the promoter nucleosome occupancy have been associated with noise. However, the relative importance of these features in noise regulation is unclear and how well these features can predict noise has not yet been assessed. Here through an integrated statistical model of gene expression noise in yeast we found that the number of regulating transcription factors (TFs) of a gene was a key predictor of noise, whereas presence of the TATA box and the promoter nucleosome occupancy had poor predictive power. With an increase in the number of regulatory TFs, there was a rise in the number of cooperatively binding TFs. In addition, an increased number of regulatory TFs meant more overlaps in TF binding sites, resulting in competition between TFs for binding to the same region of the promoter. Through modeling of TF binding to promoter and application of stochastic simulations, we demonstrated that competition and cooperation among TFs could increase noise. Thus, our work uncovers a process of noise regulation that arises out of the dynamics of gene regulation and is not dependent on any specific transcription factor or specific promoter sequence.

Expression Profile of Sphingosine Kinase 1 Isoforms in Human Cancer Tissues and Cells: Importance and Clinical Relevance of the Neglected 1b-Isoform

BACKGROUND

Overexpression of sphingosine kinase 1 (SphK1) is casually associated with many types of cancer, and inhibitors of SphK1 sensitize tumors to chemotherapy. SphK1 is expressed as two major isoforms, SphK1a and SphK1b. To date, no information has been reported on the SphK1 isoform expression profile and its clinical relevance.

OBJECTIVE

The objective is to examine the expression profile of the SphK1a and SPhK1b isoforms in human cancer and noncancer tissues and cell lines and explore their clinical relevance.

METHODS

We used PCR to qualitatively examine the expression profile of these two isoforms in breast, liver, and prostate cancer tissues plus paired adjacent tissues and in 11 cancer and normal cell lines (breast, cervical, bone, prostate, colon, brain, mesothelioma tumor and benign, and human kidney cells).

RESULTS

We found that SphK1a was ubiquitously expressed in all cancer cells and tissues tested; in contrast, SphK1b was only expressed in selective cell types in breast, prostate, and lung cancer.

CONCLUSIONS

Our data suggest that SphK1a is important for generic SphK1/S1P functions, and SphK1b mediates specialized and/or unique pathways in a specific type of tissue and could be a biomarker for cancer. This discovery is important for future SphK1-related cancer research and may have clinical implications in drug development associated with SphK1-directed cancer treatment.

Nuclear mRNA metabolism drives selective basket assembly on a subset of nuclear pore complexes in budding yeast

To determine which transcripts should reach the cytoplasm for translation, eukaryotic cells have established mechanisms to regulate selective mRNA export through the nuclear pore complex (NPC). The nuclear basket, a substructure of the NPC protruding into the nucleoplasm, is thought to function as a stable platform where mRNA-protein complexes (mRNPs) are rearranged and undergo quality control prior to export, ensuring that only mature mRNAs reach the cytoplasm. Here, we use proteomic, genetic, live-cell, and single-molecule resolution microscopy approaches in budding yeast to demonstrate that basket formation is dependent on RNA polymerase II transcription and subsequent mRNP processing. We further show that while all NPCs can bind Mlp1, baskets assemble only on a subset of nucleoplasmic NPCs, and these basket-containing NPCs associate a distinct protein and RNA interactome. Taken together, our data point toward NPC heterogeneity and an RNA-dependent mechanism for functionalization of NPCs in budding yeast through nuclear basket assembly.

3' Untranslated Regions Are Modular Entities That Determine Polyadenylation Profiles

The 3' ends of eukaryotic mRNAs are generated by cleavage of nascent transcripts followed by polyadenylation, which occurs at numerous sites within 3' untranslated regions (3' UTRs) but rarely within coding regions. An individual gene can yield many 3'-mRNA isoforms with distinct half-lives. We dissect the relative contributions of protein-coding sequences (open reading frames [ORFs]) and 3' UTRs to polyadenylation profiles in yeast. ORF-deleted derivatives often display strongly decreased mRNA levels, indicating that ORFs contribute to overall mRNA stability. Poly(A) profiles, and hence relative isoform half-lives, of most (9 of 10) ORF-deleted derivatives are very similar to their wild-type counterparts. Similarly, in-frame insertion of a large protein-coding fragment between the ORF and 3' UTR has minimal effect on the poly(A) profile in all 15 cases tested. Last, reciprocal ORF/3'-UTR chimeric genes indicate that the poly(A) profile is determined by the 3' UTR. Thus, 3' UTRs are self-contained modular entities sufficient to determine poly(A) profiles and relative 3'-isoform half-lives. In the one atypical instance, ORF deletion causes an upstream shift of poly(A) sites, likely because juxtaposition of an unusually high AT-rich stretch directs polyadenylation closely downstream. This suggests that long AT-rich stretches, which are not encountered until after coding regions, are important for restricting polyadenylation to 3' UTRs.

Mitotic checkpoint gene expression is tuned by codon usage bias

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

Global SLAM-seq for accurate mRNA decay determination and identification of NMD targets

Gene expression analysis requires accurate measurements of global RNA degradation rates, earlier problematic with methods disruptive to cell physiology. Recently, metabolic RNA labeling emerged as an efficient and minimally invasive technique applied in mammalian cells. Here, we have adapted SH-linked alkylation for the metabolic sequencing of RNA (SLAM-seq) for a global mRNA stability study in yeast using 4-thiouracil pulse-chase labeling. We assign high-confidence half-life estimates for 67.5% of expressed ORFs, and measure a median half-life of 9.4 min. For mRNAs where half-life estimates exist in the literature, their ranking order was in good agreement with previous data, indicating that SLAM-seq efficiently classifies stable and unstable transcripts. We then leveraged our yeast protocol to identify targets of the nonsense-mediated decay (NMD) pathway by measuring the change in RNA half-lives, instead of steady-state RNA level changes. With SLAM-seq, we assign 580 transcripts as putative NMD targets, based on their measured half-lives in wild-type and upf3Δ mutants. We find 225 novel targets, and observe a strong agreement with previous reports of NMD targets, 61.2% of our candidates being identified in previous studies. This indicates that SLAM-seq is a simpler and more economic method for global quantification of mRNA half-lives. Our adaptation for yeast yielded global quantitative measures of the NMD effect on transcript half-lives, high correlation with RNA half-lives measured previously with more technically challenging protocols, and identification of novel NMD regulated transcripts that escaped prior detection.

Secondary structures in RNA synthesis, splicing and translation

Even though the functional role of mRNA molecules is primarily decided by the nucleotide sequence, several properties are determined by secondary structure conformations. Examples of secondary structures include long range interactions, hairpins, R-loops and G-quadruplexes and they are formed through interactions of non-adjacent nucleotides. Here, we discuss advances in our understanding of how secondary structures can impact RNA synthesis, splicing, translation and mRNA half-life. During RNA synthesis, secondary structures determine RNA polymerase II (RNAPII) speed, thereby influencing splicing. Splicing is also determined by RNA binding proteins and their binding rates are modulated by secondary structures. For the initiation of translation, secondary structures can control the choice of translation start site. Here, we highlight the mechanisms by which secondary structures modulate these processes, discuss advances in technologies to detect and study them systematically, and consider the roles of RNA secondary structures in disease.

Modulation of RNA stability regulates gene expression in two opposite ways: through buffering of RNA levels upon global perturbations and by supporting adapted differential expression

The steady state levels of RNAs, often referred to as expression levels, result from a well-balanced combination of RNA transcription and decay. Alterations in RNA levels will therefore result from tight regulation of transcription rates, decay rates or both. Here, we explore the role of RNA stability in achieving balanced gene expression and present genome-wide RNA stabilities in Drosophila melanogaster male and female cells as well as male cells depleted of proteins essential for dosage compensation. We identify two distinct RNA-stability mediated responses involved in regulation of gene expression. The first of these responds to acute and global changes in transcription and thus counteracts potentially harmful gene mis-expression by shifting the RNA stability in the direction opposite to the transcriptional change. The second response enhances inter-individual differential gene expression by adjusting the RNA stability in the same direction as a transcriptional change. Both mechanisms are global, act on housekeeping as well as non-housekeeping genes and were observed in both flies and mammals. Additionally, we show that, in contrast to mammals, modulation of RNA stability does not detectably contribute to dosage compensation of the sex-chromosomes in D. melanogaster.

Probing RNA structures and functions by solvent accessibility: an overview from experimental and computational perspectives

Characterizing RNA structures and functions have mostly been focused on 2D, secondary and 3D, tertiary structures. Recent advances in experimental and computational techniques for probing or predicting RNA solvent accessibility make this 1D representation of tertiary structures an increasingly attractive feature to explore. Here, we provide a survey of these recent developments, which indicate the emergence of solvent accessibility as a simple 1D property, adding to secondary and tertiary structures for investigating complex structure–function relations of RNAs.

Improved alpharetrovirus-based Gag.MS2 particles for efficient and transient delivery of CRISPR-Cas9 into target cells

DNA-modifying technologies, such as the CRISPR-Cas9 system, are promising tools in the field of gene and cell therapies. However, high and prolonged expression of DNA-modifying enzymes may cause cytotoxic and genotoxic side effects and is therefore unwanted in therapeutic approaches. Consequently, development of new and potent short-term delivery methods is of utmost importance. Recently, we developed non-integrating gammaretrovirus- and MS2 bacteriophage-based Gag.MS2 (g.Gag.MS2) particles for transient transfer of non-retroviral CRISPR-Cas9 RNA into target cells. In the present study, we further improved the technique by transferring the system to the alpharetroviral vector platform (a.Gag.MS2), which significantly increased CRISPR-Cas9 delivery into target cells and allowed efficient targeted knockout of endogenous TP53/Trp53 genes in primary murine fibroblasts as well as primary human fibroblasts, hepatocytes, and cord-blood-derived CD34⁺ stem and progenitor cells. Strikingly, co-packaging of Cas9 mRNA and multiple single guide RNAs (sgRNAs) into a.Gag.MS2 chimera displayed efficient targeted knockout of up to three genes. Co-transfection of single-stranded DNA donor oligonucleotides during CRISPR-Cas9 particle production generated all-in-one particles, which mediated up to 12.5% of homology-directed repair in primary cell cultures. In summary, optimized a.Gag.MS2 particles represent a versatile tool for short-term delivery of DNA-modifying enzymes into a variety of target cells, including primary murine and human cells.

Engineered pegRNAs improve prime editing efficiency

Prime editing enables the installation of virtually any combination of point mutations, small insertions or small deletions in the DNA of living cells. A prime editing guide RNA (pegRNA) directs the prime editor protein to the targeted locus and also encodes the desired edit. Here we show that degradation of the 3' region of the pegRNA that contains the reverse transcriptase template and the primer binding site can poison the activity of prime editing systems, impeding editing efficiency. We incorporated structured RNA motifs to the 3' terminus of pegRNAs that enhance their stability and prevent degradation of the 3' extension. The resulting engineered pegRNAs (epegRNAs) improve prime editing efficiency 3-4-fold in HeLa, U2OS and K562 cells and in primary human fibroblasts without increasing off-target editing activity. We optimized the choice of 3' structural motif and developed pegLIT, a computational tool to identify non-interfering nucleotide linkers between pegRNAs and 3' motifs. Finally, we showed that epegRNAs enhance the efficiency of the installation or correction of disease-relevant mutations.

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.

A compensatory link between cleavage/polyadenylation and mRNA turnover regulates steady-state mRNA levels in yeast

Cells have compensatory mechanisms to coordinate the rates of major biological processes, thereby permitting growth in a wide variety of conditions. Here, we uncover a compensatory link between cleavage/polyadenylation in the nucleus and messenger RNA (mRNA) turnover in the cytoplasm. On a global basis, same-gene 3' mRNA isoforms with twofold or greater differences in half-lives have steady-state mRNA levels that differ by significantly less than a factor of 2. In addition, increased efficiency of cleavage/polyadenylation at a specific site is associated with reduced stability of the corresponding 3' mRNA isoform. This inverse relationship between cleavage/polyadenylation and mRNA isoform half-life reduces the variability in the steady-state levels of mRNA isoforms, and it occurs in all four growth conditions tested. These observations suggest that during cleavage/polyadenylation in the nucleus, mRNA isoforms are marked in a manner that persists upon translocation to the cytoplasm and affects the activity of mRNA degradation machinery, thus influencing mRNA stability.

An Inducible System for Silencing Establishment Reveals a Stepwise Mechanism in Which Anchoring at the Nuclear Periphery Precedes Heterochromatin Formation

In eukaryotic cells, silent chromatin is mainly found at the nuclear periphery forming subnuclear compartments that favor silencing establishment. Here, we set up an inducible system to monitor silencing establishment at an ectopic locus in relation with its subnuclear localization in budding yeast. We previously showed that introducing LacI bound lacO arrays in proximity to gene flanked by HML silencers favors the recruitment of the yeast silencing complex SIR at this locus, leading to its silencing and anchoring at the nuclear periphery. Using an inducible version of this system, we show that silencing establishment is a stepwise process occurring over several cell cycles, with the progressive recruitment of the SIR complex. In contrast, we observed a rapid, SIR-independent perinuclear anchoring, induced by the high amount of LacI binding at the lacO array leading to nucleosome eviction at this array and to the phosphorylation of H2A in the neighboring nucleosomes by Mec1 kinase. While the initial phosphorylation of H2A (H2A-P) and perinuclear anchoring are independent of the SIR complex, its latter recruitment stabilizes H2A-P and reinforces the perinuclear anchoring. Finally, we showed that Sir3 spreading stabilizes nucleosomes and limits the access of specific DNA-binding protein to DNA.

Effects of sequence motifs in the yeast 3' untranslated region determined from massively parallel assays of random sequences

BACKGROUND

The 3' untranslated region (UTR) plays critical roles in determining the level of gene expression through effects on activities such as mRNA stability and translation. Functional elements within this region have largely been identified through analyses of native genes, which contain multiple co-evolved sequence features.

RESULTS

To explore the effects of 3' UTR sequence elements outside of native sequence contexts, we analyze hundreds of thousands of random 50-mers inserted into the 3' UTR of a reporter gene in the yeast Saccharomyces cerevisiae. We determine relative protein expression levels from the fitness of transformants in a growth selection. We find that the consensus 3' UTR efficiency element significantly boosts expression, independent of sequence context; on the other hand, the consensus positioning element has only a small effect on expression. Some sequence motifs that are binding sites for Puf proteins substantially increase expression in the library, despite these proteins generally being associated with post-transcriptional downregulation of native mRNAs. Our measurements also allow a systematic examination of the effects of point mutations within efficiency element motifs across diverse sequence backgrounds. These mutational scans reveal the relative in vivo importance of individual bases in the efficiency element, which likely reflects their roles in binding the Hrp1 protein involved in cleavage and polyadenylation.

CONCLUSIONS

The regulatory effects of some 3' UTR sequence features, like the efficiency element, are consistent regardless of sequence context. In contrast, the consequences of other 3' UTR features appear to be strongly dependent on their evolved context within native genes.

Replication stress inhibits synthesis of histone mRNAs in yeast by removing Spt10p and Spt21p from the histone promoters

Proliferating cells coordinate histone and DNA synthesis to maintain correct stoichiometry for chromatin assembly. Histone mRNA levels must be repressed when DNA replication is inhibited to prevent toxicity and genome instability due to free non-chromatinized histone proteins. In mammalian cells, replication stress triggers degradation of histone mRNAs, but it is unclear if this mechanism is conserved from other species. The aim of this study was to identify the histone mRNA decay pathway in the yeast Saccharomyces cerevisiae and determine the mechanism by which DNA replication stress represses histone mRNAs. Using reverse transcription-quantitative PCR and chromatin immunoprecipitation–quantitative PCR, we show here that histone mRNAs can be degraded by both 5′ → 3′ and 3′ → 5′ pathways; however, replication stress does not trigger decay of histone mRNA in yeast. Rather, replication stress inhibits transcription of histone genes by removing the histone gene–specific transcription factors Spt10p and Spt21p from histone promoters, leading to disassembly of the preinitiation complexes and eviction of RNA Pol II from histone genes by a mechanism facilitated by checkpoint kinase Rad53p and histone chaperone Asf1p. In contrast, replication stress does not remove SCB-binding factor transcription complex, another activator of histone genes, from the histone promoters, suggesting that Spt10p and Spt21p have unique roles in the transcriptional downregulation of histone genes during replication stress. Together, our data show that, unlike in mammalian cells, replication stress in yeast does not trigger decay of histone mRNAs but inhibits histone transcription.

Quasi-essentiality of RNase Y in Bacillus subtilis is caused by its critical role in the control of mRNA homeostasis

RNA turnover is essential in all domains of life. The endonuclease RNase Y (rny) is one of the key components involved in RNA metabolism of the model organism Bacillus subtilis. Essentiality of RNase Y has been a matter of discussion, since deletion of the rny gene is possible, but leads to severe phenotypic effects. In this work, we demonstrate that the rny mutant strain rapidly evolves suppressor mutations to at least partially alleviate these defects. All suppressor mutants had acquired a duplication of an about 60 kb long genomic region encompassing genes for all three core subunits of the RNA polymerase—α, β, β′. When the duplication of the RNA polymerase genes was prevented by relocation of the rpoA gene in the B. subtilis genome, all suppressor mutants carried distinct single point mutations in evolutionary conserved regions of genes coding either for the β or β’ subunits of the RNA polymerase that were not tolerated by wild type bacteria. In vitro transcription assays with the mutated polymerase variants showed a severe decrease in transcription efficiency. Altogether, our results suggest a tight cooperation between RNase Y and the RNA polymerase to establish an optimal RNA homeostasis in B. subtilis cells.

From Yeast to Mammals, the Nonsense-Mediated mRNA Decay as a Master Regulator of Long Non-Coding RNAs Functional Trajectory

The Nonsense-Mediated mRNA Decay (NMD) has been classically viewed as a translation-dependent RNA surveillance pathway degrading aberrant mRNAs containing premature stop codons. However, it is now clear that mRNA quality control represents only one face of the multiple functions of NMD. Indeed, NMD also regulates the physiological expression of normal mRNAs, and more surprisingly, of long non-coding (lnc)RNAs. Here, we review the different mechanisms of NMD activation in yeast and mammals, and we discuss the molecular bases of the NMD sensitivity of lncRNAs, considering the functional roles of NMD and of translation in the metabolism of these transcripts. In this regard, we describe several examples of functional micropeptides produced from lncRNAs. We propose that translation and NMD provide potent means to regulate the expression of lncRNAs, which might be critical for the cell to respond to environmental changes.

Thermostability, Tunability, and Tenacity of RNA as Rubbery Anionic Polymeric Materials in Nanotechnology and Nanomedicine-Specific Cancer Targeting with Undetectable Toxicity

RNA nanotechnology is the bottom-up self-assembly of nanometer-scale architectures, resembling LEGOs, composed mainly of RNA. The ideal building material should be (1) versatile and controllable in shape and stoichiometry, (2) spontaneously self-assemble, and (3) thermodynamically, chemically, and enzymatically stable with a long shelf life. RNA building blocks exhibit each of the above. RNA is a polynucleic acid, making it a polymer, and its negative-charge prevents nonspecific binding to negatively charged cell membranes. The thermostability makes it suitable for logic gates, resistive memory, sensor set-ups, and NEM devices. RNA can be designed and manipulated with a level of simplicity of DNA while displaying versatile structure and enzyme activity of proteins. RNA can fold into single-stranded loops or bulges to serve as mounting dovetails for intermolecular or domain interactions without external linking dowels. RNA nanoparticles display rubber- and amoeba-like properties and are stretchable and shrinkable through multiple repeats, leading to enhanced tumor targeting and fast renal excretion to reduce toxicities. It was predicted in 2014 that RNA would be the third milestone in pharmaceutical drug development. The recent approval of several RNA drugs and COVID-19 mRNA vaccines by FDA suggests that this milestone is being realized. Here, we review the unique properties of RNA nanotechnology, summarize its recent advancements, describe its distinct attributes inside or outside the body and discuss potential applications in nanotechnology, medicine, and material science.

TRENDseq-A highly multiplexed high throughput RNA 3' end sequencing for mapping alternative polyadenylation

Alternative polyadenylation (APA) is a widespread and highly dynamic mechanism of gene regulation. It affects more than 70% of all genes, resulting in transcript isoforms with distinct 3' end termini. APA thereby considerably expands the diversity of the transcriptome 3' end (TREND). This leads to mRNA isoforms with profoundly different physiological effects, by affecting protein output, production of distinct protein isoforms, or modulating protein localization. APA is globally regulated in various conditions, including developmental and adaptive programs. Since perturbations of APA can disrupt biological processes, ultimately resulting in most devastating disorders, querying the APA landscape is crucial to decipher underlying mechanisms, resulting consequences and potential diagnostic and therapeutic implications. Here we provide a detailed step-by-step protocol for TRENDseq, a method for transcriptome-wide high-throughput sequencing of polyadenylated RNA 3' ends in a highly multiplexed fashion. TRENDseq exploits linear amplification of the starting material to improve sensitivity while significantly reducing the amount of input material. It thereby represents a powerful tool to study APA in numerous experimental set-ups and/or limited human samples in a highly multiplexed and reproducible manner.