Supplementary Material for: DNA methylation pattern in somatotroph pituitary neuroendocrine tumors

Introduction. Growth hormone secretion by sporadic somatotroph neuroendocrine pituitary tumors (PitNETs) is 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 was aimed to investigate difference in DNA methylation pattern 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 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 gene body and promoter. This subtype has lower methylation level of 5’ gene regions and CpG islands than the remaining tumors. Subtype 2 PitNETs are densely granulated, 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, 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.