Monitoring the 5'UTR landscape reveals isoform switches to drive translational efficiencies in cancer

GWAS of multiple neuropathology endophenotypes identifies new risk loci and provides insights into the genetic risk of dementia

Genome-wide association studies (GWAS) have identified >80 Alzheimer's disease and related dementias (ADRD)-associated genetic loci. However, the clinical outcomes used in most previous studies belie the complex nature of underlying neuropathologies. Here we performed GWAS on 11 ADRD-related neuropathology endophenotypes with participants drawn from the following three sources: the National Alzheimer's Coordinating Center, the Religious Orders Study and Rush Memory and Aging Project, and the Adult Changes in Thought study (n = 7,804 total autopsied participants). We identified eight independent significantly associated loci, of which four were new (COL4A1, PIK3R5, LZTS1 and APOC2). Separately testing known ADRD loci, 19 loci were significantly associated with at least one neuropathology after false-discovery rate adjustment. Genetic colocalization analyses identified pleiotropic effects and quantitative trait loci. Methylation in the cerebral cortex at two sites near APOC2 was associated with cerebral amyloid angiopathy. Studies that include neuropathology endophenotypes are an important step in understanding the mechanisms underlying genetic ADRD risk.

The regulatory landscape of interacting RNA and protein pools in cellular homeostasis and cancer

Biological systems encompass intricate networks governed by RNA-protein interactions that play pivotal roles in cellular functions. RNA and proteins constituting 1.1% and 18% of the mammalian cell weight, respectively, orchestrate vital processes from genome organization to translation. To date, disentangling the functional fraction of the human genome has presented a major challenge, particularly for noncoding regions, yet recent discoveries have started to unveil a host of regulatory functions for noncoding RNAs (ncRNAs). While ncRNAs exist at different sizes, structures, degrees of evolutionary conservation and abundances within the cell, they partake in diverse roles either alone or in combination. However, certain ncRNA subtypes, including those that have been described or remain to be discovered, are poorly characterized given their heterogeneous nature. RNA activity is in most cases coordinated through interactions with RNA-binding proteins (RBPs). Extensive efforts are being made to accurately reconstruct RNA-RBP regulatory networks, which have provided unprecedented insight into cellular physiology and human disease. In this review, we provide a comprehensive view of RNAs and RBPs, focusing on how their interactions generate functional signals in living cells, particularly in the context of post-transcriptional regulatory processes and cancer.

Isoform-Level Transcriptome Analysis of Peripheral Blood Mononuclear Cells from Breast Cancer Patients Identifies a Disease-Associated RASGEF1A Isoform

Background: Breast cancer (BC) comprises multiple subtypes with distinct molecular features, which differ in their interplay with host immunity, prognosis, and treatment. Non-invasive blood analyses can provide valuable insights into systemic immunity during cancer. The aim of this study was to analyze the expression of transcriptional isoforms in peripheral blood mononuclear cells (PBMCs) from BC patients and healthy women to identify potential BC immune biomarkers. Methods: RNA sequencing and isoform-level bioinformatics were performed on PBMCs from 12 triple-negative and 13 luminal A patients. Isoform expression validation by qRT-PCR and clinicopathological correlations were performed in a larger cohort (156 BC patients and 32 healthy women). Results: Transcriptional analyses showed a significant (p < 0.001) decrease in the ENST00000374459 RASGEF1A isoform in PBMCs of BC compared to healthy subjects, indicating disease-related expression changes. The decrease was associated with higher ctDNA and Ki-67 values. Conclusions: The levels of the RASGEF1A transcriptional isoform ENST00000374459 may have the potential to distinguish between BC and healthy subjects. The downregulation of ENST00000374459 in breast cancer is associated with higher proliferation and ctDNA shedding. Specialized bioinformatics analyses such as isoform analyses hold significant promise in the detection of biomarkers, since standard RNA sequencing analyses may overlook specific transcriptional changes that may be disease-associated and biologically important.

Novel Pannexin 1 isoform is increased in cancer

Pannexin 1 (PANX1) is upregulated in many cancers, where its activity and signalling promote tumorigenic properties. Here, we report a novel ∼25 kDa isoform of human PANX1 (hPANX1-25K) which lacks the N-terminus and was detected in several human cancer cell lines including melanoma, osteosarcoma, breast cancer and glioblastoma multiforme. This isoform was increased upon hPANX1 CRISPR/Cas9 deletion targeting the first exon near M1, suggesting a potential alternative translation initiation (ATI) site. hPANX1-25K was confirmed to be a hPANX1 isoform via mass spectrometry, can be N-linked glycosylated at N254, and can interact with both β-catenin and full length hPANX1. A double deletion of hPANX1 and hPANX1-25K reduces cell growth and viability in cancer cells. hPANX1-25K is prevalent throughout melanoma progression, and its levels are increased in squamous cell carcinoma cells and patient-derived tumours, compared to keratinocytes and normal skin, indicating that it may be differentially regulated in normal and cancer cells.

Alternative splicing of ALDOA confers tamoxifen resistance in breast cancer

The cancer-associated alternative splicing (AS) events generate cancer-related transcripts which are involved in uncontrolled cell proliferation and drug resistance. However, the key AS variants implicated in tamoxifen (TAM) resistance in breast cancer remain elusive. In the current study, we investigated the landscape of AS events in nine pairs of primary and relapse breast tumors from patients receiving TAM-based therapy. We unrevealed a notable association between the inclusion of exon 7.2 in the 5'untranslated region (5'UTR) of ALDOA mRNA and TAM resistance. Mechanistically, the inclusion of ALDOA exon 7.2 enhances the translation efficiency of the transcript, resulting in increased ALDOA protein expression, mTOR pathway activity, and the promotion of TAM resistance in breast cancer cells. Moreover, the inclusion of exon 7.2 in ALDOA mRNA is mediated by MSI1 via direct interaction. In addition, elevated inclusion of ALDOA exon 7.2 or expression of MSI1 is associated with an unfavorable prognosis in patients undergoing endocrine therapy. Notably, treatment with Aldometanib, an ALDOA inhibitor, effectively restrains the growth of TAM-resistant breast cancer cells in vitro and in vivo. The present study unveils the pivotal role of an AS event in ALDOA, under the regulation of MSI1, in driving TAM resistance in breast cancer. Therefore, this study provides a promising therapeutic avenue targeting ALDOA to combat TAM resistance.

Current limitations in predicting mRNA translation with deep learning models

BACKGROUND

The design of nucleotide sequences with defined properties is a long-standing problem in bioengineering. An important application is protein expression, be it in the context of research or the production of mRNA vaccines. The rate of protein synthesis depends on the 5' untranslated region (5'UTR) of the mRNAs, and recently, deep learning models were proposed to predict the translation output of mRNAs from the 5'UTR sequence. At the same time, large data sets of endogenous and reporter mRNA translation have become available.

RESULTS

In this study, we use complementary data obtained in two different cell types to assess the accuracy and generality of currently available models for predicting translational output. We find that while performing well on the data sets on which they were trained, deep learning models do not generalize well to other data sets, in particular of endogenous mRNAs, which differ in many properties from reporter constructs.

CONCLUSIONS

These differences limit the ability of deep learning models to uncover mechanisms of translation control and to predict the impact of genetic variation. We suggest directions that combine high-throughput measurements and machine learning to unravel mechanisms of translation control and improve construct design.

In vivo single-cell CRISPR uncovers distinct TNF programmes in tumour evolution

The tumour evolution model posits that malignant transformation is preceded by randomly distributed driver mutations in cancer genes, which cause clonal expansions in phenotypically normal tissues. Although clonal expansions can remodel entire tissues1-3, the mechanisms that result in only a small number of clones transforming into malignant tumours remain unknown. Here we develop an in vivo single-cell CRISPR strategy to systematically investigate tissue-wide clonal dynamics of the 150 most frequently mutated squamous cell carcinoma genes. We couple ultrasound-guided in utero lentiviral microinjections, single-cell RNA sequencing and guide capture to longitudinally monitor clonal expansions and document their underlying gene programmes at single-cell transcriptomic resolution. We uncover a tumour necrosis factor (TNF) signalling module, which is dependent on TNF receptor 1 and involving macrophages, that acts as a generalizable driver of clonal expansions in epithelial tissues. Conversely, during tumorigenesis, the TNF signalling module is downregulated. Instead, we identify a subpopulation of invasive cancer cells that switch to an autocrine TNF gene programme associated with epithelial-mesenchymal transition. Finally, we provide in vivo evidence that the autocrine TNF gene programme is sufficient to mediate invasive properties and show that the TNF signature correlates with shorter overall survival of patients with squamous cell carcinoma. Collectively, our study demonstrates the power of applying in vivo single-cell CRISPR screening to mammalian tissues, unveils distinct TNF programmes in tumour evolution and highlights the importance of understanding the relationship between clonal expansions in epithelia and tumorigenesis.

Alternative translation initiation by ribosomal leaky scanning produces multiple isoforms of the Pif1 helicase

In budding yeast, the integrity of both the nuclear and mitochondrial genomes relies on dual-targeted isoforms of the conserved Pif1 helicase, generated by alternative translation initiation (ATI) of PIF1 mRNA from two consecutive AUG codons flanking a mitochondrial targeting signal. Here, we demonstrate that ribosomal leaky scanning is the specific ATI mechanism that produces not only these, but also novel, previously uncharacterized Pif1 isoforms. Both in-frame, downstream AUGs as well as near-cognate start codons contribute to the generation of these alternative isoforms. This has crucial implications for the rational design of genuine separation-of-function alleles and provides an explanation for the suboptimal behaviour of the widely employed mitochondrial- (pif1-m1) and nuclear-deficient (pif1-m2) alleles, with mutations in the first or second AUG codon, respectively. We have taken advantage of this refined model to develop improved versions of these alleles, which will serve as valuable tools to elucidate novel functions of this helicase and to disambiguate previously described genetic interactions of PIF1 in the context of nuclear and mitochondrial genome stability.

Upstream open reading frames: new players in the landscape of cancer gene regulation

The translation of RNA by ribosomes represents a central biological process and one of the most dysregulated processes in cancer. While translation is traditionally thought to occur exclusively in the protein-coding regions of messenger RNAs (mRNAs), recent transcriptome-wide approaches have shown abundant ribosome activity across diverse stretches of RNA transcripts. The most common type of this kind of ribosome activity occurs in gene leader sequences, also known as 5' untranslated regions (UTRs) of the mRNA, that precede the main coding sequence. Translation of these upstream open reading frames (uORFs) is now known to occur in upwards of 25% of all protein-coding genes. With diverse functions from RNA regulation to microprotein generation, uORFs are rapidly igniting a new arena of cancer biology, where they are linked to cancer genetics, cancer signaling, and tumor-immune interactions. This review focuses on the contributions of uORFs and their associated 5'UTR sequences to cancer biology.

Contributions of transcriptional noise to leukaemia evolution: KAT2A as a case-study

Transcriptional noise is proposed to participate in cell fate changes, but contributions to mammalian cell differentiation systems, including cancer, remain associative. Cancer evolution is driven by genetic variability, with modulatory or contributory participation of epigenetic variants. Accumulation of epigenetic variants enhances transcriptional noise, which can facilitate cancer cell fate transitions. Acute myeloid leukaemia (AML) is an aggressive cancer with strong epigenetic dependencies, characterized by blocked differentiation. It constitutes an attractive model to probe links between transcriptional noise and malignant cell fate regulation. Gcn5/KAT2A is a classical epigenetic transcriptional noise regulator. Its loss increases transcriptional noise and modifies cell fates in stem and AML cells. By reviewing the analysis of KAT2A-depleted pre-leukaemia and leukaemia models, I discuss that the net result of transcriptional noise is diversification of cell fates secondary to alternative transcriptional programmes. Cellular diversification can enable or hinder AML progression, respectively, by differentiation of cell types responsive to mutations, or by maladaptation of leukaemia stem cells. KAT2A-dependent noise-responsive genes participate in ribosome biogenesis and KAT2A loss destabilizes translational activity. I discuss putative contributions of perturbed translation to AML biology, and propose KAT2A loss as a model for mechanistic integration of transcriptional and translational control of noise and fate decisions. This article is part of a discussion meeting issue 'Causes and consequences of stochastic processes in development and disease'.

(Alternative) transcription start sites as regulators of RNA processing

Alternative transcription start site usage (ATSS) is a widespread regulatory strategy that enables genes to choose between multiple genomic loci for initiating transcription. This mechanism is tightly controlled during development and is often altered in disease states. In this review, we examine the growing evidence highlighting a role for transcription start sites (TSSs) in the regulation of mRNA isoform selection during and after transcription. We discuss how the choice of transcription initiation sites influences RNA processing and the importance of this crosstalk for cell identity and organism function. We also speculate on possible mechanisms underlying the integration of transcriptional and post-transcriptional processes.

Quantitative profiling of human translation initiation reveals regulatory elements that potently affect endogenous and therapeutically modified mRNAs

mRNA therapeutics offer a potentially universal strategy for the efficient development and delivery of therapeutic proteins. Current mRNA vaccines include chemically modified nucleotides to reduce cellular immunogenicity. Here, we develop an efficient, high-throughput method to measure human translation initiation on therapeutically modified as well as endogenous RNAs. Using systems-level biochemistry, we quantify ribosome recruitment to tens of thousands of human 5' untranslated regions and identify sequences that mediate 250-fold effects. We observe widespread effects of coding sequences on translation initiation and identify small regulatory elements of 3-6 nucleotides that are sufficient to potently affect translational output. Incorporation of N1-methylpseudouridine (m1Ψ) selectively enhances translation by specific 5' UTRs that we demonstrate surpass those of current mRNA vaccines. Our approach is broadly applicable to dissect mechanisms of human translation initiation and engineer more potent therapeutic mRNAs.

Highlights

Measurement of >30,000 human 5' UTRs reveals a 250-fold range of translation outputSystematic mutagenesis demonstrates the causality of short (3-6nt) regulatory elementsN1-methylpseudouridine alters translation initiation in a sequence-specific mannerOptimal modified 5' UTRs outperform those in the current class of mRNA vaccines.

VAV2 orchestrates the interplay between regenerative proliferation and ribogenesis in both keratinocytes and oral squamous cell carcinoma

VAV2 is an activator of RHO GTPases that promotes and maintains regenerative proliferation-like states in normal keratinocytes and oral squamous cell carcinoma (OSCC) cells. Here, we demonstrate that VAV2 also regulates ribosome biogenesis in those cells, a program associated with poor prognosis of human papilloma virus-negative (HPV-) OSCC patients. Mechanistically, VAV2 regulates this process in a catalysis-dependent manner using a conserved pathway comprising the RAC1 and RHOA GTPases, the PAK and ROCK family kinases, and the c-MYC and YAP/TAZ transcription factors. This pathway directly promotes RNA polymerase I activity and synthesis of 47S pre-rRNA precursors. This process is further consolidated by the upregulation of ribosome biogenesis factors and the acquisition of the YAP/TAZ-dependent undifferentiated cell state. Finally, we show that RNA polymerase I is a therapeutic Achilles' heel for both keratinocytes and OSCC patient-derived cells endowed with high VAV2 catalytic activity. Collectively, these findings highlight the therapeutic potential of modulating VAV2 and the ribosome biogenesis pathways in both preneoplastic and late progression stages of OSCC.

No country for old methods: New tools for studying microproteins

Microproteins encoded by small open reading frames (sORFs) have emerged as a fascinating frontier in genomics. Traditionally overlooked due to their small size, recent technological advancements such as ribosome profiling, mass spectrometry-based strategies and advanced computational approaches have led to the annotation of more than 7000 sORFs in the human genome. Despite the vast progress, only a tiny portion of these microproteins have been characterized and an important challenge in the field lies in identifying functionally relevant microproteins and understanding their role in different cellular contexts. In this review, we explore the recent advancements in sORF research, focusing on the new methodologies and computational approaches that have facilitated their identification and functional characterization. Leveraging these new tools hold great promise for dissecting the diverse cellular roles of microproteins and will ultimately pave the way for understanding their role in the pathogenesis of diseases and identifying new therapeutic targets.

Combining a prioritization strategy and functional studies nominates 5'UTR variants underlying inherited retinal disease

BACKGROUND

5' untranslated regions (5'UTRs) are essential modulators of protein translation. Predicting the impact of 5'UTR variants is challenging and rarely performed in routine diagnostics. Here, we present a combined approach of a comprehensive prioritization strategy and functional assays to evaluate 5'UTR variation in two large cohorts of patients with inherited retinal diseases (IRDs).

METHODS

We performed an isoform-level re-analysis of retinal RNA-seq data to identify the protein-coding transcripts of 378 IRD genes with highest expression in retina. We evaluated the coverage of their 5'UTRs by different whole exome sequencing (WES) kits. The selected 5'UTRs were analyzed in whole genome sequencing (WGS) and WES data from IRD sub-cohorts from the 100,000 Genomes Project (n = 2397 WGS) and an in-house database (n = 1682 WES), respectively. Identified variants were annotated for 5'UTR-relevant features and classified into seven categories based on their predicted functional consequence. We developed a variant prioritization strategy by integrating population frequency, specific criteria for each category, and family and phenotypic data. A selection of candidate variants underwent functional validation using diverse approaches.

RESULTS

Isoform-level re-quantification of retinal gene expression revealed 76 IRD genes with a non-canonical retina-enriched isoform, of which 20 display a fully distinct 5'UTR compared to that of their canonical isoform. Depending on the probe design, 3-20% of IRD genes have 5'UTRs fully captured by WES. After analyzing these regions in both cohorts, we prioritized 11 (likely) pathogenic variants in 10 genes (ARL3, MERTK, NDP, NMNAT1, NPHP4, PAX6, PRPF31, PRPF4, RDH12, RD3), of which 7 were novel. Functional analyses further supported the pathogenicity of three variants. Mis-splicing was demonstrated for the PRPF31:c.-9+1G>T variant. The MERTK:c.-125G>A variant, overlapping a transcriptional start site, was shown to significantly reduce both luciferase mRNA levels and activity. The RDH12:c.-123C>T variant was found in cis with the hypomorphic RDH12:c.701G>A (p.Arg234His) variant in 11 patients. This 5'UTR variant, predicted to introduce an upstream open reading frame, was shown to result in reduced RDH12 protein but unaltered mRNA levels.

CONCLUSIONS

This study demonstrates the importance of 5'UTR variants implicated in IRDs and provides a systematic approach for 5'UTR annotation and validation that is applicable to other inherited diseases.

Intra-promoter switch of transcription initiation sites in proliferation signaling-dependent RNA metabolism

Global changes in transcriptional regulation and RNA metabolism are crucial features of cancer development. However, little is known about the role of the core promoter in defining transcript identity and post-transcriptional fates, a potentially crucial layer of transcriptional regulation in cancer. In this study, we use CAGE-seq analysis to uncover widespread use of dual-initiation promoters in which non-canonical, first-base-cytosine (C) transcription initiation occurs alongside first-base-purine initiation across 59 human cancers and healthy tissues. C-initiation is often followed by a 5' terminal oligopyrimidine (5'TOP) sequence, dramatically increasing the range of genes potentially subjected to 5'TOP-associated post-transcriptional regulation. We show selective, dynamic switching between purine and C-initiation site usage, indicating transcription initiation-level regulation in cancers. We additionally detail global metabolic changes in C-initiation transcripts that mark differentiation status, proliferative capacity, radiosensitivity, and response to irradiation and to PI3K-Akt-mTOR and DNA damage pathway-targeted radiosensitization therapies in colorectal cancer organoids and cancer cell lines and tissues.

The regulatory landscape of 5' UTRs in translational control during zebrafish embryogenesis

The 5' UTRs of mRNAs are critical for translation regulation, but their in vivo regulatory features are poorly characterized. Here, we report the regulatory landscape of 5' UTRs during early zebrafish embryogenesis using a massively parallel reporter assay of 18,154 sequences coupled to polysome profiling. We found that the 5' UTR is sufficient to confer temporal dynamics to translation initiation, and identified 86 motifs enriched in 5' UTRs with distinct ribosome recruitment capabilities. A quantitative deep learning model, DaniO5P, revealed a combined role for 5' UTR length, translation initiation site context, upstream AUGs and sequence motifs on in vivo ribosome recruitment. DaniO5P predicts the activities of 5' UTR isoforms and indicates that modulating 5' UTR length and motif grammar contributes to translation initiation dynamics. This study provides a first quantitative model of 5' UTR-based translation regulation in early vertebrate development and lays the foundation for identifying the underlying molecular effectors.

The Effects of Deregulated Ribosomal Biogenesis in Cancer

Ribosomes are macromolecular ribonucleoprotein complexes assembled from RNA and proteins. Functional ribosomes arise from the nucleolus, require ribosomal RNA processing and the coordinated assembly of ribosomal proteins (RPs), and are frequently hyperactivated to support the requirement for protein synthesis during the self-biosynthetic and metabolic activities of cancer cells. Studies have provided relevant information on targeted anticancer molecules involved in ribosome biogenesis (RiBi), as increased RiBi is characteristic of many types of cancer. The association between unlimited cell proliferation and alterations in specific steps of RiBi has been highlighted as a possible critical driver of tumorigenesis and metastasis. Thus, alterations in numerous regulators and actors involved in RiBi, particularly in cancer, significantly affect the rate and quality of protein synthesis and, ultimately, the transcriptome to generate the associated proteome. Alterations in RiBi in cancer cells activate nucleolar stress response-related pathways that play important roles in cancer-targeted interventions and immunotherapies. In this review, we focus on the association between alterations in RiBi and cancer. Emphasis is placed on RiBi deregulation and its secondary consequences, including changes in protein synthesis, loss of RPs, adaptive transcription and translation, nucleolar stress regulation, metabolic changes, and the impaired ribosome biogenesis checkpoint.

Translational Control of Metabolism and Cell Cycle Progression in Hepatocellular Carcinoma

The liver is a metabolic hub characterized by high levels of protein synthesis. Eukaryotic initiation factors, eIFs, control the first phase of translation, initiation. Initiation factors are essential for tumor progression and, since they regulate the translation of specific mRNAs downstream of oncogenic signaling cascades, may be druggable. In this review, we address the issue of whether the massive translational machinery of liver cells contributes to liver pathology and to the progression of hepatocellular carcinoma (HCC); it represents a valuable biomarker and druggable target. First, we observe that the common markers of HCC cells, such as phosphorylated ribosomal protein S6, belong to the ribosomal and translational apparatus. This fact is in agreement with observations that demonstrate a huge amplification of the ribosomal machinery during the progression to HCC. Some translation factors, such as eIF4E and eIF6, are then harnessed by oncogenic signaling. In particular, the action of eIF4E and eIF6 is particularly important in HCC when driven by fatty liver pathologies. Indeed, both eIF4E and eIF6 amplify at the translational level the production and accumulation of fatty acids. As it is evident that abnormal levels of these factors drive cancer, we discuss their therapeutic value.