- Open Access
Analysis of whole genome-wide microRNA transcriptome profiling in invasive pituitary adenomas and non-invasive pituitary adenomas
Chinese Neurosurgical Journal volume 5, Article number: 27 (2019)
Dysregulation of microRNAs (miRNAs) plays a critical role during the occurrence and progress of pituitary adenomas (PAs). However, the roles of miRNAs in the invasiveness of PA are poorly understood. This study aims to more comprehensively and specific define the relationship between altered miRNA and PA invasion.
The differential expression of miRNAs (DEMs) between invasive PAs (IPAs) and non-invasive PAs (NPAs) was explored by RNA sequencing and which functions were analyzed by gene ontology (GO) as well as Kyoto Encyclopedia of Genes and Genomes (KEGG). The miRNA-mRNA network was predicted with bioinformatics.
We identified 31 upregulated miRNAs and 24 downregulated miRNAs in IPAs compared with NPAs. GO analysis and KEGG pathway analysis showed the DEMs were mainly associated with cell proliferation and cell cycle pathway. In addition, on the count of predicted miRNA-mRNA network, two hub miRNAs were identified.
Our results demonstrate the miRNA-mRNA network in detail, which suggest that miRNA may be a promising target in diagnosis and therapy for IPAs.
PA is one of the most common intracranial tumors with an incidence of 10–15% . Generally, PAs are considered benign, but some of them are invasive, which invade the adjacent structures such as sphenoid sinus, cavernous sinus, and diaphragma sellae [2,3,4]. Therefore, IPAs are not only more difficult to complete surgical resection, but also more likely to recurrence after surgery.
MiRNAs are single-strand non-coding RNAs of approximately 19–23 nt, which regulate gene expression at the post-transcriptional level [5, 6], and they can also act as tumor suppressor genes or oncogenes in various tumors . For instance, miR-193b exerts tumor suppressive effects in human acute myeloid leukemia by inducing tumor cell apoptosis and G1/S arrest , while miR-210-3p plays an oncogene role in prostate cancer by promoting cancer cell epithelial–mesenchymal transition and bone metastasis via NF-κB signaling pathway . Altered expression of many miRNAs has been described in PAs, and specific miRNA signatures are related to clinical and therapeutic characteristics of the tumors . However, comprehensive and specific researches of relationships between miRNAs and invasiveness of PAs are still rare.
In order to better understand the mechanism of invasiveness in PAs, it is necessary to clarify the miRNA regulatory network in IPAs. In this study, we detected DEMs in IPAs and NPAs by RNA sequencing, and established the co-expression network contain miRNAs and predicted target genes by Cytoscape. In addition, the expression of the most upregulated miR-665 and the most downregulated miR-149-3p in IPAs was screen out. Moreover, we explained the potential functions of the two key miRNAs in invasive behavior of PAs by GO analysis and KEGG pathway analysis.
Patients and samples
Seven tumor samples were obtained from patients with PAs who underwent operation at the Department of Neurosurgery, 1st Affiliated Hospital of Kunming Medical University for identification of miRNAs by high-throughput sequencing. None of these patients has been received radiotherapy or chemotherapy before surgery. Tumor samples were divided into 2 groups according to invasive behavior proved by surgical findings and pathology: IPA and NPA. All patients were informed according to inform consent approved by the Ethics Boardof the 1st Affiliated Hospital of Kunming Medical University. Immediately following separation, the fresh tumor samples were placed in sterile, RNase-free 2.0-mL cryotubes. Then, samples were soaked in Trizol and stored at − 80 °C for following analysis.
RNA isolation, library preparation, and sequencing analysis
Total RNA was extracted from tissue samples by Trizol regent. The integrity of total RNA was detected by agarose electrophoresis and which was quantified by NanoDrop spectrophotometer. Then, the sequencing sample library was constructed by the following steps: ribosomal RNA removal, fragmentation, first-strand complementary DNA (cDNA) synthesis, second-strand cDNA synthesis, terminal repair 3′ terminal addition, ligation, and enrichment. The libraries were sequence on an Illumina Hiseq 2500/2000 platform.
MiRNA expression analysis
MiRNA expression levels were estimated by the TPM (transcript per million) through the following criteria: Normalized expression = mapped readcount/Total reads × 106 . All data were analyzed using the DESeq2 R package (1.8.3). log2FC > 1 and p < 0.05 were considered as the cutoff values for DEMs screening .
MiRNA-mRNA network construction
Based upon results of DEM analysis and target gene prediction, the miRNA-mRNA pairs were extracted to construct the miRNA-mRNA regulatory network. Then, the regulatory network was visualized using Cytoscape_v3.5.1.
Target gene prediction, gene ontology, and pathway enrichment analysis
Target genes of the DEMs were predicted using major online tools, including miRanda (http://miranda.org.uk/), PITA (http://genie.weizmann.ac.il/pubs/mir07/mir07_data.html), and RNAhybrid (https://bibiserv.cebitec.uni-bielefeld.de/rnahybrid/) . In order to analyze the main functions of the predicted target genes for the DEMs, we performed GO analysis . Moreover, KEGG  pathway enrichment analysis was used to find out the significant pathway of predicted target genes for the DEMs. A Go term or KEGG pathway with FDR < 0.05 was considered statistically significant. Top 10 enriched GO terms and pathways of DEMs were ranked by enrichment score (− log10(p value)).
Quality assessment of sequencing data
The results showed the patterns of gene expression among the samples were similar, and the Pearson correlation between samples is similar in the two groups (Fig. 1a, b). In view of the fact that the length of human miRNAs is generally 19–23 nt, the analysis of the length of reads of each sample shows that 19–23 nt accounts for a higher proportion of each of 2 groups (Fig. 1c, d). And reads of the length region were included in the analysis range. In addition, by comparing and analyzing the error rate among all samples, we found that the read error rate at 19–23 nt was far less than 0.05% (Fig. 1e, f). All above results indicate the sequencing data is reliable for further bioinformation analysis.
Apportionment and annotation of DEMs
After summarizing and classifying the sequence reads into different RNA categories, such as miRNA, sn/snoRNA, tRNA, and rRNAs, the pie chart is drawn to annotate and classify the total reads. The proportion of known miRNAs (NIN = 63.77%, INV = 56.92%) and newly discovered miRNAs (NIN = 0.04%, INV = 0.02%) can be obtained (Fig. 2a, b).
Next, Venn diagram of the DEMs in the two groups was drawn by FunRich3.1.3. A total of 136 DEMs were found in IPA, and 187 DEMs were found in NPA (Fig. 2c).
And then, the volcano map was used to infer the overall distribution of miRNAs. ▢Log2 FC > 1 and p < 0.05 were used as thresholds to screen the DEMs. A total of 31 significantly upregulated and 24 significantly downregulated miRNAs were screened out (Fig. 2d).
At last, these 55 DEMs were used to construct a hierarchical clustering analysis map (Fig. 2e).
GO and KEGG pathway analysis
In order to explore the functions of DEMs, target genes of these miRNAs were predicted by miRanda, PITA, and RNAhybrid. GO analysis and KEGG pathway analysis were used to their target pool. As results, the results of the top 10 enriched biological process (BP) terms showed that the target gene of DEMs was associated with cell proliferation, cell cycle, and apoptosis process. The top 10 KEGG pathways showed the DEMs might be involved in CDK4/6 signaling pathway, PI3K-Akt signaling pathway, and apoptosis pathways (Fig. 3).
Based on the DEMs and the prediction of target genes. The miRNA-mRNA network was generated by Cytoscape(v3.5.1). As we all know, hub nodes play important roles in biological networks. According to the degree of DEMs calculated by Cytoscape, a total of 12 miRNAs with higher values were identified (Table 1). Then the miRNA-mRNA subnetwork was generated, which included 5 upregulated miRNAs, 7 downregulated miRNAs, and 258 predicted target genes (Fig. 4).
Although PAs are classified into multiple subgroups based on the histological structure, pathological types, and hormone secretion [16, 17], the definition of clinically IPAs currently differs in the literature. In according to the 4th edition of the WHO classification of endocrine tumors published in 2017, the invasiveness of PAs can be evaluated through the tumor proliferative capacity by mitotic count and Ki-67 index , indicating that the tumor proliferation capacity is closely related to invasion of PAs.
MiRNAs, mediate post-transcriptional regulation, play an important role in epigenetic regulation. Its precursor is cleavage by Dicer enzyme, and then combined with AGO protein and other components to form RISC (RNA immuno-silencing complex) to play a key role in silencing or degradation of target mRNA [19, 20]. Some researchers have reported the regulating functions of miRNAs in PAs . Upregulated miR-34a can significantly inhibit the proliferation of PA cells GH4C1 and promote apoptosis by regulating SOX7 . Overexpression of miR-16 inhibits the proliferation of PA cells HP75 and promotes apoptosis by targeting HMGA2 expression . However, the comprehensive and specific effects of miRNAs in PA invasion behavior are still rarely reported . Thus, we recognized 55DEMs in IPAs by RNA sequencing analysis. Then, according to Cytoscape calculation, a miRNA-mRNA co-expression network including 5 upregulated miRNAs, 7 downregulated miRNAs, and 258 predicted target genes was generated in IPAs.
To further understand the potential function of miRNAs, GO and KEEG pathway analysis were applied for analyzing their possible biological functions and the signaling pathways in IPAs. The outcome of GO analysis showed that the target genes for the DEMs in IPAs were enriched for genes associated with cell proliferation, cell cycle, and apoptosis process, which is consistent with the active proliferation characteristics of IPA cells. Meanwhile, KEEG pathway analysis showed that the DEMs might be involved in CDK4/6 signaling pathway, PI3K-Akt signaling pathway, and apoptosis pathways. It is noted that the enriched terms of target genes for the DEMs from KEEG analysis were consistent with GO analysis results, and they are all closely related to cell proliferation. As mentioned above, tumor proliferation capacity is the most important feature of IPAs.
Additionally, the miRNA-mRNA network was constructed and the outcome showed that some hub miRNAs play important roles in IPAs. The hub miRNAs were identified by the degree of DEMs calculated by Cytoscape. Intriguingly, we found that these miRNAs have similar biological functions. For instance, hsa-miR-665, the most upregulated node in IPAs, can promote tumor cell proliferation or cell cycle progression in hepatocellular carcinoma . In contrast, it was found that overexpression of has-miR-149-3p, the most downregulated node in IPAs, can inhibit the proliferation and invasion of tumor cells in bladder cancer and renal epithelial cell carcinoma .
In this study, we identified the DEMs related to the invasive behavior of PAs. Further study is needed to confirm the exact relationship between DEMs and invasive behavior of PAs and clarify the molecular mechanism of miRNA affecting the invasion of PAs.
Availability of data and materials
Please contact author for data requests.
Differential expression miRNA
Invasive pituitary adenoma
Encyclopedia of Genes and Genomes
Non-invasive pituitary adenoma
Aflorei ED. Epidemiology and etiopathogenesis of pituitary adenomas. J Neuro-Oncol. 2014;117(3):379–94.
Knosp E, Steiner E, Kitz K. Pituitary adenomas with invasion of the cavernous sinus space: a magnetic resonance imaging classification compared with surgical findings. Neurosurgery. 1993,33(4):610–617; discussion 617-8.
Farrell WE. Molecular pathogenesis of pituitary tumors. Front Neuroendocrinol. 2000;21(3):174–98.
Yang Q. Molecular network basis of invasive pituitary adenoma: a review. Front Endocrinol. 2019;10:7.
Berindan-Neagoe I, Monroig Pdel C, Pasculli B. MicroRNAome genome: a treasure for cancer diagnosis and therapy. CA Cancer J Clin. 2014;64(5):311–36.
Treiber T, Treiber N. Regulation of microRNA biogenesis and its crosstalk with other cellular pathways. Nat Rev Mol Cell Biol. 2019;20(1):5–20.
Rupaimoole R, Slack FJ. MicroRNA therapeutics: towards a new era for the management of cancer and other diseases. Nat Rev Drug Discov. 2017;16(3):1–19.
Bhayadia R, Krowiorz K, Haetscher N, et al. Endogenous tumor suppressor microRNA-193b: therapeutic and prognostic value in acute myeloid leukemia. J Clin Oncol. 2018;36(10):1007–17.
Ren D, Yang Q, Dai Y, et al. Oncogenic miR-210-3p promotes prostate cancer cell EMT and bone metastasis via NF-κB signaling pathway. Mol Cancer. 2017;16(1):117.
Gentilin E, Degli Uberti E, Zatelli MC. Strategies to use microRNAs as therapeutic targets. Best Pract Res Clin Endocrinol Metab. 2016;30(5):629–39.
Zhou L, Chen J, Li Z, et al. Integrated profiling of microRNAs and mRNAs: microRNAs located on Xq273 associate with clear cell renal cell carcinoma. PloS One. 2010;5(12):e15224.
Lai Y. A statistical method for the conservative adjustment of false discovery rate (q-value). BMC Bioinformatics. 2017;18(Suppl 3):69.
Wu HJ, Ma YK, Chen T, et al. PsRobot: a web-based plant small RNA meta-analysis toolbox. Nuc Acids Res. 2012;40(Web server issue):W22–8.
Chen CM, Lu YL, Sio CP, et al. Gene ontology based housekeeping gene selection for RNA-seq normalization. Methods (San Diego, Calif.), 2014, 67(3): 354–363.
Mao X, Cai T, Olyarchuk JG. Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics (Oxford, England). 2005;21(19):3787–93.
Micko AS, Wöhrer A, Wolfsberger S. Invasion of the cavernous sinus space in pituitary adenomas: endoscopic verification and its correlation with an MRI-based classification. J Neurosurg. 2015;122(4):803–11.
Kalinin PL, Sharipov OI, Pronin IN, et al. Endoscopic transsphenoidal resection of pituitary adenomas invading the cavernous sinus. Zh Vopr Neirokhir Im N N Burdenko. 2016;80(4):63–74.
Lopes MBS. The 2017 World Health Organization classification of tumors of the pituitary gland: a summary. Acta Neuropathol. 2017;134(4):521–35.
Shukla GC, Singh J. MicroRNAs: processing, maturation, target recognition and regulatory functions. Mol Cell Pharmacol. 2011;3(3):83–92.
Meister G, Landthaler M, Patkaniowska A, et al. Human Argonaute2 mediates RNA cleavage targeted by miRNAs and siRNAs. Mol Cell. 2004;15(2):185–97.
Németh K, Darvasi O, Szücs N, et al. The role of miRNAs in the pathogenesis of pituitary adenomas. Orvosi Hetilap. 2018;159(7):252–9.
Yang Z, Zhang T, Wang Q. Overexpression of microRNA-34a attenuates proliferation and induces apoptosis in pituitary adenoma cells via SOX7. Mol Ther Oncol. 2018;10:40–7.
Niu Y, Zhou H, Liu Y, et al. miR-16 regulates proliferation and apoptosis of pituitary adenoma cells by inhibiting HMGA2. Oncol Lett. 2019;17(2):2491–7.
Feng Y, Mao ZG, Wang X, et al. MicroRNAs and target genes in pituitary adenomas. Horm Metab Res. 2018;50(3):179–92.
Hu Y, Yang C, Yang S, et al. miR-665 promotes hepatocellular carcinoma cell migration, invasion, and proliferation by decreasing Hippo signaling through targeting PTPRB. Cell Death Dis. 2018;9(10):954.
Zhan Y, Chen Z, Li Y, et al. Long non-coding RNA DANCR promotes malignant phenotypes of bladder cancer cells by modulating the miR-149/MSI2 axis as a ceRNA. J Exp Clin Cancer Res. 2018;37(1):273.
This study was supported by the Applied Basic Research Foundation of Yunnan Province, the Research Innovation Team of Yunnan Province, the Medical Academic Leader Training Program of Yunnan Province, and the Outstanding Youth Scientists and Technicians Training Program of Yunnan Province.
Ethics approval and consent to participate
The study was approved by the Ethics Board of the 1st Affiliated Hospital of Kunming Medical University (NO. 2018NS01).
Consent for publication
The authors declare that they have no competing interests.
About this article
Cite this article
Zhang, C., Qian, Y., Qiao, Y. et al. Analysis of whole genome-wide microRNA transcriptome profiling in invasive pituitary adenomas and non-invasive pituitary adenomas. Chin Neurosurg Jl 5, 27 (2019). https://doi.org/10.1186/s41016-019-0177-4
- Pituitary adenoma