Pheochromocytomas can usually be confirmed or excluded using currently available biochemical tests of catecholamine excess. Follow-up tests are, nevertheless, often required to distinguish false-positive from true-positive results. The glucagon stimulation test represents one such test; its diagnostic utility is, however, unclear.

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.

Adrenal steroid hormones, which regulate a plethora of physiological functions, are produced via tightly controlled pathways. Investigations of these pathways, based on experimental data, can be facilitated by computational modeling for calculations of metabolic rate alterations. We therefore used a model system, based on mass balance and mass reaction equations, to kinetically evaluate adrenal steroidogenesis in human adrenal cortex-derived NCI H295R cells. For this purpose a panel of 10 steroids was measured by liquid chromatographic-tandem mass spectrometry. Time-dependent changes in cell incubate concentrations of steroids - including cortisol, aldosterone, dehydroepiandrosterone and their precursors - were measured after incubation with angiotensin II, forskolin and abiraterone. Model parameters were estimated based on experimental data using weighted least square fitting. Time-dependent angiotensin II- and forskolin-induced changes were observed for incubate concentrations of precursor steroids with peaks that preceded maximal increases in aldosterone and cortisol. Inhibition of 17-alpha-hydroxylase/17,20-lyase with abiraterone resulted in increases in upstream precursor steroids and decreases in downstream products. Derived model parameters, including rate constants of enzymatic processes, appropriately quantified observed and expected changes in metabolic pathways at multiple conversion steps. Our data demonstrate limitations of single time point measurements and the importance of assessing pathway dynamics in studies of adrenal cortical cell line steroidogenesis. Our analysis provides a framework for evaluation of steroidogenesis in adrenal cortical cell culture systems and demonstrates that computational modeling-derived estimates of kinetic parameters are an effective tool for describing perturbations in associated metabolic pathways.

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.

Mutations in the AAAS gene coding for the nuclear pore complex protein ALADIN lead to the autosomal recessive disorder triple A syndrome. Triple A patients present with a characteristic phenotype including alacrima, achalasia and adrenal insufficiency. Patient fibroblasts show increased levels of oxidative stress, and several in vitro studies have demonstrated that the nucleoporin ALADIN is involved in both the cellular oxidative stress response and adrenal steroidogenesis. It is known that ALADIN knock-out mice lack a phenotype resembling human triple A syndrome. The objective of this study was to determine whether the application of chronic oxidative stress by ingestion of paraquat would generate a triple A-like phenotype in ALADIN null mice. Adult male mice were fed either a paraquat (0.25 g/kg diet) or control diet for 11 days. After application of chronic oxidative stress, ALADIN knock-out mice presented with an unexpected compensated glutathione metabolism, but lacked a phenotype resembling human triple A syndrome. We did not observe increased levels of oxidative stress and alterations in adrenal steroidogenesis in mice depleted for ALADIN. This study stresses the species-specific role of the nucleoporin ALADIN, which in mice involves a novel compensatory mechanism for regulating the cellular glutathione redox response.

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.

Mosaic or somatic EPAS1 mutations are associated with a range of phenotypes including pheochromocytoma and/or paraganglioma (PPGL), polycythemia and somatostatinoma. The pathogenic potential of germline EPAS1 variants however is not well understood. We report a number of germline EPAS1 variants occurring in patients with PPGL, including a novel variant c.739C>A (p.Arg247Ser); a previously described variant c.1121T>A (p.Phe374Tyr); several rare variants, c.581A>G (p.His194Arg), c.2353C>A (p.Pro785Thr) and c.2365A>G (p.Ile789Val); a common variant c.2296A>C (p.Thr766Pro). We performed detailed functional studies to understand their pathogenic role in PPGL. In transient transfection studies, EPAS1/HIF-2α p.Arg247Ser, p.Phe374Tyr and p.Pro785Thr were all stable in normoxia. In co-immunoprecipitation assays, only the novel variant p.Arg247Ser showed diminished interaction with pVHL. A direct interaction between HIF-2α Arg247 and pVHL was confirmed in structural models. Transactivation was assessed by means of a HRE-containing reporter gene in transiently transfected cells, and significantly higher reporter activity was only observed with EPAS1/HIF-2α p.Phe374Tyr and p.Pro785Thr. In conclusion, three germline EPAS1 variants (c.739C>A (p.Arg247Ser), c.1121T>A (p.Phe374Tyr) and c.2353C>A (p.Pro785Thr)) all have some functional features in common with somatic activating mutations. Our findings suggest that these three germline variants are hypermorphic alleles that may act as modifiers to the expression of PPGLs.

Endogenous steroid hormones, especially glucocorticoids and mineralocorticoids, derive from the adrenal cortex, and drastic or sustained changes in their circulatory levels affect multiple organ systems. Although hypoxia signaling in steroidogenesis has been suggested, knowledge on the true impact of the HIFs (Hypoxia-Inducible Factors) in the adrenocortical cells of vertebrates is scant. By creating a unique set of transgenic mouse lines, we reveal a prominent role for HIF1α in the synthesis of virtually all steroids in vivo. Specifically, mice deficient in HIF1α in adrenocortical cells displayed enhanced levels of enzymes responsible for steroidogenesis and a cognate increase in circulatory steroid levels. These changes resulted in cytokine alterations and changes in the profile of circulatory mature hematopoietic cells. Conversely, HIF1α overexpression resulted in the opposite phenotype of insufficient steroid production due to impaired transcription of necessary enzymes. Based on these results, we propose HIF1α to be a vital regulator of steroidogenesis as its modulation in adrenocortical cells dramatically impacts hormone synthesis with systemic consequences. In addition, these mice can have potential clinical significances as they may serve as essential tools to understand the pathophysiology of hormone modulations in a number of diseases associated with metabolic syndrome, auto-immunity or even cancer.

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.

We previously identified a novel syndrome in patients characterized by paraganglioma, somatostatinoma, and polycythemia. In these patients, polycythemia occurs long before any tumor develops, and tumor removal only partially corrects polycythemia, with recurrence occurring shortly after surgery. Genetic mosaicism of gain-of-function mutations of the EPAS1 gene (encoding HIF2α) located in the oxygen degradation domain (ODD), typically p.530-532, was shown as the etiology of this syndrome. The aim of the present investigation was to demonstrate that these mutations are necessary and sufficient for the development of the symptoms. We developed transgenic mice with a gain-of-function Epas1A529V mutation (corresponding to human EPAS1A530V), which demonstrated elevated levels of erythropoietin and polycythemia, a decreased urinary metanephrine-to-normetanephrine ratio, and increased expression of somatostatin in the ampullary region of duodenum. Further, inhibition of HIF2α with its specific inhibitor PT2385 significantly reduced erythropoietin levels in the mutant mice. However, polycythemia persisted after PT2385 treatment, suggesting an alternative erythropoietin-independent mechanism of polycythemia. These findings demonstrate the vital roles of EPAS1 mutations in the syndrome development and the great potential of the Epas1A529V animal model for further pathogenesis and therapeutics studies.

Catecholamines and adrenocortical steroids are important regulators of blood pressure. Bidirectional relationships between adrenal steroids and catecholamines have been established but whether this is relevant to patients with pheochromocytoma is unclear.

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.

The adrenal gland is important for many physiological and pathophysiological processes, but studies are often restricted by limited availability of sample material. Improved methods for sample preparation are needed to facilitate analyses of multiple classes of adrenal metabolites and macromolecules in a single sample. A procedure was developed for preparation of chromaffin cells, mouse adrenals, and human chromaffin tumors that allows for multi-omics analyses of different metabolites and preservation of native proteins. To evaluate the new procedure, aliquots of samples were also prepared using conventional procedures. Metabolites were analyzed by liquid-chromatography with mass spectrometry or electrochemical detection. Metabolite contents of chromaffin cells and tissues analyzed with the new procedure were similar or even higher than with conventional methods. Catecholamine contents were comparable between both procedures. The TCA cycle metabolites, cis-aconitate, isocitate, and α-ketoglutarate were detected at higher concentrations in cells, while in tumor tissue only isocitrate and potentially fumarate were measured at higher contents. In contrast, in a broad untargeted metabolomics approach, a methanol-based preparation procedure of adrenals led to a 1.3-fold higher number of detected metabolites. The established procedure also allows for simultaneous investigation of adrenal hormones and related enzyme activities as well as proteins within a single sample. This novel multi-omics approach not only minimizes the amount of sample required and overcomes problems associated with tissue heterogeneity, but also provides a more complete picture of adrenal function and intra-adrenal interactions than previously possible.

Renal cell carcinoma (RCC) is among the 10 most common cancer entities and can be categorised into distinct subtypes by differential expression of Krebs cycle genes. We investigated the predictive value of several targeted metabolites with regards to tumour stages and patient survival in an unselected cohort of 420 RCCs. Unsupervised hierarchical clustering of metabolite ratios identified two main clusters separated by α-ketoglutarate (α-KG) levels and sub-clusters with differential levels of the oncometabolite 2-hydroxyglutarate (2HG). Sub-clusters characterised by high 2HG were enriched in higher tumour stages, suggesting metabolite profiles might be suitable predictors of tumour stage or survival. Bootstrap forest models based on single metabolite signatures showed that lactate, 2HG, citrate, aspartate, asparagine, and glutamine better predicted the cancer-specific survival (CSS) of clear cell RCC patients, whereas succinate and α-ketoglutarate were better CSS predictors for papillary RCC patients. Additionally, this assay identifies rare cases of tumours with SDHx mutations, which are caused predominantly by germline mutations and which predispose to development of different neoplasms. Hence, analysis of selected metabolites should be further evaluated for potential utility in liquid biopsies, which can be obtained using less invasive methods and potentially facilitate disease monitoring for both patients and caregivers.

Hypertension is a major global health problem with high prevalence and complex associated health risks. Primary hypertension (PHT) is most common and the reasons behind primary hypertension are largely unknown. Endocrine hypertension (EHT) is another complex form of hypertension with an estimated prevalence varying from 3 to 20% depending on the population studied. It occurs due to underlying conditions associated with hormonal excess mainly related to adrenal tumours and sub-categorised: primary aldosteronism (PA), Cushing's syndrome (CS), pheochromocytoma or functional paraganglioma (PPGL). Endocrine hypertension is often misdiagnosed as primary hypertension, causing delays in treatment for the underlying condition, reduced quality of life, and costly antihypertensive treatment that is often ineffective. This study systematically used targeted metabolomics and high-throughput machine learning methods to predict the key biomarkers in classifying and distinguishing the various subtypes of endocrine and primary hypertension. The trained models successfully classified CS from PHT and EHT from PHT with 92% specificity on the test set. The most prominent targeted metabolites and metabolite ratios for hypertension identification for different disease comparisons were C18:1, C18:2, and Orn/Arg. Sex was identified as an important feature in CS vs. PHT classification.

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.