Regulation of epinephrine biosynthesis in HRAS-mutant paragangliomas.

The biochemical phenotype of paragangliomas (PGLs) is highly dependent on the underlying genetic background and tumor location. PGLs at extra-adrenal locations usually do not express phenylethanolamine N-methyltransferase (PNMT), the enzyme required for epinephrine production, which was explained by the absence of glucocorticoids. PGLs with pathogenic variants (PVs) in Harvey rat sarcoma viral oncogene homolog (HRAS) can occur in or outside of the adrenal, but always synthesize epinephrine independently of the localization. Here, we characterize the signaling pathways through which PVs in HRAS influence PNMT expression. Catecholamines, cortisol, and transcriptional features of PGL tissues with known genetic background were analyzed. Genetically modified rat pheochromocytoma cells carrying PVs in Hras were generated and analyzed for regulation of Pnmt expression. Elevated epinephrine contents in PGLs with PVs in HRAS were accompanied by enrichment in mitogen-activated protein kinase (MAPK) signaling compared to PGLs with PVs in genes that activate hypoxia pathways. In vitro, Hras PVs increased Pnmt expression and epinephrine biosynthesis through increased phosphorylation of stimulatory protein 1 via MAPK signaling. Here, we provide a molecular mechanism that explains the PV-dependent epinephrine production of PGLs.

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.

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.

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.

Succinate Mediates Tumorigenic Effects via Succinate Receptor 1: Potential for New Targeted Treatment Strategies in Succinate Dehydrogenase Deficient Paragangliomas.

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.

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.

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.

Targeting Cyclooxygenase-2 in Pheochromocytoma and Paraganglioma: Focus on Genetic Background.

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.

Using genome wide association studies to identify common QTL regions in three different genetic backgrounds based on Iberian pig breed.

One of the major limitation for the application of QTL results in pig breeding and QTN identification has been the limited number of QTL effects validated in different animal material. The aim of the current work was to validate QTL regions through joint and specific genome wide association and haplotype analyses for growth, fatness and premier cut weights in three different genetic backgrounds, backcrosses based on Iberian pigs, which has a major role in the analysis due to its high productive relevance. The results revealed nine common QTL regions, three segregating in all three backcrosses on SSC1, 0-3 Mb, for body weight, on SSC2, 3-9 Mb, for loin bone-in weight, and on SSC7, 3 Mb, for shoulder weight, and six segregating in two of the three backcrosses, on SSC2, SSC4, SSC6 and SSC10 for backfat thickness, shoulder and ham weights. Besides, 18 QTL regions were specifically identified in one of the three backcrosses, five identified only in BC_LD, seven in BC_DU and six in BC_PI. Beyond identifying and validating QTL, candidate genes and gene variants within the most interesting regions have been explored using functional annotation, gene expression data and SNP identification from RNA-Seq data. The results allowed us to propose a promising list of candidate mutations, those identified in PDE10A, DHCR7, MFN2 and CCNY genes located within the common QTL regions and those identified near ssc-mir-103-1 considered PANK3 regulators to be further analysed.

Deciphering the regulation of porcine genes influencing growth, fatness and yield-related traits through genetical genomics.

Genetical genomics approaches aim at identifying quantitative trait loci for molecular traits, also known as intermediate phenotypes, such as gene expression, that could link variation in genetic information to physiological traits. In the current study, an expression GWAS has been carried out on an experimental Iberian × Landrace backcross in order to identify the genomic regions regulating the gene expression of those genes whose expression is correlated with growth, fat deposition, and premium cut yield measures in pig. The analyses were conducted exploiting Porcine 60K SNP BeadChip genotypes and Porcine Expression Microarray data hybridized on mRNA from Longissimus dorsi muscle. In order to focus the analysis on productive traits and reduce the number of analyses, only those probesets whose expression showed significant correlation with at least one of the seven phenotypes of interest were selected for the eGWAS. A total of 63 eQTL regions were identified with effects on 36 different transcripts. Those eQTLs overlapping with phenotypic QTLs on SSC4, SSC9, SSC13, and SSC17 chromosomes previously detected in the same animal material were further analyzed. Moreover, candidate genes and SNPs were analyzed. Among the most promising results, a long non-coding RNA, ALDBSSCG0000001928, was identified, whose expression is correlated with premium cut yield. Association analysis and in silico sequence domain annotation support TXNRD3 polymorphisms as candidate to regulate ALDBSSCG0000001928 expression, which can be involved in the transcriptional regulation of surrounding genes, affecting productive and meat quality traits.