Pan-cancer analysis of mRNA stability for decoding tumour post-transcriptional programs

Differential gene expression in uterine endometrioid cancer cells and adjusted normal tissue

Endometrial cancer is a significant public health concern with rising incidence rates globally. Understanding the molecular mechanisms underlying this disease is crucial for developing effective therapeutic strategies. Our study aimed to characterize transcriptional changes in endometrial cancer tissues compared to adjusted healthy tissue. Using RNA sequencing, we identified 2483 differentially expressed genes (DEGs), including protein-coding genes, long non-coding RNAs (lncRNAs), and microRNAs (miRNAs). Notably, several known cancer-related genes were differentially expressed, such as MYC, AKT3, CCND1, and CDKN2A. Pathway analysis revealed significant alterations in cell cycle regulation, several signaling pathways, and metabolic processes. These findings provide valuable insights into the molecular pathways dysregulated in endometrial cancer. Our results may contribute to the development of novel therapeutic targets and biomarkers for this disease.

β-Ionone suppresses colorectal tumorigenesis by activating OR51E2, a potential tumor suppressor

BACKGROUND

Olfactory receptors (ORs) are present in non-olfactory tissues and contribute to diverse biological roles beyond smell perception. Among them, OR51E2 has been associated with cancer biology, and its activator, β-ionone, a natural terpenoid, is known to have anticancer effects.

PURPOSE

This study aimed to clarify the tumor-suppressive role of OR51E2 in colorectal cancer (CRC), unravel the regulatory mechanism underlying its downregulation, and evaluate the therapeutic potential of β-ionone, an OR51E2 ligand, in CRC progression.

STUDY DESIGN AND METHODS

OR51E2 expression was analyzed in human CRC tissues, matched adjacent normal tissues, and cell lines. The involvement of N6-methyladenosine (m6A) modification of OR51E2 mRNA stability was examined using METTL3/14 and YTHDF1/2/3 knockdown experiments. β-Ionone-mediated effects on intracellular calcium signaling, cell proliferation, migration, and apoptosis were evaluated in an OR51E2-dependent manner. The therapeutic efficacy of β-ionone was further evaluated in vivo using a xenograft model in nude mice.

RESULTS

OR51E2 mRNA expression and immunoreactivity were significantly reduced in CRC cells and tissues due to decreased mRNA stability. Knockdown of METTL3/14 or YTHDF1/2/3 increased OR51E2 mRNA and protein expression and inhibited CRC cell proliferation. Treatment with STM2457, an METTL3 inhibitor, restored OR51E2 expression and suppressed CRC cell proliferation. β-Ionone, a ligand of OR51E2, increased intracellular calcium levels, decreased MEK/ERK phosphorylation, and inhibited CRC cell proliferation while inducing apoptosis. These effects were abolished in OR51E2 knockdown cells. In a xenograft model, β-ionone administration (5 and 10 mg/kg body weight) significantly reduced tumor growth.

CONCLUSION

This study identifies m6A modification as a critical mechanism underlying the downregulation of OR51E2 in CRC. Activation of OR51E2 by β-ionone suppresses CRC cell proliferation and induces apoptosis by elevating intracellular calcium levels, which inhibits the MEK-ERK pathway. These findings highlight OR51E2 as a potential therapeutic target and suggest that β-ionone or m6A inhibition may represent novel strategies for CRC treatment.

Resistance to Radiation Enhances Metastasis by Altering RNA Metabolism

The cellular programs that mediate therapy resistance are often important drivers of metastasis, a phenomenon that needs to be understood better to improve screening and treatment options for cancer patients. Although this issue has been studied extensively for chemotherapy, less is known about a causal link between resistance to radiation therapy and metastasis. We investigated this problem in triple-negative breast cancer (TNBC) and established that radiation resistant tumor cells have enhanced metastatic capacity, especially to bone. Resistance to radiation increases the expression of integrin β3 (ITGβ3), which promotes enhanced migration and invasion. Bioinformatic analysis and subsequent experimentation revealed an enrichment of RNA metabolism pathways that stabilize ITGβ3 transcripts. Specifically, the RNA binding protein heterogenous nuclear ribonucleoprotein L (HNRNPL), whose expression is regulated by Nrf2, mediates the formation of circular RNAs (circRNAs) that function as competing endogenous RNAs (ceRNAs) for the family of let-7 microRNAs that target ITGβ3. Collectively, our findings identify a novel mechanism of radiation-induced metastasis that is driven by alterations in RNA metabolism.

Multiplexed assays of human disease-relevant mutations reveal UTR dinucleotide composition as a major determinant of RNA stability

Untranslated regions (UTRs) contain crucial regulatory elements for RNA stability, translation and localization, so their integrity is indispensable for gene expression. Approximately 3.7% of genetic variants associated with diseases occur in UTRs, yet a comprehensive understanding of UTR variant functions remains limited due to inefficient experimental and computational assessment methods. To systematically evaluate the effects of UTR variants on RNA stability, we established a massively parallel reporter assay on 6555 UTR variants reported in human disease databases. We examined the RNA degradation patterns mediated by the UTR library in two cell lines, and then applied LASSO regression to model the influential regulators of RNA stability. We found that UA dinucleotides and UA-rich motifs are the most prominent destabilizing element. Gain of UA dinucleotide outlined mutant UTRs with reduced stability. Studies on endogenous transcripts indicate that high UA-dinucleotide ratios in UTRs promote RNA degradation. Conversely, elevated GC content and protein binding on UA dinucleotides protect high-UA RNA from degradation. Further analysis reveals polarized roles of UA-dinucleotide-binding proteins in RNA protection and degradation. Furthermore, the UA-dinucleotide ratio of both UTRs is a common characteristic of genes in innate immune response pathways, implying a coordinated stability regulation through UTRs at the transcriptomic level. We also demonstrate that stability-altering UTRs are associated with changes in biobank-based health indices, underscoring the importance of precise UTR regulation for wellness. Our study highlights the importance of RNA stability regulation through UTR primary sequences, paving the way for further exploration of their implications in gene networks and precision medicine.

Characterization of RNF144B and PPP2R2A identified by a novel approach using TCGA data in ovarian cancer

TCGA has identified predominant somatic copy number alterations (SCNA) affecting numerous genes in HGSOC. To identify cancer-driver genes from the regions of SCNA, we have devised a scoring system that integrates information from different genetic alterations. Applying this scoring system to the TCGA-HGSOC dataset (n = 316) we have identified several well-known and novel putative cancer genes in HGSOC. We functionally validated the roles of two previously unknown genes, RNF144B and PPP2R2A. RNF144B, an E3 ubiquitin-ligase is amplified and overexpressed in 16% of HGSOC (TCGA). Overexpression of RNF144B in ovarian cancer cells increased cell proliferation, colony formation, and migration. RNF144B was significantly overexpressed in 50% of primary tumors from patients with HGSOC compared to the ovary. Further, it had significantly reduced expression in tumors after chemotherapy. PPP2R2A, the regulatory subunit of PP2A is deleted and downregulated in 38% of HGSOCs (TCGA). Overexpression of PPP2R2A inhibited cell proliferation, colony-formation, migration, and invasion in ovarian cancer cells. In OVCAR-5, which expresses low levels of PPP2R2A, Niraparib inhibited cell proliferation. PPP2R2A was not expressed in 72% of HGSOCs. This report demonstrates this approach to identifying genes from the TCGA data. Further experiments are required to conclusively prove the role of these genes in the pathogenesis of ovarian cancer.

RNA Binding of GAPDH Controls Transcript Stability and Protein Translation in Acute Myeloid Leukemia

Dysregulation of RNA binding proteins (RBPs) is a hallmark in cancerous cells. In acute myeloid leukemia (AML) RBPs are key regulators of tumor proliferation. While classical RBPs have defined RNA binding domains, RNA recognition and function in AML by non-canonical RBPs (ncRBPs) remain unclear. Given the inherent complexity of targeting AML broadly, our goal was to uncover potential ncRBP candidates critical for AML survival using a CRISPR/Cas-based screening. We identified the glycolytic enzyme glyceraldehyde-3-phosphate dehydrogenase (GAPDH) as a pro-proliferative factor in AML cells. Based on cross-linking and immunoprecipitation (CLIP), we are defining the global targetome, detecting novel RNA targets mainly located within 5'UTRs, including GAPDH, RPL13a, and PKM. The knockdown of GAPDH unveiled genetic pathways related to ribosome biogenesis, translation initiation, and regulation. Moreover, we demonstrated a stabilizing effect through GAPDH binding to target transcripts including its own mRNA. The present findings provide new insights on the RNA functions and characteristics of GAPDH in AML.

Deep generative AI models analyzing circulating orphan non-coding RNAs enable detection of early-stage lung cancer

Liquid biopsies have the potential to revolutionize cancer care through non-invasive early detection of tumors. Developing a robust liquid biopsy test requires collecting high-dimensional data from a large number of blood samples across heterogeneous groups of patients. We propose that the generative capability of variational auto-encoders enables learning a robust and generalizable signature of blood-based biomarkers. In this study, we analyze orphan non-coding RNAs (oncRNAs) from serum samples of 1050 individuals diagnosed with non-small cell lung cancer (NSCLC) at various stages, as well as sex-, age-, and BMI-matched controls. We demonstrate that our multi-task generative AI model, Orion, surpasses commonly used methods in both overall performance and generalizability to held-out datasets. Orion achieves an overall sensitivity of 94% (95% CI: 87%-98%) at 87% (95% CI: 81%-93%) specificity for cancer detection across all stages, outperforming the sensitivity of other methods on held-out validation datasets by more than  ~ 30%.

Paraspeckle-independent co-transcriptional regulation of nuclear microRNA biogenesis by SFPQ

MicroRNAs (miRNAs) play crucial roles in physiological functions and disease, but the regulation of their nuclear biogenesis remains poorly understood. Here, BioID on Drosha, the catalytic subunit of the microprocessor complex, reveals its proximity to splicing factor proline- and glutamine (Q)-rich (SFPQ), a multifunctional RNA-binding protein (RBP) involved in forming paraspeckle nuclear condensates. SFPQ depletion impacts both primary and mature miRNA expression, while other paraspeckle proteins (PSPs) or the paraspeckle scaffolding RNA NEAT1 do not, indicating a paraspeckle-independent role. Comprehensive transcriptomic analyses show that SFPQ loss broadly affects RNAs and miRNA host gene (HG) expression, influencing both their transcription and the stability of their products. Notably, SFPQ protects the oncogenic miR-17∼92 polycistron from degradation by the nuclear exosome targeting (NEXT)-exosome complex and is tightly linked with its overexpression across a broad variety of cancers. Our findings reveal a dual role for SFPQ in regulating miRNA HG transcription and stability, as well as its significance in cancers.

A unified model for interpretable latent embedding of multi-sample, multi-condition single-cell data

Single-cell analysis across multiple samples and conditions requires quantitative modeling of the interplay between the continuum of cell states and the technical and biological sources of sample-to-sample variability. We introduce GEDI, a generative model that identifies latent space variations in multi-sample, multi-condition single-cell datasets and attributes them to sample-level covariates. GEDI enables cross-sample cell state mapping on par with state-of-the-art integration methods, cluster-free differential gene expression analysis along the continuum of cell states, and machine learning-based prediction of sample characteristics from single-cell data. GEDI can also incorporate gene-level prior knowledge to infer pathway and regulatory network activities in single cells. Finally, GEDI extends all these concepts to previously unexplored modalities that require joint consideration of dual measurements, such as the joint analysis of exon inclusion/exclusion reads to model alternative cassette exon splicing, or spliced/unspliced reads to model the mRNA stability landscapes of single cells.

The E3 ubiquitin ligase Nedd4L preserves skeletal muscle stem cell quiescence by inhibiting their activation

Adult stem cells play a critical role in tissue repair and maintenance. In tissues with slow turnover, including skeletal muscle, these cells are maintained in a mitotically quiescent state yet remain poised to re-enter the cell cycle to replenish themselves and regenerate the tissue. Using a panomics approach we show that the PAX7/NEDD4L axis acts against muscle stem cell activation in homeostatic skeletal muscle. Our findings suggest that PAX7 transcriptionally activates the E3 ubiquitin ligase Nedd4L and that the conditional genetic deletion of Nedd4L impairs muscle stem cell quiescence, with an upregulation of cell cycle and myogenic differentiation genes. Loss of Nedd4L in muscle stem cells results in the expression of doublecortin (DCX), which is exclusively expressed during their in vivo activation. Together, these data establish that the ubiquitin proteasome system, mediated by Nedd4L, is a key contributor to the muscle stem cell quiescent state in adult mice.

Dynamic RNA polymerase II occupancy drives differentiation of the intestine under the direction of HNF4

Terminal differentiation requires massive restructuring of the transcriptome. During intestinal differentiation, the expression patterns of nearly 4,000 genes are altered as cells transition from progenitor cells in crypts to differentiated cells in villi. We identify dynamic occupancy of RNA polymerase II (Pol II) to gene promoters as the primary driver of transcriptomic shifts during intestinal differentiation in vivo. Changes in enhancer-promoter looping interactions accompany dynamic Pol II occupancy and are dependent upon HNF4, a pro-differentiation transcription factor. Using genetic loss-of-function, chromatin immunoprecipitation sequencing (ChIP-seq), and immunoprecipitation (IP) mass spectrometry, we demonstrate that HNF4 collaborates with chromatin remodelers and loop-stabilizing proteins and facilitates Pol II occupancy at hundreds of genes pivotal to differentiation. We also explore alternate mechanisms that drive differentiation gene expression and find that pause-release of Pol II and post-transcriptional mRNA stability regulate smaller subsets of differentially expressed genes. These studies provide insights into the mechanisms of differentiation in renewing adult tissue.

Full-spectral genome analysis of natural killer/T cell lymphoma highlights impacts of genome instability in driving its progression

BACKGROUND

Natural killer/T cell lymphoma (NKTCL) is a clinically and genetically heterogeneous disease with poor prognosis. Genome sequencing and mutation characterization provides a powerful approach for patient stratification, treatment target discovery, and etiology identification. However, previous studies mostly concentrated on base-level mutations in primary NKTCL, whereas the large-scale genomic alterations in NKTCL and the mutational landscapes in relapsed/refractory NKTCL remain largely unexplored.

METHODS

Here, we assembled whole-genome sequencing and whole-exome sequencing data from 163 patients with primary or relapsed/refractory NKTCL and compared their somatic mutational landscapes at both nucleotide and structure levels.

RESULTS

Our study not only confirmed previously reported common NKTCL mutational targets like STAT3, TP53, and DDX3X but also unveiled several novel high-frequency mutational targets such as PRDM9, DST, and RBMX. In terms of the overall mutational landscape, we observed striking differences between primary and relapsed/refractory NKTCL patient groups, with the latter exhibits higher levels of tumor mutation burden, copy number variants (CNVs), and structural variants (SVs), indicating a strong signal of genomic instability. Complex structural rearrangements such as chromothripsis and focal amplification are also significantly enriched in relapsed/refractory NKTCL patients, exerting a substantial impact on prognosis. Accordingly, we devised a novel molecular subtyping system (i.e., C0-C4) with distinct prognosis by integrating potential driver mutations at both nucleotide and structural levels, which further provides an informative guidance for novel treatments that target these specific driver mutations and genome instability as a whole.

CONCLUSIONS

The striking differences underlying the mutational landscapes between the primary and relapsed/refractory NKTCL patients highlight the importance of genomic instability in driving the progression of NKTCL. Our newly proposed molecular subtyping system is valuable in assisting patient stratification and novel treatment design towards a better prognosis in the age of precision medicine.

Dynamic RNA Polymerase II Recruitment Drives Differentiation of the Intestine under the direction of HNF4

Terminal differentiation requires a massive restructuring of the transcriptome. During intestinal differentiation, the expression patterns of nearly 4000 genes are altered as cells transition from progenitor cells in crypts to differentiated cells in villi. We identified dynamic recruitment of RNA Polymerase II (Pol II) to gene promoters as the primary driver of transcriptomic shifts during intestinal differentiation in vivo. Changes in enhancer-promoter looping interactions accompany dynamic Pol II recruitment and are dependent upon HNF4, a pro-differentiation transcription factor. Using genetic loss-of- function, ChIP-seq and IP mass spectrometry, we demonstrate that HNF4 collaborates with chromatin remodelers and loop-stabilizing proteins and facilitates Pol II recruitment at hundreds of genes pivotal to differentiation. We also explore alternate mechanisms which drive differentiation gene expression and find pause-release of Pol II and post- transcriptional mRNA stability regulate smaller subsets of differentially expressed genes. These studies provide insights into the mechanisms of differentiation in a renewing adult tissue.

Missense mutations in Myc Box I influence MYC cellular localization, mRNA partitioning and turnover to promote leukemogenesis

Somatic missense mutations in the phosphodegron domain of the MYC gene ( M YC Box I) are detected in the dominant clones of a subset of acute myeloid leukemia (AML) patients, but the mechanisms by which they contribute to AML are unknown. To unveil unique proprieties of MBI MYC mutant proteins, we systematically compared the cellular and molecular consequences of expressing similar oncogenic levels of wild type and MBI mutant MYC. We found that MBI MYC mutants can accelerate leukemia by driving unique transcriptional signatures in highly selected, myeloid progenitor subpopulations. Although these mutations increase MYC stability, they overall dampen MYC chromatin localization and lead to a cytoplasmic accumulation of the mutant proteins. This phenotype is coupled with increased translation of RNA binding proteins and nuclear export machinery, which results in altered RNA partitioning and accelerated decay of select transcripts encoding proapoptotic and proinflammatory genes. Heterozygous knockin mice harboring the germline MBI mutation Myc p.T73N exhibit cytoplasmic MYC localization, myeloid progenitors' expansion with similar transcriptional signatures to the overexpression model, and eventually develop hematological malignancies. This study uncovers that MBI MYC mutations alter MYC localization and disrupt mRNA subcellular distribution and turnover of select transcripts to accelerate tumor initiation and growth.