Abstract
Microarray data of non-functioning pituitary adenomas (NFPAs) were analyzed to disclose novel genes and pathways involved in NFPA tumorigenesis. Raw microarray data were downloaded from Gene Expression Omnibus. Data pre-treatment and differential analysis were conducted using packages in R. Functional and pathway enrichment analyses were performed using package GOs-tats. A protein-protein interaction (PPI) network was constructed using server STRING and Cytoscape. Known genes involved in pituitary adenomas (PAs), were obtained from the Comparative Toxicogenomics Database. A total of 604 differentially expressed genes (DEGs) were identifed between NFPAs and controls, including 177 up- and 427 down-regulated genes. Jak-STAT and p53 signaling pathways were significantly enriched by DEGs. The PPI network of DEGs was constructed, containing 99 up- and 288 down-regulated known disease genes (e.g. EGFR and ESR1) as well as 16 up- and 17 down-regulated potential novel NFPAs-related genes (e.g. COL4A5, LHX3, MSN, and GHSR). Genes like COL4A5, LHX3, MSN, and GHSR and pathways such as p53 signaling and Jak-STAT signaling, might participate in NFPA development. Although further validations are required, these findings might provide guidance for future basic and therapy researches.