In order to identify genetic factors related to thyroid cancer susceptibility, we adopted a candidate gene approach. We studied tag- and putative functional SNPs in genes involved in thyroid cell differentiation and proliferation, and in genes found to be differentially expressed in thyroid carcinoma. A total of 768 SNPs in 97 genes were genotyped in a Spanish series of 615 cases and 525 controls, the former comprising the largest collection of patients with this pathology from a single population studied to date. SNPs in an LD block spanning the entire FOXE1 gene showed the strongest evidence of association with papillary thyroid carcinoma susceptibility. This association was validated in a second stage of the study that included an independent Italian series of 482 patients and 532 controls. The strongest association results were observed for rs1867277 (OR[per-allele] = 1.49; 95%CI = 1.30-1.70; P = 5.9x10(-9)). Functional assays of rs1867277 (NM_004473.3:c.-283G>A) within the FOXE1 5' UTR suggested that this variant affects FOXE1 transcription. DNA-binding assays demonstrated that, exclusively, the sequence containing the A allele recruited the USF1/USF2 transcription factors, while both alleles formed a complex in which DREAM/CREB/alphaCREM participated. Transfection studies showed an allele-dependent transcriptional regulation of FOXE1. We propose a FOXE1 regulation model dependent on the rs1867277 genotype, indicating that this SNP is a causal variant in thyroid cancer susceptibility. Our results constitute the first functional explanation for an association identified by a GWAS and thereby elucidate a mechanism of thyroid cancer susceptibility. They also attest to the efficacy of candidate gene approaches in the GWAS era.

The molecular genetics of inherited phaeochromocytoma have received considerable attention, but the somatic genetic and epigenetic events that characterise tumourigenesis in sporadic phaeochromocytomas are less well defined. Previously, we found considerable overlap between patterns of promoter region tumour suppressor gene (TSG) hypermethylation in two neural crest tumours, neuroblastoma and phaeochromocytoma. In order to identify candidate biomarkers and epigenetically inactivated TSGs in phaeochromocytoma and neuroblastoma, we characterised changes in gene expression in three neuroblastoma cell lines after treatment with the demethylating agent 5-azacytidine. Promoter region methylation status was then determined for 28 genes that demonstrated increased expression after demethylation. Three genes HSP47, homeobox A9 (HOXA9) and opioid binding protein (OPCML) were methylated in >10% of phaeochromocytomas (52, 17 and 12% respectively). Two of the genes, epithelial membrane protein 3 (EMP3) and HSP47, demonstrated significantly more frequent methylation in neuroblastoma than phaeochromocytoma. These findings extend epigenotype of phaeochromocytoma and identify candidate genes implicated in sporadic phaeochromocytoma tumourigenesis.

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.

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.

Hereditary leiomyomatosis and renal cell cancer syndrome (HLRCC) is a very rare hereditary disorder characterized by cutaneous leiomyomas (CLMs), uterine leiomyomas (ULMs), renal cysts (RCys) and renal cell cancers (RCCs). We aimed to describe the genetics, clinical features and potential genotype-phenotype associations in the largest cohort of fumarate hydratase enzyme mutation carriers known from Spain using a multicentre, retrospective study of individuals with a genetic or clinical diagnosis of HLRCC. We collected clinical information from medical records, analysed genetic variants and looked for genotype-phenotype associations. Analyses were performed using R 3.6.0. software. We included 197 individuals: 74 index cases and 123 relatives. CLMs were diagnosed in 65% of patients, ULMs in 90% of women, RCys in 37% and RCC in 10.9%. Twenty-seven different pathogenic variants were detected, 12 (44%) of them not reported previously. Patients with missense pathogenic variants showed higher frequencies of CLMs, ULMs and RCys, than those with loss-of-function variants (p = 0.0380, p = 0.0015 and p = 0.024, respectively). This is the first report of patients with HLRCC from Spain. The frequency of RCCs was lower than those reported in the previously published series. Individuals with missense pathogenic variants had higher frequencies of CLMs, ULMs and RCys.

Metastatic breast cancer (MBC) progressing after endocrine therapy frequently activates PI3K/AKT/mTOR pathway. The BOLERO-2 trial showed that everolimus-exemestane achieves increased progression free survival (PFS) compared with exemestane. However, there is great inter-patient variability in toxicity and response to exemestane-everolimus treatment. The objective of this study was to perform an exploratory study analyzing the implication of single nucleotide polymorphisms (SNPs) on outcomes from this treatment through a pharmacogenetic analysis.

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.

Pheochromocytomas/paragangliomas (PPGLs) with pathogenic mutations in the succinate dehydrogenase subunit B (SDHB) are associated with a high metastatic risk. Somatostatin receptor 2 (SSTR2)-dependent imaging is the most sensitive imaging modality for SDHB-related PPGLs, suggesting that SSTR2 expression is a significant cell surface therapeutic biomarker of such tumors.

Pheochromocytomas and paragangliomas (PPGLs) with activated pseudohypoxic pathways are associated with an immature catecholamine phenotype and carry a higher risk for metastasis. For improved understanding of the underlying mechanisms we investigated the impact of hypoxia and pseudohypoxia on catecholamine biosynthesis in pheochromocytoma cells naturally lacking Hif2α (MPC and MTT) or expressing both Hif1α and Hif2α (PC12). Cultivation under extrinsic hypoxia or in spheroid culture (intrinsic hypoxia) increased cellular dopamine and norepinephrine contents in all cell lines. To distinguish further between Hif1α- and Hif2α-driven effects we expressed Hif2α in MTT and MPC-mCherry cells (naturally lacking Hif2α). Presence of Hif2α resulted in similarly increased cellular dopamine and norepinephrine under hypoxia as in the control cells. Furthermore, hypoxia resulted in enhanced phosphorylation of tyrosine hydroxylase (TH). A specific knockdown of Hif1α in PC12 diminished these effects. Pseudohypoxic conditions, simulated by expression of Hif2α under normoxia resulted in increased TH phosphorylation, further stimulated by extrinsic hypoxia. Correlations with PPGL tissue data led us to conclude that catecholamine biosynthesis under hypoxia is mainly mediated through increased phosphorylation of TH, regulated as a short-term response (24-48 h) by HIF1α. Continuous activation of hypoxia-related genes under pseudohypoxia leads to a HIF2α-mediated phosphorylation of TH (permanent status).

Pheochromocytomas and paragangliomas (PPGL) are rare neuroendocrine tumors with a strong hereditary background and a large genetic heterogeneity. Identification of the underlying genetic cause is crucial for the management of patients and their families as it aids differentiation between hereditary and sporadic cases. To improve diagnostics and clinical management we tailored an enrichment based comprehensive multi-gene next generation sequencing panel applicable to both analyses of tumor tissue and blood samples. We applied this panel to tumor samples and compared its performance to our current routine diagnostic approach. Routine diagnostic sequencing of 11 PPGL susceptibility genes was applied to blood samples of 65 unselected PPGL patients at a single center in Dresden, Germany. Predisposing germline mutations were identified in 19 (29.2%) patients. Analyses of 28 PPGL tumor tissues using the dedicated PPGL panel revealed pathogenic or likely pathogenic variants in known PPGL susceptibility genes in 21 (75%) cases, including mutations in IDH2, ATRX and HRAS. These mutations suggest sporadic tumor development. Our results imply a diagnostic benefit from extended molecular tumor testing of PPGLs and consequent improvement of patient management. The approach is promising for determination of prognostic biomarkers that support therapeutic decision-making.

Papillary Thyroid Cancer (PTC) is a heterogeneous and complex disease; susceptibility to PTC is influenced by the joint effects of multiple common, low-penetrance genes, although relatively few have been identified to date. Here we applied a rigorous combined approach to assess both the individual and epistatic contributions of genetic factors to PTC susceptibility, based on one of the largest series of thyroid cancer cases described to date. In addition to identifying the involvement of TSHR variation in classic PTC, our pioneer study of epistasis revealed a significant interaction between variants in STK17B and PAX8. The interaction was detected by MD-MBR (p = 0.00010) and confirmed by other methods, and then replicated in a second independent series of patients (MD-MBR p = 0.017). Furthermore, we demonstrated an inverse correlation between expression of PAX8 and STK17B in a set of cell lines derived from human thyroid carcinomas. Overall, our work sheds additional light on the genetic basis of thyroid cancer susceptibility, and suggests a new direction for the exploration of the inherited genetic contribution to disease using association studies.

Hypomethylation of DNA is a hallmark of cancer and its analysis as tumor biomarker has been proposed, but its determination in clinical settings is hampered by lack of standardized methodologies. Here, we present QUAlu (Quantification of Unmethylated Alu), a new technique to estimate the Percentage of UnMethylated Alu (PUMA) as a surrogate for global hypomethylation. QUAlu consists in the measurement by qPCR of Alu repeats after digestion of genomic DNA with isoschizomers with differential sensitivity to DNA methylation. QUAlu performance has been evaluated for reproducibility, trueness and specificity, and validated by deep sequencing. As a proof of use, QUAlu has been applied to a broad variety of pathological examination specimens covering five cancer types. Major findings of the preliminary application of QUAlu to clinical samples include: (1) all normal tissues displayed similar PUMA; (2) tumors showed variable PUMA with the highest levels in lung and colon and the lowest in thyroid cancer; (3) stools from colon cancer patients presented higher PUMA than those from control individuals; (4) lung squamous cell carcinomas showed higher PUMA than lung adenocarcinomas, and an increasing hypomethylation trend associated with smoking habits. In conclusion, QUAlu is a simple and robust method to determine Alu hypomethylation in human biospecimens and may be easily implemented in research and clinical settings.

Using an animal model for the study of murine thymic lymphomas, frequent losses of heterozygosity (LOHs) on six chromosomal regions were reported. To determine the existence of LOH loci in human lymphomas, we screened 43 non-Hodgkin's lymphomas by using polymorphisms mapped in each of the orthologous human region. All LOH regions (1p32-p36, 9p21, 9p22-p23, 9q22-q34 and 10q23-q24) were observed in B-cell lymphomas at different frequencies. In addition, analysis of the tumor suppressor genes, p16(INK4a), p73 and PTEN located in 9p21, 1p36 and 10q23, respectively, revealed the participation of p16(INK4a) and p73 but not of PTEN. Importantly, two subregions of LOH were observed in 9q22-q34, one of them not previously described. Our results support the usefulness of murine models to study human lymphoid neoplasias and reveal novel regions on 9q22-q34 candidate for containing tumor suppressor genes involved in human lymphomas.

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.

Comprehensive molecular studies on tumours are needed to delineate immortalization process steps and identify sensitive prognostic biomarkers in thyroid cancer.

The limited knowledge about the PTEN hamartoma tumor syndrome (PHTS) makes its diagnosis a challenging task. We aimed to define the clinical and genetic characteristics of this syndrome in the Spanish population and to identify new genes potentially associated with the disease.

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.

Hypoxia activates pathways associated with tumor progression, metastatic spread, and alterations in the immune microenvironment leading to an immunosuppressive phenotype. In particular, the upregulation of PD-L1, a target for therapy with checkpoint inhibitors, is well-studied in several tumors. However, the relationship between hypoxia and PD-L1 regulation in pheochromocytomas and paragangliomas (PPGL), and especially in paragangliomas treated with embolization, is still largely unexplored. We investigated the expression of the hypoxia-marker HIF-2α and of PD-L1 in a PPGL-cohort with and without embolization as potential biomarkers that may predict the response to treatment with HIF-2α and checkpoint inhibitors. A total of 29 tumor samples from 25 patients who were operated at a single center were included and analyzed utilizing immunohistochemistry (IHC) for PD-L1 and HIF-2α. Embolization prior to surgery was performed in seven (24%) tumors. PD-L1 expression in tumor cells of head and neck paragangliomas (HNPGLs) receiving prior embolization (median PD-L1 positivity: 15%) was significantly higher as compared to PD-L1 expression in HNPGLs without prior embolization (median PD-L1 positivity: 0%) (p = 0.008). Consistently, significantly more HNPGLs with prior embolization were positive for HIF-2α (median nuclear HIF-2α positivity: 40%) as compared to HNPGLs without prior embolization (median nuclear HIF-2α positivity: 0%) (p = 0.016). Our results support the hypothesis that embolization with subsequent hypoxia leads to the upregulation of both PD-L1 and HIF-2α in HNPGLs, and could thus facilitate targeted treatment with HIF-2α and checkpoint inhibitors in the case of inoperable, locally advanced, or metastatic disease.

Medullary thyroid carcinoma accounts for 2% to 5% of thyroid malignancies, of which 75% are sporadic and the remaining 25% are hereditary and related to multiple endocrine neoplasia type 2 syndrome. Despite a genotype-phenotype correlation with specific germline RET mutations, knowledge of pathways specifically associated with each mutation and with non-RET-mutated sporadic MTC remains lacking. Gene expression patterns have provided a tool for identifying molecular events related to specific tumor types and to different clinical features that could help identify novel therapeutic targets. Using transcriptional profiling of 49 frozen MTC specimens classified as RET mutation, we identified PROM1, LOXL2, GFRA1, and DKK4 as related to RET(M918T) and GAL as related to RET(634) mutation. An independent series of 19 frozen and 23 formalin-fixed, paraffin-embedded (FFPE) MTCs was used for validation by RT-qPCR. Two tissue microarrays containing 69 MTCs were available for IHC assays. According to pathway enrichment analysis and gene ontology biological processes, genes associated with the MTC(M918T) group were involved mainly in proliferative, cell adhesion, and general malignant metastatic effects and with Wnt, Notch, NFκB, JAK/Stat, and MAPK signaling pathways. Assays based on silencing of PROM1 by siRNAs performed in the MZ-CRC-1 cell line, harboring RET(M918T), caused an increase in apoptotic nuclei, suggesting that PROM1 is necessary for survival of these cells. This is the first report of PROM1 overexpression among primary tumors.

Liquid biopsies from cancer patients have the potential to improve diagnosis and prognosis. The assessment of surrogate markers of tumor progression in circulating extracellular vesicles could be a powerful non-invasive approach in this setting. We have characterized extracellular vesicles purified from the lymphatic drainage also known as exudative seroma (ES) of stage III melanoma patients obtained after lymphadenectomy. Proteomic analysis showed that seroma-derived exosomes are enriched in proteins resembling melanoma progression. In addition, we found that the BRAFV600E mutation can be detected in ES-derived extracellular vesicles and its detection correlated with patients at risk of relapse.