Skip to main content

Identification and gene expression profiling of human gonadotrophic pituitary adenoma stem cells



Gonadotrophic pituitary adenoma is a major subtype of pituitary adenoma in the sellar region, but it is rarely involved in the hypersecretion of hormones into blood; thus, it is commonly regarded as “non-functioning.” Its tumorigenic mechanisms remain unknown. The aim of this study was to identify human gonadotrophic pituitary adenoma stem cells (hPASCs) and explore the underlying gene expression profiles. In addition, the potential candidate genes involved in the invasive properties of pituitary adenoma were examined.


The hPASCs from 14 human gonadotrophic pituitary adenoma clinical samples were cultured and verified via immunohistochemistry. Genetic profiling of hPASCs and the matched tumor cells was performed through RNA-sequencing and subjected to enrichment analysis. By aligning the results with public databases, the candidate genes were screened and examined in invasive and non-invasive gonadotrophic pituitary adenomas using Real-time polymerase chain reaction.


The hPASCs were successfully isolated and cultured from gonadotrophic pituitary adenoma in vitro, which were identified as positive for generic stem cell markers (Sox2, Oct4, Nestin and CD133) via immunohistochemical staining. The hPASCs could differentiate into the tumor cells expressing follicle-stimulating hormone in the presence of fetal bovine serum in the culture medium. Through RNA-sequencing, 1352 differentially expressed genes were screened and identified significantly enriched in various gene ontologies and important pathways. The expression levels of ANXA2, PMAIP1, SPRY2, C2CD4A, APOD, FGF14 and FKBP10 were significantly upregulated while FNDC5 and MAP3K4 were downregulated in the invasive gonadotrophic pituitary adenomas compared to the non-invasive ones.


Genetic profiling of hPASCs may explain the tumorigenesis and invasiveness of gonadotrophic pituitary adenoma. ANXA2 may serve as a potential therapeutic target for the treatment of gonadotrophic pituitary adenoma.


Pituitary adenoma accounts for approximately 15% of primary intracranial neoplasms and is generally considered to be benign [1]. Gonadotrophic lineage adenoma, the second most common subtype of pituitary adenoma following lactotroph lineage [2], can express follicle-stimulating hormone (FSH) or leucine hormone (LH) that may not be hypersecreted into the blood. The most common clinical symptom of gonadotrophic pituitary adenoma is compression to the adjacent structures such as the optic chiasma or invasion of the cavernous sinus [3]. Therefore, this subtype of pituitary adenoma is often considered “non-functioning.” The mechanisms of gonadotrophic pituitary adenoma tumorigenesis are still undefined, which poses challenges to the development of effective adjuvant therapy following surgical resection.

Tumor stem cells (TSC) are identified as a specific subpopulation of cells possessing stem-cell properties in certain types of bulky tumors [4, 5]. TSC have been shown to manifest a certain degree of resistance to radio- or chemotherapy, which is one possible reason for tumor relapse and metastasis [6]. Presumably, they arise from the normal stem/progenitor cells in the postnatal organs [7]. In humans, single-cell transcriptomic analysis discovered that the developing pituitary gland has a major stem cell population [8]. Accumulating evidence has pointed to a specific subpopulation of pituitary stem or progenitor cells that may give rise to pituitary tumors [9]. Previously, Xu isolated human pituitary adenoma stem cells (hPASCs) and verified their cellular characteristics, including sphere-forming, self-renewal, and multipotency [10]. Compared to differentiated daughter cells, these adenoma stem-like cells could express increased levels of stemness-associated genes, anti-apoptotic proteins, and pituitary progenitor cell markers and form tumors in the mouse model [10]. Later, Wurth characterized more cellular properties of hPASCs in somatotropinoma and non-functioning pituitary adenoma [11].

To date, research on pituitary adenoma stem cells has been conducted mainly in somatotropinoma and non-functioning adenoma [12], while very few studies have been performed on other subtypes. Here, we isolated and cultured hPASCs from human gonadotrophic pituitary adenoma samples resected during neurosurgery. The aim was to elucidate the genetic profiling of hPASCs via RNA-sequencing and to further test potential candidate genes in invasive gonadotrophic pituitary adenoma. Using this in vitro hPASC model, we hope to achieve a better understanding of the tumorigenesis and invasive behavior of gonadotrophic pituitary adenoma.

Materials and methods

Patients and samples

Fourteen human pituitary adenoma samples obtained during neurosurgery in our department were used for hPASC culture and sequencing. The detailed clinical information of each patient is shown in Additional file 1: Table S1.

Clinically, invasive pituitary adenoma can be defined using the Knosp grade on magnetic resonance imaging [13]. If the edge of the tumor is extended beyond the lateral tangent of the intra- and supra-cavernous internal carotid artery, it is defined invasive and the grade will be 3 or 4. Wherease, non-invasive tumors were located within the tangent of the internal carotid artery and categorized as grades 0, 1 and 2 [13]. Seven non-invasive and nine invasive gonadotrophic pituitary adenoma samples were used for the verification of candidate gene expression. The clinical information is summarized in Additional file 1: Table S2. All the samples were pathologically diagnosed as gonadotrophic pituitary adenoma by two pathologists. This study was approved by the Ethics Committee of Beijing Tiantan Hospital, and written informed consent was available for all patients.

Culture and differentiation of hPASCs

During surgery, the pituitary adenoma specimen was resected and immediately divided into two portions. One portion was “snap-frozen” and stored in liquid nitrogen for pathological immunostaining and sequencing analysis, and the other portion was transferred to the laboratory for hPASC culture. The tumor specimen was thoroughly washed with 1X phosphate buffer saline (PBS) and minced with a sharp scalpel into small pieces under sterile conditions. The cells were collected in 10 ml Dulbecco’s modified Eagle’s Medium /F-12 (DMEM/F-12) and centrifuged at 1000 rpm for 2 min. The cell pellet was dissolved in 2 ml 1X Accutase enzyme (Stemcell Technology, USA) and incubated at 37 °C for 5 min. The dissociated cells were centrifuged at 1,000 rpm for 3 min, and the cell pellet was incubated with red blood cell lysis buffer (154 mM NH4Cl, 10 mM KHCO3, 0.1 mM EDTA, pH 7.4) for 5 min to eradicate contamination from red blood cells. After washing, the cells were dissolved and cultured in the stem cell medium composed of DMEM/F-12 supplemented with 2 mM L-glutamine, 1% penicillin–streptomycin, 1X B27 (50X, Life Technologies, USA), 20 ng/ml basic fibroblast growth factor (bFGF, Peprotech) and 20 ng/ml epidermal growth factor (EGF, Peprotech) at 37 °C with 5% CO2 in an incubator. The stem cell medium was refreshed every 5 days.

For differentiation, the cultured hPASCs were seeded in a 24-well plate, and the culture medium was replaced with the differentiation medium composed of DMEM/F-12 supplemented with 15% fetal bovine serum (FBS), 2 mM L-glutamine, and 1% penicillin–streptomycin. The differentiation medium was refreshed every 5 days.


The hPASCs were cultured in a poly-D-lysine-coated 6-well plate for attachment, then washed with 1X PBS and fixed with 4% paraformaldehyde for 30 min at room temperature. The fixed cells were incubated in 0.3% Triton™ X-100 (Sigma-Aldrich, USA) for 15 min and washed with 1X PBS three times. The cells were blocked with 5% bovine serum albumin (BSA) for 1 h and incubated with primary antibody at 4 °C overnight. The cells were then washed with 1X PBS and incubated with secondary antibody conjugated to Alexa-Fluor 488 or Alexa-Fluor 647 (abcam, USA) at 1:1000 dilution for 1 h at room temperature. The images were captured using a Zeiss microscopic imaging system with fluorescence emission system at different magnifications. For tissue staining, the samples were sectioned at a thickness of 10 μm per slice and subjected to the above procedures.

All antibodies were purchased from Abcam unless otherwise stated. The antibodies used in the experiments were as follows: rabbit anti-Sox2 (1:800), rabbit anti-Oct4 (1:500), rabbit anti-Nestin (1:200), rabbit anti-CD133 (1:500), rabbit anti-ANXA2(1:1000), and mouse anti-FSHβ (1:800, Santa Cruz Biotechnology).

Cell count

The cells in each culture well were washed carefully with 1X PBS before harvesting. Then, the cells were collected in a 1.5-ml tube and centrifuged at 1000 rpm for 20 s. The cell pellet was resuspended in 1-ml DMEM. For counting, 10-μl of cell suspension was added to a hemacytometer, and the number of cells in the central gridded area was manually counted under the microscope. For each sample, the total number was calculated with four separate wells, each performed in duplicate.

Enzyme-linked immunosorbent assay (ELISA)

ELISA was performed to quantify the level of FSH in the culture medium according to manufacturer’s protocol (Elabscience, USA). Briefly, 100 μl standard or culture medium sample was added to each well and incubated for 90 min at 37 °C. The liquid was then replaced with 100 μl biotinylated detection Ab for 1 h at 37 °C. After washing, 100 μl HRP conjugate was added to the well and incubated for 30 min at 37 °C. Each well was washed thoroughly before 90 μl substrate reagent was added and incubated for 15 min at 37 °C. The reaction was stopped by adding 50 μl stop solution, and the mixture was evaluated at 450 nm using a spectrophotometer (Bio-Tek, USA).

RNA extraction

Total RNA from the hPASCs or tumor specimen was extracted by TRIzol reagent (Invitrogen, USA) according to manufacturer’s instructions. RNA integrity was assessed using the RNA Nano 6000 Assay Kit of the Bioanalyzer 2100 system (Agilent Technologies, USA). Total RNA was stored at − 80 °C before use.

Preparation of cDNA library and RNA-sequencing

RNA sequencing was performed by the commercial provider Novogene Co. Ltd. Briefly, first-strand cDNA was synthesized using a random hexamer primer and M-MuLV Reverse Transcriptase (RNase H). Second strand cDNA synthesis was performed using DNA Polymerase I and RNase H. After the preferential selection of 250–300-bp cDNA fragments, polymerase chain reaction (PCR) was performed using Phusion High-Fidelity DNA polymerase, universal PCR primers, and Index (X) Primer. PCR products were purified (AMPure XP system), and the quality of cDNA library was assessed on the Agilent Bioanalyzer 2100 system.

Clustering of the index-coded samples was performed on a cBot Cluster Generation System using TruSeq PE Cluster Kit v3-cBot-HS (Illumina) according to the manufacturer’s instructions. After cluster generation, the library preparations were sequenced on an Illumina Novaseq platform.

Abundance of differentially expressed genes (DEG)

After read mapping, the counts for differentially expressed genes from each sequenced library were adjusted and analyzed using the “edgeR” package from R software (3.22.5). The P-values were adjusted by the Benjamini & Hochberg method. An adjusted P-value of 0.05 and fold change of Log2FC were set as the threshold for significant differential expression.

Enrichment analysis of differentially expressed genes

Gene Ontology (GO) enrichment analysis of differentially expressed genes was implemented using the cluster Profiler R package. GO terms with corrected p < 0.05 were considered significantly enriched by differentially expressed genes.

The Cluster Profiler R package was used to test the statistical enrichment of differential expression genes in the pathways.

Protein–protein interaction (PPI) network

GSE26966 dataset from GEO (Gene Expression Omnibus) contains the expression profiles of 14 gonadotrophic pituitary adenomas in comparison to 3 normal pituitary tissues in the published research [14]. The differentially expressed genes of hPASCs were selected according to the following criteria: Log2FC > 1.5 or <  − 1.5, and adjusted p < 0.05. The genes presented in both datasets were illustrated using a Venn diagram, which were further subjected to PPI analysis using the STRING database ( The network was outlined by Cytoscape software (3.9.1).

Real-time PCR

First-strand cDNA was synthesized by combining total RNA (2 μg) with a reverse transcription mixture (RevertAid™ First Strand cDNA Synthesis Kit, Thermo-fisher, USA) in a total reaction volume of 20 μl, according to the manufacturer’s instructions. The mixture was incubated at 42 °C for 60 min and 70 °C for 10 min and then returned to ice. One microliter of the mixture was extracted for Real-time PCR quantification together with 25-μl SYBR Green PCR Master Mix (TaKaRa, Japan) in a TaKaRa 7500 Real-time system (TaKaRa, Japan). β-Actin served as the internal control. The relative expression level was quantified using the 2−∆∆ct method. The primers used for qPCR are mentioned in Additional file 1: Table S3.

Statistical analysis

Statistical analyses were performed using GraphPad Prism 6.0 software (GraphPad Software, Inc., USA). Data are shown as the mean ± standard error of the mean (SEM) or the mean ± standard deviation (SD). A t-test was used to compare the statistical differences. A P-value < 0.05 was considered to indicate a significant difference.


Culture and verification of hPASCs

Human gonadotrophic pituitary adenoma stem cells were cultured and verified in each sample of 14 patients. The clinical information regarding the included patients is listed in Table 1. On day 0, the tumor cells were completely dissociated into single cells and cultured in serum-free medium supplemented with EGF and bFGF. Small spheres appeared to form after 3 days and became more apparent on day 14. At this point, no other cells such as erythrocytes or fibroblasts were obviously present in the culture by microscopic observation. Figure 1A shows the representative morphological appearance of hPASCs from one patient at different time points. With the addition of a fresh stem culture medium, the spheres continued growing for 21 days. Moreover, the size of the hPASC spheres varied between samples. Figure 1B shows the morphology of hPASC from the tumor samples of three different patients on day 14.

Table 1 Summary of clinical information for hPASC culture
Fig. 1
figure 1

The morphology of human gonadotrophic adenoma stem cells. A Representative morphological changes of hPASC on day 3, day 14 and day 21 by microscopic observation from Case 3. B Representative morphology of hPASC culture from patient Case 1, Case 4 and Case 8 on day 14. C Immunohistochemistry of Sox2, Oct4, CD133 and Nestin in hPASCs cultured on day 14. All the markers were positively stained in the hPASCs. Scale bar is 100um for all panels

Immunostaining was performed to examine the expression of the generic stem cell markers Sox2, Oct4, CD133 and Nestin in the hPASCs cultured in the stem cell medium on day 14. As shown in Fig. 1C, positive staining was observed for all four markers in the cultured hPASCs.

Differentiation of hPASCs

The differentiation capacity of the cultured hPASCs was examined by supplementing FBS at a final concentration of 15% to the medium. After 5 days of differentiation, the expression level of Sox2 was reduced in differentiated daughter cells compared to the hPASCs. FSHβ, a marker for FSH, was positively stained in the differentiated cells but not the hPASCs (Fig. 2A). After the 5-day differentiation, the number of differentiated daughter cells was significantly increased by 2.19 ± 0.29 (p < 0.005) compared to the hPASC (Fig. 2B, C). However, the level of FSH secreted by the differentiated cells into the medium was not significantly changed (Fig. 2D).

Fig. 2
figure 2

Human gonadotrophic adenoma stem cells possessed the capacity of differentiation. A Immunostaining of FSH and Sox2 in hPASCs and the differentiated daughter cells. Compared to the hPASCs, the differentiated cells expressed a much lower level of Sox2 but strong FSH. B Representative images of hPASCs and the differentiated cells. C Quantification of cell number between hPASCs and the differentiated cells. After 5-day differentiation, the number of cells in the differentiated group was significantly higher than the hPASC group. (**p < 0.005; n = 4 independent experiments). D Measurement of FSH level in the hPASC and the differentiated cells by ELISA. The level of FSH in the differentiated cells was not significantly changed compared to the hPASC group (p > 0.05; n = 4 independent experiments). d-cells: differentiated cells. Scale bar is 100um for all panels

DEG profiles between hPASC and matched tumor samples

To explore the genetic profiling of gonadotrophic pituitary adenoma stem cells, RNA-sequencing was performed in the hPASCs and the matched tumor samples. Principal component analysis (PCA) showed that both groups manifested a significant distinction regarding the overall expression pattern (Fig. 3A). There were 1352 differentially expressed genes (Log2FC > 1 or <  − 1, adjusted p < 0.05) between the hPASCs and matched tumor samples. In detail, 843 genes were significantly upregulated (Log2FC > 1, adjusted p < 0.05) and 509 genes were downregulated (Log2FC <  − 1, adjusted p < 0.05) in the hPASC group. The volcano plot depicting the expression pattern of these genes is shown in Fig. 3B.

Fig. 3
figure 3

Illustration of differentially expressed genes in hPASC and the matched tumor cells. A Principal Components Analysis (PCA) showed a distinctive expression pattern between hPASCs (blue) and the matched tumor cells (red). B Volcano plot showed the distribution of the upregulated and downregulated genes presented in hPASCs and the tumor cells. Those upregulated genes (Log2FC > 1, p < 0.05) were represented as red dots while downregulated ones were in blue (Log2FC < -1, p < 0.05). C KEGG enrichment of differentially expressed genes between hPASC and the matched tumor cells. (D-E) GO analysis of differentially expressed genes between hPASC and the matched tumor cells

Based on the differentially expressed genes, enrichment analysis was performed to reveal key factors and signaling pathways. The top-rated KEGG pathways were focal adhesion, ECM–receptor interaction and proteoglycans in cancer (Fig. 3C). GO analysis identified key events, including the intrinsic apoptotic signaling pathway, the sterol metabolic process, Rho protein signal transduction, the regulation of angiogenesis, and cell adhesion (Fig. 3D). These molecules were mainly enriched in the basement membrane, transmembrane transporter complex, secretory granule lumen, and endoplasmic reticulum lumen (Fig. 3E).

Core protein–protein interaction (PPI) network

GSE26966 dataset contained the expression profiles of 14 gonadotrophic pituitary adenomas compared with 3 normal pituitary tissues. To screen the essential genes altered during the tumorigenesis of hPASCs, 512 of 1352 candidate genes were selected according to specific criteria (Log2FC ≥ 1.5 or ≤  − 1.5, p < 0.05) through RNA-sequencing and aligned with the genes listed in GSE26966. One hundred and thirteen genes were shared in both datasets (Fig. 4A), and their expression pattern was shown in Fig. 4B between hPASCs and the matched tumors. The PPI network among those genes was predicted using the string database (Fig. 4C).

Fig. 4
figure 4

Analysis of shared differentially expressed genes in both GSE26966 dataset and hPASCs. A Venn diagram showed 113 differentially expressed genes were exclusively contained in both GSE26966 dataset and RNA-sequencing results of hPASCs. B Heatmap of 113 genes expressed in hPASCs and the matched tumor cells. C Protein–protein Interaction (PPI) of 113 genes was analyzed by and the network was depicted using cytoscape

Expression of candidate genes in non-invasive and invasive gonadotrophic pituitary adenomas

The invasiveness of pituitary adenoma into the surrounding structure indicated the aggressiveness of the tumor. In clinical practice, pituitary adenomas at Knosp grades 3 or 4 are generally considered invasive, while the grade is 0–2 normally for non-invasive tumors (Fig. 5A). Nine genes were significantly differentially expressed between seven non-invasive and nine invasive gonadotrophic pituitary adenoma samples (Fig. 5B). In the invasive gonadotrophic pituitary adenoma (Knosp grade 3–4), seven genes were significantly upregulated (ANXA2, PMAIP1, SPRY2, C2CD4A, APOD, FGF14 and FKBP10) and two genes were downregulated (FNDC5 and MAP3K4).

Fig. 5
figure 5

Measurement of gene expression between non-invasive and invasive gonadotrophic pituitary adenoma. A MRI images of pituitary adenoma in Knosp grade 1 and Knosp grade 4. B The expression of nine candidate genes was quantified and found significantly altered between 7 non-invasive (Knosp grade 1–2) and 9 invasive (Knosp grade 3–4) human gonadotrophic pituitary adenoma clinical samples. *p < 0.05; **p < 0.01; ****p < 0.001

Immunohistochemistry staining

The protein level of ANXA2 was examined via immunohistochemical staining of the clinical samples. As shown in Fig. 6, the immunoreactivity of ANXA2 was significantly stronger in the invasive gonadotrophic pituitary adenoma (Knosp grade 4) than in the non-invasive one (Knosp grade 1).

Fig. 6
figure 6

Immunohistochemistry of ANXA2 in non-invasive (Knosp grade 1) and invasive gonadotrophic pituitary adenoma (Knosp grade 4). Magnification scale × 200


Whole-exome sequencing of silent gonadotrophic pituitary adenoma has detected some putative driver genes such as platelet-derived growth factor-D (PDGFD), N-myc downregulated gene family member 4 (NDRG4), and Zipper sterile-α-motif kinase (ZAK); however, further analysis indicated that these genes do not significantly contribute to tumorigenesis [15]. Other epigenetic factors, such as microRNAs, DNA methylation, and environmental factors were also proposed to explain the molecular mechanism [16,17,18]. Tumor stem cells have been proposed to explain the tumorigenesis and biological properties of pituitary adenoma for many years. Wurth’s review has summarized convincing evidence of the existence of hPASCs in pituitary adenoma and their potential roles [12]. In this study, hPASCs from gonadotrophic pituitary adenoma were successfully established according to a published protocol with some modifications. To the best of our knowledge, no previous studies have evaluated the genetic profiling of hPASCs in this tumor subtype.

The morphology of the hPASCs forming pituispheres was similar to that reported in the previous literature [10, 11]. The culture time-frame of hPASCs in this study was 21 days, but in our observation, the pituispheres could be expanded continuously to a maximum of 60 days (data not shown). This suggested the self-renewal capacity of these hPASCs. Immunohistochemical staining was commonly used for verification of tumor stem cells in pituitary adenoma; four generic markers, namely Sox2 [11, 19,20,21], CD133 [22, 23], Nestin [10, 23] and Oct4 [23, 24] were selected and showed immunopositivity in the cultured hPASCs, thus validating the properties of their stemness.

Next, we tested whether the hPASC pituispheres had the capacity to differentiate into tumor cells following the previous protocol [10, 11]. After adding serum to the hPASC culture for 5 days, the total number of cells was significantly increased, indicating that hPASCs could generate some new tumor cells. This is similar to some other types of tumor stem cells like glioblastoma, possibly due to intrinsic alteration of the cell cycle [25]. Moreover, serum supplementation of the gonadotrophic adenoma stem cells not only drove the proliferation but also facilitated the expression of the respective hormone FSH. This was consistent with the previous observation on somatotropinoma stem cells [11]. Although the FSH was expressed in the differentiated cells, ELISA results showed that the level of hormone in the medium from the differentiated cells was not significantly changed. We speculated that the secretive mechanism of the hormone from the cells was not altered during this process, or else the intrinsic molecular properties of the tumor remained unchanged.

Peculis demonstrated that somatic mutations of primary pituitary tumor samples could still be detected in the exomes of respective tumor pituispheres [26], indicating that gene expression, rather than mutations, play a primary role during tumorigenesis. Our RNA-sequencing data from both hPASCs and matched tumor samples have identified the differentially expressed genes during the process. By enrichment analysis, the putative top signaling pathways from KEGG analysis were focal adhesion, ECM–receptor interaction, and proteoglycans in cancer. GO analysis predicted hemostasis, angiogenesis, MAP kinase activity, and cell adhesion involved in the tumorigenesis. Focal adhesion involves many pro-survival signaling molecules, including growth factor receptors, intracellular molecules, and integrins, and plays important roles in the regulation of cell mobility, proliferation, and tumor cell survival [27]. The MAPK pathway has also been widely reported in various studies on pituitary adenoma tumorigenesis [28]. More work is needed to elucidate other potential key pathways during the process of differentiation of hPASCs.

To narrow down the scope of candidate genes, we aligned our RNA-sequencing data with public databases. The rationale was to determine whether the genes that played an essential role in the tumorigenesis co-existed in both databases. The alignment indicated 113 genes, and their interactions can shed a light on the underlying molecular mechanisms of tumorigenesis. Since TSC has been shown to be associated with the aggressiveness of tumors in previous studies [29], we tested the expression of some candidate genes in the invasive and non-invasive gonadotrophic pituitary adenomas. In clinical practice, the Knosp grading system has been widely used in the pituitary adenoma surgery assessment [13]. Some previous studies have already used this system to screen the differences in genetic profiles between invasive and non-invasive pituitary adenomas [30]. Nine genes were verified as being significantly differentially expressed between non-invasive and invasive gonadotrophic pituitary adenomas, and further staining indicates that the gene that produces the ANXA2 protein was one of the key genes. ANXA2 is an important member of the annexin protein family and involved in various cellular functions, including cell division, calcium signaling, and cell growth [31]. Recent studies have identified ANXA2 expression in several cancers where it promotes tumor progression via proliferation, migration, epithelial-mesenchymal transformation (EMT), invasion and stem cell formation [32]. Further analysis of the molecular mechanisms of ANXA2 in the tumorigenesis and invasion of pituitary adenoma is thus needed.


In this study, we provide insight into the molecular biology of hPASC cultured from human gonadotrophic pituitary adenoma. The presence and differentiation capacity of hPASC were validated in the culture. The genetic profiles of hPASCs and their matched tumor samples were compared. Enrichment analysis identified focal adhesion, ECM–receptor interaction, proteoglycans in cancer, intrinsic apoptotic signaling pathway, sterol metabolic process, Rho protein signal transduction, and regulation of angiogenesis as pivotal pathways involved in the tumorigenesis. There were seven candidate genes significantly increased in the invasive human gonadotrophic pituitary adenoma compared to the non-invasive one. ANXA2 was verified to be one of the key candidate genes in the process of tumor invasiveness. Further studies should be conducted to provide more details on the molecular mechanism of ANXA2 in tumorigenesis and invasiveness. In addition, more therapeutic strategies for treating hPASCs should be developed and assessed, which may reduce tumor invasion and relapse.

Availability of data and materials

The datasets generated and/or analysed during the current study are not publicly available because it is used for further analysis, but are available from the corresponding author on reasonable request.



Human pituitary adenoma stem cells


Follicle-stimulating hormone


Leucine hormone


Tumor stem cells


Phosphate buffer saline


Dulbecco's modified eagle’s medium


Basic fibroblast growth factor


Epidermal growth factor


Fetal bovine serum


Enzyme-linked immunosorbent assay


Gene ontology (GO)


Kyoto encyclopedia of genes and genomes (KEGG)


Gene expression omnibus


Principal component analysis


Protein–protein interaction


  1. Sivakumar W, Chamoun R, Nguyen V, Couldwell WT (2011) Incidental pituitary adenomas. Neurosurg Focus 31:E18

    Article  Google Scholar 

  2. Neou M, Villa C, Armignacco R, Jouinot A, Raffin-Sanson ML, Septier A et al (2020) Pangenomic classification of pituitary neuroendocrine tumors. Cancer Cell 37:123-134.e5

    Article  CAS  Google Scholar 

  3. Drummond J, Roncaroli F, Grossman AB, Korbonits M (2019) Clinical and pathological aspects of silent pituitary adenomas. J Clin Endocrinol Metab 104:2473–2489

    Article  Google Scholar 

  4. Bonnet D, Dick JE (1997) Human acute myeloid leukemia is organized as a hierarchy that originates from a primitive hematopoietic cell. Nat Med 3:730–737

    Article  CAS  Google Scholar 

  5. Wang JC (2010) Good cells gone bad: the cellular origins of cancer. Trends Mol Med 16:145–151

    Article  CAS  Google Scholar 

  6. Abbaszadegan MR, Bagheri V, Razavi MS, Momtazi AA, Sahebkar A, Gholamin M (2017) Isolation, identification, and characterization of cancer stem cells: a review. J Cell Physiol 232:2008–2018

    Article  CAS  Google Scholar 

  7. Abdullah LN, Chow EK-H (2013) Mechanisms of chemoresistance in cancer stem cells. Clin Transl Med 2:3

    Article  Google Scholar 

  8. Zhang S, Cui Y, Ma X, Yong J, Yan L, Yang M et al (2020) Single-cell transcriptomics identifies divergent developmental lineage trajectories during human pituitary development. Nat Commun 11:5275

    Article  CAS  Google Scholar 

  9. Andoniadou CL, Matsushima D, Mousavy Gharavy SN, Signore M, Mackintosh AI, Schaeffer M et al (2013) Sox2(+) stem/progenitor cells in the adult mouse pituitary support organ homeostasis and have tumor-inducing potential. Cell Stem Cell 13:433–445

    Article  CAS  Google Scholar 

  10. Xu Q, Yuan X, Tunici P, Liu G, Fan X, Xu M et al (2009) Isolation of tumour stem-like cells from benign tumours. Br J Cancer 101:303–311

    Article  CAS  Google Scholar 

  11. Wurth R, Barbieri F, Pattarozzi A, Gaudenzi G, Gatto F, Fiaschi P et al (2017) Phenotypical and pharmacological characterization of stem-like cells in human pituitary adenomas. Mol Neurobiol 54:4879–4895

    Article  CAS  Google Scholar 

  12. Würth R, Thellung S, Corsaro A, Barbieri F, Florio T (2020) Experimental evidence and clinical implications of pituitary adenoma stem cells. Front Endocrinol (Lausanne) 11:54

    Article  Google Scholar 

  13. Knosp E, Steiner E, Kitz K, 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.

    Article  CAS  Google Scholar 

  14. Michaelis KA, Knox AJ, Xu M, Kiseljak-Vassiliades K, Edwards MG, Geraci M et al (2011) Identification of growth arrest and DNA-damage-inducible gene beta (GADD45beta) as a novel tumor suppressor in pituitary gonadotrope tumors. Endocrinology 152:3603–3613

    Article  CAS  Google Scholar 

  15. Newey PJ, Nesbit MA, Rimmer AJ, Head RA, Gorvin CM, Attar M et al (2013) Whole-exome sequencing studies of nonfunctioning pituitary adenomas. J Clin Endocrinol Metab 98:E796-800

    Article  CAS  Google Scholar 

  16. Pease M, Ling C, Mack WJ, Wang K, Zada G (2013) The role of epigenetic modification in tumorigenesis and progression of pituitary adenomas: a systematic review of the literature. PLoS ONE 8:e82619

    Article  Google Scholar 

  17. Kober P, Boresowicz J, Rusetska N, Maksymowicz M, Goryca K, Kunicki J et al (2018) DNA methylation profiling in nonfunctioning pituitary adenomas. Mol Cell Endocrinol 473:194–204

    Article  CAS  Google Scholar 

  18. Cheng S, Xie W, Miao Y, Guo J, Wang J, Li C et al (2019) Identification of key genes in invasive clinically non-functioning pituitary adenoma by integrating analysis of DNA methylation and mRNA expression profiles. J Transl Med 17:407

    Article  CAS  Google Scholar 

  19. Mertens F, Gremeaux L, Chen J, Fu Q, Willems C, Roose H et al (2015) Pituitary tumors contain a side population with tumor stem cell-associated characteristics. Endocr Relat Cancer 22:481–504

    Article  CAS  Google Scholar 

  20. Peverelli E, Giardino E, Treppiedi D, Meregalli M, Belicchi M, Vaira V et al (2017) Dopamine receptor type 2 (DRD2) and somatostatin receptor type 2 (SSTR2) agonists are effective in inhibiting proliferation of progenitor/stem-like cells isolated from nonfunctioning pituitary tumors. Int J Canc J Int du Cancer 140:1870–1880

    Article  CAS  Google Scholar 

  21. Capatina C, Cimpean AM, Raica M, Coculescu M, Poiana C (2019) SOX 2 expression in human pituitary adenomas-correlations with pituitary function. In vivo (Athens, Greece) 33:79–83

    CAS  Google Scholar 

  22. Yunoue S, Arita K, Kawano H, Uchida H, Tokimura H, Hirano H (2011) Identification of CD133+ cells in pituitary adenomas. Neuroendocrinology 94:302–312

    Article  CAS  Google Scholar 

  23. Chen L, Ye H, Wang X, Tang X, Mao Y, Zhao Y et al (2014) Evidence of brain tumor stem progenitor-like cells with low proliferative capacity in human benign pituitary adenoma. Cancer Lett 349:61–66

    Article  CAS  Google Scholar 

  24. Orciani M, Davis S, Appolloni G, Lazzarini R, Mattioli-Belmonte M, Ricciuti RA et al (2015) Isolation and characterization of progenitor mesenchymal cells in human pituitary tumors. Cancer Gene Ther 22:9–16

    Article  CAS  Google Scholar 

  25. Koyama-Nasu R, Nasu-Nishimura Y, Todo T, Ino Y, Saito N, Aburatani H et al (2013) The critical role of cyclin D2 in cell cycle progression and tumorigenicity of glioblastoma stem cells. Oncogene 32:3840–3845

    Article  CAS  Google Scholar 

  26. Peculis R, Mandrika I, Petrovska R, Dortane R, Megnis K, Nazarovs J et al (2020) Pituispheres contain genetic variants characteristic to pituitary adenoma tumor tissue. Front Endocrinol 11:313

    Article  Google Scholar 

  27. Mitra SK, Hanson DA, Schlaepfer DD (2005) Focal adhesion kinase: in command and control of cell motility. Nat Rev Mol Cell Biol 6:56–68

    Article  CAS  Google Scholar 

  28. Cakir M, Grossman AB (2009) Targeting MAPK (Ras/ERK) and PI3K/Akt pathways in pituitary tumorigenesis. Expert Opin Ther Targets 13:1121–1134

    Article  CAS  Google Scholar 

  29. Lytle NK, Barber AG, Reya T (2018) Stem cell fate in cancer growth, progression and therapy resistance. Nat Rev Cancer 18:669–680

    Article  CAS  Google Scholar 

  30. Shen X, Liu Q, Xu J, Wang Y (2021) Correlation between the expression of interleukin-6, STAT3, e-cadherin and n-cadherin protein and invasiveness in nonfunctional pituitary adenomas. J Neurol Surg B Skull Base 82:e59-69

    Article  Google Scholar 

  31. Mussunoor S, Murray GI (2008) The role of annexins in tumour development and progression. J Pathol 216:131–140

    Article  CAS  Google Scholar 

  32. Sharma MC (2019) Annexin A2 (ANX A2): an emerging biomarker and potential therapeutic target for aggressive cancers. Int J Canc J Int du Canc 144:2074–2081

    Article  CAS  Google Scholar 

Download references


We thank Dr. Qi Li and Dr. Wenyan He in helping the cell culture and molecular biological experiments in China National Clinical Research Center for Neurological Diseases, Beijing Tiantan Hospital,Capital Medical University.


This work was supported by Beijing Natural Science Foundation (7212008). Peng Kang is supported by Beijing Hospitals Authority Youth Programme (QML20180504) and Beijing Program of Excellent Young Investigator (2015000021469G221).

Author information

Authors and Affiliations



Conception of the work: PK, WJ; Acquisition of clinical samples: PK, LHY, JL, GJJ, WJ; Cell culture and molecular analysis: PK, LHY, PLL, LJ, JYP, JLZ; Bioinformatical analysis: PK, LHY; Data treatment and statistical analysis: PK, LHY, JL, SCM; Draft manuscript: PK, LHY, WJ; All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Wang Jia or Peng Kang.

Ethics declarations

Ethics approval and consent to participate

The procedure was in accordance with the ethical standards and approved by the ethics committee of Beijing Tiantan Hospital (KY202215501).

Consent for publication

Not available.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

. Table S1. Clinical demographics of gonadotrophic pituitary adenoma patients in hPASC culture and sequencing; Table S2. Clinical demographics of invasive and non-invasive gonadotrophic pituitary adenoma patients for Real-time PCR quantification; Table S3. PCR primers for verification of 9 candidate genes in invasive and non-invasive gonadotrophic pituitary adenoma.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Yuan, L., Li, P., Li, J. et al. Identification and gene expression profiling of human gonadotrophic pituitary adenoma stem cells. acta neuropathol commun 11, 24 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Gonadotrophic pituitary adenoma
  • Tumor stem cells
  • RNA-sequencing