Abstract
Pituitary adenomas (PA) are one of the most common primary brain tumors and comprise 15% of brain neoplasms. Most PAs are histologically benign but can cause significant morbidity. The genetic profile of PAs is poorly understood. We used single-cell RNA sequencing using the 10X genomic platform to investigate cellular heterogeneity in twelve non-functioning pituitary adenoma samples from nine patients including site-specific (core vs edge) samples from three patients. Our analysis identified discrete clusters of cells associated with activation of specific functional pathways including lipid metabolism, angiogenic, and antigen presentation and processing pathways regardless of location within the tumor. MALT1, a lncRNA associated with increased proliferation and metastasis was ubiquitously expressed amongst these samples. Analysis of the core vs edge samples showed two specific clusters with activated invasion-promoting pathways including PI3k/AKT signaling, Wnt signaling (Wnt6 and FZD4), and epithelial-mesenchymal transition (TGFB1, SMAD1, ZEB1, and SNAI2) in the edge of the tumors. The activated Wnt signaling cascade drove a proinflammatory tumor microenvironment induced by the expression of IL-1, IL-17, and Toll-like receptors (TLR6 and TLR7/8) resulting in suppression of Tregs. Copy number variation analysis using the CONICS-CNV algorithm highlighted distinct chromosomal alterations within our samples that led to insight into clonal variations within each tumor with loss of chromosome 2 an early event in tumorigenesis and gain/loss of chromosome 19 as late events. Mapping the copy number variation analysis with the somatic variant analysis using the Vartrix algorithm identified novel driver mutations within these tumors. These findings help define the molecular fingerprint of pituitary adenomas and provide insights which could be utilized for better management of these tumors.
Validation of copy number variation analysis for next-generation sequencing diagnostics
