Deep and accurate detection of m6A RNA modifications using miCLIP2 and m6Aboost machine learning.

N6-methyladenosine (m6A) is the most abundant internal RNA modification in eukaryotic mRNAs and influences many aspects of RNA processing. miCLIP (m6A individual-nucleotide resolution UV crosslinking and immunoprecipitation) is an antibody-based approach to map m6A sites with single-nucleotide resolution. However, due to broad antibody reactivity, reliable identification of m6A sites from miCLIP data remains challenging. Here, we present miCLIP2 in combination with machine learning to significantly improve m6A detection. The optimized miCLIP2 results in high-complexity libraries from less input material. Importantly, we established a robust computational pipeline to tackle the inherent issue of false positives in antibody-based m6A detection. The analyses were calibrated with Mettl3 knockout cells to learn the characteristics of m6A deposition, including m6A sites outside of DRACH motifs. To make our results universally applicable, we trained a machine learning model, m6Aboost, based on the experimental and RNA sequence features. Importantly, m6Aboost allows prediction of genuine m6A sites in miCLIP2 data without filtering for DRACH motifs or the need for Mettl3 depletion. Using m6Aboost, we identify thousands of high-confidence m6A sites in different murine and human cell lines, which provide a rich resource for future analysis. Collectively, our combined experimental and computational methodology greatly improves m6A identification.

The oncofetal RNA-binding protein IGF2BP1 is a druggable, post-transcriptional super-enhancer of E2F-driven gene expression in cancer.

The IGF2 mRNA-binding protein 1 (IGF2BP1) is a non-catalytic post-transcriptional enhancer of tumor growth upregulated and associated with adverse prognosis in solid cancers. However, conserved effector pathway(s) and the feasibility of targeting IGF2BP1 in cancer remained elusive. We reveal that IGF2BP1 is a post-transcriptional enhancer of the E2F-driven hallmark in solid cancers. IGF2BP1 promotes G1/S cell cycle transition by stabilizing mRNAs encoding positive regulators of this checkpoint like E2F1. This IGF2BP1-driven shortening of the G1 cell cycle phase relies on 3'UTR-, miRNA- and m6A-dependent regulation and suggests enhancement of cell cycle progression by m6A-modifications across cancers. In addition to E2F transcription factors, IGF2BP1 also stabilizes E2F-driven transcripts directly indicating post-transcriptional 'super'-enhancer role of the protein in E2F-driven gene expression in cancer. The small molecule BTYNB disrupts this enhancer function by impairing IGF2BP1-RNA association. Consistently, BTYNB interferes with E2F-driven gene expression and tumor growth in experimental mouse tumor models.

Comprehensive analysis of translation from overexpressed circular RNAs reveals pervasive translation from linear transcripts.

Circular RNAs (circRNAs) encompass a widespread and conserved class of RNAs, which are generated by back-splicing of downstream 5' to upstream 3' splice sites. CircRNAs are tissue-specific and have been implicated in diseases including cancer. They can function as sponges for microRNAs (miRNAs) or RNA binding proteins (RBPs), for example. Moreover, some contain open reading frames (ORFs) and might be translated. The functional relevance of such peptides, however, remains largely elusive. Here, we report that the ORF of circZNF609 is efficiently translated when expressed from a circZNF609 overexpression construct. However, endogenous proteins could not be detected. Moreover, initiation of circZNF609 translation is independent of m6A-generating enzyme METTL3 or RNA sequence elements such as internal ribosome entry sites (IRESs). Surprisingly, a comprehensive mutational analysis revealed that deletion constructs, which are deficient in producing circZNF609, still generate the observed protein products. This suggests that the apparent circZNF609 translation originates from trans-splicing by-products of the overexpression plasmids and underline that circRNA overexpression constructs need to be evaluated carefully, particularly when functional studies are performed.

IGF2BP1 is the first positive marker for anaplastic thyroid carcinoma diagnosis.

Anaplastic thyroid carcinomas (ATC) are rare, but represent the most lethal malignancy of the thyroid. Selective molecular markers and drivers distinguishing ATC from other thyroid carcinomas of follicular origin remain largely unknown, limiting advances in diagnosis and treatment. In a retrospective study, we analyzed gene expression in 36 ATC, 18 poorly differentiated, 132 papillary, and 55 follicular thyroid carcinoma, as well as 124 paired and unpaired normal thyroid tissues in three independent cohorts by RNA-sequencing and immunohistochemistry. RNA-sequencing data in the test cohort suggested selective ATC protein biomarkers. Evaluation of these revealed that ATCs are characterized by the de novo expression of various testis antigens, including melanoma-associated antigen A3 (MAGEA3), but most importantly the oncofetal IGF2 mRNA binding protein 1 (IGF2BP1). Shallow whole genome sequencing essentially excluded that IGF2BP1 upregulation results from gene copy number alterations. Immunohistochemical analyses in all three tumor cohorts confirmed the selective de novo expression of IGF2BP1 protein in ATC. In sum, 75% (27/36) of all tested ATC and 0.5% (1/204) of poorly and well-differentiated thyroid carcinoma tissue samples were positive for IGF2BP1 protein. This indicates that IGF2BP1 protein expression identifies ATC with a diagnostic odds ratio of 612 (95% CI: 74.6-5021). In addition, we found that MAGEA3 is exclusively, although less consistently upregulated in ATC, presenting with an odds ratio of 411 (95% CI: 23.8-7098.7). Importantly, we provide confirmatory evidence that IGF2BP1 and MAGEA3 expression distinguishes ATC from poorly differentiated thyroid carcinoma. IGF2BP1 furthermore identified ATC foci within low-grade follicular thyroid carcinoma. In conclusion, IGF2BP1 represents the most promising single-gene marker available for ATC, followed by MAGEA3, improving on current techniques. Robust markers are essential to help distinguish this high-grade malignancy from other thyroid carcinomas, to guide surgical decision making, therapy and post-resection/therapy monitoring strategies.

IGF2BP1 promotes SRF-dependent transcription in cancer in a m6A- and miRNA-dependent manner.

The oncofetal mRNA-binding protein IGF2BP1 and the transcriptional regulator SRF modulate gene expression in cancer. In cancer cells, we demonstrate that IGF2BP1 promotes the expression of SRF in a conserved and N6-methyladenosine (m6A)-dependent manner by impairing the miRNA-directed decay of the SRF mRNA. This results in enhanced SRF-dependent transcriptional activity and promotes tumor cell growth and invasion. At the post-transcriptional level, IGF2BP1 sustains the expression of various SRF-target genes. The majority of these SRF/IGF2BP1-enhanced genes, including PDLIM7 and FOXK1, show conserved upregulation with SRF and IGF2BP1 synthesis in cancer. PDLIM7 and FOXK1 promote tumor cell growth and were reported to enhance cell invasion. Consistently, 35 SRF/IGF2BP1-dependent genes showing conserved association with SRF and IGF2BP1 expression indicate a poor overall survival probability in ovarian, liver and lung cancer. In conclusion, these findings identify the SRF/IGF2BP1-, miRNome- and m6A-dependent control of gene expression as a conserved oncogenic driver network in cancer.

Telomerase reactivation induces progression of mouse Braf V600E -driven thyroid cancers without telomere lengthening.

Mutations in the promoter of the telomerase reverse transcriptase ( TERT ) gene are the paradigm of a cross-cancer alteration in a non-coding region. TERT promoter mutations (TPMs) are biomarkers of poor prognosis in several tumors, including thyroid cancers. TPMs enhance TERT transcription, which is otherwise silenced in adult tissues, thus reactivating a bona fide oncoprotein. To study TERT deregulation and its downstream consequences, we generated a Tert mutant promoter mouse model via CRISPR/Cas9 engineering of the murine equivalent locus (Tert -123C>T ) and crossed it with thyroid-specific Braf V600E -mutant mice. We also employed an alternative model of Tert overexpression (K5-Tert). Whereas all Braf V600E animals developed well-differentiated papillary thyroid tumors, 29% and 36% of Braf V600E +Tert -123C>T and Braf V600E +K5-Tert mice progressed to poorly differentiated thyroid cancers at week 20, respectively. Braf+Tert tumors showed increased mitosis and necrosis in areas of solid growth, and older animals from these cohorts displayed anaplastic-like features, i.e., spindle cells and macrophage infiltration. Murine Tert promoter mutation increased Tert transcription in vitro and in vivo , but temporal and intra-tumoral heterogeneity was observed. RNA-sequencing of thyroid tumor cells showed that processes other than the canonical Tert-mediated telomere maintenance role operate in these specimens. Pathway analysis showed that MAPK and PI3K/AKT signaling, as well as processes not previously associated with this tumor etiology, involving cytokine and chemokine signaling, were overactivated. Braf+Tert animals remained responsive to MAPK pathway inhibitors. These models constitute useful pre-clinical tools to understand the cell-autonomous and microenvironment-related consequences of Tert-mediated progression in advanced thyroid cancers and other aggressive tumors carrying TPMs.

Papillary Thyroid Carcinomas with Tall Cell Features: An Intermediate Entity Between Classic and Tall Cell Subtypes.

Background: While the diagnosis of papillary thyroid carcinomas (PTCs) with tall cell features (PTCtcf) is often made for carcinomas with histological features intermediate between classic and tall cell subtypes of PTC (tcPTC), its comparative signature to that of either tcPTC or classic PTC is less clear. The objective of this study was to perform an integrative clinicopathologic and genomic analysis elucidating the spectrum of tcPTC, PTCtcf, and classic PTC. Methods: We analyzed all consecutive patients with tcPTC and PTCtcf evaluated at a tertiary academic referral center between 2005 and 2020, as well as a comparative cohort of classic PTC, in a retrospective observational cohort analysis. Clinicopathologic data were compared among the three groups, including progression-free survival (PFS), recurrent/persistent disease, and a negative composite outcome of death, progression, or need for advanced therapy. To specifically understand differences between tcPTC and PTCtcf, targeted next-generation sequencing was performed in a subset of these cohorts. Results: A total of 292 patients were analyzed (81 tcPTC, 65 PTCtcf, 146 classic PTC). Thirteen percent of tcPTC versus 8% of PTCtcf versus 1% of classic PTC had the advanced American Joint Committee on Cancer stage (p = 0.002). Similarly, macroscopic extrathyroidal extension was observed in 38% of tcPTC, 14% of PTCtcf, and 12% of classic PTC (p < 0.001). The 5-year PFS was 76.5%, 81.5%, and 88.3% for tcPTC, PTCtcf, and classic PTC, respectively, while the rates of the negative composite outcome 40.2% for tcPTC, 20.7% for PTCtcf, and 11.2% for classic PTC (p < 0.001). In a multivariable Cox regression analysis, the negative composite outcome was independently associated with tcPTC (HR 4.3 [confidence interval 1.1-16.1], p = 0.03). tcPTC demonstrated substantially more hotspot TERT promoter mutations than PTCtcf (44% vs. 6%, p = 0.012). Conclusions: Our study demonstrates a continuum of disease-specific risk of PTC, pointing at PTCtcf as an intermediate entity between tcPTC and classic PTC. These data provide a more refined understanding of risk at time of presentation, while better elucidating the diversity of genomic drivers.

Telomerase Upregulation Induces Progression of Mouse BrafV600E-Driven Thyroid Cancers and Triggers Nontelomeric Effects.

Mutations in the promoter of the telomerase reverse transcriptase (TERT) gene are the paradigm of a cross-cancer alteration in a noncoding region. TERT promoter mutations (TPM) are biomarkers of poor prognosis in cancer, including thyroid tumors. TPMs enhance TERT transcription, which is otherwise silenced in adult tissues, thus reactivating a bona fide oncoprotein. To study TERT deregulation and its downstream consequences, we generated a Tert mutant promoter mouse model via CRISPR/Cas9 engineering of the murine equivalent locus (Tert-123C>T) and crossed it with thyroid-specific BrafV600E-mutant mice. We also employed an alternative model of Tert overexpression (K5-Tert). Whereas all BrafV600E animals developed well-differentiated papillary thyroid tumors, 29% and 36% of BrafV600E+Tert-123C>T and BrafV600E+K5-Tert mice progressed to poorly differentiated cancers at week 20, respectively. Tert-upregulated tumors showed increased mitosis and necrosis in areas of solid growth, and older animals displayed anaplastic-like features, that is, spindle cells and macrophage infiltration. Murine TPM increased Tert transcription in vitro and in vivo, but temporal and intratumoral heterogeneity was observed. RNA-sequencing of thyroid tumor cells showed that processes other than the canonical Tert-mediated telomere maintenance role operate in these specimens. Pathway analysis showed that MAPK and PI3K/AKT signaling, as well as processes not previously associated with this tumor etiology, involving cytokine, and chemokine signaling, were overactivated. These models constitute useful preclinical tools to understand the cell-autonomous and microenvironment-related consequences of Tert-mediated progression in advanced thyroid cancers and other aggressive tumors carrying TPMs. IMPLICATIONS: Telomerase-driven cancer progression activates pathways that can be dissected and perhaps therapeutically exploited.

Additional Oncogenic Alterations in RAS-Driven Differentiated Thyroid Cancers Associate with Worse Clinicopathologic Outcomes.

PURPOSE: RAS mutations occur across the spectrum of thyroid neoplasms, and more tools are needed for better prognostication. The objective of this study was to evaluate how additional genetic events affecting key genes modify prognosis in patients with RAS-mutant thyroid cancers, and specifically differentiated thyroid cancers (DTC). EXPERIMENTAL DESIGN: We performed a clinical-genomic analysis of consecutive patients with DTC, poorly differentiated (PDTC), or anaplastic thyroid cancer (ATC) between January 2014 and December 2021, in whom a custom-targeted next-generation sequencing assay was performed. Patients harboring RAS mutations were included, and we compared their clinical features and outcomes based upon the presence of additional oncogenic alterations. RESULTS: Seventy-eight patients were identified, with 22% (17/78) harboring a driver RAS mutation plus an additional oncogenic alteration. All six (100%) ATCs had an additional mutation. Compared with DTCs harboring a solitary RAS mutation, patients with DTC with RAS and additional mutation(s) were more likely to be classified as American Thyroid Association high-risk of recurrence (77% vs. 12%; P < 0.001) and to have larger primary tumors (4.7 vs. 2.5 cm; P = 0.002) and advanced stage (III or IV) at presentation (67% vs. 3%; P < 0.001). Importantly, over an average 65-month follow-up, DTC-specific-mortality was more than 10-fold higher (20% vs. 1.8%; P = 0.011) when additional mutations were identified. CONCLUSIONS: Identification of key additional mutations in patients with RAS-mutant thyroid cancers confers a more aggressive phenotype, increases mortality risk in DTC, and can explain the diversity of RAS-mutated thyroid neoplasia. These data support genomic profiling of DTCs to inform prognosis and clinical decision-making.

MiRNA Deregulation Distinguishes Anaplastic Thyroid Carcinoma (ATC) and Supports Upregulation of Oncogene Expression.

Anaplastic thyroid carcinoma (ATC) is the most fatal and rapidly evolving endocrine malignancy invading the head and neck region and accounts for up to 50% of thyroid cancer-associated deaths. Deregulation of the microRNA (miRNA) expression promotes thyroid carcinoma progression by modulating the reorganization of the ATC transcriptome. Here, we applied comparative miRNA-mRNA sequencing on a cohort of 28 thyroid carcinomas to unravel the association of deregulated miRNA and mRNA expression. This identified 85 miRNAs significantly deregulated in ATC. By establishing a new analysis pipeline, we unraveled 85 prime miRNA-mRNA interactions supporting the downregulation of candidate tumor suppressors and the upregulation of bona fide oncogenes such as survivin (BIRC5) in ATC. This miRNA-dependent reprogramming of the ATC transcriptome provided an mRNA signature comprising 65 genes sharply distinguishing ATC from other thyroid carcinomas. The validation of the deregulated protein expression in an independent thyroid carcinoma cohort demonstrates that miRNA-dependent oncogenes comprised in this signature, the transferrin receptor TFRC (CD71) and the E3-ubiquitin ligase DTL, are sharply upregulated in ATC. This upregulation is sufficient to distinguish ATC even from poorly differentiated thyroid carcinomas (PDTC). In sum, these findings provide new diagnostic tools and a robust resource to explore the key miRNA-mRNA regulation underlying the progression of thyroid carcinoma.

MSI1 Promotes the Expression of the GBM Stem Cell Marker CD44 by Impairing miRNA-Dependent Degradation.

The stem cell marker Musashi1 (MSI1) is highly expressed during neurogenesis and in glioblastoma (GBM). MSI1 promotes self-renewal and impairs differentiation in cancer and non-malignant progenitor cells. However, a comprehensive understanding of its role in promoting GBM-driving networks remains to be deciphered. We demonstrate that MSI1 is highly expressed in GBM recurrences, an oncologist's major defiance. For the first time, we provide evidence that MSI1 promotes the expression of stem cell markers like CD44, co-expressed with MSI1 within recurrence-promoting cells at the migrating front of primary GBM samples. With GBM cell models of pediatric and adult origin, including isolated primary tumorspheres, we show that MSI1 promotes stem cell-like characteristics. Importantly, it impairs CD44 downregulation in a 3'UTR- and miRNA-dependent manner by controlling mRNA turnover. This regulation is disturbed by the previously reported MSI1 inhibitor luteolin, providing further evidence for a therapeutic target potential of MSI1 in GBM treatment.

Synthetic circular miR-21 RNA decoys enhance tumor suppressor expression and impair tumor growth in mice.

Naturally occurring circular RNAs efficiently impair miRNA functions. Synthetic circular RNAs may thus serve as potent agents for miRNA inhibition. Their therapeutic effect critically relies on (i) the identification of optimal miRNA targets, (ii) the optimization of decoy structures and (iii) the development of efficient formulations for their use as drugs. In this study, we extensively explored the functional relevance of miR-21-5p in cancer cells. Analyses of cancer transcriptomes reveal that miR-21-5p is the by far most abundant miRNA in human cancers. Deletion of the MIR21 locus in cancer-derived cells identifies several direct and indirect miR-21-5p targets, including major tumor suppressors with prognostic value across cancers. To impair miR-21-5p activities, we evaluate synthetic, circular RNA decoys containing four repetitive binding elements. In cancer cells, these decoys efficiently elevate tumor suppressor expression and impair tumor cell vitality. For their in vivo delivery, we for the first time evaluate the formulation of decoys in polyethylenimine (PEI)-based nanoparticles. We demonstrate that PEI/decoy nanoparticles lead to a significant inhibition of tumor growth in a lung adenocarcinoma xenograft mouse model via the upregulation of tumor suppressor expression. These findings introduce nanoparticle-delivered circular miRNA decoys as a powerful potential therapeutic strategy in cancer treatment.

The Critical Role of RNA m6A Methylation in Cancer.

Since the identification of the first RNA demethylase and the establishment of methylated RNA immunoprecipitation-sequencing methodology 6 to 7 years ago, RNA methylation has emerged as a widespread phenomenon and a critical regulator of transcript expression. This new layer of regulation is termed "epitranscriptomics." The most prevalent RNA methylation, N 6-methyladenosine (m6A), occurs in approximately 25% of transcripts at the genome-wide level and is enriched around stop codons, in 5'- and 3'-untranslated regions, and within long internal exons. RNA m6A modification regulates RNA splicing, translocation, stability, and translation into protein. m6A is catalyzed by the RNA methyltransferases METTL3, METTL14, and METTL16 (writers), is removed by the demethylases FTO and ALKBH5 (erasers), and interacts with m6A-binding proteins, such as YTHDF1 and IGF2BP1 (readers). RNA methyltransferases, demethylases, and m6A-binding proteins are frequently upregulated in human cancer tissues from a variety of organ origins, increasing onco-transcript and oncoprotein expression, cancer cell proliferation, survival, tumor initiation, progression, and metastasis. Although RNA methyltransferase inhibitors are not available yet, FTO inhibitors have shown promising anticancer effects in vitro and in animal models of cancer. Further screening for selective and potent RNA methyltransferase, demethylase, or m6A-binding protein inhibitors may lead to compounds suitable for future clinical trials in cancer patients.