Summary We report a comprehensive molecular characterization of pheochromocytomas and paragangliomas (PCC/PGLs), a rare tumor type. Multi-platform integration revealed that PCC/PGLs are driven by diverse alterations affecting multiple genes and pathways. Pathogenic germline mutations occurred in eight PCC/PGL susceptibility genes. We identified CSDE1 as a somatically-mutated driver gene, complementing four known drivers (HRAS, RET, EPAS1, NF1). We also discovered fusion genes in PCC/PGL, involving MAML3, BRAF, NGFR and NF1. Integrated analysis classified PCC/PGLs into four molecularly-defined groups: a kinase signaling subtype, a pseudohypoxia subtype, a Wnt-altered subtype, driven by MAML3 and CSDE1, and a cortical admixture subtype. Correlates of metastatic PCC/PGL included the MAML3 fusion gene. This integrated molecular characterization provides a comprehensive foundation for developing PCC/PGL precision medicine.

Our findings suggest that 11KT is the dominant bioactive androgen in children during adrenarche and PremA. Its androgenic capacity suggests that it may be responsible for the phenotypic changes seen in these phenomena.

Adrenocortical carcinoma (ACC) is a rare neoplasm with a heterogeneous outcome and limited treatment options. To understand its molecular and genomic landscape as a part of The Cancer Genome Atlas (TCGA) project, we performed the genomic, transcriptomic, epigenomic and proteomic profiling of 91 ACCs. We identified potential driving alterations including amplifications (TERT, TERF2 and CDK4), deletions (ZNRF3, CDKN2A and RB1) and point mutations in genes unknown to participate in adrenal disease (RPL22) and genes known to initiate familial syndromes that occasionally include adrenocortical neoplasms (TP53, CTNNB1, PRKAR1A, MEN1). The finding of PRKAR1A expands the catalogue of pathogenic pathways underlying ACC, suggesting of the protein kinase alpha signaling pathway as a potential target for molecular interventions. Novel transcript fusions potentially leading to overactive kinases included EXOSC10-MTOR and PPP1CB-BRE. DNA copy number analysis unveiled prevalent DNA losses leading to hypodiploidy as well as whole genome doubling (WGD) in 51% of ACC. The similar penetrance of loss of heterozygosity before and after WGD suggests a sequential development from hypodiploidy to polyploidy along the doubling in a subset of ACCs, which was endorsed by the worse outcome for WGD samples relative to nonWGD ACCs. An association between TERT expression and WGD was observed, suggesting a role for telomere regulation. These findings present ACC as a model disease for studies of WGD which is a frequent event in many tumor types. Unsupervised clustering of DNA methylation, copy number, gene expression, miRNA expression and protein abundance converged into three classes with specific biological characteristics and a respective median event free survival of 8, 38 and >100 months (p-value 1.7e-13). Comparison of the subtypes suggested additional drivers such as protein kinase C (PKC) phosphorylation and upregulation of a miRNA cluster at chromosome Xq27.3, which complemented the genomic alterations identified in these subtypes. To gain more insights into this rare cancer type, we placed ACC in a broader context of cancer genomic profiles including an array of other cancer types. These analyses revealed interesting shared features, including beta-catenin activation with a subset of endometroid cancer, DNA mismatch repair gene mutational signature with gastrointestinal cancers and a smoking signature with lung cancer. These findings highlight the commonalities between ACC and other lineages of cancer. Taken together, we found Wnt signaling pathway and p53/Rb signaling pathway were the most frequently altered pathways in ACC. Meanwhile, new players surfaced from our analyses including the PKA and PKC pathways. Our results present a comprehensive genomic landscape and refined molecular classification of ACC improve our understanding of its pathogenesis, and will ultimately improve the care of patients. Citation Format: Siyuan Zheng, Andrew D. Cherniack, Ninad Dewal, Richard A. Moffitt, Ludmila Danilova, Bradley A. Murray, Antonio M. Lerario, Tobias Else, Theo A. Knijnenburg, Giovanni Ciriello, Seungchan Kim, Guillaume Assie, Olena Morozova, Rehan Akbani, Juliann Shih, Katherine A. Hoadley, Toni K. Choueiri, Jens Waldmann, Ozgur Mete, Gordon A. Robertson, Matthew Meyerson, Michael J. Demeure, Felix Beuschlein, Anthony Gill, Ana C. Latronico, Maria C. Fragosa, Leslie Cope, Electron Kebebew, Mouhammed A. Habra, Timothy G. Whitsett, Kimberly J. Bussey, William E. Rainey, Sylvia Asa, Jérôme Bertherat, Martin Fassnacht, David A. Wheeler, The Cancer Genome Atlas Research Network, Gary D. Hammer, Thomas J. Giordano, Roel Verhaak. Comprehensive Pan-Genomic characterization of adrenocortical carcinoma. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 2976. doi:10.1158/1538-7445.AM2015-2976

Molecular diagnosis is rarely established in 46,XX testicular (T) disorder of sex development (DSD) individuals with atypical genitalia. The Wilms' tumour factor‐1 (WT1) gene is involved in early gonadal development in both sexes. Classically, WT1 deleterious variants are associated with 46,XY disorders of sex development (DSD) because of gonadal dysgenesis. We report a novel frameshift WT1 variant identified in an SRY‐negative 46,XX testicular DSD girl born with atypical genitalia. Target massively parallel sequencing involving DSD‐related genes identified a novel heterozygous WT1 c.1453_1456del; p.Arg485Glyfs*14 variant located in the fourth zinc finger of the protein which is absent in the population databases. Segregation analysis and microsatellite analysis confirmed the de novo status of the variant that is predicted to be deleterious by in silico tools and to increase WT1 target activation in crystallographic model. This novel and predicted activating frameshift WT1 variant leading to the 46,XX testicular DSD phenotype includes the fourth zinc‐finger DNA‐binding domain defects in the genetic aetiology of 46,XX DSD.

Lysyl oxidase-like 3 (LOXL3), belonging to the lysyl oxidase family, is responsible for the crosslinking in collagen or elastin. The cellular localization of LOXL3 is in the extracellular space by reason of its canonical function. In tumors, the presence of LOXL3 has been associated with genomic stability, cell proliferation, and metastasis. In silico analysis has shown that glioblastoma was among tumors with the highest LOXL3 expression levels. LOXL3 silencing of U87MG cells by siRNA led to the spreading of the tumor cell surface, and the transcriptome analysis of these cells revealed an upregulation of genes coding for extracellular matrix, cell adhesion, and cytoskeleton components, convergent to an increase in cell adhesion and a decrease in cell invasion observed in functional assays. Significant correlations of LOXL3 expression with genes coding for tubulins were observed in the mesenchymal subtype in the TCGA RNA-seq dataset of glioblastoma (GBM). Conversely, genes involved in endocytosis and lysosome formation, along with MAPK-binding proteins related to focal adhesion turnover, were downregulated, which may corroborate the observed decrease in cell viability and increase in the rate of cell death. Invasiveness is a major determinant of the recurrence and poor outcome of GBM patients, and downregulation of LOXL3 may contribute to halting the tumor cell invasion.

Background Microglia are the tissue-resident macrophages of the CNS. They originate in the yolk sac, colonize the CNS during embryonic development and form a self-sustaining population with limited turnover. A consequence of their relative slow turnover is that microglia can serve as a long-term memory for inflammatory or neurodegenerative events. Methods Using ATAC-, ChIP- and RNA-sequencing, we characterized the epigenomes and transcriptomes of FACS-purified microglia from mice exposed to different stimuli. A repeated endotoxin challenge (LPS) was used to induce tolerance in microglia, while genotoxic stress (DNA repair deficiency-induced accelerated aging through Ercc1 deficiency) resulted in primed (hypersensitive) microglia. Results Whereas the enrichment of permissive epigenetic marks at enhancer regions could explain training (hyper-responsiveness) of primed microglia to an LPS challenge, the tolerized response of microglia seems to be regulated by loss of permissive epigenetic marks. We identify that inflammatory stimuli and accelerated aging as a result of genotoxic stress activate distinct gene networks. These gene networks and associated biological processes are partially overlapping, which is likely driven by specific transcription factor networks, resulting in altered epigenetic signatures and distinct functional (desensitized vs. primed) microglia phenotypes. Conclusion This study provides insight into epigenetic profiles and transcription factor networks associated with transcriptional signatures of tolerized and trained microglia in vivo, leading to a better understanding of innate immune memory of microglia.

Familial primary aldosteronism (PA) is rare and mostly diagnosed in early-onset hypertension (HT). However, 'sporadic' bilateral adrenal hyperplasia (BAH) is the most frequent cause of PA and remains without genetic etiology in most cases. Our aim was to investigate new genetic defects associated with BAH and PA. We performed whole-exome sequencing (paired blood and adrenal tissue) in 6 patients with PA caused by BAH that underwent unilateral adrenalectomy. Additionally, we conducted functional studies in adrenal hyperplastic tissue and transfected cells to confirm the pathogenicity of the identified genetic variants. Rare germline variants in phosphodiesterase 2A (PDE2A) and 3B (PDE3B) genes were identified in 3 patients. The PDE2A heterozygous variant (p.Ile629Val) was identified in a patient with BAH and early-onset HT at 13 yrs of age. Two PDE3B heterozygous variants (p.Arg217Gln and p.Gly392Val) were identified in patients with BAH and HT diagnosed at 18 and 33 yrs of age, respectively. A strong PDE2A staining was found in all cases of BAH in zona glomerulosa and/or micronodules (that were also positive for CYP11B2). PKA activity in frozen tissue was significantly higher in BAH from patients harboring PDE2A and PDE3B variants. PDE2A and PDE3B variants significantly reduced protein expression in mutant transfected cells compared to WT. Interestingly, PDE2A and PDE3B variants increased SGK1 and SCNN1G/ENaCg at mRNA or protein levels. In conclusion, PDE2A and PDE3B variants were associated with PA caused by BAH. These novel genetic findings expand the spectrum of genetic etiologies of PA.

Lack of routine fresh or frozen tissue is a barrier to widespread transcriptomic analysis of adrenal cortical tumors and an impediment to translational research in endocrinology and endocrine oncology. Our group has previously pioneered the use of targeted amplicon-based next-generation sequencing for archival formalin-fixed paraffin-embedded (FFPE) adrenal tissue specimens to characterize the spectrum of somatic mutations in various forms of primary aldosteronism. Herein, we developed and validated a novel 194-amplicon targeted next-generation RNA sequencing (RNAseq) assay for transcriptomic analysis of adrenal tumors using clinical-grade FFPE specimens. Targeted RNAseq-derived expression values for 27 adrenal cortical tumors, including aldosterone-producing adenomas (APA; n=8), cortisol-producing adenomas (CPA; n=11), and adrenal cortical carcinomas (ACC; n=8), highlighted known differentially-expressed genes (DEGs; i. e., CYP11B2, IGF2, etc.) and tumor type-specific transcriptional modules (i. e., high cell cycle/proliferation transcript expression in ACC, etc.), and a subset of DEGs was validated orthogonally using quantitative reverse transcription PCR (qRT-PCR). Finally, unsupervised hierarchical clustering using a subset of high-confidence DEGs revealed three discrete clusters representing APA, CPA, and ACC tumors with corresponding unique gene expression signatures, suggesting potential clinical utility for a transcriptomic-based approach to tumor classification. Overall, these data support the use of targeted amplicon-based RNAseq for comprehensive transcriptomic profiling of archival FFPE adrenal tumor material and indicate that this approach may facilitate important translational research opportunities for the study of these tumors.