DNA Methylation Pattern in Somatotroph Pituitary Neuroendocrine Tumors.

INTRODUCTION: Growth hormone secretion by sporadic somatotroph neuroendocrine pituitary tumors (PitNETs) is a major cause of acromegaly. These tumors are relatively heterogenous in terms of histopathological and molecular features. Our previous transcriptomic profiling of somatotroph tumors revealed three distinct molecular subtypes. This study aimed to investigate the difference in DNA methylation patterns in subtypes of somatotroph PitNETs and its role in distinctive gene expression. METHODS: Genome-wide DNA methylation was investigated in 48 somatotroph PitNETs with EPIC microarrays. Gene expression was assessed with RNAseq. Bisulfite pyrosequencing and qRT-PCR were used for verifying the results of DNA methylation and gene expression. RESULTS: Clustering tumor samples based on methylation data reflected the transcriptome-related classification. Subtype 1 tumors are densely granulated without GNAS mutation, characterized by high expression of NR5A1 (SF-1) and GIPR. The expression of both genes is correlated with specific methylation of the gene body and promoter. This subtype has a lower methylation level of 5' gene regions and CpG islands than the remaining tumors. Subtype 2 PitNETs are densely granulated and frequently GNAS-mutated, while those in subtype 3 are mainly sparsely granulated. Methylation/expression analysis indicates that ∼50% genes located in differentially methylated regions are those differentially expressed between tumor subtypes. Correlation analysis revealed DNA methylation-controlled genes, including CDKN1B, CCND2, EBF3, CDH4, CDH12, MGMT, STAT5A, PLXND1, PTPRE, and MMP16, and genes encoding ion channels and semaphorins. CONCLUSION: DNA methylation profiling confirmed the existence of three molecular subtypes of somatotroph PitNETs. High expression of NR5A1 and GIPR in subtype 1 tumors is correlated with specific methylation of both genes.

Fast clustering and cell-type annotation of scATAC data using pre-trained embeddings.

Data from the single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) are now widely available. One major computational challenge is dealing with high dimensionality and inherent sparsity, which is typically addressed by producing lower dimensional representations of single cells for downstream clustering tasks. Current approaches produce such individual cell embeddings directly through a one-step learning process. Here, we propose an alternative approach by building embedding models pre-trained on reference data. We argue that this provides a more flexible analysis workflow that also has computational performance advantages through transfer learning. We implemented our approach in scEmbed, an unsupervised machine-learning framework that learns low-dimensional embeddings of genomic regulatory regions to represent and analyze scATAC-seq data. scEmbed performs well in terms of clustering ability and has the key advantage of learning patterns of region co-occurrence that can be transferred to other, unseen datasets. Moreover, models pre-trained on reference data can be exploited to build fast and accurate cell-type annotation systems without the need for other data modalities. scEmbed is implemented in Python and it is available to download from GitHub. We also make our pre-trained models available on huggingface for public use. scEmbed is open source and available at https://github.com/databio/geniml. Pre-trained models from this work can be obtained on huggingface: https://huggingface.co/databio.

Transcriptomic Classification of Pituitary Neuroendocrine Tumors Causing Acromegaly.

Acromegaly results from growth hormone hypersecretion, predominantly caused by a somatotroph pituitary neuroendocrine tumor (PitNET). Acromegaly-causing tumors are histologically diverse. Our aim was to determine transcriptomic profiles of various somatotroph PitNETs and to evaluate clinical implication of differential gene expression. A total of 48 tumors were subjected to RNA sequencing, while expression of selected genes was assessed in 134 tumors with qRT-PCR. Whole-transcriptome analysis revealed three transcriptomic groups of somatotroph PitNETs. They differ in expression of numerous genes including those involved in growth hormone secretion and known prognostic genes. Transcriptomic subgroups can be distinguished by determining the expression of marker genes. Analysis of the entire cohort of patients confirmed differences between molecular subtypes of tumors. Transcriptomic group 1 includes ~20% of acromegaly patients with GNAS mutations-negative, mainly densely granulated tumors that co-express GIPR and NR5A1 (SF-1). SF-1 expression was verified with immunohistochemistry. Transcriptomic group 2 tumors are the most common (46%) and include mainly GNAS-mutated, densely granulated somatotroph and mixed PitNETs. They have a smaller size and express favorable prognosis-related genes. Transcriptomic group 3 includes predominantly sparsely granulated somatotroph PitNETs with low GNAS mutations frequency causing ~35% of acromegaly. Ghrelin signaling is implicated in their pathogenesis. They have an unfavorable gene expression profile and higher invasive growth rate.