Quantifying 3'UTR length from scRNA-seq data reveals changes independent of gene expression.

Although more than half of all genes generate transcripts that differ in 3'UTR length, current analysis pipelines only quantify the amount but not the length of mRNA transcripts. 3'UTR length is determined by 3' end cleavage sites (CS). We map CS in more than 200 primary human and mouse cell types and increase CS annotations relative to the GENCODE database by 40%. Approximately half of all CS are used in few cell types, revealing that most genes only have one or two major 3' ends. We incorporate the CS annotations into a computational pipeline, called scUTRquant, for rapid, accurate, and simultaneous quantification of gene and 3'UTR isoform expression from single-cell RNA sequencing (scRNA-seq) data. When applying scUTRquant to data from 474 cell types and 2134 perturbations, we discover extensive 3'UTR length changes across cell types that are as widespread and coordinately regulated as gene expression changes but affect mostly different genes. Our data indicate that mRNA abundance and mRNA length are two largely independent axes of gene regulation that together determine the amount and spatial organization of protein synthesis.

Subcytoplasmic location of translation controls protein output.

The cytoplasm is highly compartmentalized, but the extent and consequences of subcytoplasmic mRNA localization in non-polarized cells are largely unknown. We determined mRNA enrichment in TIS granules (TGs) and the rough endoplasmic reticulum (ER) through particle sorting and isolated cytosolic mRNAs by digitonin extraction. When focusing on genes that encode non-membrane proteins, we observed that 52% have transcripts enriched in specific compartments. Compartment enrichment correlates with a combinatorial code based on mRNA length, exon length, and 3' UTR-bound RNA-binding proteins. Compartment-biased mRNAs differ in the functional classes of their encoded proteins: TG-enriched mRNAs encode low-abundance proteins with strong enrichment of transcription factors, whereas ER-enriched mRNAs encode large and highly expressed proteins. Compartment localization is an important determinant of mRNA and protein abundance, which is supported by reporter experiments showing that redirecting cytosolic mRNAs to the ER increases their protein expression. In summary, the cytoplasm is functionally compartmentalized by local translation environments.

The FXR1 network acts as signaling scaffold for actomyosin remodeling.

It is currently not known that mRNAs fulfill structural roles in the cytoplasm. Here, we report the FXR1 network, an mRNA-protein (mRNP) network present throughout the cytoplasm: FXR1 packages exceptionally long mRNAs that serve as an underlying network scaffold and concentrate FXR1 molecules, which have multiple protein binding sites. The proximity of FXR1 molecules makes the FXR1 network a hub for transient interactions of proteins lacking RNA-binding domains. We show that the FXR1 network is necessary for RhoA signaling-induced actomyosin reorganization to provide spatial proximity between kinases and their substrates. A point mutation in FXR1, which is found in its FMR1 homolog and causes Fragile X syndrome, disrupts the network. FXR1 network disruption prevents actomyosin remodeling-an essential and ubiquitous process for the regulation of cell shape, migration, and synaptic function. These findings uncover a structural role for cytoplasmic mRNA and show how the FXR1 RNA-binding protein as part of the FXR1 network acts as organizer of signaling reactions.

mRNA interactions with disordered regions control protein activity.

The cytoplasm is compartmentalized into different translation environments. mRNAs use their 3'UTRs to localize to distinct cytoplasmic compartments, including TIS granules (TGs). Many transcription factors, including MYC, are translated in TGs. It was shown that translation of proteins in TGs enables the formation of protein complexes that cannot be established when these proteins are translated in the cytosol, but the mechanism is poorly understood. Here we show that MYC protein complexes that involve binding to the intrinsically disordered region (IDR) of MYC are only formed when MYC is translated in TGs. TG-dependent protein complexes require TG-enriched mRNAs for assembly. These mRNAs bind to a new and widespread RNA-binding domain in neutral or negatively charged IDRs in several transcription factors, including MYC. RNA-IDR interaction changes the conformational ensemble of the IDR, enabling the formation of MYC protein complexes that act in the nucleus and control functions that cannot be accomplished by cytosolically-translated MYC. We propose that certain mRNAs have IDR chaperone activity as they control IDR conformations. In addition to post-translational modifications, we found a novel mode of protein activity regulation. Since RNA-IDR interactions are prevalent, we suggest that mRNA-dependent control of protein functional states is widespread.

Context-specific regulation and function of mRNA alternative polyadenylation.

Alternative cleavage and polyadenylation (APA) is a widespread mechanism to generate mRNA isoforms with alternative 3' untranslated regions (UTRs). The expression of alternative 3' UTR isoforms is highly cell type specific and is further controlled in a gene-specific manner by environmental cues. In this Review, we discuss how the dynamic, fine-grained regulation of APA is accomplished by several mechanisms, including cis-regulatory elements in RNA and DNA and factors that control transcription, pre-mRNA cleavage and post-transcriptional processes. Furthermore, signalling pathways modulate the activity of these factors and integrate APA into gene regulatory programmes. Dysregulation of APA can reprogramme the outcome of signalling pathways and thus can control cellular responses to environmental changes. In addition to the regulation of protein abundance, APA has emerged as a major regulator of mRNA localization and the spatial organization of protein synthesis. This role enables the regulation of protein function through the addition of post-translational modifications or the formation of protein-protein interactions. We further discuss recent transformative advances in single-cell RNA sequencing and CRISPR-Cas technologies, which enable the mapping and functional characterization of alternative 3' UTRs in any biological context. Finally, we discuss new APA-based RNA therapeutics, including compounds that target APA in cancer and therapeutic genome editing of degenerative diseases.

Enhancers regulate 3' end processing activity to control expression of alternative 3'UTR isoforms.

Multi-UTR genes are widely transcribed and express their alternative 3'UTR isoforms in a cell type-specific manner. As transcriptional enhancers regulate mRNA expression, we investigated if they also regulate 3'UTR isoform expression. Endogenous enhancer deletion of the multi-UTR gene PTEN did not impair transcript production but prevented 3'UTR isoform switching which was recapitulated by silencing of an enhancer-bound transcription factor. In reporter assays, enhancers increase transcript production when paired with single-UTR gene promoters. However, when combined with multi-UTR gene promoters, they change 3'UTR isoform expression by increasing 3' end processing activity of polyadenylation sites. Processing activity of polyadenylation sites is affected by transcription factors, including NF-κB and MYC, transcription elongation factors, chromatin remodelers, and histone acetyltransferases. As endogenous cell type-specific enhancers are associated with genes that increase their short 3'UTRs in a cell type-specific manner, our data suggest that transcriptional enhancers integrate cellular signals to regulate cell type-and condition-specific 3'UTR isoform expression.

A working model for condensate RNA-binding proteins as matchmakers for protein complex assembly.

Most cellular processes are carried out by protein complexes, but it is still largely unknown how the subunits of lowly expressed complexes find each other in the crowded cellular environment. Here, we will describe a working model where RNA-binding proteins in cytoplasmic condensates act as matchmakers between their bound proteins (called protein targets) and newly translated proteins of their RNA targets to promote their assembly into complexes. Different RNA-binding proteins act as scaffolds for various cytoplasmic condensates with several of them supporting translation. mRNAs and proteins are recruited into the cytoplasmic condensates through binding to specific domains in the RNA-binding proteins. Scaffold RNA-binding proteins have a high valency. In our model, they use homotypic interactions to assemble condensates and they use heterotypic interactions to recruit protein targets into the condensates. We propose that unoccupied binding sites in the scaffold RNA-binding proteins transiently retain recruited and newly translated proteins in the condensates, thus promoting their assembly into complexes. Taken together, we propose that lowly expressed subunits of protein complexes combine information in their mRNAs and proteins to colocalize in the cytoplasm. The efficiency of protein complex assembly is increased by transient entrapment accomplished by multivalent RNA-binding proteins within cytoplasmic condensates.

Generation of 3'UTR knockout cell lines by CRISPR/Cas9-mediated genome editing.

In addition to the protein code, messenger RNAs (mRNAs) also contain untranslated regions (UTRs). 3'UTRs span the region between the translational stop codon and the poly(A) tail. Sequence elements located in 3'UTRs are essential for pre-mRNA processing. 3'UTRs also contain elements that can regulate protein abundance, localization, and function. At least half of all human genes use alternative cleavage and polyadenylation (APA) to further diversify the regulatory potential of protein functions. Traditional gene editing approaches are designed to disrupt the production of functional proteins. Here, we describe a method that allows investigators to manipulate 3'UTR sequences of endogenous genes for both single- 3'UTR and multi-3'UTR genes. As 3'UTRs can regulate individual functions of proteins, techniques to manipulate 3'UTRs at endogenous gene loci will help to disentangle multi-functionality of proteins. Furthermore, the ability to directly examine the impact of gene regulatory elements in 3'UTRs will provide further insights into their functional significance.

Endogenous p53 expression in human and mouse is not regulated by its 3'UTR.

The TP53 gene encodes the tumor suppressor p53 which is functionally inactivated in many human cancers. Numerous studies suggested that 3'UTR-mediated p53 expression regulation plays a role in tumorigenesis and could be exploited for therapeutic purposes. However, these studies did not investigate post-transcriptional regulation of the native TP53 gene. Here, we used CRISPR/Cas9 to delete the human and mouse TP53/Trp53 3'UTRs while preserving endogenous mRNA processing. This revealed that the endogenous 3'UTR is not involved in regulating p53 mRNA or protein expression neither in steady state nor after genotoxic stress. Using reporter assays, we confirmed the previously observed repressive effects of the isolated 3'UTR. However, addition of the TP53 coding region to the reporter had a dominant negative impact on expression as its repressive effect was stronger and abrogated the contribution of the 3'UTR. Our data highlight the importance of genetic models in the validation of post-transcriptional gene regulatory effects.

In vivo reconstitution finds multivalent RNA-RNA interactions as drivers of mesh-like condensates.

Liquid-like condensates have been thought to be sphere-like. Recently, various condensates with filamentous morphology have been observed in cells. One such condensate is the TIS granule network that shares a large surface area with the rough endoplasmic reticulum and is important for membrane protein trafficking. It has been unclear how condensates with mesh-like shapes but dynamic protein components are formed. In vitro and in vivo reconstitution experiments revealed that the minimal components are a multivalent RNA-binding protein that concentrates RNAs that are able to form extensive intermolecular mRNA-mRNA interactions. mRNAs with large unstructured regions have a high propensity to form a pervasive intermolecular interaction network that acts as condensate skeleton. The underlying RNA matrix prevents full fusion of spherical liquid-like condensates, thus driving the formation of irregularly shaped membraneless organelles. The resulting large surface area may promote interactions at the condensate surface and at the interface with other organelles.

Gain of Additional BIRC3 Protein Functions through 3'-UTR-Mediated Protein Complex Formation.

Alternative 3' untranslated regions (3' UTRs) are widespread, but their functional roles are largely unknown. We investigated the function of the long BIRC3 3' UTR, which is upregulated in leukemia. The 3' UTR does not regulate BIRC3 protein localization or abundance but is required for CXCR4-mediated B cell migration. We established an experimental pipeline to study the mechanism of regulation and used mass spectrometry to identify BIRC3 protein interactors. In addition to 3'-UTR-independent interactors involved in known BIRC3 functions, we detected interactors that bind only to BIRC3 protein encoded from the mRNA with the long 3' UTR. They regulate several functions, including CXCR4 trafficking. We further identified RNA-binding proteins differentially bound to the alternative 3' UTRs and found that cooperative binding of Staufen and HuR mediates 3'-UTR-dependent complex formation. We show that the long 3' UTR is required for the formation of specific protein complexes that enable additional functions of BIRC3 protein beyond its 3'-UTR-independent functions.

What Are 3' UTRs Doing?

3' untranslated regions (3' UTRs) of messenger RNAs (mRNAs) are best known to regulate mRNA-based processes, such as mRNA localization, mRNA stability, and translation. In addition, 3' UTRs can establish 3' UTR-mediated protein-protein interactions (PPIs), and thus can transmit genetic information encoded in 3' UTRs to proteins. This function has been shown to regulate diverse protein features, including protein complex formation or posttranslational modifications, but is also expected to alter protein conformations. Therefore, 3' UTR-mediated information transfer can regulate protein features that are not encoded in the amino acid sequence. This review summarizes both 3' UTR functions-the regulation of mRNA and protein-based processes-and highlights how each 3' UTR function was discovered with a focus on experimental approaches used and the concepts that were learned. This review also discusses novel approaches to study 3' UTR functions in the future by taking advantage of recent advances in technology.

3' UTRs Regulate Protein Functions by Providing a Nurturing Niche during Protein Synthesis.

Messenger RNAs (mRNAs) are the templates for protein synthesis as the coding region is translated into the amino acid sequence. mRNAs also contain 3' untranslated regions (3' UTRs) that harbor additional elements for the regulation of protein function. If the amino acid sequence of a protein is necessary and sufficient for its function, we call it 3' UTR-independent. In contrast, functions that are accomplished by protein complexes whose formation requires the presence of a specific 3' UTR are 3' UTR-dependent protein functions. We showed that 3' UTRs can regulate protein activity without affecting protein abundance, and alternative 3' UTRs can diversify protein functions. We currently think that the regulation of protein function by 3' UTRs is facilitated by the local environment at the site of protein synthesis, which we call the nurturing niche for nascent proteins. This niche is composed of the mRNA and the bound proteins that consist of RNA-binding proteins and recruited proteins. It enables the formation of specific protein complexes, as was shown for TIS granules, a recently discovered cytoplasmic membraneless organelle. This finding suggests that changing the niche for nascent proteins will alter protein activity and function, implying that cytoplasmic membraneless organelles can regulate protein function in a manner that is independent of protein abundance.

A Membraneless Organelle Associated with the Endoplasmic Reticulum Enables 3'UTR-Mediated Protein-Protein Interactions.

Approximately half of human genes generate mRNAs with alternative 3' untranslated regions (3'UTRs). Through 3'UTR-mediated protein-protein interactions, alternative 3'UTRs enable multi-functionality of proteins with identical amino acid sequence. While studying how information on protein features is transferred from 3'UTRs to proteins, we discovered that the broadly expressed RNA-binding protein TIS11B forms a membraneless organelle, called TIS granule, that enriches membrane protein-encoding mRNAs with multiple AU-rich elements. TIS granules form a reticular meshwork intertwined with the endoplasmic reticulum (ER). The association between TIS granules and the ER creates a subcellular compartment-the TIGER domain-with a biophysically and biochemically distinct environment from the cytoplasm. This compartment promotes 3'UTR-mediated interaction of SET with membrane proteins, thus allowing increased surface expression and functional diversity of proteins, including CD47 and PD-L1. The TIGER domain is a subcellular compartment that enables formation of specific and functionally relevant protein-protein interactions that cannot be established outside.

Two high-risk susceptibility loci at 6p25.3 and 14q32.13 for Waldenström macroglobulinemia.

Waldenström macroglobulinemia (WM)/lymphoplasmacytic lymphoma (LPL) is a rare, chronic B-cell lymphoma with high heritability. We conduct a two-stage genome-wide association study of WM/LPL in 530 unrelated cases and 4362 controls of European ancestry and identify two high-risk loci associated with WM/LPL at 6p25.3 (rs116446171, near EXOC2 and IRF4; OR = 21.14, 95% CI: 14.40-31.03, P = 1.36 × 10-54) and 14q32.13 (rs117410836, near TCL1; OR = 4.90, 95% CI: 3.45-6.96, P = 8.75 × 10-19). Both risk alleles are observed at a low frequency among controls (~2-3%) and occur in excess in affected cases within families. In silico data suggest that rs116446171 may have functional importance, and in functional studies, we demonstrate increased reporter transcription and proliferation in cells transduced with the 6p25.3 risk allele. Although further studies are needed to fully elucidate underlying biological mechanisms, together these loci explain 4% of the familial risk and provide insights into genetic susceptibility to this malignancy.

Widespread intronic polyadenylation inactivates tumour suppressor genes in leukaemia.

DNA mutations are known cancer drivers. Here we investigated whether mRNA events that are upregulated in cancer can functionally mimic the outcome of genetic alterations. RNA sequencing or 3'-end sequencing techniques were applied to normal and malignant B cells from 59 patients with chronic lymphocytic leukaemia (CLL)1-3. We discovered widespread upregulation of truncated mRNAs and proteins in primary CLL cells that were not generated by genetic alterations but instead occurred by intronic polyadenylation. Truncated mRNAs caused by intronic polyadenylation were recurrent (n = 330) and predominantly affected genes with tumour-suppressive functions. The truncated proteins generated by intronic polyadenylation often lack the tumour-suppressive functions of the corresponding full-length proteins (such as DICER and FOXN3), and several even acted in an oncogenic manner (such as CARD11, MGA and CHST11). In CLL, the inactivation of tumour-suppressor genes by aberrant mRNA processing is substantially more prevalent than the functional loss of such genes through genetic events. We further identified new candidate tumour-suppressor genes that are inactivated by intronic polyadenylation in leukaemia and by truncating DNA mutations in solid tumours4,5. These genes are understudied in cancer, as their overall mutation rates are lower than those of well-known tumour-suppressor genes. Our findings show the need to go beyond genomic analyses in cancer diagnostics, as mRNA events that are silent at the DNA level are widespread contributors to cancer pathogenesis through the inactivation of tumour-suppressor genes.

Systems-based approaches enable identification of gene targets which improve the flavour profile of low-ethanol wine yeast strains.

Metabolic engineering has been vital to the development of industrial microbes such as the yeast Saccharomyces cerevisiae. However, sequential rounds of modification are often needed to achieve particular industrial design targets. Systems biology approaches can aid in identifying genetic targets for modification through providing an integrated view of cellular physiology. Recently, research into the generation of commercial yeasts that can produce reduced-ethanol wines has resulted in metabolically-engineered strains of S. cerevisiae that are less efficient at producing ethanol from sugar. However, these modifications led to the concomitant production of off-flavour by-products. A combination of transcriptomics, proteomics and metabolomics was therefore used to investigate the physiological changes occurring in an engineered low-ethanol yeast strain during alcoholic fermentation. Integration of 'omics data identified several metabolic reactions, including those related to the pyruvate node and redox homeostasis, as being significantly affected by the low-ethanol engineering methodology, and highlighted acetaldehyde and 2,4,5-trimethyl-1,3-dioxolane as the main off-flavour compounds. Gene remediation strategies were then successfully applied to decrease the formation of these by-products, while maintaining the 'low-alcohol' phenotype. The data generated from this comprehensive systems-based study will inform wine yeast strain development programmes, which, in turn, could potentially play an important role in assisting winemakers in their endeavour to produce low-alcohol wines with desirable flavour profiles.

Combining liquid chromatography and tandem mass spectrometry approaches to the study of monoterpene glycosides (aroma precursors) in wine grape.

Monoterpene-glycosides are important aroma precursors that, undergo hydrolysis, confer intense floral notes to the wines. Therefore, the knowledge of the nature of the sugar residues and the structure of these molecules is of great interest. In present study, liquid chromatography (LC) separation coupled with different mass spectrometry (MS) experiments for the characterization of these compounds were explored. The LC parameters were tuned to optimize the resolution between the analytes present in grape sample extracts. Twenty principal peaks with a relative abundance >1% were selected and divided in 4 classes characterized by different molecular weight. In general, positive ionization of the studied compounds displayed the [M + NH4 ]+ ion as base peak. On the contrary, a distribution between [M + Cl]- and [M + HCOO]- species was observed in negative ion mode. However, a clear differentiation between the studied compounds was only possible by combining both LC and tandem MS (MS/MS). Indeed, by applying a series of energy resolved MS/MS experiments and monitoring both positive and negative ions, a structural characterization of the analytes was achieved. The proposed LC-MS/MS approach provided the profile of monoterpenol-diglycosides and allowed the identification of a series of isobaric terpene-diglycosides in grape. The study of their MS/MS spectra indicated the structure of geranic and/or nerolic acid aglycones. To verify the interest of studied compounds, a preliminary evaluation of the intensity of signals of these glycosides were carried out. The obtained results showed a significant difference between the grape samples collected in two different vintages.