Identification of the enhancer RNAs related to tumorgenesis of pituitary neuroendocrine tumors

Abstract

Background: Pituitary neuroendocrine tumors (PitNETs), which originate from the pituitary gland, account for 10%–15% of all intracranial neoplasms. Recent studies have indicated that enhancer RNAs (eRNAs) exert regulatory effects on tumor growth. However, the mechanisms underlying the eRNA-mediated tumorigenesis of PitNETs have not been elucidated. Methods: Normal pituitary and PitNETs tissues were used to identify the differentially expressed eRNAs (DEEs). Immune gene sets and hallmarks of cancer gene sets were quantified based on single sample gene set enrichment analysis (ssGSEA) algorithm using GSVA. The perspective of immune cells among all samples was calculated by the CIBERSORT algorithm. Moreover, the regulatory network composed of key DEEs, target genes of eRNAs, hallmarks of cancer gene sets, differentially expressed TF, immune cells and immune gene sets were constructed by Pearson correlation analysis. Small molecular anti-PitNETs drugs were explored by CMap analysis and the accuracy of the study was verified by in vitro and in vivo experiments, ATAC-seq and ChIP-seq. Results: In this study, data of 134 PitNETs and 107 non-tumorous pituitary samples were retrieved from a public database to identify differentially expressed genes. In total, 1128 differentially expressed eRNAs (DEEs) (494 upregulated eRNAs and 634 downregulated eRNAs) were identified. Next, the correlation of DEEs with cancer-related and immune-related gene signatures was examined to establish a co-expression regulatory network comprising 18 DEEs, 50 potential target genes of DEEs, 5 cancer hallmark gene sets, 2 differentially expressed transcription factors, 4 immune cell types, and 4 immune gene sets. Based on this network, the following four therapeutics for PitNETs were identified using Connectivity Map analysis: ciclopirox, bepridil, clomipramine, and alexidine. The growth-inhibitory effects of these therapeutics were validated using in vitro experiments. Ciclopirox exerted potential growth-inhibitory effects on PitNETs. Among the DEEs, GNLY, HOXB7, MRPL33, PRDM16, TCF7, and ZNF26 were determined to be potential diagnostic and therapeutic biomarkers for PitNETs. Conclusion: This study illustrated the significant influence of eRNAs on the occurrence and development of PitNETs. By constructing the co-expression regulation network, GNLY, HOXB6, MRPL33, PRDM16, TCF7, and ZNF26 were identified as relatively significant DEEs which were considered as the novel biomarkers of diagnosis and treatment of PitNETs. This study demonstrated the roles of eRNAs in the occurrence and development of PitNETs and revealed that ciclopirox was a potential therapeutic for pituitary adenomas.