Identification of Functional Genes in Pterygium Based on Bioinformatics Analysis

Xu, Yuting; Qiao, Chen; He, Siying; Lu, Chen; Dong, Shiqi; Wu, Xiying; Yan, Ming; Zheng, Fang

doi:https://doi.org/10.1155/2020/2383516

BioMed Research International

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusion Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2020 | Article ID 2383516 | https://doi.org/10.1155/2020/2383516

Identification of Functional Genes in Pterygium Based on Bioinformatics Analysis

Yuting Xu,¹Chen Qiao,²Siying He,³Chen Lu,⁴Shiqi Dong,¹Xiying Wu,¹Ming Yan,¹and Fang Zheng²

Academic Editor: Mitsuru Nakazawa

Received20 Jul 2020

Revised12 Oct 2020

Accepted24 Oct 2020

Published23 Nov 2020

Abstract

Purpose. The competing endogenous RNA (ceRNA) network regulatory has been investigated in the occurrence and development of many diseases. This research aimed at identifying the key RNAs of ceRNA network in pterygium and exploring the underlying molecular mechanism. Methods. Differentially expressed long noncoding RNAs (lncRNAs), microRNAs (miRNAs), and mRNAs were obtained from the Gene Expression Omnibus (GEO) database and analyzed with the R programming language. LncRNA and miRNA expressions were extracted and pooled by the GEO database and compared with those in published literature. The lncRNA-miRNA-mRNA network was constructed of selected lncRNAs, miRNAs, and mRNAs. Metascape was used to perform Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses on mRNAs of the ceRNA network and to perform Protein-Protein Interaction (PPI) Network analysis on the String website to find candidate hub genes. The Comparative Toxicogenomic Database (CTD) was used to find hub genes closely related to pterygium. The differential expressions of hub genes were verified using the reverse transcription-real-time fluorescent quantitative PCR (RT-qPCR). Result. There were 8 lncRNAs, 12 miRNAs, and 94 mRNAs filtered to construct the primary ceRNA network. A key lncRNA LIN00472 ranking the top 1 node degree was selected to reconstruct the LIN00472 network. The GO and KEGG pathway enrichment showed the mRNAs in ceRNA networks mainly involved in homophilic cell adhesion via plasma membrane adhesion molecules, developmental growth, regulation of neuron projection development, cell maturation, synapse assembly, central nervous system neuron differentiation, and PID FOXM1 PATHWAY. According to the Protein-Protein Interaction Network (PPI) analysis on mRNAs in LINC00472 network, 10 candidate hub genes were identified according to node degree ranking. Using the CTD database, we identified 8 hub genes closely related to pterygium; RT-qPCR verified 6 of them were highly expressed in pterygium. Conclusion. Our research found LINC00472 might regulate 8 hub miRNAs (miR-29b-3p, miR-183-5p, miR-138-5p, miR-211-5p, miR-221-3p, miR-218-5p, miR-642a-5p, miR-5000-3p) and 6 hub genes (CDH2, MYC, CCNB1, RELN, ERBB4, RB1) in the ceRNA network through mainly PID FOXM1 PATHWAY and play an important role in the development of pterygium.

1. Introduction

Pterygium is a common ocular surface disease characterized by a triangular-shaped growth consisting of fibrotic subconjunctival connective tissue and hypertrophy of the overlying conjunctival epithelium [1]. The population affected by pterygium reached 200 million globally, and the prevalence in China was 108.65 million [2, 3]. Pterygium has long been considered as a chronic degenerative condition; however, after abnormal expression of the p53 protein was found in the epithelium [4], pterygium is now considered ultraviolet-related uncontrolled cell proliferation, similar to tumors [5]. The development of pterygium is a complicated process involving cell proliferation, migration, inflammatory infiltrates, fibrosis, angiogenesis, and extracellular matrix breakdown [2]. However, the mechanism of pterygium formation and development is still not completely understood. Studies have demonstrated that noncoding RNAs served crucial roles in numerous diseases [6–8]. They are classified into two classes: small noncoding RNAs (microRNAs (miRNAs), small interfering RNAs, and transfer RNAs) and long noncoding RNAs (lncRNAs). Salmena et al. proposed the competing endogenous RNA (ceRNA) hypothesis, wherein lncRNAs harboring miRNA response elements competed with one another to bind to a common miRNA and thereby acted as molecular “sponges” and depressed the target genes of the miRNAs [9]. The ceRNA network has been demonstrated in numerous diseases, particularly in cancer. By establishing the lncRNA-miRNA-mRNA network, Wang et al. identified functional genes in heart failure and provided further insights into the important roles of the ceRNA network in heart failure [10]. Yang et al. discovered that lncRNA RP11-169F17.1 and RP11-669N7.2 could become novel biomarkers of stomach adenocarcinoma assessed by the construction of the ceRNA network [11]. Zhu et al. reported that 4 lncRNAs might act as potential therapeutic targets or candidate prognostic biomarkers in clear cell kidney carcinoma by reconstruction and comprehensive analysis of the ceRNA regulatory network [12].

However, there is only 1 research on constructing ceRNA networks for pterygium [13]. The hub genes of the network and their roles in the development of pterygium have not been fully understood. Therefore, we collected the lncRNA, miRNA, and mRNA datasets and published literature; filtered out differentially expressed long noncoding RNAs (DELs), microRNAs (DEMis), and mRNAs (DEMs) to constructed ceRNA networks; and identified hub genes closely related to pterygium, based on the functional and pathway analysis on mRNAs in the ceRNA network.

2. Materials and Methods

2.1. Data Collection

The Gene Expression Omnibus (GEO) is an international, public functional genomics data repository for high-throughput microarray and next-generation sequences [14]. The published pterygium gene expression profiles (GSE83627, GSE21346, GSE51995, and GSE2513) were downloaded [15–19]. The GSE83627 lncRNA data set was collected using GPL14550 platforms (Agilent-028004 SurePrint G3 Human GE 8 × 60 K Microarray, Agilent Technologies, Inc.) and included samples from 4 pterygium and 4 healthy conjunctiva tissues. The miRNA expression data of GSE21346 were based on the GPL7723 platforms (miRCURYLNA microRNA Array, v.11.0-hsa, mmu & rno, Exiqon A/S) and consisted of samples from 3 pterygium and 3 healthy conjunctiva tissues. The mRNA expression data of GSE51995 were from 4 pterygium samples and 4 healthy conjunctiva tissues and were based on GPL14550 platforms (Agilent-028004 SurePrint G3 Human GE 8 × 60 K Microarray, Agilent Technologies, Inc.). Another mRNA expression data of GSE2513 were from 8 pterygium samples and 4 healthy conjunctiva tissues and were based on GPL96 platforms (Affymetrix Human Genome U133A Array, Affymetrix). In order to make the datasets more abundant, we executed a systemic search in PubMed, CNKI, and Web of Science and found 2 microarray studies on DELs and 7 studies on DEMis, including 4 microarray assays studies and 3 experimental verification studies.

2.2. Data Preprocessing

All downloaded data sets were standardized through the limma package of R 3.6.3, and the standardized expression matrix was used for the differential expression analysis. The DELs, DEMis, and DEMs between pterygium samples and healthy control conjunctiva tissues were selected according to and values <0.05 [20]. Sequentially, DELs and DEMis were pooled based on the microarray data in the searched studies, respectively, with the Draw Venn Diagram website (http://bioinformatics.psb.ugent.be/webtools/Venn/). And the pooled results were used as candidate DELs and DEMis for the following analysis. The volcano plots and the heatmaps of DELs, DEMis, and DEMs expressions were set up using the ggpolt 2 package and pheatmap package in R 3.6.3 separately.

2.3. Prediction of Target lncRNAs and mRNAs of miRNAs

The lncRNA targets of the candidate miRNAs were predicted using DIANA-LncBase v3 (https://diana.e-ce.uth.gr/lncbasev3/interactions), which identifies miRNA and lncRNA interactions derived from manually curated publications and the analysis of 153 AGO CLIP-Seq libraries [21], and starbase 2.0 (http://starbase.sysu.edu.cn/) [22]. The predicted DELs could be selected for ceRNA network construction. Data for miRNA and mRNA interactions were downloaded from miRDB (http://www.mirdb.org/), TargetScanHuman7.1 (http://www.targetscan.org/), and miWalk2.0 (http://zmf.umm.uni-heidelberg.de/apps/zmf/mirwalk2/miRretsys-self.html), and DEMs which appeared in all three databases were used for the following analysis.

2.4. Construction of the lncRNA-miRNA-mRNA Network

Based on the associations between lncRNAs-miRNAs and miRNAs-mRNAs, the ceRNA network was constructed. First, the upregulated and downregulated RNAs (lncRNAs, miRNAs, and mRNAs) were assigned with values <0.05; then, the lncRNA-miRNA-mRNA network was visualized by using the Cytoscape 3.7.2 software [23], and all node degrees of the ceRNA network were calculated using the CytoHubba plugin [24] of Cytoscape 3.7.2. The significant lncRNA was selected according to the node degree ranking, and a new ceRNA network was reconstructed centered on the top lncRNA.

2.5. Functional Analysis on mRNAs in ceRNA Network for Discovering Hub Genes

To analyze the function of mRNAs in ceRNA networks, GO and KEGG analyses were performed by Metascape (http://metascape.org/) [25]. The protein-protein interaction (PPI) network of mRNAs in the ceRNA network was constructed by the STRING database with the confidence as the cutoff level, providing the critical assessment and integration of protein-protein interactions of genes, including direct (physical) as well as indirect (function) associations [26, 27]. Then, we used the CytoHubba plugin of Cytoscape 3.7.2 to analyze the PPI network and to choose the candidate hub genes with top node degrees. Sequentially, we identified hub genes closely related to pterygium using the Comparative Toxicogenomic Database (CTD), which provides interactions between diseases-genes and diseases-drugs.

2.6. Validation of Hub Genes

2.6.1. Patient Sample Preparation

There were 20 primary pterygium specimens obtained in Zhongnan Hospital of Wuhan University during pterygium surgery. The control conjunctival tissues were derived from healthy conjunctival tissues on the temporal side of the surgical eye of the same patient, with a size of . There were 9 males aged 40-84 () years old and 11 females aged 40-79 () years old. All participants are of Han nationality and have no blood relationship. This study was reviewed and approved by the Ethics Committee of Zhongnan Hospital of Wuhan University, and the written informed consent of the patients has been obtained before the study. Patients with recurrent pterygium, history of severe eye trauma, eye infection, cataract, glaucoma, acute dacryocystitis, and severe mental illness were excluded.

2.6.2. Tissue RNA Extraction and RT-qPCR

All specimens were collected during the operation and immediately placed in RNA protective agent (RNA later™, SIGMA, USA), stored at -20°C for usage. RNA was extracted from 80 mg specimen tissue using TRIZOL reagent. The relative expression of hub genes was measured using a reverse kit (HiScript® III RT SuperMix for qPCR (+gDNA wiper)) and ChamQ™ Universal SYBR® qPCR Master Mix (Vazyme company, Nanjing, China) on CFX96 fluorescence qPCR instrument (Bio-Rad, USA). And the PCR reaction was carried as follows: predenaturation at 95°C for 30 seconds, denaturation at 95°C for 3-10 seconds, and annealing extension at 60°C for 10-30 seconds, for 40 cycles. Using GAPDH as internal references, the relative expression of genes was calculated by the 2^-△Ct method. The primer sequences are listed in Supplementary Table 1.

3. Results

3.1. Filtered DELs, DEMis, and DEMs

There are four recruited pterygium datasets from the GEO search, including the lncRNA data set (GSE83627), miRNA data set (GSE21346), and mRNA data sets (GSE51995 and GSE2513). In Figure 1, the expression of DELs, DEMis, and DEMs is shown by volcano plots and heat maps (GSE83627 (Figure 1(e)), GSE21346 (Figure 1(f)), GSE51995 (Figure 1(g)), GSE2513 (Figure 1(h)). After the limma package differential analysis, 214 upregulated lncRNAs were selected (Figure 1(a)), 4 upregulated miRNAs (miR-138-5p, miR-184, miR-642a-5p, miR-1298-5p) were found (Figure 1(b)), and there are 360 downregulated and 537 upregulated mRNAs in GSE51995 (Figure 1(c)) and GSE2513 (Figure 1(d)) datasets. Using online search keywords (lncRNA and pterygium) in PubMed, CNKI, and Web of Science, we found two publications about lncRNA microarray data of pterygium [28, 29]. There were 10 upregulated lncRNAs and 10 downregulated lncRNAs in Liu’s study and 10 upregulated lncRNAs and 7 downregulated lncRNAs in Zheng’s research (Supplementary Table 2). Then, we used the Venn diagram to make an intersection analysis with lncRNAs in the two publications and previously selected DELs from datasets (Figure 2(a)). So, we merged the lncRNAs from two studies with DELs; then, finally, 251 lncRNAs were obtained to be DELs in the following study, of which 17 are downregulated lncRNAs and 234 upregulated lncRNAs. Moreover, there were seven references including 4 microarray articles and 3 experimental validation (Supplementary Table 2). In the 4 microarray articles, Engelsvold et al. [30] found 14 upregulated and 12 downregulated miRNAs; Cui et al. [31] found 10 differently expressed miRNAs, including 1 upregulated and 9 downregulated; Lee et al. [32] found 4 upregulated miRNAs (miR-143-3p, miR-181a-2-3p, miR-377-5p, miR-411-5p); and Lan et al. found 2 upregulated miRNAs in pterygium (miR-138-5p, miRPlus-E1233) [33]. Three experimental validation studies confirmed that miR-218-5p was downregulated, and miR-182-5p, miR-183-5p, miR-184, and miR-221-3p were upregulated [34–36]. The Venn diagram was used to remove the duplicate miRNAs and combine the reported miRNA with 4 DEMis filtered in GSE21346 (Figure 2(b)). Then, totally, 46 differentially expressed miRNAs were obtained, among which 21 miRNAs were downregulated and 25 miRNAs were upregulated. Eventually, a total of 360 downregulated mRNAs and 537 upregulated mRNAs of two datasets in GEO were filtered. In a word, we got 251 lncRNAs, 46 miRNAs, and 897 mRNAs as candidate RNAs for subsequent study.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 1

The differentially expressional pattern of RNAs from pterygium and control conjunctiva. (a–d) Volcano plots of DELs, DEMis, and DEMs: (a) DELs (GSE83627); (b) DEMis (GSE21346); (c) DEMs (GSE51995); (d) DEMs (GSE2513). Each point indicates an RNA. Green dots denote the downregulated RNAs, red points represent the upregulated RNAs under the same thresholds, and gray points indicate RNAs that did not change significantly. The criteria are based on and values <0.05. (e–h) Heat maps of DERNAs (e) DELs (GSE83627), (f) DEMis (GSE21346), (g) DEMs (GSE51995), and (h) DEMs (GSE2513). Each column represents one sample, and each row refers to an RNA. The color legend is on the top-right of the figure. Green indicates RNAs with a lower expression relative to the geometrical means; red indicates genes with a greater expression relative to the geometrical means. DELs represent differentially expressed long noncoding RNAs; DEMis represent differentially expressed microRNAs; DEMs represent differentially expressed mRNAs.

(a)

(b)

3.2. miRNA Predicted lncRNA and mRNA Targets

Firstly, the associations between the DELs and DEMis were assessed to reveal the 12 miRNAs targeting 8 lncRNAs. Subsequently, 94 mRNAs targeted by 12 miRNAs were identified. Finally, the 8 coexpressed lncRNAs, 12 coexpressed miRNAs, and 94 coexpressed mRNAs were selected to construct the primary ceRNA network (Table 1).

3.3. lncRNA-miRNA-mRNA Networks

The lncRNA-miRNA-mRNA ceRNA network was constructed by Cytoscape 3.7.2 (Figure 3(a)), which consisted of above 8 lncRNA nodes, 12 miRNA nodes, 94 mRNA nodes, and 136 edges. According to the node degree ranking, we selected the highest node degree RNA (LINC00472) as the center to construct the LINC00472 network (Figure 3(b)). The new ceRNA network centered by LINC00472 was constructed of 1 lncRNA nodes, 8 miRNA nodes, 84 mRNA nodes, and 106 edges (Figure 4(a)). This suggested that LINC00472 might be a key lncRNA in pterygium functioning through reacting with the 8 miRNAs (miR-29b-3p, miR-183-5p, miR-138-5p, miR-211-5p, miR-221-3p, miR-218-5p, miR-642a-5p, miR-5000-3p). Among them, miRNA-221-3p [34], miR-183-5p [35], and miR-218-5p [36] have been reported to have important function in pterygium. The other 5 miRNAs have been reported to play an important role in cancers. For example, Drummond et al. found miR-29b-3p was involved in a novel mechanism wherein cigarette smoke promoted accelerated cardiac and renal tissue injury in chronic kidney diseases [37]. And miR-138-5p can inhibit the malignant progression of prostate cancer [38] and lung cancer growth, through the miR-138-5p/FOXC1 pathway [39]. miR-211-5p has been confirmed to be related to a variety of cancers, such as cervical cancer [40], the nonsmall cell lung cancer cell [41], breast cancer [42], oral squamous cell carcinoma [43], and papillary thyroid cancer [44]. miR-642a-5p and miR-5000-3p were found have important effects on colon cancer [45, 46]. As we know, pterygium is considered as uncontrolled cell proliferation and is similar to tumor [5], which supports that these miRNAs could also play an important role in the development of pterygium.

(a)

(b)

(c)

Figure 3

Construction of the primary lncRNA-miRNA-mRNA network and enrichment analysis of the network. (a) The overview of the whole ceRNA network. This interaction network contains 8 lncRNAs, 12 miRNAs, and 94 mRNAs. The regulatory relationship among lncRNA, miRNA, and mRNA was visualized using Cytoscape 3.7.2. (b) The top 20 mRNA node degrees of the primary ceRNA network visualize with a bar chart colored. (c) The top 19 enriched GO and KEGG pathway terms visualize with a bar chart colored. Purple denoted lncRNAs, red denoted miRNAs, and green denoted mRNAs. Triangles denoted upregulated RNAs, and polygons denoted downregulated RNAs.

(a)

(b)

3.4. Functional Analysis of mRNAs Related to Pterygium

The 94 mRNAs in the primary ceRNA network was investigated by the GO and KEGG pathway enrichment analysis with Metascape. Metascape results were dominated by functional categories, including regulation of homophilic cell adhesion via plasma membrane adhesion molecules, developmental growth, regulation of neuron projection development, adherens junction interactions, cell maturation, synapse assembly, central nervous system neuron differentiation, positive regulation of mRNA metabolic process, negative regulation of cell-substrate adhesion, tissue morphogenesis, and blood vessel morphogenesis. Two pathways were enriched, PID FOXM1 PATHWAY and PID ATF2 PATHWAY. The occurrence of pterygium is closely related to the epithelial-mesenchymal transitions (EMT) and angiogenesis under the tissue [47]. FOXM1 is a member of the forkhead box (FOX) transcription factor family. Previous reports have demonstrated that FOXM1 is oncogenic, in breast cancer [48], hepatocellular carcinoma [49], prostate cancer [50], lung cancer [51], and colorectal cancer [52], and plays important roles in angiogenesis, invasion, and metastasis [53]. It may indicate the mRNAs in the ceRNA network are the biogenesis of the transformation of conjunctival cells into fibroblasts and the formation of subconjunctival neovascularization according to these biological processes and pathways, especially by regulating cell maturation, tissue morphogenesis, blood vessel morphogenesis, and PID FOXM1 PATHWAY (Figure 3(c)).

In order to explore the role of LINC00472 in pterygium, we performed GO and KEGG enrichment analysis again on the 84 mRNAs in the LINC00472 ceRNA network. We found that these genes were enriched in the same biology as the primary ceRNA network, but also in the cell part morphogenesis, regulation of synaptic transmission, glutamatergic, cellular response to lipid, negative regulation of cell differentiation, nervous system development, temperature homeostasis, NABA CORE MATRISOME, regulation of behavior, and response to growth factor. In addition, the mRNAs in both networks are enriched in PID FOXM1 PATHWAY, indicating that LINC00472 is likely to regulate the growth of pterygium through PID FOXM1 PATHWAY (Figure 4(b)).

3.5. PPI Network Construction and Analysis of Hub Genes

The PPI network of mRNAs in the LINC00472 network was obtained by using the STRING database, including 36 nodes and 48 edges, and the nodes contain 24 downregulated genes and 12 upregulated genes (Figure 5(a)). The top 10 genes evaluated by connectivity degree in the PPI network were present in CytoHubba, and the result showed that CDH2 was the most significant genes with a connectivity degree of 10, followed by MYC, CCNB1, PCDHA10, RELN, CCNA2, ERBB4, CDH11, PCDHA, and RB1 (Figure 5(b)). We defined the 10 genes as candidate hub genes in the PPI network. It means LINC00472 may regulate pterygium through these 10 key genes.

(a)

(b)

3.6. Identified Hub Genes Closely Related to Pterygium

We searched for genes related to pterygium in the CTD and found 6163 mRNAs, these genes are associated with pterygium or its descendants. A gene has either a curated association to the disease (M marker/mechanism and/or T therapeutic) or an inferred association via a curated chemical interaction. Then, we found 8 overlapped genes in the database, CCNB1, CDH11, MYC, CCNA2, ERBB4, RELN, RB1, and CDH2 (Figure 6), as selected hub genes. It implied that the 8 genes in the LINC00472 network may play a key role in pterygium.

3.7. Hub Genes Were Verified

Through RT-qPCR, we used 20 groups of pterygium and control conjunctival tissues to verify the above 8 hub genes. There are 6 genes that are highly expressed in pterygium, MYC, CDH2, CCNB1, RELN, RB1, and ERBB4. Compared with the control conjunctiva, only the expression of CDH11 and CCNA2 showed no significant differences with pterygium (Figure 7). The results indicate that the prediction of the hub genes is reasonably reliable and could be used to discover more hub genes in pterygium.

4. Discussion

Pterygium is a common disease in ophthalmology, and its occurrence is related to the proliferation of subconjunctival fibrous tissue and the formation of new blood vessels [15]. The treatment of pterygium is mainly surgical treatment, but the high postoperative recurrence rate is still a problem that is currently difficult to solve [54]. Therefore, effective drug therapy combined with the operation to reduce the recurrence rate of pterygium has become a research hot spot. Recently, with the development of high-throughput genomic platforms and bioinformatic analysis, more and more noncoding RNAs, especially miRNAs and lncRNAs, have been discovered. They serve an important role in the regulation of multiple biological processes, including the development of pterygium [55, 56].

In this study, we mined pterygium-associated RNA expression datasets from the GEO to construct the ceRNA network. Through the analysis of the 4 data sets and 9 publications about lncRNAs and miRNAs, we obtained 251 DELs, 46 DEMis, and 897 DEMs. By the link of miRNA, the primary ceRNA network of pterygium was constructed with 8 lncRNA nodes, 12 miRNA nodes, 94 mRNA nodes, and 136 edges. Since LINC00472 has the highest node degree value of all lncRNAs in the ceRNA network (Figure 3(b)), we took LINC00472 as the center and built the LINC00472 network. And we found genes in the primary ceRNA network mainly enriched in homophilic cell adhesion via plasma membrane adhesion molecules, developmental growth, regulation of neuron projection development, adherens junctions interactions, cell maturation, synapse assembly, central nervous system neuron differentiation, positive regulation of mRNA metabolic process, negative regulation of cell-substrate adhesion, tissue morphogenesis, blood vessel morphogenesis, and PID FOXM1 PATHWAY, while genes in the LINC00472 network enriched in homophilic cell adhesion via plasma membrane adhesion molecules, developmental growth, cell part morphogenesis, regulation of synaptic transmission, glutamatergic, cellular response to lipid, negative regulation of cell differentiation, nervous system development, and temperature homeostasis. In addition, the genes in both networks enriched in PID FOXM1 PAYHWAY, which represents cell proliferation [57], angiogenesis, invasion, and metastasis [53]. This indicated that LINC00472 might play an important role in pterygium through FOXM1 PATHWAY. After the PPI network analysis on mRNAs in the LINC00472 network was visualized by Cytoscape 3.7.2, 10 candidate hub genes were found out. We searched the 10 hub genes in the CTD; 8 hub genes were identified to be involved in pterygium. They are CCNB1, CDH11, MYC, CCNA2, ERBB4, RELN, RB1, and CDH2, and 6 of them were confirmed highly expressed in pterygium by wet experiments (CCNB1, MYC, ERBB4, RELN, RB1, and CDH2). We presume that the LINC00472 network might function in pterygium via the FOXM1 PATHWAY, especially through the regulation of CCNB1, MYC, ERBB4, RELN, RB1, and CDH2 in the ceRNA network. And the data mining based on the public data and publications in the present study is useful and reliable for digging emerging therapy targets of pterygium.

5. Conclusion

The present study increased the understanding of the ceRNA-associated regulatory mechanism in the pterygium and identified LINC00472 as a key lncRNA in the whole ceRNA network, which might be highly related to pterygium via the FOXM1 PATHWAY in the LINC00472 network, especially by regulating its related 6 hub genes (CCNB1, MYC, ERBB4, RELN, RB1, and CDH2).

Data Availability

The microarray data (GSE83627, GSE21346, GSE51995, GSE2513) are obtained from the GEO database.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Authors’ Contributions

M. Yan designed the experiments. F. Zheng, S. He, C. Lu, C. Qiao, and Y. Xu wrote the manuscript. M. Yan and F. Zheng conducted the expression analysis. Y. Xu, C. Qiao, C. Lu, and S. He carried out the study and collected important background information. S. Dong and X. Wu helped to do some supplements. All authors read and approved the final manuscript. Yuting Xu and Chen Qiao contributed equally to this work.

Acknowledgments

This work was supported by the Grant of Nation Natural Science Foundation of China (grant no. 81770898) and the research project of Wuhan Health and Family Planning Commission (grant no. WX18Q46).

Supplementary Materials

Supplementary Table 1: Primers used for RT-qPCR of hub genes. Supplementary Table 2: Differentially expressed lncRNA and miRNAs in pterygium from 9 publications [28–36]. (Supplementary Materials)

References

J. Chui, N. D. Girolamo, D. Wakefield, and M. T. Coroneo, “The pathogenesis of pterygium: current concepts and their therapeutic implications,” The Ocular Surface, vol. 6, no. 1, pp. 24–43, 2008.
View at: Publisher Site | Google Scholar
M. Fuest, Y. C. Liu, M. T. Coroneo, and J. S. Mehta, “Femtosecond laser assisted pterygium surgery,” Cornea, vol. 36, no. 7, pp. 889–892, 2017.
View at: Publisher Site | Google Scholar
P. Song, X. Chang, M. Wang, L. An, and C. Pan, “Variations of pterygium prevalence by age, gender and geographic characteristics in China: A systematic review and meta-analysis,” PLoS One, vol. 12, no. 3, p. e0174587, 2017.
View at: Publisher Site | Google Scholar
K. W. Kim, S. H. Park, and J. C. Kim, “Fibroblast biology in pterygia,” Experimental Eye Research, vol. 142, pp. 32–39, 2016.
View at: Publisher Site | Google Scholar
D. T. Tan, W. Y. Tang, Y. P. Liu, H. S. Goh, and D. R. Smith, “Apoptosis and apoptosis related gene expression in normal conjunctiva and pterygium,” The British Journal of Ophthalmology, vol. 84, no. 2, pp. 212–216, 2000.
View at: Publisher Site | Google Scholar
S. Xu, D. Kamato, P. J. Little, S. Nakagawa, J. Pelisek, and Z. G. Jin, “Targeting epigenetics and non-coding RNAs in atherosclerosis: from mechanisms to therapeutics,” Pharmacology & Therapeutics, vol. 196, pp. 15–43, 2019.
View at: Publisher Site | Google Scholar
T. Sallam, J. Sandhu, and P. Tontonoz, “Long noncoding RNA discovery in cardiovascular disease: decoding form to function,” Circulation Research, vol. 122, no. 1, pp. 155–166, 2018.
View at: Publisher Site | Google Scholar
C. M. Wong, F. H. C. Tsang, and I. O. L. Ng, “Non-coding RNAs in hepatocellular carcinoma: molecular functions and pathological implications,” Nature Reviews Gastroenterology & Hepatology, vol. 15, no. 3, pp. 137–151, 2018.
View at: Publisher Site | Google Scholar
L. Salmena, L. Poliseno, Y. Tay, L. Kats, and P. P. Pandolfi, “A ceRNA hypothesis: the Rosetta Stone of a hidden RNA language?” Cell, vol. 146, no. 3, pp. 353–358, 2011.
View at: Publisher Site | Google Scholar
G. Wang, X. Zheng, Y. Zheng et al., “Construction and analysis of the lncRNA-miRNA-mRNA network based on competitive endogenous RNA reveals functional genes in heart failure,” Molecular Medicine Reports, vol. 19, no. 2, pp. 994–1003, 2019.
View at: Publisher Site | Google Scholar
J. Yang and H. Song, “Identification of long noncoding RNA RP11-169F17.1 and RP11-669N7.2 as novel prognostic biomarkers of stomach adenocarcinoma based on integrated bioinformatics analysis,” Epigenomics, vol. 11, no. 11, pp. 1307–1321, 2019.
View at: Publisher Site | Google Scholar
H. Zhu, J. Lu, H. Zhao et al., “Functional long noncoding RNAs (lncRNAs) in clear cell kidney carcinoma revealed by reconstruction and comprehensive analysis of the lncRNA-miRNA-mRNA regulatory network,” Medical Science Monitor, vol. 24, pp. 8250–8263, 2018.
View at: Publisher Site | Google Scholar
N. Xu, Y. Cui, J. Dong, and L. Huang, “Exploring the molecular mechanisms of pterygium by constructing lncRNA-miRNA-mRNA regulatory network,” Investigative Opthalmology & Visual Science, vol. 61, no. 8, p. 12, 2020.
View at: Publisher Site | Google Scholar
T. Barrett, S. E. Wilhite, P. Ledoux et al., “NCBI GEO: archive for functional genomics data sets--update,” Nucleic Acids Research, vol. 41, Database issue, pp. D991–D995, 2013.
View at: Publisher Site | Google Scholar
W. Lan, A. Hou, R. Lakshminarayanan, Y. P. Lim, and L. Tong, “Linc-9432 is a novel pterygium lincRNA which regulates differentiation of fibroblasts,” FEBS Letters, vol. 592, no. 7, pp. 1173–1184, 2018.
View at: Publisher Site | Google Scholar
A. Hou, W. Lan, K. P. Law et al., “Evaluation of global differential gene and protein expression in primary pterygium: S100A8 and S100A9 as possible drivers of a signaling network,” PLoS One, vol. 9, no. 5, article e97402, 2014.
View at: Publisher Site | Google Scholar
A. K. Riau, T. T. Wong, W. Lan et al., “Aberrant DNA methylation of matrix remodeling and cell adhesion related genes in pterygium,” PLoS One, vol. 6, no. 2, article e14687, 2011.
View at: Publisher Site | Google Scholar
Y. W. Wong, J. Chew, H. Yang, D. T. Tan, and R. Beuerman, “Expression of insulin-like growth factor binding protein-3 in pterygium tissue,” The British Journal of Ophthalmology, vol. 90, no. 6, pp. 769–772, 2006.
View at: Publisher Site | Google Scholar
S. He, H. Sun, Y. Huang et al., “Identification and interaction analysis of significant genes and MicroRNAs in pterygium,” BioMed Research International, vol. 2019, Article ID 2767512, 12 pages, 2019.
View at: Publisher Site | Google Scholar
M. E. Ritchie, B. Phipson, D. Wu et al., “limma powers differential expression analyses for RNA-sequencing and microarray studies,” Nucleic Acids Research, vol. 43, no. 7, article e47, 2015.
View at: Publisher Site | Google Scholar
M. D. Paraskevopoulou, I. S. Vlachos, D. Karagkouni et al., “DIANA-LncBase v2: indexing microRNA targets on non-coding transcripts,” Nucleic Acids Research, vol. 44, no. D1, pp. D231–D238, 2016.
View at: Publisher Site | Google Scholar
J. H. Li, S. Liu, H. Zhou, L. H. Qu, and J. H. Yang, “starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein–RNA interaction networks from large-scale CLIP-Seq data,” Nucleic Acids Research, vol. 42, no. D1, pp. D92–D97, 2013.
View at: Publisher Site | Google Scholar
P. Shannon, A. Markiel, O. Ozier et al., “Cytoscape: a software environment for integrated models of biomolecular interaction networks,” Genome Research, vol. 13, no. 11, pp. 2498–2504, 2003.
View at: Publisher Site | Google Scholar
C. H. Chin, S. H. Chen, H. H. Wu, C. W. Ho, M. T. Ko, and C. Y. Lin, “cytoHubba: identifying hub objects and sub-networks from complex interactome,” BMC Systems Biology, vol. 8, Supplement 4, p. S11, 2014.
View at: Publisher Site | Google Scholar
Y. Zhou, B. Zhou, L. Pache et al., “Metascape provides a biologist-oriented resource for the analysis of systems-level datasets,” Nature Communications, vol. 10, no. 1, article 1523, 2019.
View at: Publisher Site | Google Scholar
D. Szklarczyk, A. Franceschini, S. Wyder et al., “STRING v10: protein-protein interaction networks, integrated over the tree of life,” Nucleic Acids Research, vol. 43, no. D1, pp. D447–D452, 2015.
View at: Publisher Site | Google Scholar
M. X. Li, L. T. Jin, T. J. Wang et al., “Identification of potential core genes in triple negative breast cancer using bioinformatics analysis,” Oncotargets and Therapy, vol. 11, pp. 4105–4112, 2018.
View at: Publisher Site | Google Scholar
J. Liu, X. Ding, L. Yuan, and X. Zhang, “Identification of pterygium-related long non-coding RNAs and expression profiling by microarray analysis,” International Journal of Molecular Medicine, vol. 38, no. 2, pp. 529–536, 2016.
View at: Publisher Site | Google Scholar
J. Zheng, Study on the differential expression profile of lncRNA in pterygium, Master, 2018.
D. H. Engelsvold, T. P. Utheim, O. K. Olstad et al., “miRNA and mRNA expression profiling identifies members of the miR-200 family as potential regulators of epithelial–mesenchymal transition in pterygium,” Experimental Eye Research, vol. 115, pp. 189–198, 2013.
View at: Publisher Site | Google Scholar
Y. H. Cui, H. Y. Li, Z. X. Gao et al., “Regulation of apoptosis by miR-122 in pterygium via targeting Bcl-w,” Investigative Ophthalmology & Visual Science, vol. 57, no. 8, pp. 3723–3730, 2016.
View at: Publisher Site | Google Scholar
J. H. Lee, S. A. Jung, Y. A. Kwon, J. L. Chung, and U. S. Kim, “Expression of microRNAs in fibroblast of pterygium,” International Journal of Ophthalmology, vol. 9, no. 7, pp. 967–972, 2016.
View at: Publisher Site | Google Scholar
W. Lan, S. Chen, and L. Tong, “MicroRNA-215 regulates fibroblast function: insights from a human fibrotic disease,” Cell Cycle, vol. 14, no. 12, pp. 1973–1984, 2015.
View at: Publisher Site | Google Scholar
C. W. Wu, Y. W. Cheng, N. Y. Hsu et al., “MiRNA-221 negatively regulated downstream p27Kip1 gene expression involvement in pterygium pathogenesis,” Molecular Vision, vol. 20, pp. 1048–1056, 2014.
View at: Google Scholar
G. İçme, A. Yilmaz, E. Dinç, A. Görür, Ş. B. Fidanci, and L. Tamer, “Assessment of miR-182, miR-183, miR-184, and miR-221 expressions in primary pterygium and comparison with the normal conjunctiva,” Eye & Contact Lens: Science & Clinical Practice, vol. 45, no. 3, pp. 208–211, 2019.
View at: Publisher Site | Google Scholar
S. Han, Y. Chen, Y. Gao, B. Sun, and Y. Kong, “MicroRNA-218-5p inhibit the migration and proliferation of pterygium epithelial cells by targeting EGFR via PI3K/Akt/mTOR signaling pathway,” Experimental Eye Research, vol. 178, pp. 37–45, 2019.
View at: Publisher Site | Google Scholar
C. A. Drummond, L. E. Crotty Alexander, S. T. Haller et al., “Cigarette smoking causes epigenetic changes associated with cardiorenal fibrosis,” Physiological Genomics, vol. 48, no. 12, pp. 950–960, 2016.
View at: Publisher Site | Google Scholar
D. Zhang, X. Liu, Q. Zhang, and X. Chen, “miR-138-5p inhibits the malignant progression of prostate cancer by targeting FOXC1,” Cancer Cell International, vol. 20, no. 1, p. 297, 2020.
View at: Publisher Site | Google Scholar
X. Bai, J. Shao, S. Zhou et al., “Inhibition of lung cancer growth and metastasis by DHA and its metabolite, RvD1, through miR-138-5p/FOXC1 pathway,” Journal of Experimental & Clinical Cancer Research, vol. 38, no. 1, p. 479, 2019.
View at: Publisher Site | Google Scholar
Y. Bai and X. Li, “hsa_circ_0008285 facilitates the progression of cervical cancer by targeting miR-211-5p/SOX4 axis,” Cancer Management and Research, vol. 12, pp. 3927–3936, 2020.
View at: Publisher Site | Google Scholar
M. Kang, J. Shi, B. Li, M. Luo, S. Xu, and X. Liu, “LncRNA DGCR5 regulates the non-small cell lung cancer cell growth, migration, and invasion through regulating miR-211-5p/EPHB6 axis,” BioFactors, vol. 45, no. 5, pp. 788–794, 2019.
View at: Publisher Site | Google Scholar
S. Yarahmadi, Z. Abdolvahabi, Z. Hesari et al., “Inhibition of sirtuin 1 deacetylase by miR-211-5p provides a mechanism for the induction of cell death in breast cancer cells,” Gene, vol. 711, p. 143939, 2019.
View at: Publisher Site | Google Scholar
Y. Guo, Y. Chen, H. Liu, and W. Yan, “Alpinetin inhibits oral squamous cell carcinoma proliferation via miR-211-5p upregulation and notch pathway deactivation,” Nutrition and Cancer, vol. 72, no. 5, pp. 757–767, 2020.
View at: Publisher Site | Google Scholar
M. Liang, J. Jia, L. Chen et al., “LncRNA MCM3AP-AS1 promotes proliferation and invasion through regulating miR-211-5p/SPARC axis in papillary thyroid cancer,” Endocrine, vol. 65, no. 2, pp. 318–326, 2019.
View at: Publisher Site | Google Scholar
C. Lin, Y. Zhang, Y. Chen, Y. Bai, and Y. Zhang, “Long noncoding RNA LINC01234 promotes serine hydroxymethyltransferase 2 expression and proliferation by competitively binding miR-642a-5p in colon cancer,” Cell Death & Disease, vol. 10, no. 2, p. 137, 2019.
View at: Publisher Site | Google Scholar
Z. S. Xia, L. Wang, T. Yu et al., “MiR-5000-3p, miR-5009-3P and miR-552: potential microRNA biomarkers of side population cells in colon cancer,” Oncology Reports, vol. 32, no. 2, pp. 589–596, 2014.
View at: Publisher Site | Google Scholar
C. W. Wu, M. L. Peng, K. T. Yeh, Y. Y. Tsai, C. C. Chiang, and Y. W. Cheng, “Corrigendum to “Inactivation of p53 in pterygium influence miR-200a expression resulting in ZEB1/ZEB2 up-regulation and EMT processing” [Exp. Eye Res. 146 (2016) 206-211],” Experimental Eye Research, vol. 151, p. 256, 2016.
View at: Publisher Site | Google Scholar
Z. Hamurcu, A. Ashour, N. Kahraman, and B. Ozpolat, “FOXM1 regulates expression of eukaryotic elongation factor 2 kinase and promotes proliferation, invasion and tumorgenesis of human triple negative breast cancer cells,” Oncotarget, vol. 7, no. 13, pp. 16619–16635, 2016.
View at: Publisher Site | Google Scholar
M. Egawa, Y. Yoshida, S. Ogura et al., “Increased expression of Forkhead box M1 transcription factor is associated with clinicopathological features and confers a poor prognosis in human hepatocellular carcinoma,” Hepatology Research, vol. 47, no. 11, pp. 1196–1205, 2017.
View at: Publisher Site | Google Scholar
Y. Liu, Y. Liu, B. Yuan et al., “FOXM1 promotes the progression of prostate cancer by regulating PSA gene transcription,” Oncotarget, vol. 8, no. 10, pp. 17027–17037, 2017.
View at: Publisher Site | Google Scholar
Y. Wang, W. Zhang, L. Wen et al., “FOXM1 confers resistance to gefitinib in lung adenocarcinoma via a MET/AKT-dependent positive feedback loop,” Oncotarget, vol. 7, no. 37, pp. 59245–59259, 2016.
View at: Publisher Site | Google Scholar
D. Wang, G. Hu, Y. du et al., “Aberrant activation of hedgehog signaling promotes cell proliferation via the transcriptional activation of forkhead Box M1 in colorectal cancer cells,” Journal of Experimental & Clinical Cancer Research, vol. 36, no. 1, p. 23, 2017.
View at: Publisher Site | Google Scholar
Y. Xie, D. Cui, L. Sui et al., “Induction of forkhead box M1 (FoxM1) by EGF through ERK signaling pathway promotes trophoblast cell invasion,” Cell and Tissue Research, vol. 362, no. 2, pp. 421–430, 2015.
View at: Publisher Site | Google Scholar
D. Hacioglu and H. Erdol, “Developments and current approaches in the treatment of pterygium,” International Ophthalmology, vol. 37, no. 4, pp. 1073–1081, 2017.
View at: Publisher Site | Google Scholar
B. Yu and G. Shan, “Functions of long noncoding RNAs in the nucleus,” Nucleus, vol. 7, no. 2, pp. 155–166, 2016.
View at: Publisher Site | Google Scholar
T. Derrien, R. Johnson, G. Bussotti et al., “The GENCODE v7 catalog of human long noncoding RNAs: analysis of their gene structure, evolution, and expression,” Genome Research, vol. 22, no. 9, pp. 1775–1789, 2012.
View at: Publisher Site | Google Scholar
B. Huang, P. Mu, Y. Yu et al., “Inhibition of EZH2 and activation of ERRγ synergistically suppresses gastric cancer by inhibiting FOXM1 signaling pathway,” Gastric Cancer, 2020.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2020 Yuting Xu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

575

Downloads

1218

Citations