Inhibition of EZH2 and activation of ERRγ synergistically suppresses gastric cancer by inhibiting FOXM1 signaling pathway

Exploring the impact of estrogen-related receptor gamma on metabolism and disease

Estrogen-related receptor gamma (ERRγ) is a member of the ERR orphan nuclear receptor family which possesses three subtypes, α, β, and γ. ERRγ is reportedly predominantly expressed in metabolically active tissues and cells, which promotes positive and negative effects in different tissues. ERRγ overexpression in the liver, pancreas, and thyroid cells is related to liver cancer, oxidative stress, reactive oxygen species (ROS) regulation, and carcinoma. Reduced ERRγ expression in the brain, immune cells, tumor cells, and energy metabolism causes neurological dysfunction, gastric cancer, and obesity. ERRγ is a constitutive receptor; however, its transcriptional activity also depends on co-regulators, agonists, and antagonists, which, when after forming a complex, can play a role in targeting and treating diseases. Moreover, ERRγ has proven crucial in regulating cellular and metabolic activity. However, many functions mediated via ERRγ remain unknown and require further exploration. Hence, considering the importance of ERRγ, this review focuses on the critical findings and interactions between ERRγ and co-regulators, agonists, and antagonists alongside its relationship with downstream and upstream signaling pathways and diseases. This review highlights new findings and provides a path to understanding the current ideas and future studies on ERRγ-mediated cellular activity.

DY131 activates ERRγ/TFAM axis to protect against metabolic disorders and acute kidney injury

Renal tubular injury is considered as the main pathological feature of acute kidney injury (AKI), and mitochondrial dysfunction in renal tubular cells is implicated in the pathogenesis of AKI. The estrogen-related receptor γ (ERRγ) is a member of orphan nuclear receptors which plays a regulatory role in mitochondrial biosynthesis, energy metabolism and many metabolic pathways. Online datasets showed a dominant expression of ERRγ in renal tubules, but the role of ERRγ in AKI is still unknown. In the present study, we investigated the role of ERRγ in the pathogenesis of AKI and the therapeutic efficacy of ERRγ agonist DY131 in several murine models of AKI. ERRγ expression was reduced in kidneys of AKI patients and AKI murine models along with a negative correlation to the severity of AKI. Consistently, silencing ERRγ in vitro enhanced cisplatin-induced tubular cells apoptosis, while ERRγ overexpression in vivo utilizing hydrodynamic-based tail vein plasmid delivery approach alleviated cisplatin-induced AKI. ERRγ agonist DY131 could enhance the transcriptional activity of ERRγ and ameliorate AKI in various murine models. Moreover, DY131 attenuated the mitochondrial dysfunction of renal tubular cells and metabolic disorders of kidneys in AKI, and promoted the expression of the mitochondrial transcriptional factor A (TFAM). Further investigation showed that TFAM could be a target gene of ERRγ and DY131 might ameliorate AKI by enhancing ERRγ-mediated TFAM expression protecting mitochondria. These findings highlighted the protective effect of DY131 on AKI, thus providing a promising therapeutic strategy for AKI.

Metabolic regulator ERRγ governs gastric stem cell differentiation into acid-secreting parietal cells

Parietal cells (PCs) produce gastric acid to kill pathogens and aid digestion. Dysregulated PC census is common in disease, yet how PCs differentiate is unclear. Here, we identify the PC progenitors arising from isthmal stem cells, using mouse models and human gastric cells, and show that they preferentially express cell-metabolism regulator and orphan nuclear receptor Estrogen-related receptor gamma (Esrrg, encoding ERRγ). Esrrg expression facilitated the tracking of stepwise molecular, cellular, and ultrastructural stages of PC differentiation. EsrrgP2ACreERT2 lineage tracing revealed that Esrrg expression commits progenitors to differentiate into mature PCs. scRNA-seq indicated the earliest Esrrg+ PC progenitors preferentially express SMAD4 and SP1 transcriptional targets and the GTPases regulating acid-secretion signal transduction. As progenitors matured, ERRγ-dependent metabolic transcripts predominated. Organoid and mouse studies validated the requirement of ERRγ for PC differentiation. Our work chronicles stem cell differentiation along a single lineage in vivo and suggests ERRγ as a therapeutic target for PC-related disorders.

Research Progress in Estrogen-related Receptor Gamma (ERRγ) Agonists and Inverse Agonists

Estrogen-related receptor gamma (ERRγ), one of three members of the ERR family, is an inducible transcription factor. ERRγ has dual functions in different tissues. The decreased expression of ERRγ in the brain, stomach, prostate, and fat cells can cause neuropsychological dysfunction, gastric cancer, prostate cancer, and obesity. However, when ERRγ is present in the liver, pancreas, and thyroid follicular cells, ERRγ overexpression is related to liver cancer, type II diabetes, oxidative liver injury, and anaplastic thyroid carcinoma. Signaling pathway studies have confirmed that ERRγ agonists or inverse agonists can regulate ERRγ expression to treat related diseases. The collision between residue Phe435 and the modulator is a key factor determining the activation or inhibition of ERRγ. Although more than 20 agonists and inverse agonists of ERRγ have been reported, no clinical studies have been found in the literature. This review summarizes the important relationship between ERRγ-related signaling pathways and diseases, research progress, and the structure-activity relationship of modulators. These findings provide guidance for further study on new ERRγ modulators.

Melatonin inhibited the progression of gastric cancer induced by Bisphenol S via regulating the estrogen receptor 1

In recent years, Bisphenol S (BPS) has increasingly been used as an alternative to Bisphenol A (BPA) in food, paper, and personal care products. It is imperative to clarify the relationship between BPS and tumors in order to treat and prevent diseases. This study discovered a new method for predicting tumor correlations between BPS interactive genes. According to analyses conducted by Gene Ontology and Kyoto Encyclopedia of Genes and Genomes, interactive genes were primarily found in gastric cancer. Based on gene-targeted prediction and molecular docking, BPS appears to exert potential gastric cancer-causing effects through estrogen receptor 1 (ESR1). In addition, gastric cancer patients' prognosis could be accurately predicted by a bisphenol-based prognostic prediction model. Subsequently, the proliferation and migration abilities of gastric cancer cells were further demonstrated to be significantly enhanced by BPS. Similarly, molecular docking analysis revealed that melatonin is also highly correlated with gastric cancer and BPS. In cell proliferation and migration assays, melatonin and BPS exposure inhibited the invasion abilities of gastric cancer cells compared to BPS-exposure. Our research provided a new direction for the exploration the correlation between cancer and environmental toxicity.

Role of Estrogen-Related Receptor γ and PGC-1α/SIRT3 Pathway in Early Brain Injury After Subarachnoid Hemorrhage

Estrogen-related receptors (ERRs) were shown to play an important role in the regulation of free radical-mediated pathology. This study aimed to investigate the neuroprotective effect of ERRγ activation against early brain injury (EBI) after subarachnoid hemorrhage (SAH) and the potential underlying mechanisms. In a rat model of SAH, the time course of ERRs and SIRT3 and the effects of ERRγ activation were investigated. ERRγ agonist DY131, selective inhibitor GSK5182, or SIRT3 selective inhibitor 3-TYP were administered intracerebroventricularly (icv) in the rat model of SAH. The use of 3-TYP was for validating SIRT3 as the downstream signaling of ERRγ activation. Post-SAH assessments included SAH grade, neurological score, Western blot, Nissl staining, and immunofluorescence staining in rats. In an vitro study, the ERRγ agonist DY131 and ERRγ siRNA were administered to primary cortical neurons stimulated by Hb, after which cell viability and neuronal deaths were accessed. Lastly, the brain ERRγ levels and neuronal death were accessed in SAH patients. We found that brain ERRγ expressions were significantly increased, but the expression of SIRT3 dramatically decreased after SAH in rats. In the brains of SAH rats, ERRγ was expressed primarily in neurons, astrocytes, and microglia. The activation of ERRγ with DY131 significantly improved the short-term and long-term neurological deficits, accompanied by reductions in oxidative stress and neuronal apoptosis at 24 h after SAH in rats. DY131 treatment significantly increased the expressions of PGC-1α, SIRT3, and Bcl-2 while downregulating the expressions of 4-HNE and Bax. ERRγ antagonist GSK5182 and SIRT3 inhibitor 3-TYP abolished the neuroprotective effects of ERRγ activation in the SAH rats. An in vitro study showed that Hb stimulation significantly increased intracellular oxidative stress in primary cortical neurons, and DY131 reduced such elevations. Primary cortical neurons transfected with the ERRγ siRNA exhibited notable apoptosis and abolished the protective effect of DY131. The examination of SAH patients' brain samples revealed increases in ERRγ expressions and neuronal apoptosis marker CC3. We concluded that ERRγ activation with DY131 ameliorated oxidative stress and neuronal apoptosis after the experimental SAH. The effects were, at least in part, through the ERRγ/PGC-1α/SIRT3 signaling pathway. ERRγ may serve as a novel therapeutic target to ameliorate EBI after SAH.

Multifaceted Transcriptional Network of Estrogen-Related Receptor Alpha in Health and Disease

Estrogen-related receptors (ERRα, β and γ in mammals) are orphan members of the nuclear receptor superfamily acting as transcription factors. ERRs are expressed in several cell types and they display various functions in normal and pathological contexts. Amongst others, they are notably involved in bone homeostasis, energy metabolism and cancer progression. In contrast to other nuclear receptors, the activities of the ERRs are apparently not controlled by a natural ligand but they rely on other means such as the availability of transcriptional co-regulators. Here we focus on ERRα and review the variety of co-regulators that have been identified by various means for this receptor and their reported target genes. ERRα cooperates with distinct co-regulators to control the expression of distinct sets of target genes. This exemplifies the combinatorial specificity of transcriptional regulation that induces discrete cellular phenotypes depending on the selected coregulator. We finally propose an integrated view of the ERRα transcriptional network.

The BRD4 inhibitor JQ1 augments the antitumor efficacy of abemaciclib in preclinical models of gastric carcinoma

BACKGROUND

Advanced gastric cancer (GC) is a lethal malignancy, harboring recurrent alterations in cell cycle pathway, especially the CDKN2A-CDK4/CDK6/CCND1 axis. However, monotherapy of CDK4/6 inhibitors has shown limited antitumor effects for GC, and combination treatments were urgently needed for CDK4/6 inhibitors.

METHODS

Here, we performed a comprehensive analysis, including drug screening, pan-cancer genomic dependency analysis, and epigenetic sequencing to identify the candidate combination with CDK4/6 inhibitors. Mechanisms were investigated by bulk RNA-sequencing and experimental validation was conducted on diverse in vitro or in vivo preclinical GC models.

RESULTS

We found that the BRD4 inhibitor JQ1 augments the antitumor efficacy of the CDK4/6 inhibitor abemaciclib (ABE). Diverse in vitro and in vivo preclinical GC models are examined and synergistic benefits from the combination therapy are obtained consistently. Mechanistically, the combination of ABE and JQ1 enhances the cell cycle arrest of GC cells and induces unique characteristics of cellular senescence through the induction of DNA damage, which is revealed by transcriptomic profiling and further validated by substantial in vitro and in vivo GC models.

CONCLUSION

This study thus proposes a candidate combination therapy of ABE and JQ1 to improve the therapeutic efficacy and worth further investigation in clinical trials for GC.

Chromatin and noncoding RNA-mediated mechanisms of gastric tumorigenesis

Gastric cancer (GC) is one of the most common and deadly cancers in the world. It is a multifactorial disease highly influenced by environmental factors, which include radiation, smoking, diet, and infectious pathogens. Accumulating evidence suggests that epigenetic regulators are frequently altered in GC, playing critical roles in gastric tumorigenesis. Epigenetic regulation involves DNA methylation, histone modification, and noncoding RNAs. While it is known that environmental factors cause widespread alterations in DNA methylation, promoting carcinogenesis, the chromatin- and noncoding RNA-mediated mechanisms of gastric tumorigenesis are still poorly understood. In this review, we focus on discussing recent discoveries addressing the roles of histone modifiers and noncoding RNAs and the mechanisms of their interactions in gastric tumorigenesis. A better understanding of epigenetic regulation would likely facilitate the development of novel therapeutic approaches targeting specific epigenetic regulators in GC.

Cancer-derived exosomal circ_0038138 enhances glycolysis, growth, and metastasis of gastric adenocarcinoma via the miR-198/EZH2 axis

This study aims to decipher the impact and downstream mechanisms of the bioinformatically identified circ_0038138 delivered by cancer-derived exosomes in gastric adenocarcinoma (GAC). Expression of circ_0038138 in clinical GAC tissues and exosomes (Exos) from clinical plasma samples (plasma-Exos) was predicted by bioinformatics analysis and validated by RT-qPCR. The binding affinity between circ_0038138, miR-198 and EZH2 was identified using luciferase activity, RIP, and RNA pull-down assays. GAC cells (AGS) were co-cultured with Exos isolated from GAC cell supernatant (GC9811-P). After co-culture, the behaviors of GAC cells including proliferation and glycolysis were assessed to identify the biological effect of exosomal circ_0038138. Also, in vivo effects of exosomal circ_0038138 on the tumorigenesis and lung metastasis of GAC cells were evaluated by developing nude mouse xenograft and metastatic models. circ_0038138 upregulation was detected in GAC tissues and plasma-Exos. Exos delivered circ_0038138 to GAC cells and potentiated the proliferative, migratory, invasive, and glycolytic potentials of GAC cells. Mechanistically, circ_0038138 competitively bound to miR-198, which in turn targeted EZH2 by binding to its 3'-UTR. Silencing of EZH2 promoted CXXC4 expression and inhibited Wnt/β-catenin pathway activation, thus repressing the malignancy and glycolysis of GAC cells. In vivo assay confirmed that exosomal circ_0038138 induced tumorigenesis and lung metastasis by regulating the miR-198/EZH2 axis. Collectively, our work suggests that the Exo-mediated transfer of circ_0038138 potentially facilitates the glycolysis, growth and metastasis of GAC cells via miR-198/EZH2 axis, which offers a potential prognostic marker and a therapeutic target for GAC.

Enhancer of zeste homolog 2 exerts functions in gastric cancer development via modulating microRNA-222-3p methylation and WEE1 expression

OBJECTIVE

Enhancer of zeste homolog 2 (EZH2) has been studied in gastric cancer (GC), while the role of EZH2 in GC via binding to microRNA (miR)-222-3p remains obscure. This research aims to unravel the regulatory mechanism of EZH2 in GC progression via the modulation of miR-222-3p/WEE1 axis.

METHODS

Initially, EZH2, miR-222-3p and WEE1 levels in GC cells and tissues were examined. Thereafter, vectors altering EZH2, miR-222-3p or WEE1 expression were transfected into GC cells to determine the malignant behaviors involved in tumorigenesis of GC cells. Finally, the targeting relations among EZH2, miR-222-3p and WEE1 were validated.

RESULTS

EZH2 and WEE1 were up-regulated while miR-222-3p was down-regulated in GC tissues and cells. The decreased EZH2, silenced WEE1, or restored miR-222-3p constrained the malignant behaviors involved in tumorigenesis of GC cells. Deletion of miR-222-3p could reverse the effect of silenced EZH2 on suppressing the biological functions of GC cells. EZH2 could bind to the promoter of miR-222-3p, and there was a targeting relationship between miR-222-3p and WEE1.

CONCLUSION

Our study demonstrates that EZH2 promotes GC development via the modulation of miR-222-3p/WEE1 axis, thus providing promising therapeutic targets for GC therapy.

Nuclear receptor estrogen-related receptor gamma suppresses colorectal cancer aggressiveness by regulating Wnt/β-catenin signalling

Colorectal cancer is the predominant cause of cancer-related death worldwide, because of lack of effective therapeutic targets. Estrogen-related receptor gamma (ESRRG), which belongs to the family of nuclear receptors, functions as an important element regulating gene transcription. In our report, we identified ESRRG as a potential tumor suppressor. The decreased level of ESRRG was initially observed in CRC and was highly associated with poor prognosis. ESRRG overexpression abrogated cell growth and metastasis in vitro and in vivo. Mechanistically, ESRRG repressed the epithelial-to-mesenchymal transition (EMT) process and antagonized Wnt signaling by regulating β-catenin degradation. In addition, significant ESRRG hypermethylation was found in CRC and inversely correlated with its expression. Consistently, the expression of ESRRG was induced after treatment with DNA demethylating agent 5-AZA. Taken together, these findings define a tumor-suppressive role of ESRRG in CRC, providing a potential novel therapeutic approach for this cancer.

High Expression of EZH2 Mediated by ncRNAs Correlates with Poor Prognosis and Tumor Immune Infiltration of Hepatocellular Carcinoma

Background: Hepatocellular carcinoma (HCC) is the predominant form of liver cancer and is accompanied by a complex regulatory network. Increasing evidence suggests that an abnormal gene expression of EZH2 is associated with HCC progression. However, the molecular mechanism by which non-coding RNAs (ncRNAs) regulate EZH2 remains elusive. Methods: The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) data were used to perform differential expression analysis and prognostic analysis. We used the Encyclopedia of RNA Interactomes (ENCORI) database to predict candidate miRNAs and lncRNAs that may bind to EZH2. Subsequently, the comprehensive analysis (including expression analysis, correlation analysis, and survival analysis) identified ncRNAs that contribute to EZH2 overexpression. Results: EZH2 was found to be upregulated in the majority of tumor types and associated with a poor prognosis. Hsa-miR-101-3p was identified as a target miRNA of EZH2. Additionally, SNHG6 and MALAT1 were identified as upstream lncRNAs of hsa-miR-101-3p. Meanwhile, correlation analysis revealed that EZH2 expression was significantly associated with the infiltration of several immune cell types in HCC. Conclusion: SNHG6 or MALAT1/hsa-miR-101-3p/EZH2 axis were identified as potential regulatory pathways in the progression of HCC.

The long and short non-coding RNAs modulating EZH2 signaling in cancer

Non-coding RNAs (ncRNAs) are a large family of RNA molecules with no capability in encoding proteins. However, they participate in developmental and biological processes and their abnormal expression affects cancer progression. These RNA molecules can function as upstream mediators of different signaling pathways and enhancer of zeste homolog 2 (EZH2) is among them. Briefly, EZH2 belongs to PRCs family and can exert functional roles in cells due to its methyltransferase activity. EZH2 affects gene expression via inducing H3K27me3. In the present review, our aim is to provide a mechanistic discussion of ncRNAs role in regulating EZH2 expression in different cancers. MiRNAs can dually induce/inhibit EZH2 in cancer cells to affect downstream targets such as Wnt, STAT3 and EMT. Furthermore, miRNAs can regulate therapy response of cancer cells via affecting EZH2 signaling. It is noteworthy that EZH2 can reduce miRNA expression by binding to promoter and exerting its methyltransferase activity. Small-interfering RNA (siRNA) and short-hairpin RNA (shRNA) are synthetic, short ncRNAs capable of reducing EZH2 expression and suppressing cancer progression. LncRNAs mainly regulate EZH2 expression via targeting miRNAs. Furthermore, lncRNAs induce EZH2 by modulating miRNA expression. Circular RNAs (CircRNAs), like lncRNAs, affect EZH2 expression via targeting miRNAs. These areas are discussed in the present review with a focus on molecular pathways leading to clinical translation.

CircGSK3B promotes RORA expression and suppresses gastric cancer progression through the prevention of EZH2 trans-inhibition

BACKGROUND

Circular RNAs (circRNAs) are a class of non-coding RNA that play critical roles in the development and pathogenesis of various cancers. The circRNA circGSK3B (hsa_circ_0003763) has been shown to enhance cell proliferation, migration, and invasion in hepatocellular carcinoma. However, the specific functions and underlying mechanistic involvement of circGSK3B in gastric cancer (GC) have not yet been explored. Our study aimed to investigate the effect of circGSK3B on the progression of GC and to identify any potential mechanisms underlying this process.

METHODS

CircRNA datasets associated with GC were obtained from the PubMed, GEO, and ArrayExpress databases, and circRNAs were validated via RT-qPCR and Sanger sequencing. Biotin-labeled RNA pull-down, mass spectrometry, RNA immunoprecipitation, and in vitro binding assays were employed to determine proteins demonstrating interactions with circGSK3B. Gene expression regulation was assessed through RT-qPCR, chromatin immunoprecipitation, and western blot assays. Gain- and loss-of-function assays were used to analyze any effects of circGSK3B and its partner regulatory molecule (EZH2) on the proliferation, invasion, and migration abilities of GC cells both in vitro and in vivo.

RESULTS

CircGSK3B was mainly identified in the nucleus. This circRNA was present at a reduced concentration in GC tissues and cell lines. Overexpression of circGSK3B was shown to inhibit the growth, invasion, and metastasis of GC cells both in vitro and in vivo. Mechanistically, circGSK3B directly interacted with EZH2, acting to suppress the binding of EZH2 and H3K27me3 to the RORA promoter, and leading to an elevation in RORA expression and ultimately the suppression of GC progression.

CONCLUSIONS

CircGSK3B acts as a tumor suppressor, reducing EZH2 trans-inhibition and GC progression. This demonstrates the potential use of this RNA as a therapeutic target for GC.

Estrogen-Related Receptor γ Agonist DY131 Ameliorates Lipopolysaccharide-Induced Acute Liver Injury

Sepsis-associated liver dysfunction remains a challenge in clinical practice with high mortality and limited specific therapies. DY131 is a pharmacological agonist of the orphan receptor estrogen-related receptor (ERR) γ which plays a crucial role in regulating energy generation, oxidative metabolism, cell apoptosis, inflammatory responses, etc. However, its role in acute liver injury is unknown. In this study, we evaluated the effect of DY131 on lipopolysaccharide (LPS)-induced liver injury. Mice were pretreated with DY131 through intraperitoneal injection at a dose of 5 mg/kg/day for 3 days prior to LPS challenge (10 mg/kg). 24 h later, they were anesthetized and sacrificed. Blood and liver tissues were collected for further studies. In a separate experiment, mice were treated with saline (vehicle) or DY131 for 3 days to evaluate the toxicity of DY131. We found that ERRγ was downregulated in the liver tissues from LPS-treated mice. Pretreatment with DY131 ameliorated LPS-induced liver injury as demonstrated by reduced liver enzyme release (ALT, AST, and LDH), improved liver morphological damage, and attenuated oxidative stress, inflammation and apoptosis. Meanwhile, DY131 had no significant side effects on hepatic and renal functions in mice. Finally, transcriptomics analysis revealed that the dysregulated pathways associated with inflammation and metabolism were significantly reversed by DY131 in LPS-treated mice, providing more evidence in favor of the protective effect of DY131 against LPS-induced liver injury. Altogether, these findings highlighted the protective effect of DY131 on LPS-induced hepatotoxicity possibly via suppressing oxidative stress, inflammation, and apoptosis.

Identification of Functional Genes in Pterygium Based on Bioinformatics Analysis

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.