MAML3-fusions modulate vascular and immune tumour microenvironment and confer high metastatic risk in pheochromocytoma and paraganglioma.

Pheochromocytomas and paragangliomas are rare neuroendocrine tumours. Around 20-25 % of patients develop metastases, for which there is an urgent need of prognostic markers and therapeutic stratification strategies. The presence of a MAML3-fusion is associated with increased metastatic risk, but neither the processes underlying disease progression, nor targetable vulnerabilities have been addressed. We have compiled a cohort of 850 patients, which has shown a 3.65 % fusion prevalence and represents the largest MAML3-positive series reported to date. While MAML3-fusions mainly cause single pheochromocytomas, we also observed somatic post-zygotic events, resulting in multiple tumours in the same patient. MAML3-tumours show increased expression of neuroendocrine-to-mesenchymal transition markers, MYC-targets, and angiogenesis-related genes, leading to a distinct tumour microenvironment with unique vascular and immune profiles. Importantly, our findings have identified MAML3-tumours specific vulnerabilities beyond Wnt-pathway dysregulation, such as a rich vascular network, and overexpression of PD-L1 and CD40, suggesting potential therapeutic targets.

SpadaHC: a database to improve the classification of variants in hereditary cancer genes in the Spanish population.

Accurate classification of genetic variants is crucial for clinical decision-making in hereditary cancer. In Spain, genetic diagnostic laboratories have traditionally approached this task independently due to the lack of a dedicated resource. Here we present SpadaHC, a web-based database for sharing variants in hereditary cancer genes in the Spanish population. SpadaHC is implemented using a three-tier architecture consisting of a relational database, a web tool and a bioinformatics pipeline. Contributing laboratories can share variant classifications and variants from individuals in Variant Calling Format (VCF) format. The platform supports open and restricted access, flexible dataset submissions, automatic pseudo-anonymization, VCF quality control, variant normalization and liftover between genome builds. Users can flexibly explore and search data, receive automatic discrepancy notifications and access SpadaHC population frequencies based on many criteria. In February 2024, SpadaHC included 18 laboratory members, storing 1.17 million variants from 4306 patients and 16 343 laboratory classifications. In the first analysis of the shared data, we identified 84 genetic variants with clinically relevant discrepancies in their classifications and addressed them through a three-phase resolution strategy. This work highlights the importance of data sharing to promote consistency in variant classifications among laboratories, so patients and family members can benefit from more accurate clinical management. Database URL: https://spadahc.ciberisciii.es/.

Clinical and molecular markers guide the genetics of pheochromocytoma and paraganglioma.

Over the past two decades, research into the genetic susceptibility behind pheochromocytoma and paraganglioma (PPGL) has surged, ranking them among the most heritable tumors. Massive sequencing combined with careful patient selection has so far identified more than twenty susceptibility genes, leading to an over-detection of variants of unknown significance (VUS) that require precise molecular markers to determine their pathogenic role. Moreover, some PPGL patients remain undiagnosed, possibly due to mutations in regulatory regions of already known genes or mutations in undiscovered genes. Accurate classification of VUS and identification of new genes require well-defined clinical and molecular markers that allow effective genetic diagnosis of most PPGLs.

Defective function of α-ketoglutarate dehydrogenase exacerbates mitochondrial ATP deficits during complex I deficiency.

The NDUFS4 knockout (KO) mouse phenotype resembles the human Complex I deficiency Leigh Syndrome. The irreversible succination of protein thiols by fumarate is increased in select regions of the NDUFS4 KO brain affected by neurodegeneration. We report that dihydrolipoyllysine-residue succinyltransferase (DLST), a component of the α-ketoglutarate dehydrogenase complex (KGDHC) of the tricarboxylic acid (TCA) cycle, is succinated in the affected regions of the NDUFS4 KO brain. Succination of DLST reduced KGDHC activity in the brainstem (BS) and olfactory bulb (OB) of KO mice. The defective production of KGDHC derived succinyl-CoA resulted in decreased mitochondrial substrate level phosphorylation (SLP), further aggravating the existing oxidative phosphorylation (OXPHOS) ATP deficit. Protein succinylation, an acylation modification that requires succinyl-CoA, was reduced in the KO mice. Modeling succination of a cysteine in the spatial vicinity of the DLST active site or introduction of succinomimetic mutations recapitulates these metabolic deficits. Our data demonstrate that the biochemical deficit extends beyond impaired Complex I assembly and OXPHOS deficiency, functionally impairing select components of the TCA cycle to drive metabolic perturbations in affected neurons.

Deubiquitinase USP9X loss sensitizes renal cancer cells to mTOR inhibition.

Mammalian target of rapamycin (mTOR) is a central regulator of mammalian metabolism and physiology. Aberrant hyperactivation of the mTOR pathway promotes tumor growth and metastasis, and can also promote tumor resistance to chemotherapy and cancer drugs; this makes mTOR an attractive cancer therapeutic target. mTOR inhibitors have been approved to treat cancer; however, the mechanisms underlying drug sensitivity remain poorly understood. Here, whole exome sequencing of three chromophobe renal cell carcinoma (chRCC) patients with exceptional mTOR inhibitor sensitivity revealed that all three patients shared somatic mutations in the deubiquitinase gene USP9X. The clonal characteristics of the mutations, which were amassed by studying multiple patients' primary and metastatic samples from various years, together with the low USP9X mutation rate in unselected chRCC series, reinforced a causal link between USP9X and mTOR inhibitor sensitivity. Rapamycin treatment of USP9X-depleted HeLa and renal cancer 786-O cells, along with the pharmacological inhibition of USP9X, confirmed that this protein plays a role in patients' sensitivity to mTOR inhibitors. USP9X was not found to exert a direct effect on mTORC1, but subsequent ubiquitylome analyses identified p62 as a direct USP9X target. Increased p62 ubiquitination and the augmented rapamycin effect upon bortezomib treatment, together with the results of p62 and LC3 immunofluorescence assays, suggested that dysregulated autophagy in USP9X-depleted cells can have a synergistic effect with mTOR inhibitors. In summary, we show that USP9X constitutes a potential novel marker of sensitivity to mTOR inhibitors in chRCC patients, and represents a clinical strategy for increasing the sensitivity to these drugs.

Genomic and immune landscape Of metastatic pheochromocytoma and paraganglioma.

The mechanisms triggering metastasis in pheochromocytoma/paraganglioma are unknown, hindering therapeutic options for patients with metastatic tumors (mPPGL). Herein we show by genomic profiling of a large cohort of mPPGLs that high mutational load, microsatellite instability and somatic copy-number alteration burden are associated with ATRX/TERT alterations and are suitable prognostic markers. Transcriptomic analysis defines the signaling networks involved in the acquisition of metastatic competence and establishes a gene signature related to mPPGLs, highlighting CDK1 as an additional mPPGL marker. Immunogenomics accompanied by immunohistochemistry identifies a heterogeneous ecosystem at the tumor microenvironment level, linked to the genomic subtype and tumor behavior. Specifically, we define a general immunosuppressive microenvironment in mPPGLs, the exception being PD-L1 expressing MAML3-related tumors. Our study reveals canonical markers for risk of metastasis, and suggests the usefulness of including immune parameters in clinical management for PPGL prognostication and identification of patients who might benefit from immunotherapy.

Genetic bases of pheochromocytoma and paraganglioma.

The genetics of pheochromocytoma and paraganglioma (PPGL) has become increasingly complex over the last two decades. The list of genes involved in the development of these tumors has grown steadily, and there are currently more than 20 driver genes implicated in either the hereditary or the sporadic nature of the disease. Although genetic diagnosis is achieved in about 75-80% of patients, genetic etiology remains unexplained in a significant percentage of cases. Patients lacking a genetic diagnosis include not only those with apparently sporadic PPGL but also patients with a family history of the disease or with multiple tumors, that meet the criteria to be considered as candidates for carrying germline mutations in yet undiscovered genes. Mutations in known PPGL genes deregulate three main signaling pathways (hypoxia, kinase signaling, and Wnt-signaling pathways), which could be the starting point for the development of personalized treatment for PPGL patients. Furthermore, the integration of results from several genomic high-throughput platforms enables the discovery of regulatory mechanisms that cannot be identified by analyzing each piece of information separately. These strategies are powerful tools for elucidating optimal therapeutic options based on molecular biomarkers in PPGL and represent an important step toward the achievement of precision medicine for patients with metastatic PPGL.

Comprehensive molecular analysis of immortalization hallmarks in thyroid cancer reveals new prognostic markers.

BACKGROUND: Comprehensive molecular studies on tumours are needed to delineate immortalization process steps and identify sensitive prognostic biomarkers in thyroid cancer. METHODS AND RESULTS: In this study, we extensively characterize telomere-related alterations in a series of 106 thyroid tumours with heterogeneous clinical outcomes. Using a custom-designed RNA-seq panel, we identified five telomerase holoenzyme-complex genes upregulated in clinically aggressive tumours compared to tumours from long-term disease-free patients, being TERT and TERC denoted as independent prognostic markers by multivariate regression model analysis. Characterization of alterations related to TERT re-expression revealed that promoter mutations, methylation and/or copy gains exclusively co-occurred in clinically aggressive tumours. Quantitative-FISH (fluorescence in situ hybridization) analysis of telomere lengths showed a significant shortening in these carcinomas, which matched with a high proliferative rate measured by Ki-67 immunohistochemistry. RNA-seq data analysis indicated that short-telomere tumours exhibit an increased transcriptional activity in the 5-Mb-subtelomeric regions, site of several telomerase-complex genes. Gene upregulation enrichment was significant for specific chromosome-ends such as the 5p, where TERT is located. Co-FISH analysis of 5p-end and TERT loci showed a more relaxed chromatin configuration in short telomere-length tumours compared to normal telomere-length tumours. CONCLUSIONS: Overall, our findings support that telomere shortening leads to a 5p subtelomeric region reorganization, facilitating the transcription and accumulation of alterations at TERT-locus.

International initiative for a curated SDHB variant database improving the diagnosis of hereditary paraganglioma and pheochromocytoma.

BACKGROUND: SDHB is one of the major genes predisposing to paraganglioma/pheochromocytoma (PPGL). Identifying pathogenic SDHB variants in patients with PPGL is essential to the management of patients and relatives due to the increased risk of recurrences, metastases and the emergence of non-PPGL tumours. In this context, the 'NGS and PPGL (NGSnPPGL) Study Group' initiated an international effort to collect, annotate and classify SDHB variants and to provide an accurate, expert-curated and freely available SDHB variant database. METHODS: A total of 223 distinct SDHB variants from 737 patients were collected worldwide. Using multiple criteria, each variant was first classified according to a 5-tier grouping based on American College of Medical Genetics and NGSnPPGL standardised recommendations and was then manually reviewed by a panel of experts in the field. RESULTS: This multistep process resulted in 23 benign/likely benign, 149 pathogenic/likely pathogenic variants and 51 variants of unknown significance (VUS). Expert curation reduced by half the number of variants initially classified as VUS. Variant classifications are publicly accessible via the Leiden Open Variation Database system (https://databases.lovd.nl/shared/genes/SDHB). CONCLUSION: This international initiative by a panel of experts allowed us to establish a consensus classification for 223 SDHB variants that should be used as a routine tool by geneticists in charge of PPGL laboratory diagnosis. This accurate classification of SDHB genetic variants will help to clarify the diagnosis of hereditary PPGL and to improve the clinical care of patients and relatives with PPGL.

Co-occurrence of mutations in NF1 and other susceptibility genes in pheochromocytoma and paraganglioma.

Introduction: The percentage of patients diagnosed with pheochromocytoma and paraganglioma (altogether PPGL) carrying known germline mutations in one of the over fifteen susceptibility genes identified to date has dramatically increased during the last two decades, accounting for up to 35-40% of PPGL patients. Moreover, the application of NGS to the diagnosis of PPGL detects unexpected co-occurrences of pathogenic allelic variants in different susceptibility genes. Methods: Herein we uncover several cases with dual mutations in NF1 and other PPGL genes by targeted sequencing. We studied the molecular characteristics of the tumours with co-occurrent mutations, using omic tools to gain insight into the role of these events in tumour development. Results: Amongst 23 patients carrying germline NF1 mutations, targeted sequencing revealed additional pathogenic germline variants in DLST (n=1) and MDH2 (n=2), and two somatic mutations in H3-3A and PRKAR1A. Three additional patients, with somatic mutations in NF1 were found carrying germline pathogenic mutations in SDHB or DLST, and a somatic truncating mutation in ATRX. Two of the cases with dual germline mutations showed multiple pheochromocytomas or extra-adrenal paragangliomas - an extremely rare clinical finding in NF1 patients. Transcriptional and methylation profiling and metabolite assessment showed an "intermediate signature" to suggest that both variants had a pathological role in tumour development. Discussion: In conclusion, mutations affecting genes involved in different pathways (pseudohypoxic and receptor tyrosine kinase signalling) co-occurring in the same patient could provide a selective advantage for the development of PPGL, and explain the variable expressivity and incomplete penetrance observed in some patients.

Analysis of Telomere Maintenance Related Genes Reveals NOP10 as a New Metastatic-Risk Marker in Pheochromocytoma/Paraganglioma.

One of the main problems we face with PPGL is the lack of molecular markers capable of predicting the development of metastases in patients. Telomere-related genes, such as TERT and ATRX, have been recently described in PPGL, supporting the association between the activation of immortalization mechanisms and disease progression. However, the contribution of other genes involving telomere preservation machinery has not been previously investigated. In this work, we aimed to analyze the prognostic value of a comprehensive set of genes involved in telomere maintenance. For this study, we collected 165 PPGL samples (97 non-metastatic/63 metastatic), genetically characterized, in which the expression of 29 genes of interest was studied by NGS. Three of the 29 genes studied, TERT, ATRX and NOP10, showed differential expression between metastatic and non-metastatic cases, and alterations in these genes were associated with a shorter time to progression, independent of SDHB-status. We studied telomere length by Q-FISH in patient samples and in an in vitro model. NOP10 overexpressing tumors displayed an intermediate-length telomere phenotype without ALT, and in vitro results suggest that NOP10 has a role in telomerase-dependent telomere maintenance. We also propose the implementation of NOP10 IHC to better stratify PPGL patients.

Prevalence of pathogenic germline variants in patients with metastatic renal cell carcinoma.

PURPOSE: Germline pathogenic variants are estimated to affect 3-5% of renal cell carcinoma (RCC) patients. However, higher mutational prevalence in non-clear cell RCC (non-ccRCC) and advanced disease has been suggested. METHODS: To clarify the prevalence of pathogenic germline variants in metastatic RCC, we sequenced 29 cancer susceptibility genes in 294 unselected metastatic RCC cases plus 21 patients with clinical hereditary features. In 145 tumors, genes frequently mutated in RCC were sequenced and methylation was assessed in selected cases. RESULTS: Germline variants in RCC predisposition genes (FH, VHL) were detected in 1.4% of the unselected metastatic patients, with higher frequency in non-ccRCC versus ccRCC (6.4% and 0.4%; P = 0.0025) and in younger patients (P = 0.036). Among the 315 studied patients, 14% of non-type 1 papillary cases (4 of 28), all metastatic <1 year after diagnosis, carried a FH germline variant with loss of heterozygosity and tumor genome hypermethylation. Variants in other cancer-associated genes (e.g., MUTYH, BRCA2, CHEK2) occurred in 5.1% of the unselected series, with unclear significance for RCC. CONCLUSION: Our findings confirm a high prevalence of pathogenic germline variants in RCC predisposition genes in metastatic non-ccRCC, and highlight that metastatic patients with papillary type 2 or unconventional histologies compatible with FH would benefit from genetic screening.

Novel DNMT3A Germline Variant in a Patient with Multiple Paragangliomas and Papillary Thyroid Carcinoma.

Over the past few years, next generation technologies have been applied to unravel the genetics of rare inherited diseases, facilitating the discovery of new susceptibility genes. We recently found germline DNMT3A gain-of-function variants in two patients with head and neck paragangliomas causing a characteristic hypermethylated DNA profile. Here, whole-exome sequencing identifies a novel germline DNMT3A variant (p.Gly332Arg) in a patient with bilateral carotid paragangliomas, papillary thyroid carcinoma and idiopathic intellectual disability. The variant, located in the Pro-Trp-Trp-Pro (PWWP) domain of the protein involved in chromatin targeting, affects a residue mutated in papillary thyroid tumors and located between the two residues found mutated in microcephalic dwarfism patients. Structural modelling of the variant in the DNMT3A PWWP domain predicts that the interaction with H3K36me3 will be altered. An increased methylation of DNMT3A target genes, compatible with a gain-of-function effect of the alteration, was observed in saliva DNA from the proband and in one independent acute myeloid leukemia sample carrying the same p.Gly332Arg variant. Although further studies are needed to support a causal role of DNMT3A variants in paraganglioma, the description of a new DNMT3A alteration in a patient with multiple clinical features suggests a heterogeneous phenotypic spectrum related to DNMT3A germline variants.

CD133 Expression in Medullary Thyroid Cancer Cells Identifies Patients with Poor Prognosis.

CONTEXT: The identification of markers able to determine medullary thyroid cancer (MTC) patients at high-risk of disease progression is critical to improve their clinical management and outcome. Previous studies have suggested that expression of the stem cell marker CD133 is associated with MTC aggressiveness. OBJECTIVE: To evaluate CD133 impact on disease progression in MTC and explore the regulatory mechanisms leading to the upregulation of this protein in aggressive tumors. PATIENTS: We compiled a series of 74 MTCs with associated clinical data and characterized them for mutations in RET and RAS proto-oncogenes, presumed to be related with disease clinical behavior. RESULTS: We found that CD133 immunohistochemical expression was associated with adverse clinicopathological features and predicted a reduction in time to disease progression even when only RET-mutated cases were considered in the analysis (log-rank test P < 0.003). Univariate analysis for progression-free survival revealed CD133 expression and presence of tumor emboli in peritumoral blood vessels as the most significant prognostic covariates among others such as age, gender, and prognostic stage. Multivariate analysis identified both variables as independent factors of poor prognosis (hazard ratio = 16.6 and 2; P = 0.001 and 0.010, respectively). Finally, we defined hsa-miR-30a-5p, a miRNA downregulated in aggressive MTCs, as a CD133 expression regulator. Ectopic expression of hsa-miR-30a-5p in MZ-CRC-1 (RETM918T) cells significantly reduced CD133 mRNA expression. CONCLUSIONS: Our results suggest that CD133 expression may be a useful tool to identify MTC patients with poor prognosis, who may benefit from a more extensive primary surgical management and follow-up.

Molecular characterization of chromophobe renal cell carcinoma reveals mTOR pathway alterations in patients with poor outcome.

Chromophobe renal cell carcinoma (chRCC) is a histologically and molecularly distinct class of rare renal tumor. TCGA studies revealed low mutational burden, with only TP53 and PTEN recurrently mutated, and discovered alterations in TERT promoter and in the electron transport chain Complex I genes. However, knowledge on drug targetable genes is limited and treatments at metastatic stage do not follow a molecular rationale. In a large series of 92 chRCC enriched with metastatic cases, we performed an in-depth characterization of mTOR pathway alterations through targeted NGS and immunohistochemistry (IHC) of phospho-S6, tuberin, and PTEN. Mutations in mitochondria, telomere maintenance and other renal cancer related genes and p53 IHC, were also assessed. The impact on metastasis development and disease specific survival was determined, using TCGA-KICH series (n = 65) for validation. mTOR pathway mutations (MTOR, TSC1, TSC2) were present in 17% of primary tumors, most of them being classified as pathogenic. Mutations were associated with positive IHC staining of phospho-S6 and PTEN (P = 0.009 and P = 0.001, respectively) and with chRCC eosinophilic variant (P = 0.039), supporting a biological relevance of the pathway. mTOR pathway mutations were associated with worse clinical outcomes. Survival analysis gave a hazard ratio of 5.5 (P = 0.027), and this association was confirmed in TCGA-KICH (HR = 10.3, P = 0.006). TP53 mutations were enriched in metastatic cases (P = 0.018), and mutations in telomere maintenance genes showed a trend in the same direction. p53 IHC staining pattern was associated with the underlying TP53 defect, and negative PTEN IHC staining (82% of cases) suggested PTEN loss as a chRCC hallmark. In conclusion, our study provides with novel genomic knowledge in chRCC and identifies novel markers of poor survival. Furthermore, this is the first study showing that mTOR pathway mutations correlate with poor prognosis, and may help to identify patients with increased sensitivity to mTOR inhibitors.

A Novel Approach for the Identification of Pharmacogenetic Variants in MT-RNR1 through Next-Generation Sequencing Off-Target Data.

Specific genetic variants in the mitochondrially encoded 12S ribosomal RNA gene (MT-RNR1) cause aminoglycoside-induced irreversible hearing loss. Mitochondrial DNA is usually not included in targeted sequencing experiments; however, off-target data may deliver this information. Here, we extract MT-RNR1 genetic variation, including the most relevant ototoxicity variant m.1555A>G, using the off-target reads of 473 research samples, sequenced through a capture-based, custom-targeted panel and whole exome sequencing (WES), and of 1245 diagnostic samples with clinical WES. Sanger sequencing and fluorescence-based genotyping were used for genotype validation. There was a correlation between off-target reads and mitochondrial coverage (rcustomPanel = 0.39, p = 2 × 10-13 and rWES = 0.67, p = 7 × 10-21). The median read depth of MT-RNR1 m.1555 was similar to the average mitochondrial genome coverage, with saliva and blood samples giving comparable results. The genotypes from 415 samples, including three m.1555G carriers, were concordant with fluorescence-based genotyping data. In clinical WES, median MT-RNR1 coverage was 56×, with 90% of samples having ≥20 reads at m.1555 position, and one m.1494T and three m.1555G carriers were identified with no evidence for heteroplasmy. Altogether, this study shows that obtaining MT-RNR1 genotypes through off-target reads is an efficient strategy that can impulse preemptive pharmacogenetic screening of this mitochondrial gene.

PTEN expression and mutations in TSC1, TSC2 and MTOR are associated with response to rapalogs in patients with renal cell carcinoma.

The mammalian target of rapamycin (mTOR) pathway inhibitors are key drugs for the treatment of many tumor types, however, there are no predictive biomarkers in clinical use. Here, we performed a molecular and immunohistochemical characterization of key mTOR pathway components in a series of 105 renal cell carcinoma patients treated with rapalogs, aimed at identifying markers of treatment response. Mutational analysis in MTOR, TSC1 and TSC2 was performed through targeted next-generation sequencing (NGS), and immunohistochemistry (IHC) was performed for PTEN, pAKT, pS6K1, pS6 and p21. Among patients with NGS data, 11 of 87 (13%) had mTOR pathway mutations (8 in MTOR, 1 in TSC1 and 2 in TSC2). When comparing the molecular data to the response of the patients, we found that partial response was more frequent in cases with mTOR pathway mutations than in those without mutations (odds ratio [OR] = 0.08, 95% confidence interval [CI] = 0.008-0.79, p = 0.030 univariate; p = 0.038 multivariable). Regarding IHC, negative PTEN staining was detected in 58% of the tumors, and it was more frequent in rapalog responder patients (OR = 0.24, 95% CI = 0.065-0.86, p = 0.029 univariate; p = 0.029 multivariable). Mutations and PTEN IHC were not mutually exclusive events and its combination improved response prediction (OR = 0.16, 95% CI = 0.04-0.62, p = 0.008 univariate; p = 0.013 multivariable). The staining of other proteins did not show and association with response and no association with PFS was observed in unselected patients. In conclusion, our findings suggest that mTOR pathway mutations, negative PTEN IHC and their combination are potential markers of rapalog response.

Hsa-miR-139-5p is a prognostic thyroid cancer marker involved in HNRNPF-mediated alternative splicing.

It is critical to identify biomarkers and functional networks associated with aggressive thyroid cancer to anticipate disease progression and facilitate personalized patient management. We performed miRNome sequencing of 46 thyroid tumors enriched with advanced disease patients with a median follow-up of 96 months. MiRNome profiles correlated with tumor-specific histopathological and molecular features, such as stromal cell infiltration and tumor driver mutation. Differential expression analysis revealed a consistent hsa-miR-139-5p downexpression in primary carcinomas from patients with recurrent/metastatic disease compared to disease-free patients, sustained in paired local metastases and validated in publicly available thyroid cancer series. Exogenous expression of hsa-miR-139-5p significantly reduced migration and proliferation of anaplastic thyroid cancer cells. Proteomic analysis indicated RICTOR, SMAD2/3 and HNRNPF as putative hsa-miR-139-5p targets in our cell system. Abundance of HNRNPF mRNA, encoding an alternative splicing factor involved in cryptic exon inclusion/exclusion, inversely correlated with hsa-miR-139-5p expression in human tumors. RNA sequencing analysis revealed 174 splicing events differentially regulated upon HNRNPF repression in our cell system, affecting genes involved in RTK/RAS/MAPK and PI3K/AKT/MTOR signaling cascades among others. These results point at the hsa-miR-139-5p/HNRNPF axis as a novel regulatory mechanism associated with the modulation of major thyroid cancer signaling pathways and tumor virulence.

Integrative multi-omics analysis identifies a prognostic miRNA signature and a targetable miR-21-3p/TSC2/mTOR axis in metastatic pheochromocytoma/paraganglioma.

Rationale: Pheochromocytomas and paragangliomas (PPGLs) are rare neuroendocrine tumors that present variable outcomes. To date, no effective therapies or reliable prognostic markers are available for patients who develop metastatic PPGL (mPPGL). Our aim was to discover robust prognostic markers validated through in vitro models, and define specific therapeutic options according to tumor genomic features. Methods: We analyzed three PPGL miRNome datasets (n=443), validated candidate markers and assessed them in serum samples (n=36) to find a metastatic miRNA signature. An integrative study of miRNome, transcriptome and proteome was performed to find miRNA targets, which were further characterized in vitro. Results: A signature of six miRNAs (miR-21-3p, miR-183-5p, miR-182-5p, miR-96-5p, miR-551b-3p, and miR-202-5p) was associated with metastatic risk and time to progression. A higher expression of five of these miRNAs was also detected in PPGL patients' liquid biopsies compared with controls. The combined expression of miR-21-3p/miR-183-5p showed the best power to predict metastasis (AUC=0.804, P=4.67·10-18), and was found associated in vitro with pro-metastatic features, such as neuroendocrine-mesenchymal transition phenotype, and increased cell migration rate. A pan-cancer multi-omic integrative study correlated miR-21-3p levels with TSC2 expression, mTOR pathway activation, and a predictive signature for mTOR inhibitor-sensitivity in PPGLs and other cancers. Likewise, we demonstrated in vitro a TSC2 repression and an enhanced rapamycin sensitivity upon miR-21-3p expression. Conclusions: Our findings support the assessment of miR-21-3p/miR-183-5p, in tumors and liquid biopsies, as biomarkers for risk stratification to improve the PPGL patients' management. We propose miR-21-3p to select mPPGL patients who may benefit from mTOR inhibitors.