Jing Guo
Qiuyue Fang1
Chuzhong Li
Yazhuo Zhang- 1Beijing Neurosurgical Institute, Capital Medical University, Beijing, China
- 2Department of Neurosurgery, Beijing Tiantan Hospital Affiliated to Capital Medical University, Beijing, China
- 3Beijing Institute for Brain Disorders Brain Tumor Center, Beijing, China
- 4China National Clinical Research Center for Neurological Diseases, Beijing, China
Purpose: Non-functioning pituitary adenoma (NFPA) is a very common type of intracranial tumor, which can be locally invasive and can have a high recurrence rate. The tumor microenvironment (TME) shows a high correlation with tumor pathogenesis and prognosis. The current study aimed to identify microenvironment-related genes in NFPAs and assess their prognostic value.
Methods: 73 NFPA tumor samples were collected from Beijing Tiantan Hospital and transcriptional expression profiles were obtained through microarray analysis. The immune and stromal scores of each sample were calculated through the ESTIMATE algorithm, and the patients were divided into high and low immune/stromal score groups. Intersection differentially expressed genes (DEGs) were then obtained to construct a protein–protein interaction (PPI) network. Potential functions and pathways of intersection DEGs were then analyzed through Gene Ontology and the Kyoto Encyclopedia of Genes and Genomes. The prognostic value of these genes was evaluated. The quantitative real-time polymerase chain reaction in another set of NFPA samples was used to confirm the credibility of the bioinformatics analysis.
Results: The immune/stromal scores were significantly correlated with cavernous sinus (CS) invasion. The Kaplan–Meier curve indicated that the high immune score group was significantly related to poor recurrence-free survival. We identified 497 intersection DEGs based on the high vs. low immune/stromal score groups. Function enrichment analyses of 497 DEGs and hub genes from the PPI network showed that these genes are mainly involved in the immune/inflammatory response, T cell activation, and the phosphatidylinositol 3 kinase-protein kinase B signaling pathway. Among the intersection DEGs, 88 genes were further verified as significantly expressed between the CS invasive group and the non-invasive group, and five genes were highly associated with NFPA prognosis.
Conclusion: We screened out a series of critical genes associated with the TME in NFPAs. These genes may play a fundamental role in the development and prognosis of NFPA and may yield new therapeutic targets.
Introduction
Non-functioning pituitary adenomas (NFPAs) are benign pituitary neoplasms. They account for 14–54% of pituitary adenomas, which are the second most common primary intracranial tumor (Fontana and Gaillard, 2009; Fernandez et al., 2010; Ostrom et al., 2015; Day et al., 2016; Ntali and Wass, 2018). The delay in the diagnosis of NFPAs is related to the clinical manifestations of the lack of excessive hormone secretion. Apart from pituitary incidentalomas, NFPAs are usually detected when the tumor is large enough to cause pressure effects on surrounding structures, with symptoms including headache, visual field defects, decreased libido and hydrocephalus (Lindholm et al., 2006; Sam et al., 2015). Surgery is preferred when treatment is needed and can be combined with adjuvant radiotherapy (Wass and Karavitaki, 2009). Although NFPAs are benign neoplasms, most demonstrate invasive growth at the time of diagnosis, manifesting as local infiltration, and have an increased risk of post-operative recurrence (Landeiro et al., 2015). However, the current understanding of the molecular mechanism that causes NFPAs and the mechanism of their invasion is still limited. Therefore, exploring the pathogenesis of NFPAs is essential for early detection of the disease and early prevention of invasive NFPAs.
Recently, increasing numbers of studies have shown that the tumor microenvironment (TME) plays a vital role in tumor progression and treatment (Tahmasebi Birgani and Carloni, 2017; Wu and Dai, 2017; Ren et al., 2018; Bian et al., 2019). The TME is the environment in which a tumor exists, including its blood vessels, lymphatic vessels, extracellular matrix, stromal cells, immune/inflammatory cells, and secreted proteins (Balkwill et al., 2012). Among the components of the TME, immune cells (lymphocytes and macrophages) and stromal cells (such as fibroblasts) are the main non-tumor cells and they could have an impact on cancer progression and may represent a new therapeutic target (Liu et al., 2018; Xiong et al., 2018; Ocana et al., 2019). For example, the poor prognosis of renal cell carcinoma is related to the infiltration of baseline CD8 T-cells (Giraldo et al., 2015; Drake and Stein, 2018). Some studies have found that macrophages can promote the invasiveness of pituitary tumors (Marques et al., 2019; Sato et al., 2019; Yagnik et al., 2019), as well as prostate (Maolake et al., 2017), ovarian (Colvin, 2014) and breast cancer (Carron et al., 2017). Kashima et al. (2019) have shown that the accumulation of tumor-associated fibroblasts (TAFs) in esophageal squamous cell carcinoma can enhance lymph node metastasis and promote tumor progression. Moreover, Lv et al. (2018) found that the expression of α-smooth muscle actin and vascular endothelial growth factor in invasive pituitary adenoma TAFs was higher than that in non-invasive TAFs and normal fibroblasts, and TAFs in invasive pituitary adenoma promoted the proliferation of GH3 cells in vitro and tumor growth in vivo. The above studies suggest that a comprehensive understanding of the TME and its associated genes may be important in improving therapeutic strategies for NFPAs. Based on the specific gene expression signatures of immune and stromal cells, Yoshihara et al. (2013) proposed a new algorithm called Estimation of Stromal and Immune cells in Malignant Tumor tissues using Expression data (ESTIMATE), which infers the proportion of stromal and immune cells in tumor samples by calculating immune and stromal scores. Some researchers have used this method to identify diagnostic and prognostic markers and immune related therapeutic targets in their studies (Jia et al., 2018; Luo et al., 2019; Pan et al., 2019).
In the current study, we used the ESTIMATE algorithm to calculate the stromal score and immune score in the TME of NFPAs to determine prognosis and target genes for NFPA treatment. Based on a high and low immune/stromal score, we obtained a series of differentially expressed genes (DEGs), and their relationship with the invasiveness and prognosis of NFPAs was investigated.
Materials and Methods
Patients and Samples
We collected 73 specimens from patients diagnosed with NFPA who underwent surgery and did not receive pre-operative radiochemotherapy at Beijing Tiantan Hospital from October 2007 to July 2014. These patients lacked clinical and biochemical evidence of adenohypophysis hormone overproduction. Post-operative recurrence of NFPA was defined as an increase of more than 2 mm in the diameter of the largest tumor measured through magnetic resonance imaging (MRI) in any direction from the day of surgery to the end of follow-up. Based on immunohistochemistry results, the NFPAs were classified into gonadotroph adenomas (GAs), silent adenomas (SAs), and null cell adenomas (NCAs). Cavernous sinus (CS) invasion was defined using the Knosp grading scale (grade 3 or 4) based on pre-operative MRI results (Knosp et al., 1993). Based on MRI data, the NFPAs were classified into microadenomas (diameter <10 mm), macroadenomas (≥10 mm) and giant adenomas (≥40 mm). The clinical data for the adenomas are described in Table 1. In addition, 15 CS invasive and 11 non-invasive NFPA clinical samples from the same hospital were collected as an independent validation cohort to verify the expressive levels using quantitative real-time polymerase chain reaction (qRT-PCR) (Supplementary Table 1). Tumor specimens were stored in liquid nitrogen within 30 min after resection until RNA was extracted. This study was approved by the local Ethics Committee and informed consent was obtained from each subject.
Table 1. Summary of patient demographics and clinical characteristics.
Total RNA Extraction and RNA Microarrays
A phenol-free mirVanaTM miRNA isolation kit (Cat # AM1561, Ambion, Austin, TX, United States) was used to extract and purify total RNA for the generation of fluorescently labeled SBC human ceRNA array V1.0 cRNA targets (4 × 180 K). The labeled cRNA target was then hybridized with the slide. After hybridization, slides were scanned using an Agilent microarray scanner (Agilent Technologies, Santa Clara, CA, United States). After extracting the data using feature extraction software 10.7 (Agilent Technologies, Santa Clara, CA, United States), a quantile algorithm was used to normalize the raw data using the Limma software package for the R program.
Differentially Expressed Genes (DEGs) Based on Immune/Stromal Scores
We applied the ESTIMATE algorithm in the R program (version 3.5.3) to calculate the immune and stromal scores for each sample (Yoshihara et al., 2013). To test the correlations between prognoses and immune/stromal score, NFPAs were classified into a high-immune/stromal score group and a low-stromal/immune score group. The R package “maxstat” was used to identify the optimal cutoff for grouping patients. CIBERSORT (Newman et al., 2015) is a deconvolution algorithm1 used to calculate the abundance of 22 types of human infiltrated immune cells in each sample. DEGs were identified by using the R program to compare the high-immune/stromal score group and the low-stromal/immune score group, and further study was considered with a | fold change| >2 and a p-value < 0.05. Venn diagrams were used to depict the common DEGs in both groups.
Construction of a Protein–Protein Interaction (PPI) Network
The protein–protein Interaction (PPI) network of DEGs was constructed using an online STRING database2 (Version 11.0) (Szklarczyk et al., 2019), and the interaction with a combined score >0.9 was considered statistically significant. The visualization and analysis of PPIs were completed using Cytoscape (version 3.3.0) software (Shannon et al., 2003). The degree distribution was analyzed using the Cytoscape plug-in and Network Analyzer.
Functional Enrichment Analysis
All of the intersection DEGs were enrichment analyzed through Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) to predict their biological function, which was performed using the ClueGo (Bindea et al., 2009) Cytoscape plugin (version 3.3.0). GO terms analysis included biological process, cellular components (CCs), and molecular functions (MFs) terms and a KEGG pathway with an adjusted p-value < 0.05 was considered to be significant.
Evaluation of Intersection DEGs
Box-plots were used to depict the intersection gene expression levels between the CS invasive and non-invasive groups. The statistical significance of the difference between the groups was calculated by using a t-test. Kaplan–Meier plots and log-rank test p-values were used to clarify the correlation between each intersection gene and prognostic performance. A forest plot was used to visualize the independent risk factors of the intersection genes for recurrence-free survival (RFS) through Cox regression analysis.
Validation of Gene Expression Through qRT-PCR
The total RNA of the validated samples was extracted and purified as described above. Then, 1 μg of total RNA was reverse transcribed into complementary DNA (cDNA) using a High Capacity cDNA Reverse Transcription Kit (0049472, Thermo Fisher) in accordance with the manufacturer’s instructions. Power SYBRTM Green PCR Master Mix (4367659, Thermo Fisher) was used for qRT-PCR with a total reaction volume of 20 μL. The housekeeping gene glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was used as an internal control gene. All primers were synthesized by Sangon Biotech (Shanghai, China). The sequences of the primers are shown in Table 2. The level of mRNA was determined using QuantStudio 3 and 5 Systems (Applied Biosystems). For relative quantitation, expression levels were calculated using the 2–ΔΔCT method. All qRT-PCR analyses were performed in triplicate.
Table 2. Primers used for qRT-PCR.
Results
The Correlation of Immune/Stromal Scores With Clinical Characteristics in NFPA
The RNA-seq data consisted of 73 samples, including 34 male and 39 female subjects with the median age at diagnosis was 52 year (range 25–73 year). The median follow-up was 60 months (range 4–98 months). We obtained complete gene expression profiles through high throughput sequencing and identified 18827 mRNAs expressed in NFPA with an expression value > 0 in the 73 samples. The stromal scores ranged from −1648.8 to 1273.7, the immune scores ranged from −1429.1 to 2376.3, and the ESTIMATE scores ranged from −2866 to 3650, as calculated using the ESTIMATE algorithm. The distribution of the immune and stromal scores did not vary by histological subtype (p = 0.84/0.16) and tumor size (p = 0.68/0.75) but was significantly different from that of the CS invasive and CS non-invasive groups (p = 0.0057/0.0013) (Figures 1A,B). The degrees of immune cell infiltration of each sample were determined by the CIBERSORT algorithm (Supplementary Figure 1). We found that the infiltration degree of B cells naive, T cells CD4 naive, and Macrophages M1 were significantly different between CS invasion and CS non-invasion group (Figure 1C). The data were analyzed using the Kruskal–Wallis test and the Mann Whitney test.
Figure 1. Immune and stromal scores were not correlated with histological subtype, tumor size, but were correlated with cavernous sinus (CS) invasion. (A) Distribution of immune scores plotted against histological subtype (p = 0.84), CS invasive classification (p = 0.0057), tumor size (p = 0.68). (B) Distribution of stromal scores plotted against histological subtype (p = 0.16), CS invasive classification (p = 0.0013), and tumor size (p = 0.75). (C) CIBERSORT analysis of the proportions of 22 tumor infiltrating immune cell types between CS invasion and CS-non-invasion group (1000 permutations). *p < 0.05; **p < 0.01; ns, not significant.
To investigate the relationship between prognosis and the immune/stromal score, 73 NFPAs were classified into high-score groups and low-score groups (see methods) and analyzed through Kaplan–Meier survival analysis and the log-rank test. The results revealed that the low immune scores were significantly correlated with better RFS (p = 0.01, Figure 2A). The low stromal scores were consistently correlated with better RFS, although this correlation was not statistically significant (p = 0.17, Figure 2B).
Figure 2. Immune conditions are associated with NFPA recurrence-free survival (RFS) and identification of DEGs based on immune scores and stromal scores. (A) Low immune scores were significantly correlated with an improved RFS (p = 0.01). (B) Stromal scores were not significantly associated with RFS (p = 0.17). (C,D) Volcano plot of DEGs from the low vs. high immune score/stromal score groups. Genes with p < 0.05 are shown in red (fold change >2) and blue (fold change <–2). Black plots represent the remaining genes (those with no significant difference). (E,F) Venn diagrams showing the number of commonly upregulated or downregulated DEGs in the stromal and immune score groups.
DEGs for the Immune/Stromal Scores in NFPA
By comparing the gene expression profiles of the high and low immune score groups, we identified 879 immune-related DEGs (cut-off criteria: p < 0.05 and | fold change| >2), including 818 upregulated genes and 61 downregulated genes (Figure 2C). Similarly, based on the high- and low-stromal score group, we identified 694 upregulated and 58 downregulated stromal-related genes (Figure 2D). In addition, Venn diagram analysis revealed that there were 478 intersection upregulated genes and 19 intersection downregulated genes in both the immune and stromal groups (Figures 2E,F). These intersection DEGs were selected for subsequent analysis.
Functional Enrichment Analysis of Intersection DEGs
Based on the ClueGO gene annotation tool, we conducted GO analysis to further evaluate the biological functions of the 497 intersection DEGs (Figure 3). Three subontologies of these DEGs were analyzed (Figure 3A). In terms of biological processes (BP), these intersection DEGs were mainly enriched in terms of T cell activation, inflammatory response, cell migration, and cytokine production; for the CC, these DEGs were primarily clustered in the extracellular matrix, major histocompatibility complex (MHC) protein complex and adherens junction; with regard to MF, they were mainly associated with cytokine activity, cell adhesion molecule (CAM) binding, and chemokine receptor binding. Moreover, enrichment analysis of the KEGG pathway suggested that they are mainly related to the phosphatidylinositol 3 kinase-protein kinase B (PI3K-Akt) signaling pathway, CAMs, and cytokine-cytokine receptor interaction (Figure 3B).
Figure 3. Function analysis of DEGs by Gene Ontology (GO) (A) and Kyoto Encyclopedia of Genes and Genomes (KEGG) (B).
PPI Network Construction of Intersection DEGs
The PPI network was constructed based on 497 intersection genes to explore hub genes and develop a thorough picture of these genes at the systems level. It contains 220 nodes and 154 interactions (Figure 4A). Figure 4B shows the top 20 DEGs based on degree distribution in the PPI network. The topological analysis suggested that HLA-E (degree = 28) was the key gene in the PPI network, which could interact with most genes as it had the highest degree. Moreover, functional enrichment analysis of these top 20 DEGs by ClueGo analysis also showed that they were significantly associated with 63 GO terms and 22 KEGG pathways, such as regulation of lymphocyte migration, regulation of T cell migration, NF-kappa B signaling pathway and TNF signaling pathway (Figure 4C).
Figure 4. Topological features of DEGs PPI network. (A) The view of the PPI network. Red indicates upregulation and green indicates downregulation. (B) The key DEGs in the PPI network with the top 20 degree distributions. (C) Function analysis of the 20 DEGs in the PPI network.
The Expression of Stromal-Immune Score Related Intersection DEGs in NFPA
Among the 497 intersection genes (478 upregulated genes and 19 downregulated genes), there are 88 genes which are significantly differentially expressed in the CS invasive group and non-invasive group (p < 0.05, | Log2(fold change)| >1). Among these genes, there are 11 genes (POSTN, LUM, GPC3, CGNL1, ABCA8, TFPI2, THBS2, BCAT1, MFAP4, HTRA3, and ADH1A) with adjust p < 0.05 and | Log2(fold change)| >1.5 (Figures 5A–K).
Figure 5. Box plots showing the different expression level of 11 genes with adjust p < 0.05 and | Log2(fold change)| >1.5 in NFPA CS invasion status [(A–K), FC: fold change]. (L) Validation of differentially expressed genes in the NFPAs through qRT-PCR. Relative expression of GPC3, ABCA8, TFPI2, and MFAP4 between CS invasive NFPAs and non-invasive NFPAs. Data were analyzed using unpaired Student’s t-test. *p < 0.05; ***p < 0.001.
Because of the limited sample size, GPC3, ABCA8, TFPI2, and MFAP4 were randomly selected from the above genes to verify their expression levels in our clinical samples. We performed qRT-PCR to measure the expression level of these genes between CS invasive NFPAs and non-invasive NFPAs (Figure 5L). Consistent with the results of the microarray data analysis, GPC3, ABCA8, TFPI2, and MFAP4 were significantly downregulated in the invasive group (n = 15) compared with the non-invasive group (n = 11).
Survival Analysis of Intersection DEGs
The correlation between the intersection DEGs and RFS was evaluated through Kaplan–Meier survival curves using the 73 NFPAs. Among the 497 intersection DEGs, five were predictive of RFS (p < 0.05), as assessed through the log-rank test (Figures 6A–E). Among these five prognostic intersection DEGs, patients with upregulated ADGRG6, CD52, GPR183, and NNMT expression demonstrated shorter RFS time than those with downregulated expression, while upregulated TBX3 was associated with favorable outcomes. Moreover, a forest plot demonstrated the prognostic effects of these five DEGs through univariable Cox regression analysis (Figure 6F).
Figure 6. Kaplan–Meier survival curves showing the impact of the expression level of five genes in terms of NFPA recurrence-free survival (RFS). (A–E) Comparison of RFS in the high (red line) and low (blue line) gene expression groups. p < 0.05 was used to assess differences in the Log rank test. (F) Forest plot of hazard ratios for these five prognostic DEGs. Hazard ratios and corresponding 95% confidence intervals were assessed through univariable Cox analysis.
Discussion
Despite advances in endoscopic techniques, most patients with NFPA are at increased risk of post-operative recurrence and a poor therapeutic response because of its invasive growth in surrounding normal tissues at the time of diagnosis (Landeiro et al., 2015; Fernandez-Miranda et al., 2018). Therefore, it is important to explore the pathogenesis and invasiveness of NFPAs, to help us understand their recurrence and provide systematic treatments. The TME is implicated in the initiation and development of various tumors (Balkwill et al., 2012; Hui and Chen, 2015), including pituitary adenomas (Lu et al., 2015; Lv et al., 2018; Sato et al., 2019; Yagnik et al., 2019). Therefore, it is important to identify TME-related genes that are involved in NFPA recurrence and invasiveness.
First, the immune score and stromal score were obtained using the ESTIMATE algorithm. Our results indicate that the immune/stromal score was not significantly correlated with histological subtype and tumor size, but was significantly related to CS invasion. We also investigate the abundance of immune cells in NFPA using the CIBERSORT algorithm, which showed that B cells naive, T cells CD4 naive, and Macrophages M1 were significantly different between the CS invasion and CS non-invasion group. Previous studies have shown that TME-related factors are associated with invasive pituitary adenomas, such as the number of M2-macrophages, the number of inflammatory cells, and increased expression of PD-L1 (Qiu et al., 2013; Sato et al., 2019). In addition, Kaplan–Meier survival analysis showed that a lower immune/stromal score was associated with better RFS. Although the stromal score was not significant in this context, possibly because of the small sample size, it also suggested that the TME may have a deep impact on tumor progression. This is consistent with previous studies showing that tumor-associated macrophages, T regulatory cells and mast cells are associated with poor prognosis in many types of cancer (Bingle et al., 2002; Balkwill et al., 2012; Hui and Chen, 2015; Fridman et al., 2017; Kemeny et al., 2020).
We analyzed 497 intersection DEGs (478 upregulated and 19 downregulated) which were obtained from the high vs. low immune/stromal score groups, as the intersection genes are the most highly conserved. Among the 497 intersection genes analyzed using GO terms, we found that most of them were enrichened in the immune/inflammatory response (BP), cytokine activity (MF) and extracellular matrix (CC). KEGG pathway enrichment analysis showed that these intersection DEGs were mainly clustered in the PI3K-Akt signaling pathway, CAMs, and cytokine-cytokine receptor interaction. The PI3K-Akt signaling pathway is a signal transduction cascade involved in cell growth and metabolism (Leslie et al., 2001). Long et al. (2019) revealed that COL6A6 could block the PI3K-Akt-pathway to suppress the growth and metastasis of pituitary adenoma. Moreover, we were able to construct a PPI module based on these 497 intersection genes. These hub genes are related to the regulation of lymphocyte migration, the regulation of T cell migration, the NF-kappa B signaling pathway and the TNF signaling pathway, which indicated that these genes are related to immune/inflammation and stromal response. We, therefore, assessed them through function enrichment analysis.
We analyzed 497 interaction genes and identified 88 genes that were significantly differentially expressed between the invasion group and the non-invasion group. To further narrow the scope, we showed 11 genes with adjust p < 0.05 and |Log2(fold change)| >1.5. As the sample size is limited, GPC3, ABCA8, TFPI2, and MFAP4 were randomly selected to verify the expression level in CS invasive NFPAs and non-invasive NFPAs through qRT-PCR to confirm the microarray analysis data. We found that the qRT-PCR results matched well with the bioinformation analysis of NFPAs. Several identified genes in our results have previously been reported to play important roles in several types of cancer. Wang et al. (2019) found that the expression of the ABC transporter gene (ABCA8) was significantly downregulated in prostate cancer compared with non-cancerous prostate tissues (Demidenko et al., 2015). This is consistent with our qRT-PCR and microarray results showing that the expression level of ABCA8 was significantly decreased in the CS invasive group vs. the non-invasive group. Guo et al. (2020) found that the increased expression of ABCA8 is correlated with poor prognosis and is associated with immune infiltration in stomach adenocarcinoma. Glypican-3 (GPC3) is a cell surface proteoglycan anchored by glycosyl-phosphatidylinositol (Pilia et al., 1996). Some studies have developed various immune-related treatment strategies against GPC3 (Sayem et al., 2016), particularly in hepatocellular carcinoma (Sawada et al., 2016; Wu et al., 2017). Tissue factor pathway inhibitor 2 (TFPI2) is a proteolytic enzyme inhibitor with a Kunitz domain that is down-regulated in several types of cancer (Konduri et al., 2001; Rollin et al., 2005; Sato et al., 2005). Gaud et al. (2011) found that TFPI2 produced by non-small lung cancer cells may contribute to the synthesis of matrix metallopeptidase by fibroblasts in the microenvironment. Microfibril associated protein 4 (MFAP4) is frequently downregulated in most human cancers and its high mRNA levels are significantly associated with better outcomes in the early stage in several cancers (Yang et al., 2019). In addition, Niu et al. (2011) have shown that MFAP4 promotes inflammatory cell recruitment and assists immunological cancer surveillance to limit cancer cell survival in early stage cancers.
Finally, five TME-related genes were identified with a significant relationship to the RFS. Among them, the four upregulated genes ADGRG6, CD52, GPR183 and NNMT correlate with unfavorable NFPA outcomes. ADGRG6 is a member of the adhesion G protein-coupled receptor family and the depletion of ADGRG6 expression in urothelial bladder carcinoma cells compromised their ability to recruit endothelial cells and induce tube formation (Wu et al., 2019). Upregulated CD52 is correlated with poor survival in patients with lung adenocarcinoma (Li et al., 2018). CD52 is a glycosylphosphatidylinositol-anchored protein expressed on the surface of normal T and B lymphocytes and is also known as the CAMPATH-1 antigen (Hale et al., 1996). In addition, nicotinamide N-methyltransferase (NNMT) was found to be overexpressed in gastric carcinoma tissue compared with adjacent tissues and is related to poor prognosis (Chen et al., 2016). NNMT is a phase II metabolizing enzyme, and mainly catalyzes the methylation of nicotinamide and other pyridines into pyridinium ions (Xie et al., 2014). Moreover, Wang et al. (2019) found that overexpressed NNMT is correlated with poor survival and chemotherapy response in breast cancer patients who received chemotherapy.
As explained above, a series of genes that we identified through bioinformatics methods are related to TME, and these genes are also related to patient prognosis and treatment. Although the function and role of some of these genes have been reported in multiple types of tumor and several genes have been validated in another set of NFPAs, further verification in large and diverse samples is still needed. Moreover, the mechanisms and potential functions of these TME-related genes in NFPAs remain to be investigated.
In summary, TME-related genes might be correlated with tumor invasion and prognosis. Our study may provide a novel signature and therapeutic target for the immune/stromal component in NFPAs.
Data Availability Statement
The data presented in the study are deposited in the Gene Expression Omnibus (https://www.ncbi.nlm.nih.gov/geo/) with accession numbers: GSE169498, further inquiries can be directed to the corresponding author/s.
Ethics Statement
The studies involving human participants were reviewed and approved by Ethics Committees of Beijing Tiantan Hospital. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
WX, CL, and YZ worked on conception, designed the research, and approved the manuscript. QF and YL contributed to collecting and analyzing the clinical data of patients. JG was dedicated to data analysis, interpretation, and drafting. All authors read and approved the final manuscript.
Funding
This study was supported by the National Natural Science Foundation of China (81771489) and the Beijing Municipal Science and Technology Commission (Z171100000117002).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.627117/full#supplementary-material
Supplementary Figure 1 | Stacked bar chart of the fraction of immune cells of the tumor microenvironment in 73 NFPA patients.
Supplementary Table 1 | Summary of validated patients demographics and characteristics.
Footnotes
References
Balkwill, F. R., Capasso, M., and Hagemann, T. (2012). The tumor microenvironment at a glance. J. Cell Sci. 125, 5591–5596. doi: 10.1242/jcs.116392
Bian, X., Xiao, Y. T., Wu, T., Yao, M., Du, L., Ren, S., et al. (2019). Microvesicles and chemokines in tumor microenvironment: mediators of intercellular communications in tumor progression. Mol. Cancer 18:50. doi: 10.1186/s12943-019-0973-7
Bindea, G., Mlecnik, B., Hackl, H., Charoentong, P., Tosolini, M., Kirilovsky, A., et al. (2009). ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics 25, 1091–1093. doi: 10.1093/bioinformatics/btp101
Bingle, L., Brown, N. J., and Lewis, C. E. (2002). The role of tumour-associated macrophages in tumour progression: implications for new anticancer therapies. J. Pathol. 196, 254–265. doi: 10.1002/path.1027
Carron, E. C., Homra, S., Rosenberg, J., Coffelt, S. B., Kittrell, F., Zhang, Y., et al. (2017). Macrophages promote the progression of premalignant mammary lesions to invasive cancer. Oncotarget 8, 50731–50746. doi: 10.18632/oncotarget.14913
Chen, C., Wang, X., Huang, X., Yong, H., Shen, J., Tang, Q., et al. (2016). Nicotinamide N-methyltransferase: a potential biomarker for worse prognosis in gastric carcinoma. Am. J. Cancer Res. 6, 649–663.
Colvin, E. K. (2014). Tumor-associated macrophages contribute to tumor progression in ovarian cancer. Front. Oncol. 4:137. doi: 10.3389/fonc.2014.00137
Day, P. F., Loto, M. G., Glerean, M., Picasso, M. F., Lovazzano, S., and Giunta, D. H. (2016). Incidence and prevalence of clinically relevant pituitary adenomas: retrospective cohort study in a Health Management Organization in Buenos Aires, Argentina. Arch. Endocrinol. Metab. 60, 554–561. doi: 10.1590/2359-3997000000195
Demidenko, R., Razanauskas, D., Daniunaite, K., Lazutka, J. R., Jankevicius, F., and Jarmalaite, S. (2015). Frequent down-regulation of ABC transporter genes in prostate cancer. BMC Cancer 15:683. doi: 10.1186/s12885-015-1689-8
Drake, C. G., and Stein, M. N. (2018). The immunobiology of kidney cancer. J. Clin. Oncol. 36, 3547–3552. doi: 10.1200/JCO.2018.79.2648
Fernandez, A., Karavitaki, N., and Wass, J. A. (2010). Prevalence of pituitary adenomas: a community-based, cross-sectional study in Banbury (Oxfordshire, UK). Clin. Endocrinol. (Oxf) 72, 377–382. doi: 10.1111/j.1365-2265.2009.03667.x
Fernandez-Miranda, J. C., Zwagerman, N. T., Abhinav, K., Lieber, S., Wang, E. W., Snyderman, C. H., et al. (2018). Cavernous sinus compartments from the endoscopic endonasal approach: anatomical considerations and surgical relevance to adenoma surgery. J. Neurosurg. 129, 430–441. doi: 10.3171/2017.2.JNS162214
Fontana, E., and Gaillard, R. (2009). [Epidemiology of pituitary adenoma: results of the first Swiss study]. Rev. Med. Suisse 5, 2172–2174.
Fridman, W. H., Zitvogel, L., Sautes-Fridman, C., and Kroemer, G. (2017). The immune contexture in cancer prognosis and treatment. Nat. Rev. Clin. Oncol. 14, 717–734. doi: 10.1038/nrclinonc.2017.101
Gaud, G., Iochmann, S., Guillon-Munos, A., Brillet, B., Petiot, S., Seigneuret, F., et al. (2011). TFPI-2 silencing increases tumour progression and promotes metalloproteinase 1 and 3 induction through tumour-stromal cell interactions. J. Cell. Mol. Med. 15, 196–208. doi: 10.1111/j.1582-4934.2009.00989.x
Giraldo, N. A., Becht, E., Pages, F., Skliris, G., Verkarre, V., Vano, Y., et al. (2015). Orchestration and prognostic significance of immune checkpoints in the microenvironment of primary and metastatic renal cell cancer. Clin. Cancer Res. 21, 3031–3040. doi: 10.1158/1078-0432.CCR-14-2926
Guo, Y., Wang, Z. W., Su, W. H., Chen, J., and Wang, Y. L. (2020). Prognostic value and immune infiltrates of ABCA8 and FABP4 in stomach adenocarcinoma. Biomed Res. Int. 2020:4145164. doi: 10.1155/2020/4145164
Hale, C., Bartholomew, M., Taylor, V., Stables, J., Topley, P., and Tite, J. (1996). Recognition of CD52 allelic gene products by CAMPATH-1H antibodies. Immunology 88, 183–190. doi: 10.1111/j.1365-2567.1996.tb00003.x
Hui, L., and Chen, Y. (2015). Tumor microenvironment: sanctuary of the devil. Cancer Lett. 368, 7–13. doi: 10.1016/j.canlet.2015.07.039
Jia, D., Li, S., Li, D., Xue, H., Yang, D., and Liu, Y. (2018). Mining TCGA database for genes of prognostic value in glioblastoma microenvironment. Aging 10, 592–605. doi: 10.18632/aging.101415
Kashima, H., Noma, K., Ohara, T., Kato, T., Katsura, Y., Komoto, S., et al. (2019). Cancer-associated fibroblasts (CAFs) promote the lymph node metastasis of esophageal squamous cell carcinoma. Int. J. Cancer 144, 828–840. doi: 10.1002/ijc.31953
Kemeny, H. R., Elsamadicy, A. A., Farber, S. H., Champion, C. D., Lorrey, S. J., Chongsathidkiet, P., et al. (2020). Targeting PD-L1 initiates effective antitumor immunity in a murine model of cushing disease. Clin. Cancer Res. 26, 1141–1151. doi: 10.1158/1078-0432.CCR-18-3486
Knosp, E., Steiner, E., Kitz, K., and Matula, C. (1993). Pituitary adenomas with invasion of the cavernous sinus space: a magnetic resonance imaging classification compared with surgical findings. Neurosurgery 33, 610–617. doi: 10.1227/00006123-199310000-00008
Konduri, S. D., Tasiou, A., Chandrasekar, N., and Rao, J. S. (2001). Overexpression of tissue factor pathway inhibitor-2 (TFPI-2), decreases the invasiveness of prostate cancer cells in vitro. Int. J. Oncol. 18, 127–131. doi: 10.3892/ijo.18.1.127
Landeiro, J. A., Fonseca, E. O., Monnerat, A. L., Taboada, G. F., Cabral, G. A., and Antunes, F. (2015). Nonfunctioning giant pituitary adenomas: invasiveness and recurrence. Surg. Neurol. Int. 6:179. doi: 10.4103/2152-7806.170536
Leslie, N. R., Biondi, R. M., and Alessi, D. R. (2001). Phosphoinositide-regulated kinases and phosphoinositide phosphatases. Chem. Rev. 101, 2365–2380. doi: 10.1021/cr000091i
Li, L., Peng, M., Xue, W., Fan, Z., Wang, T., Lian, J., et al. (2018). Integrated analysis of dysregulated long non-coding RNAs/microRNAs/mRNAs in metastasis of lung adenocarcinoma. J. Transl. Med. 16:372. doi: 10.1186/s12967-018-1732-z
Lindholm, J., Nielsen, E. H., Bjerre, P., Christiansen, J. S., Hagen, C., Juul, S., et al. (2006). Hypopituitarism and mortality in pituitary adenoma. Clin. Endocrinol. 65, 51–58. doi: 10.1111/j.1365-2265.2006.02545.x
Liu, Z., Zhu, Y., Xu, L., Zhang, J., Xie, H., Fu, H., et al. (2018). Tumor stroma-infiltrating mast cells predict prognosis and adjuvant chemotherapeutic benefits in patients with muscle invasive bladder cancer. Oncoimmunology 7:e1474317. doi: 10.1080/2162402X.2018.1474317
Long, R., Liu, Z., Li, J., and Yu, H. (2019). COL6A6 interacted with P4HA3 to suppress the growth and metastasis of pituitary adenoma via blocking PI3K-Akt pathway. Aging 11, 8845–8859. doi: 10.18632/aging.102300
Lu, J. Q., Adam, B., Jack, A. S., Lam, A., Broad, R. W., and Chik, C. L. (2015). Immune cell infiltrates in pituitary adenomas: more macrophages in larger adenomas and more t cells in growth hormone adenomas. Endocr. Pathol. 26, 263–272. doi: 10.1007/s12022-015-9383-6
Luo, Y., Zeng, G., and Wu, S. (2019). Identification of microenvironment-related prognostic genes in bladder cancer based on gene expression profile. Front. Genet. 10:1187. doi: 10.3389/fgene.2019.01187
Lv, L., Zhang, S., Hu, Y., Zhou, P., Gao, L., Wang, M., et al. (2018). Invasive pituitary adenoma-derived tumor-associated fibroblasts promote tumor progression both in vitro and in vivo. Exp. Clin. Endocrinol. Diabetes 126, 213–221. doi: 10.1055/s-0043-119636
Maolake, A., Izumi, K., Shigehara, K., Natsagdorj, A., Iwamoto, H., Kadomoto, S., et al. (2017). Tumor-associated macrophages promote prostate cancer migration through activation of the CCL22-CCR4 axis. Oncotarget 8, 9739–9751. doi: 10.18632/oncotarget.14185
Marques, P., Barry, S., Carlsen, E., Collier, D., Ronaldson, A., Awad, S., et al. (2019). Chemokines modulate the tumour microenvironment in pituitary neuroendocrine tumours. Acta Neuropathol. Commun. 7:172. doi: 10.1186/s40478-019-0830-3
Newman, A. M., Liu, C. L., Green, M. R., Gentles, A. J., Feng, W., Xu, Y., et al. (2015). Robust enumeration of cell subsets from tissue expression profiles. Nat. Methods 12, 453–457. doi: 10.1038/nmeth.3337
Niu, D., Peatman, E., Liu, H., Lu, J., Kucuktas, H., Liu, S., et al. (2011). Microfibrillar-associated protein 4 (MFAP4) genes in catfish play a novel role in innate immune responses. Dev. Comp. Immunol. 35, 568–579. doi: 10.1016/j.dci.2011.01.002
Ntali, G., and Wass, J. A. (2018). Epidemiology, clinical presentation and diagnosis of non-functioning pituitary adenomas. Pituitary 21, 111–118. doi: 10.1007/s11102-018-0869-3
Ocana, M. C., Martinez-Poveda, B., Quesada, A. R., and Medina, M. A. (2019). Metabolism within the tumor microenvironment and its implication on cancer progression: an ongoing therapeutic target. Med. Res. Rev. 39, 70–113. doi: 10.1002/med.21511
Ostrom, Q. T., Gittleman, H., Fulop, J., Liu, M., Blanda, R., Kromer, C., et al. (2015). CBTRUS Statistical Report: Primary Brain and Central Nervous System Tumors Diagnosed in the United States in 2008-2012. Neuro. Oncol. 17(Suppl. 4), iv1–iv62. doi: 10.1093/neuonc/nov189
Pan, X. B., Lu, Y., Huang, J. L., Long, Y., and Yao, D. S. (2019). Prognostic genes in the tumor microenvironment in cervical squamous cell carcinoma. Aging 11, 10154–10166. doi: 10.18632/aging.102429
Pilia, G., Hughes-Benzie, R. M., Mackenzie, A., Baybayan, P., Chen, E. Y., Huber, R., et al. (1996). Mutations in GPC3, a glypican gene, cause the Simpson-Golabi-Behmel overgrowth syndrome. Nat. Genet. 12, 241–247. doi: 10.1038/ng0396-241
Qiu, L., Yang, J., Wang, H., Zhu, Y., Wang, Y., and Wu, Q. (2013). Expression of T-helper-associated cytokines in the serum of pituitary adenoma patients preoperatively and postperatively. Med. Hypotheses 80, 781–786. doi: 10.1016/j.mehy.2013.03.011
Ren, B., Cui, M., Yang, G., Wang, H., Feng, M., You, L., et al. (2018). Tumor microenvironment participates in metastasis of pancreatic cancer. Mol. Cancer 17:108. doi: 10.1186/s12943-018-0858-1
Rollin, J., Iochmann, S., Bléchet, C., Hubé, F., Régina, S., Guyétant, S., et al. (2005). Expression and methylation status of tissue factor pathway inhibitor-2 gene in non-small-cell lung cancer. Br. J. Cancer 92, 775–783. doi: 10.1038/sj.bjc.6602298
Sam, A. H., Shah, S., Saleh, K., Joshi, J., Roncaroli, F., Robinson, S., et al. (2015). Clinical outcomes in patients with nonfunctioning pituitary adenomas managed conservatively. Clin. Endocrinol. (Oxf) 83, 861–865. doi: 10.1111/cen.12860
Sato, M., Tamura, R., Tamura, H., Mase, T., Kosugi, K., Morimoto, Y., et al. (2019). Analysis of tumor angiogenesis and immune microenvironment in non-functional pituitary endocrine tumors. J. Clin. Med. 8:695. doi: 10.3390/jcm8050695
Sato, N., Parker, A. R., Fukushima, N., Miyagi, Y., Iacobuzio-Donahue, C. A., Eshleman, J. R., et al. (2005). Epigenetic inactivation of TFPI-2 as a common mechanism associated with growth and invasion of pancreatic ductal adenocarcinoma. Oncogene 24, 850–858. doi: 10.1038/sj.onc.1208050
Sawada, Y., Yoshikawa, T., Ofuji, K., Yoshimura, M., Tsuchiya, N., Takahashi, M., et al. (2016). Phase II study of the GPC3-derived peptide vaccine as an adjuvant therapy for hepatocellular carcinoma patients. Oncoimmunology 5, e1129483. doi: 10.1080/2162402X.2015.1129483
Sayem, M. A., Tomita, Y., Yuno, A., Hirayama, M., Irie, A., Tsukamoto, H., et al. (2016). Identification of glypican-3-derived long peptides activating both CD8(+) and CD4(+) T cells; prolonged overall survival in cancer patients with Th cell response. Oncoimmunology 5:e1062209. doi: 10.1080/2162402X.2015.1062209
Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504. doi: 10.1101/gr.1239303
Szklarczyk, D., Gable, A. L., Lyon, D., Junge, A., Wyder, S., Huerta-Cepas, J., et al. (2019). STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 47, D607–D613. doi: 10.1093/nar/gky1131
Tahmasebi Birgani, M., and Carloni, V. (2017). Tumor microenvironment, a paradigm in hepatocellular carcinoma progression and therapy. Int. J. Mol. Sci. 18:405. doi: 10.3390/ijms18020405
Wang, Y., Zeng, J., Wu, W., Xie, S., Yu, H., Li, G., et al. (2019). Nicotinamide N-methyltransferase enhances chemoresistance in breast cancer through SIRT1 protein stabilization. Breast Cancer Res. 21:64. doi: 10.1186/s13058-019-1150-z
Wass, J. A., and Karavitaki, N. (2009). Nonfunctioning pituitary adenomas: the Oxford experience. Nat. Rev. Endocrinol. 5, 519–522. doi: 10.1038/nrendo.2009.147
Wu, Q., Pi, L., Le Trinh, T., Zuo, C., Xia, M., Jiao, Y., et al. (2017). A novel vaccine targeting glypican-3 as a treatment for hepatocellular carcinoma. Mol. Ther. 25, 2299–2308. doi: 10.1016/j.ymthe.2017.08.005
Wu, S., Ou, T., Xing, N., Lu, J., Wan, S., Wang, C., et al. (2019). Whole-genome sequencing identifies ADGRG6 enhancer mutations and FRS2 duplications as angiogenesis-related drivers in bladder cancer. Nat. Commun. 10:720. doi: 10.1038/s41467-019-08576-5
Wu, T., and Dai, Y. (2017). Tumor microenvironment and therapeutic response. Cancer Lett. 387, 61–68. doi: 10.1016/j.canlet.2016.01.043
Xie, X., Yu, H., Wang, Y., Zhou, Y., Li, G., Ruan, Z., et al. (2014). Nicotinamide N-methyltransferase enhances the capacity of tumorigenesis associated with the promotion of cell cycle progression in human colorectal cancer cells. Arch. Biochem. Biophys. 564, 52–66. doi: 10.1016/j.abb.2014.08.017
Xiong, Y., Wang, K., Zhou, H., Peng, L., You, W., and Fu, Z. (2018). Profiles of immune infiltration in colorectal cancer and their clinical significant: a gene expression-based study. Cancer Med. 7, 4496–4508. doi: 10.1002/cam4.1745
Yagnik, G., Rutowski, M. J., Shah, S. S., and Aghi, M. K. (2019). Stratifying nonfunctional pituitary adenomas into two groups distinguished by macrophage subtypes. Oncotarget 10, 2212–2223. doi: 10.18632/oncotarget.26775
Yang, J., Song, H., Chen, L., Cao, K., Zhang, Y., Li, Y., et al. (2019). Integrated analysis of microfibrillar-associated proteins reveals MFAP4 as a novel biomarker in human cancers. Epigenomics 11, 1635–1651. doi: 10.2217/epi-2018-0080
Keywords: non-functioning pituitary adenoma, tumor microenvironment, immune, invasive, recurrence
Citation: Guo J, Fang Q, Liu Y, Xie W, Li C and Zhang Y (2021) Screening and Identification of Key Microenvironment-Related Genes in Non-functioning Pituitary Adenoma. Front. Genet. 12:627117. doi: 10.3389/fgene.2021.627117
Received: 08 November 2020; Accepted: 08 March 2021;
Published: 27 April 2021.
Edited by:
Lixin Cheng, Jinan University, ChinaCopyright © 2021 Guo, Fang, Liu, Xie, Li and Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yazhuo Zhang, zyz2004520@yeah.net