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.

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.

Cyclooxygenase 2 (COX-2) is a key enzyme of the tumorigenesis-inflammation interface and can be induced by hypoxia. A pseudohypoxic transcriptional signature characterizes pheochromocytomas and paragangliomas (PPGLs) of the cluster I, mainly represented by tumors with mutations in von Hippel-Lindau (VHL), endothelial PAS domain-containing protein 1 (EPAS1), or succinate dehydrogenase (SDH) subunit genes. The aim of this study was to investigate a possible association between underlying tumor driver mutations and COX-2 in PPGLs. COX-2 gene expression and immunoreactivity were examined in clinical specimens with documented mutations, as well as in spheroids and allografts derived from mouse pheochromocytoma (MPC) cells. COX-2 in vivo imaging was performed in allograft mice. We observed significantly higher COX-2 expression in cluster I, especially in VHL-mutant PPGLs, however, no specific association between COX-2 mRNA levels and a hypoxia-related transcriptional signature was found. COX-2 immunoreactivity was present in about 60% of clinical specimens as well as in MPC spheroids and allografts. A selective COX-2 tracer specifically accumulated in MPC allografts. This study demonstrates that, although pseudohypoxia is not the major determinant for high COX-2 levels in PPGLs, COX-2 is a relevant molecular target. This potentially allows for employing selective COX-2 inhibitors as targeted chemotherapeutic agents and radiosensitizers. Moreover, available models are suitable for preclinical testing of these treatments.

Pheochromocytomas (PCCs) and paragangliomas (PGLs) are neuroendocrine tumors that are mostly benign. Metastatic disease does occur in about 10% of cases of PCC and up to 25% of PGL, and for these patients no effective therapies are available. Patients with mutations in the succinate dehydrogenase subunit B (SDHB) gene tend to have metastatic disease. We hypothesized that a down-regulation in the active succinate dehydrogenase B subunit should result in notable changes in cellular metabolic profile and could present a vulnerability point for successful pharmacological targeting.

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.

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.

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.

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.

Abemaciclib is an ATP-competitive, reversible kinase inhibitor selective for CDK4 and CDK6 that has shown antitumor activity as a single agent in hormone receptor positive (HR+) metastatic breast cancer in clinical trials. Here, we examined the mechanistic effects of abemaciclib treatment using in vitro and in vivo breast cancer models. Treatment of estrogen receptor positive (ER+) breast cancer cells with abemaciclib alone led to a decrease in phosphorylation of Rb, arrest at G1, and a decrease in cell proliferation. Moreover, abemaciclib exposure led to durable inhibition of pRb, TopoIIα expression and DNA synthesis, which were maintained after drug removal. Treatment of ER+ breast cancer cells also led to a senescence response as indicated by accumulation of β-galactosidase, formation of senescence-associated heterochromatin foci, and a decrease in FOXM1 positive cells. Continuous exposure to abemaciclib altered breast cancer cell metabolism and induced apoptosis. In a xenograft model of ER+ breast cancer, abemaciclib monotherapy caused regression of tumor growth. Overall these data indicate that abemaciclib is a CDK4 and CDK6 inhibitor that, as a single agent, blocks breast cancer cell progression, and upon longer treatment can lead to sustained antitumor effects through the induction of senescence, apoptosis, and alteration of cellular metabolism.

Pheochromocytomas and paragangliomas (PPGLs) provide some of the clearest genetic evidence for the critical role of metabolism in the tumorigenesis process. Approximately 40% of PPGLs are caused by driver germline mutations in 16 known susceptibility genes, and approximately half of these genes encode members of the tricarboxylic acid (TCA) cycle. Taking as a starting point the involvement of the TCA cycle in PPGL development, we aimed to identify unreported mutations that occurred in genes involved in this key metabolic pathway and that could explain the phenotypes of additional individuals who lack mutations in known susceptibility genes. To accomplish this, we applied a targeted sequencing of 37 TCA-cycle-related genes to DNA from 104 PPGL-affected individuals with no mutations in the major known predisposing genes. We also performed omics-based analyses, TCA-related metabolite determination, and 13C5-glutamate labeling assays. We identified five germline variants affecting DLST in eight unrelated individuals (∼7%); all except one were diagnosed with multiple PPGLs. A recurrent variant, c.1121G>A (p.Gly374Glu), found in four of the eight individuals triggered accumulation of 2-hydroxyglutarate, both in tumors and in a heterologous cell-based assay designed to functionally evaluate DLST variants. p.Gly374Glu-DLST tumors exhibited loss of heterozygosity, and their methylation and expression profiles are similar to those of EPAS1-mutated PPGLs; this similarity suggests a link between DLST disruption and pseudohypoxia. Moreover, we found positive DLST immunostaining exclusively in tumors carrying TCA-cycle or EPAS1 mutations. In summary, this study reveals DLST as a PPGL-susceptibility gene and further strengthens the relevance of the TCA cycle in PPGL development.

Paragangliomas and pheochromocytomas (PPGLs) are chromaffin tumors associated with severe catecholamine-induced morbidities. Surgical removal is often curative. However, complete resection may not be an option for patients with succinate dehydrogenase subunit A-D (SDHx) mutations. SDHx mutations are associated with a high risk for multiple recurrent, and metastatic PPGLs. Treatment options in these cases are limited and prognosis is dismal once metastases are present. Identification of new therapeutic targets and candidate drugs is thus urgently needed. Previously, we showed elevated expression of succinate receptor 1 (SUCNR1) in SDHB PPGLs and SDHD head and neck paragangliomas. Its ligand succinate has been reported to accumulate due to SDHx mutations. We thus hypothesize that autocrine stimulation of SUCNR1 plays a role in the pathogenesis of SDHx mutation-derived PPGLs. We confirmed elevated SUCNR1 expression in SDHx PPGLs and after SDHB knockout in progenitor cells derived from a human pheochromocytoma (hPheo1). Succinate significantly increased viability of SUCNR1-transfected PC12 and ERK pathway signaling compared to control cells. Candidate SUCNR1 inhibitors successfully reversed proliferative effects of succinate. Our data reveal an unrecognized oncometabolic function of succinate in SDHx PPGLs, providing a growth advantage via SUCNR1.

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.