RNA-biology ruling cancer progression? Focus on 3'UTRs and splicing

Exploring the frequency of a TP53 polyadenylation signal variant in tumor DNA from patients diagnosed with lung adenocarcinomas, sarcomas and uterine leiomyomas

The TP53 3'UTR variant rs78378222 A>C has been detected in different tumor types as a somatic alteration that reduces p53 expression through modification of polyadenylation and miRNA regulation. Its prevalence is not yet known in all tumors. Herein, we examine tumor tissue prevalence of rs7837822 in Brazilian cohorts of patients from south and southeast regions diagnosed with lung adenocarcinoma (LUAD, n=586), sarcoma (SARC, n=188) and uterine leiomyoma (ULM, n=41). The minor allele (C) was identified in heterozygosity in 6/586 LUAD tumors (prevalence = 1.02 %) and none of the SARC and ULM samples. Additionally, next generation sequencing analysis revealed that all variant-positive tumors (n=4) with sample availability had additional pathogenic or likely pathogenic somatic variants in the TP53 coding regions. Among them, 3/4 (75 %) had the same pathogenic or likely pathogenic sequence variant (allele frequency <0.05 in tumor DNA) namely c.751A>C (p.Ile251Leu). Our results indicate a low somatic prevalence of rs78378222 in LUAD, ULM and SARC tumors from Brazilian patients, which suggests that no further analysis of this variant in the specific studied regions of Brazil is warranted. However, these findings should not exclude tumor molecular testing of this TP53 3'UTR functional variant for different populations.

Versatile RNA: overlooked gems of the transcriptome

The critical role of RNA, its use and targetability concerning different aspects of human health are gaining more attention because our understanding of the versatility of RNA has dramatically evolved over the last decades. We now appreciate that RNA is far more critical than a messenger molecule and possesses many complicated functions. As a multifunctional molecule with its sequence, flexible structures and enzymatic abilities, RNA is genuinely powerful. Mammalian transcriptomes consist of a dynamically regulated plethora of coding and noncoding RNA types. However, some aspects of RNA metabolism remain to be explored. In this Viewpoint, we focus on the transcriptome's unconventional and possibly overlooked aspects to emphasize the importance of RNA in mammalian systems.

Differential expression of mRNA 3'-end isoforms in cervical and ovarian cancers

Early diagnosis and therapeutic targeting are continuing challenges for gynecological cancers. Here, we focus on cancer transcriptomes and describe the differential expression of 3'UTR isoforms in patients using an algorithm to detect differential poly(A) site usage. We find primarily 3'UTR shortening cases in cervical cancers compared with the normal cervix. We show differential expression of alternate 3'-end isoforms of FOXP1, VPS4B, and OGT in HPV16-positive patients who develop high-grade cervical lesions compared with the infected but non-progressing group. In contrast, in ovarian cancers, 3'UTR lengthening is more evident compared with normal ovary tissue. Nevertheless, highly malignant ovarian tumors have unique 3'UTR shortening events (e.g., CHRAC1, SLC16A1, and TOP2A), some of which correlate with upregulated protein levels in tumors. Overall, our study shows isoform level deregulation in gynecological cancers and highlights the complexity of the transcriptome. This transcript diversity could help identify novel cancer genes and provide new possibilities for diagnosis and therapy.

The 3' Non-Coding Sequence Negatively Regulates PD-L1 Expression, and Its Regulators Are Systematically Identified in Pan-Cancer

The 3'-untranslated region (3'-UTR) of PD-L1 is significantly longer than the coding sequences (CDSs). However, its role and regulators have been little studied. We deleted whole 3'-UTR region by CRISPR-Cas9. Prognostic analysis was performed using online tools. Immune infiltration analysis was performed using the Timer and Xcell packages. Immunotherapy response prediction and Cox regression was performed using the R software. MicroRNA network analysis was conducted by the Cytoscape software. The level of PD-L1 was significantly and dramatically up-regulated in cells after deleting the 3'-UTR. Additionally, we discovered a panel of 43 RNA-binding proteins (RBPs) whose expression correlates with PD-L1 in the majority of cancer cell lines and tumor tissues. Among these RBPs, PARP14 is widely associated with immune checkpoints, the tumor microenvironment, and immune-infiltrating cells in various cancer types. We also identified 38 microRNAs whose individual expressions are associated with PD-L1 across different cancers. Notably, miR-3139, miR-4761, and miR-15a-5p showed significant associations with PD-L1 in most cancer types. Furthermore, we revealed 21 m6A regulators that strongly correlate with PD-L1. Importantly, by combining the identified RBP and m6A regulators, we established an immune signature consisting of RBMS1, QKI, ZC3HAV1, and RBM38. This signature can be used to predict the responsiveness of cancer patients to immune checkpoint blockade treatment. We demonstrated the critical role of the 3'-UTR in the regulation of PD-L1 and identified a significant number of potential PD-L1 regulators across various types of cancer. The biomarker signature generated from our findings shows promise in predicting patient prognosis. However, further biological investigation is necessary to explore the potential of these PD-L1 regulators.

Differential impact of BRAFV600E isoforms on tumorigenesis in a zebrafish model of melanoma

BRAFV600E comes as two main splicing variants. The well-studied ref isoform and the recently discovered X1 isoform are co-expressed in cancer cells and differ in terms of 3'UTR length and sequence, as well as C-term protein sequence. Here, we use a melanoma model in zebrafish to study the role played by each isoform in larval pigmentation, nevi formation, and their progression into melanoma tumours. We show that both BRAFV600E-ref and BRAFV600E-X1 proteins promote larval pigmentation and nevi formation, while melanoma-free survival curves performed in adult fish indicate that BRAFV600E-ref protein is a much stronger melanoma driver that BRAFV600E-X1 protein. Crucially, we also show that the presence of the 3'UTR suppresses the effect of ref protein. Our data highlight the necessity to undertake a systematic study of BRAFV600E isoforms, in order to uncover the full spectrum of their kinase-(in)dependent and coding-(in)dependent functions, hence to develop more informed strategies for therapeutic targeting.

PARP1 negatively regulates MAPK signaling by impairing BRAF-X1 translation

In human cells BRAF oncogene is invariably expressed as a mix of two coding transcripts: BRAF-ref and BRAF-X1. These two mRNA isoforms, remarkably different in the sequence and length of their 3'UTRs, are potentially involved in distinct post-transcriptional regulatory circuits. Herein, we identify PARP1 among the mRNA Binding Proteins that specifically target the X1 3'UTR in melanoma cells. Mechanistically, PARP1 Zinc Finger domain down-regulates BRAF expression at the translational level. As a consequence, it exerts a negative impact on MAPK pathway, and sensitizes melanoma cells to BRAF and MEK inhibitors, both in vitro and in vivo. In summary, our study unveils PARP1 as a negative regulator of the highly oncogenic MAPK pathway in melanoma, through the modulation of BRAF-X1 expression.

Loss-of-function mutation in Pcsk1 increases serum APOA1 level and LCAT activity in mice

Background

The convertase subtilisin/kexin family 1 gene (PCSK1) has been associated in various human genetics studies with a wide spectrum of metabolic phenotypes, including early-onset obesity, hyperphagia, diabetes insipidus, and others. Despite the evident influence of PCSK1 on obesity and the known functions of other PCSKs in lipid metabolism, the role of PCSK1 specifically in lipid and cholesterol metabolism remains unclear. This study evaluated the effect of loss of PCSK1 function on high-density lipoprotein (HDL) metabolism in mice.

Results

HDL cholesterol, apolipoprotein A1 (APOA1) levels in serum and liver, and the activities of two enzymes (lecithin-cholesterol acyltransferase, LCAT and phospholipid transfer protein, PLTP) were evaluated in 8-week-old mice with a non-synonymous single nucleotide mutation leading to an amino acid substitution in PCSK1, which results in a loss of protein’s function. Mutant mice had similar serum HDL cholesterol concentration but increased levels of serum total and mature APOA1, and LCAT activity in comparison to controls.

Conclusions

This study presents the first evaluation of the role of PCSK1 in HDL metabolism using a loss-of-function mutant mouse model. Further investigations will be needed to determine the underlying molecular mechanism.

Identification of an mRNA isoform switch for HNRNPA1 in breast cancers

Roles of HNRNPA1 are beginning to emerge in cancers; however, mechanisms causing deregulation of HNRNPA1 function remain elusive. Here, we describe an isoform switch between the 3'-UTR isoforms of HNRNPA1 in breast cancers. We show that the dominantly expressed isoform in mammary tissue has a short half-life. In breast cancers, this isoform is downregulated in favor of a stable isoform. The stable isoform is expressed more in breast cancers, and more HNRNPA1 protein is synthesized from this isoform. High HNRNPA1 protein levels correlate with poor survival in patients. In support of this, silencing of HNRNPA1 causes a reversal in neoplastic phenotypes, including proliferation, clonogenic potential, migration, and invasion. In addition, silencing of HNRNPA1 results in the downregulation of microRNAs that map to intragenic regions. Among these miRNAs, miR-21 is known for its transcriptional upregulation in breast and numerous other cancers. Altogether, the cancer-specific isoform switch we describe here for HNRNPA1 emphasizes the need to study gene expression at the isoform level in cancers to identify novel cases of oncogene activation.

Double-Stranded RNA Structural Elements Holding the Key to Translational Regulation in Cancer: The Case of Editing in RNA-Binding Motif Protein 8A

Mesothelioma is an aggressive cancer associated with asbestos exposure. RNA-binding motif protein 8a (RBM8A) mRNA editing increases in mouse tissues upon asbestos exposure. The aim of this study was to further characterize the role of RBM8A in mesothelioma and the consequences of its mRNA editing. RBM8A protein expression was higher in mesothelioma compared to mesothelial cells. Silencing RBM8A changed splicing patterns in mesothelial and mesothelioma cells but drastically reduced viability only in mesothelioma cells. In the tissues of asbestos-exposed mice, editing of Rbm8a mRNA was associated with increased protein immunoreactivity, with no change in mRNA levels. Increased adenosine deaminase acting on dsRNA (ADAR)-dependent editing of Alu elements in the RBM8A 3′UTR was observed in mesothelioma cells compared to mesothelial cells. Editing stabilized protein expression. The unedited RBM8A 3′UTR had a stronger interaction with Musashi (MSI) compared to the edited form. The silencing of MSI2 in mesothelioma or overexpression of Adar2 in mesothelial cells resulted in increased RBM8A protein levels. Therefore, ADAR-dependent editing contributes to maintaining elevated RBM8A protein levels in mesothelioma by counteracting MSI2-driven downregulation. A wider implication of this mechanism for the translational control of protein expression is suggested by the editing of similarly structured Alu elements in several other transcripts.

A C-term truncated EIF2Bγ protein encoded by an intronically polyadenylated isoform introduces unfavorable EIF2Bγ-EIF2γ interactions

Eukaryotic translation initiates upon recruitment of the EIF2-GTP·Met-tRNA_i ternary complex (TC) to the ribosomes. EIF2 (α, β, γ subunits) is a GTPase. The GDP to GTP exchange within EIF2 is facilitated by the guanine nucleotide exchange factor EIF2B (α-ε subunits). During stress-induced conditions, phosphorylation of the α-subunit of EIF2 turns EIF2 into an inhibitor of EIF2B. In turn, inhibition of EIF2B decreases TC formation and triggers the internal stress response (ISR), which determines the cell fate. Deregulated ISR has been linked to neurodegenerative disorders and cancer, positioning EIF2B as a promising therapeutic target. Hence, a better understanding of the mechanisms/factors that regulate EIF2B activity is required. Here, combining transcript and protein level analyses, we describe an intronically polyadenylated (IPA) transcript of EIF2B's γ-subunit. We show that the IPA mRNA isoform is translated into a C-terminus truncated protein. Using structural modeling, we predict that the truncated EIF2Bγ protein has unfavorable interactions with EIF2γ, leading to a potential decrease in the stability of the nonproductive EIF2:EIF2B complex. While we discovered and confirmed the IPA mRNA isoform in breast cancer cells, the expression of this isoform is not cancer-specific and is widely present in normal tissues. Overall, our data show that a truncated EIF2Bγ protein co-exists with the canonical protein and is an additional player to regulate the equilibrium between productive and nonproductive states of the EIF2:EIF2B complex. These results may have implications in stress-induced translation control in normal and disease states. Our combinatorial approach demonstrates the need to study noncanonical mRNA and protein isoforms to understand protein interactions and intricate molecular mechanisms.

Non-Coding Variants in Cancer: Mechanistic Insights and Clinical Potential for Personalized Medicine

The cancer genome is characterized by extensive variability, in the form of Single Nucleotide Polymorphisms (SNPs) or structural variations such as Copy Number Alterations (CNAs) across wider genomic areas. At the molecular level, most SNPs and/or CNAs reside in non-coding sequences, ultimately affecting the regulation of oncogenes and/or tumor-suppressors in a cancer-specific manner. Notably, inherited non-coding variants can predispose for cancer decades prior to disease onset. Furthermore, accumulation of additional non-coding driver mutations during progression of the disease, gives rise to genomic instability, acting as the driving force of neoplastic development and malignant evolution. Therefore, detection and characterization of such mutations can improve risk assessment for healthy carriers and expand the diagnostic and therapeutic toolbox for the patient. This review focuses on functional variants that reside in transcribed or not transcribed non-coding regions of the cancer genome and presents a collection of appropriate state-of-the-art methodologies to study them.