Logo Medical Science Monitor

Call: +1.631.470.9640
Mon - Fri 10:00 am - 02:00 pm EST

Contact Us

Logo Medical Science Monitor Logo Medical Science Monitor Logo Medical Science Monitor

20 September 2021: Database Analysis  

Integrated Bioinformatic Analysis of Competing Endogenous RNA Network of Choriocarcinoma

Qianxia Tan1ABE, Zhihui Tan1BC, Junliang Liu1CD, Yanqun Mo1CF, Huining Liu1AG*

DOI: 10.12659/MSM.931475

Med Sci Monit 2021; 27:e931475

0 Comments

Abstract

BACKGROUND: Numerous studies have demonstrated that noncoding RNAs are involved in choriocarcinoma (CC). The competing endogenous RNA (ceRNA) network plays an important role in the occurrence and development of carcinoma. However, the involvement of the ceRNA network in CC remains unclear. The current study aimed to investigate the regulatory mechanism of ceRNA in CC.

MATERIAL AND METHODS: We downloaded the messenger RNAs (mRNAs) expression profiles (GSE20510 and GSE65654) and microRNAs (miRNAs) expression profiles (GSE32346 and GSE130489) from GEO datasets. The limma package of R software was used to identify differentially expressed RNAs (DERNAs). Then, we performed functional annotation of the differentially expressed mRNAs (DEmRNAs). TargetScan, miRDB, miRWalk, and Starbase were used to construct a CC-specific ceRNA network and select key molecules.

RESULTS: The results identified a total of 177 DEmRNAs and 189 differentially expressed miRNAs (DEmiRNAs) between the trophoblast and CC cell line samples. Ten differentially expressed lncRNAs (DElncRNAs) were obtained based on experimental studies. The DEmRNAs were mainly enriched in cell proliferation, positive regulation of the apoptotic process, and cell death. A total of 10 genes were ascertained as hub genes. Based on DEmRNAs, DEmiRNAs, and DElncRNAs, a CC-specific ceRNA network was established. Five DElncRNAs, 15 DEmiRNAs, and 45 DEmRNAs were identified. In addition, LINC00261, MEG3, MALAT1, H19, and OGFRP1 were identified as 5 key lncRNAs in choriocarcinoma.

CONCLUSIONS: This study provides novel insights into CC mechanisms and identified potential therapeutic targets for CC.

Keywords: Carcinoma, Choriocarcinoma, Gene Expression Profiling, Computational Biology, Datasets as Topic, Female, Humans, Pregnancy, RNA, Messenger, Uterine Neoplasms

Background

Choriocarcinoma (CC) is an aggressive, malignant trophoblastic neoplasm, either gestational or nongestational in origin [1]. Although CC is sensitive to chemotherapy, about 20% of patients with drug-resistant CC experience treatment failure and are at increased risk of metastatic lesions, hysterectomy, and decreased fertility [2,3]. Early metastasis and drug resistance in chemotherapy remain a clinical challenge in CC. Therefore, potential biomarkers and therapeutic targets in CC individualized treatment strategies need to be investigated.

Historically, most of the non-protein-coding parts of the human genome have been considered junk DNA. Although less than 2% of the human genome encodes proteins, the majority of all nucleotides can be detected under certain conditions [4]. MicroRNAs (miRNAs) are small endogenous RNAs with a regulatory function, consisting of 20–25 nucleotides. miRNA research in cancer is expanding because miRNAs are promising candidates for cancer biomarker development [5]. Functional studies have demonstrated that miRNA disorders are causal in many cancer cases, and miRNA mimics and miRNA inhibitors have shown promise as novel therapeutic agents [6]. Long noncoding RNAs (lncRNAs), a class of novel noncoding RNAs, are defined as transcripts of more than 200 nucleotides [7]. lncRNA-dysregulated expression is widespread in the pathogenesis of various human diseases, including cancer [8]. Studying lncRNA functions will provide researchers with tools and opportunities to develop lncRNA-based therapeutics for cancer [9].

Salmena and colleagues first proposed the concept of competitive endogenous RNA (ceRNA). Evidence has shown that pseudogene RNAs can act as a ‘sponge’ by competitively binding common miRNAs, releasing or weakening inhibitory effects through isolating miRNAs away from parental mRNA [10]. Accumulating evidence revealed the regulatory network of miRNAs and their target genes in carcinoma. Studies have shown that LINC01133 can competitively bind to miR-106a-3p and regulate adenomatous polyposis coli expression in gastric cancer [11]. Besides acting as a ceRNA sponge to absorb miR-570, LINC00612 can elevate the expression of PHD finger protein 14 in bladder cancer [12]. Previous studies suggested that lncRNA may act as ceRNA in CC. For example, lncRNA MALAT1 promoted the proliferation of CC by miR-218-mediated Fbxw8 regulation [13]. lncRNA RNA PCA3 contributes to the progression of CC by acting as a ceRNA against miR-106b [14]. A recent study reported targeted regulating of miR-515-5p by lncRNA LOXL1-AS1 in the proliferation and migration of CC [15].

How the ceRNA network is regulated in CC remains unclear. In the present study, we built a CC-specific ceRNA regulatory network using the RNA expression profiles collected from Gene Expression Omnibus (GEO) datasets. The research aimed to initially discover lncRNA-miRNA-mRNA ceRNA-mediated regulatory mechanisms in CC. This study may provide novel potential therapeutic targets for CC.

Material and Methods

SELECTION OF THE DATASETS:

We selected datasets including trophoblast cells and CC cell lines from the GEO database (http://www.ncbi.nlm.nih.gov/geo/) of microarray data on gene expression obtained from mRNA and miRNA. Two expression profiling datasets, GSE20510 and GSE65654, which identify mRNA involved in CC, were obtained from the GEO (GPL96 platform, [HG-U133A] Affymetrix Human Genome U133A Array, and GPL 10558 platform, Illumina HumanHT-12 V4.0 expression beadchip). Trophoblast cells included first-trimester cytotrophoblasts and immortalized extravillous trophoblast cell line HTR-8, as well as CC cell lines including JEG-3 and BeWo. The GSE20510 dataset contained 3 trophoblast samples and 7 CC cell line samples. The GSE65654 dataset contained 2 trophoblast samples and 4 CC cell line samples. We also downloaded the expression profiles of miRNA GSE32346 (2 trophoblast samples and 4 CC cell line samples) and GSE130489 (3 trophoblast samples and 3 CC cell line samples).

DATA QUALITY CONTROL:

Intra-group data repeatability in the was assessed by Pearson’s correlation test. We used the software and operating environment in the R programming language for drawing graphs and statistical analysis. The heat map in R software was used to visualize correlations between all samples. Principal component analysis (PCA) is an unsupervised statistical analysis technique and is a widely used sample clustering method. It is an efficient information compression method, which can reduce the dimensionality in the dataset and discover latent variables among samples. Sample clustering analysis was performed to assess the intra-group data repeatability.

SCREENING DEMRNAS AND DEMIRNAS:

We combined trophoblast samples and CC cell line samples data. The “Limma” package in R software was used to screen out DEmRNAs and DEmiRNAs. DEmRNAs were identified if | Log2 Fold Change (FC)| >1.5 and P value <0.05. Meanwhile, DEmiRNAs were defined if |Log2 FC| >2 and P value <0.05. Volcano maps were drawn using “ggplot2 package” operated by R software. The Circos was conducted by Omistudio, an online tool (http://www.Omistudio.cn/). The Circos delineated expression level and the overlaps between differently expressed RNA data visualized in the heatmap simultaneously. The Online Venn tool (http://bioinformatics.psb.ugent.be/webtools/Venn/) was used to create Venn diagrams. Venn diagrams exhibit common DEmRNAs and DEmiRNAs shared between the 2 datasets.

FUNCTIONAL ENRICHMENT OF DEMRNAS:

The online tool DAVID (https://david.ncifcrf.gov/home.jsp/) (version 6.8) was used for enrichment analysis. DAVID was used to reveal Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis enriched by DEmRNAs. GO is an ontology database widely used in bioinformatics, including cellular components, molecular functions, and biological processes. KEGG is a database resource for understanding high-level functions and utilities of the biological system. The “ggplot2 package” in R software was used to visualize the enrichment analysis. Cytoscape plug-ins ClueGO and CluePedia were used to perform enrichment analysis for DEmRNAs again. The hypergeometric significance P value<0.05 was regarded as significant.

PROTEIN–PROTEIN INTERACTION (PPI) NETWORK CONSTRUCTION:

Search Tool for the Retrieval of Interacting Genes (STRING; http://string-db.org/) was used to map the common DEmRNAs. The PPI network was constructed by Cytoscape 3.7.1. The app Molecular Complex Detection (MCODE) from the Cytoscape software suite was used to perform the gene network clustering analysis. MCODE identified the most significant PPI network modules; the submodule of PPI network was a degree MCODE scores >5, cut-off=2, Maximum depth=100, node score cut-off=0.2, and k-score=2. In CytoHubba, a plug-in of Cytoscape, the key genes were explored under the degrees set (degrees ≥10). The expression level of hub genes was drawn using the “pheatmap package” in R.

CONSTRUCTION OF A CC-SPECIFIC CERNA NETWORK:

To improve the accuracy of DERNAs, we excluded the DEmRNAs and DEmiRNAs that were not paired in miRDB (http://www.mirdb.org/), TargetScan (http://www.targetscan.org/), and miRWalk (http://mirwalk.umm.uni-heidelberg.de/). The common DERNAs were collected to construct a choriocarcinoma-specific ceRNA network. The prediction of lncRNA-miRNA interaction was investigated by Starbase (http://starbase.sysu.edu.cn/). We searched PubMed, Web of science, and Medline for DElncRNAs in CC based on experimental research. Finally, Cytoscape software was used to visualize the CC-specific ceRNA network.

STATISTICAL ANALYSIS:

The statistical significance of differences between the 2 groups was assessed using the t test. All analyses were conducted using R software 4.0.3. A P value <0.05 was considered statistically significant.

Results

QUALITY OF THE DATASETS:

Principal component analysis (PCA) and Pearson’s correlation test were used to test the intra-group data repeatability. As shown in Figure 1A1, 1A2, the samples in the trophoblast group had strong correlations in the GSE20510 and GSE65654 datasets, and the samples in the CC group were similar to the trophoblast group. Similarly, PCA results showed that the intra-group data repeatability for GSE20510 and GSE65654 was acceptable, and the trophoblast group and the CC group were separated from each other (Figure 1B1, 1B2). In addition, the samples in the trophoblast group had strong correlations in the GSE32346 and GSE130489, and the samples in the CC group were similar to the trophoblast group (Figure 1A3, 1A4). Based on the PCA, the intra-group data repeatability for GSE32346 and GSE130489 was acceptable, and the trophoblast group and the CC group were separated from each other (Figure 1B3, 1B4).

DEMRNAS AND DEMIRNAS IDENTIFIED BETWEEN TROPHOBLAST AND CHORIOCARCINOMA:

After the analyses of GSE20510, GSE65654, GSE32346, and GSE130489 datasets with the “Limma package” in R software, the DEmRNA between trophoblast cells and CC cell lines were presented in volcano plots (Figure 1C1, 1C2). Meanwhile, the DEmiRNAs were shown in volcano plots (Figure 1C3, 1C4). The overlap between differentially expressed mRNA datasets of GSE20510 and GSE65654 and the expression level visualized in the heatmap simultaneously were presented in Circos (Figure 1D1). The DEmiRNA lists of GSE32346 and GSE130489 were also performed in Circos (Figure 1D2). A total of 1205 DEmRNAs in the GSE20510 dataset, 854 DEmRNAs in the GSE65654 dataset, 316 DEmiRNAs in the GSE32346, and 324 DEmiRNAs in the GSE130489 were identified. Venn diagrams were constructed and revealed that 177 DEmRNAs (Supplementary Table 1) and 189 DEmiRNAs (Supplementary Table 2) were contained within the 2 examined datasets simultaneously (Figure 1E1, 1E2).

GO AND KEGG ENRICHMENT ANALYSIS:

Functional enrichment analysis was conducted to explore the biological functions of DEmRNAs. The top 5 biological processes (BP) terms were enriched in positive regulation of transcription from RNA polymerase II promoter, cell proliferation, heart development, angiogenesis, and positive regulation of the apoptotic process (Figure 2A, Table 1). The top 5 cell components (CC) terms of DEmRNAs were plasma membrane, extracellular exosome, cytosol, extracellular region, and extracellular space (Figure 2B). The top 3 molecular function (MF) terms were protein binding, sequence-specific DNA binding, and protein homodimerization activity (Figure 2C). The KEGG pathway analysis disclosed that the main enriched pathways included pathways in cancer, signaling pathways regulating pluripotency of stem cells, Wnt signaling pathway, proteoglycans in cancer, and Rap1 signaling pathway (Figure 2D). We used Cytoscape plug-ins ClueGO and CluePedia to analyze DEmRNAs and found that positive regulation of stress-activated MAPK cascade, interleukin-1 beta production, response to hydrogen peroxide, and connective tissue development were involved (Figure 2E).

PPI NETWORK CONSTRUCTION AND MODULE ANALYSIS:

A PPI network with 127 nodes and 280 edges was acquired; the network had an interaction score >0.4 according to the STRING online database (Figure 3A). The nodes represent genes, and the edges correspond to the links between genes. Blue represents downregulated genes and red represents upregulated genes. We used MCODE and CytoHubba in Cytoscape to perform network gene clustering to identify the critical module and hub genes. As shown in Figures 3B, 4 key modules with 9 upregulated genes and 15 downregulated genes were identified. The list of key modules genes is shown in Table 2. Furthermore, 10 genes (CYR61, FSTL1, CTGF, FBN1, APOA2, MXRA8, CKAP4, MGAT4A, IGFBP3, and IGFSP7) were determined as hub genes with the greatest degree of network connection (Figure 4A). Hierarchical cluster analysis revealed that the hub genes could effectively differentiate the CC cells from the trophoblast cells (Figure 4B). Among these genes, MGAT4A, APOA2, and MXRA8 were highly expressed in the CC cell line samples, while FSTL1, IGFBP3, IGFBP7, CYR61, FBN1, CKAP4, and CTGF were downregulated in the CC cell line samples. The hub genes may play a vital role in the occurrence or development of CC.

CC-SPECIFIC CERNA NETWORK:

Among the 15 obtained DEmiRNAs (Table 3), miR-10b, miR-155, miR-598, miR-411, miR-370, miR-134, miR-539, miR-10a, miR-376c, miR-495, and miR-152 served as the key upregulation miRNAs; miR-372, miR-520f, miR-934, and miR-1270 served as the key downregulation miRNAs. Due to the lack of high-throughput sequencing lncRNA in CC, we searched the database for DElncRNAs based on experimental research. Five DElncRNAs (LINC00261, MEG3, MALAT1, H19, and OGFRP1) were considered key lncRNAs regulating their target mRNAs’ expression by competing with miRNAs in CC. In total, 5 DElncRNAs (4 upregulated and 1 downregulated), 15 interacting DEmiRNAs (11 upregulated and 4 downregulated), and 45 DEmRNAs (18 upregulated and 27 downregulated) participated in the formation of CC-specific ceRNA regulatory network (Figure 5).

Discussion

Choriocarcinoma is a malignant tumor that occurs frequently in women of childbearing age. Emerging studies have indicated the regulatory roles of miRNAs and lncRNAs in the occurrence and development of CC. CeRNAs are transcripts that cross-regulate each other through competition with shared miRNAs. According to the ceRNA hypothesis, lncRNAs could serve as a sponge of target miRNAs [16]. CeRNA provides a novel way to explore cancer diagnosis, prognosis, and targeted treatment, but the ceRNA regulatory mechanism in CC has been unclear. Therefore, we constructed a CC-specific ceRNA network to investigate the regulatory mechanisms in CC.

In this research, we analyzed the differentially expressed RNAs to construct a ceRNA network. Herein, a train of bioinformatic analyses was performed on 4 independent high-throughput datasets (from trophoblast and CC cell line). By the analysis of GSE20510 and GSE65654 datasets, 177 common DEmRNAs were identified, of which 82 and 95 were upregulated and downregulated, respectively. It is well known that hub genes play critical roles in biological networks. Among the DEmRNAs, 21 hub nodes were identified by MCODE and 10 potential hub genes (CYR61, FSTL1, CTGF, FBN1, APOA2, MXRA8, CKAP4, MGAT4A, IGFBP3, and IGFSP7) were obtained using the CytoHubba plug-in of Cytoscape. Alpha-1,3-mannosyl-glycoprotein 4-beta-N-acetylglucosaminyltransferase A (MGAT4A) is a key enzyme in the pathway of synthesizing complex N-glycans. Several previous studies showed that MGAT4A is involved in the progression of choriocarcinoma, hepatocarcinoma, and colorectal carcinoma [17–19]. Chemokine receptors belong to the G-protein-coupled receptors (GPCRs) family, which appears to be involved in inflammatory diseases, tumor growth, and metastasis. A positive correlation between chemokine receptor expression and worse prognosis has been found in most but not all cancers. CC-chemokine receptor 7 (CCR7) is significantly overexpressed in colorectal tumor tissues compared with paired normal tissues. Significant differences were also found in overall survival and disease-free survival [20]. In breast cancer, it has been reported that the chemokine receptors CCR7 and their ligands are associated with metastasis [21]. This is consistent with our finding that CCR7 is a key node and has relatively high expression in CC.

The GO analyses showed that the DEmRNAs are mainly enriched in cell proliferation, positive regulation of the apoptotic process, and cell death. The KEGG pathway analysis also showed that the DEmRNAs mostly participate in pathways in cancer, signaling pathways regulating pluripotency of stem cells, Rap1 signaling pathway, and Wnt signaling pathway. The Wnt signaling pathway promotes decidualization, endometrial function, and trophoblast differentiation, and changes in Wnt signaling components were found in cancers of reproductive tissues, endometriosis, and gestational diseases [23]. Wnt/β-catenin signaling regulates StarD7 transcription in CC cell line JEG-3 cells through the canonical pathway [23].

Noncoding RNAs (ncRNAs) are of interest because of their versatile roles in biological and pathological processes, particularly the widely studied small miRNAs and long noncoding RNAs (lncRNAs) [24]. Through the analysis of GSE32346 and GSE130489 datasets, 189 common DEmiRNAs were obtained, of which 147 were downregulated and 42 were upregulated. Fifteen DEmiRNAs (11 upregulated and 4 downregulated) serve as key target miRNAs in the ceRNA network. It has been reported that miR-1270 is a negative regulator in liver carcinoma and is involved in post-transcriptional regulation of centromere protein M, and hepatitis B virus X protein can promote liver carcinoma by inhibiting miR1270 [25]. miR-1270 is aberrantly upregulated in papillary thyroid cancer. Reduced miR-1270 expression suppresses cancer cell proliferation and migration, which was reversed by suppression of cancer cell invasion downregulation [26]. miR-411 acts as an oncogene in osteosarcoma and esophageal squamous cell carcinoma [27,28]. On the contrary, miR-411 inhibits the growth and metastasis of bladder cancer, cervical cancer, and gastric cancer [29–31].

Five DElncRNAs, 15 DEmiRNAs, and 45 DEmRNAs were identified as abnormally expressed RNAs and we built a CC-specific network. This study found MALAT1 and LINC00261 can increase CCR7 via competing with miR-1270. The metastasis-associated lung adenocarcinoma transcript 1(MALAT1) was originally identified as a prognostic parameter for lung adenocarcinoma or squamous cell carcinoma [32]. A series of data has been collected linking MALAT1 with cancer types and provided new insights into its biogenesis, and molecular functions [33]. Several studies suggested that MALAT1 can be aberrantly expressed in choriocarcinoma, ovarian, cervical, colorectal cancers, and breast cancer [13,34–37]. MALAT1 may be an oncogene that promotes CC proliferation and tumor growth through miR-218-mediated F-box/WD repeat-containing protein 8 regulation [13]. Low expression of lncRNA LINC00261 was reported in multiple cancers, including breast, colon, and lung cancer. These results indicate LINC00261 is a tumor suppressor that halts cellular proliferation, invasion, and tumor growth [38–40]. LINC00261 was significantly underexpressed in CC tissues and CC cell lines. Overexpression of LINC00261 resulted in decreased cell proliferation and cell cycle arrest in G0/G1 phase, suppressed cell migration and invasion, and promoted cell apoptosis and the relative activities of caspase 3 and caspase 9 in CC cells [41]. We speculated that upregulating MALAT1 indirectly suppresses the expression of TCF4 via competitively binding to miR-155. miR-155 is involved in the syncytialization of trophoblast cells induced by lipopolysaccharides [42]. It has been reported that LIN28A has a functional role in regulating trophoblast differentiation and function, and that loss of LIN28A in human trophoblast is sufficient to induce differentiation [43]. In the CC ceRNA network, LIN28A was upregulated in CC and might be regulated by miR-598.

Conclusions

In the present study, we constructed a CC-specific ceRNA regulatory network by bioinformatics analysis. A CC-specific ceRNA regulatory network might provide novel insights into CC mechanisms and identify potential therapeutic targets for CC.

Figures

(A) A1–A4 represents Pearson’s correlation analysis of GSE20510, GSE65654, GSE32346, and GSE130489 datasets, respectively. The color reflects the intensity of the correlation. The greater the absolute value of a number, the higher the correlation. (B1–B4) PCA of samples from the GSE20510, GSE65654, GSE32346, and GSE130489 datasets. In the figure, the x-axis is principal component 1 (PC1) and the y-axis is principal component 2 (PC2). (C1–C4) The volcano plot of the GSE20510, GSE65654, GSE32346, and GSE130489 datasets. (D1) Circos of GSE20510 and GSE65654 was performed using the online tool Omistudio. The outer circle represents gene lists, the inner circle represents the gene expression level, and the red curve represents the overlap between differentially expressed gene lists. (D2) Circos of GSE32346 and GSE130489. (E1, E2) 177 sharing genes between the GSE20510 and GSE65654 datasets are shown in Venn diagram, and 189 differentially expressed miRNAs were contained in the GSE32346 and GSE130489 simultaneously.Figure 1. (A) A1–A4 represents Pearson’s correlation analysis of GSE20510, GSE65654, GSE32346, and GSE130489 datasets, respectively. The color reflects the intensity of the correlation. The greater the absolute value of a number, the higher the correlation. (B1–B4) PCA of samples from the GSE20510, GSE65654, GSE32346, and GSE130489 datasets. In the figure, the x-axis is principal component 1 (PC1) and the y-axis is principal component 2 (PC2). (C1–C4) The volcano plot of the GSE20510, GSE65654, GSE32346, and GSE130489 datasets. (D1) Circos of GSE20510 and GSE65654 was performed using the online tool Omistudio. The outer circle represents gene lists, the inner circle represents the gene expression level, and the red curve represents the overlap between differentially expressed gene lists. (D2) Circos of GSE32346 and GSE130489. (E1, E2) 177 sharing genes between the GSE20510 and GSE65654 datasets are shown in Venn diagram, and 189 differentially expressed miRNAs were contained in the GSE32346 and GSE130489 simultaneously. (A) The top 14 changes in DEmRNAs’ biological processes (BP) in trophoblast and choriocarcinoma samples through GO enrichment analyses. (B) Cellular components (CC) analysis of DEmRNAs. (C) Molecular functions (MF) analysis of DEmRNAs. (D) The KEGG pathway analysis of DEmRNAs. (E) Using Cytoscape plug-in ClueGO and CluePedia to display the network of enriched terms colored by cluster identity, where nodes that share the same cluster identity are typically close to each other. DEmRNAs – differentially expressed mRNAs; GO – Gene Ontology; KEGG – Kyoto Encyclopedia of Genes and Genomes.Figure 2. (A) The top 14 changes in DEmRNAs’ biological processes (BP) in trophoblast and choriocarcinoma samples through GO enrichment analyses. (B) Cellular components (CC) analysis of DEmRNAs. (C) Molecular functions (MF) analysis of DEmRNAs. (D) The KEGG pathway analysis of DEmRNAs. (E) Using Cytoscape plug-in ClueGO and CluePedia to display the network of enriched terms colored by cluster identity, where nodes that share the same cluster identity are typically close to each other. DEmRNAs – differentially expressed mRNAs; GO – Gene Ontology; KEGG – Kyoto Encyclopedia of Genes and Genomes. (A) The protein–protein interaction (PPI) analysis in choriocarcinoma. The blue and red nodes represent downregulated and upregulated genes, respectively. The degree value is displayed by the size of the nodes. (B) The most significant modules in the PPI network.Figure 3. (A) The protein–protein interaction (PPI) analysis in choriocarcinoma. The blue and red nodes represent downregulated and upregulated genes, respectively. The degree value is displayed by the size of the nodes. (B) The most significant modules in the PPI network. (A) Cytohubba a plug-in of Cytoscape identified the hub genes from the PPI network. (B) Hierarchical cluster analysis suggested that the hub genes could effectively differentiate the choriocarcinoma samples from the trophoblast samples in the GSE20510 and GSE65654 datasets. Blue represents downregulated genes, red represents upregulated genes.Figure 4. (A) Cytohubba a plug-in of Cytoscape identified the hub genes from the PPI network. (B) Hierarchical cluster analysis suggested that the hub genes could effectively differentiate the choriocarcinoma samples from the trophoblast samples in the GSE20510 and GSE65654 datasets. Blue represents downregulated genes, red represents upregulated genes. The lncRNA-miRNA-mRNA ceRNA network constructed by Cytoscape. The diamond, square, and circle represent lncRNA, miRNA, and mRNA, respectively. Blue stands for upregulated RNAs, red stands for downregulated RNAs. lncRNAs – long noncoding RNAs; miRNAs – microRNAs; mRNA – messenger RNAs; ceRNA – competing endogenous RNA.Figure 5. The lncRNA-miRNA-mRNA ceRNA network constructed by Cytoscape. The diamond, square, and circle represent lncRNA, miRNA, and mRNA, respectively. Blue stands for upregulated RNAs, red stands for downregulated RNAs. lncRNAs – long noncoding RNAs; miRNAs – microRNAs; mRNA – messenger RNAs; ceRNA – competing endogenous RNA.

References

1. Savage J, Adams E, Veras E, Choriocarcinoma in women: Analysis of a case series with genotyping: Am J Surg Pathol, 2017; 41(12); 1593-606

2. Yu S, Wu C, Tan Q, Long noncoding RNA H19 promotes chemotherapy resistance in choriocarcinoma cells: J Cell Biochem, 2019; 120(9); 15131-44

3. Wu C, Yu S, Tan Q, Role of AhR in regulating cancer stem cell-like characteristics in choriocarcinoma: Cell Cycle, 2018; 17(18); 2309-20

4. Djebali S, Davis CA, Merkel A, Landscape of transcription in human cells: Nature, 2012; 489(7414); 101-8

5. Lu TX, Rothenberg ME, MicroRNA: J Allergy Clin Immunol, 2018; 141(4); 1202-7

6. 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); 203-22

7. Kopp F, Mendell JT, Functional classification and experimental dissection of long noncoding RNAs: Cell, 2018; 172(3); 393-407

8. Ma Y, Zhang J, Wen L, Membrane-lipid associated lncRNA: A new regulator in cancer signaling: Cancer Lett, 2018; 419; 27-29

9. Choudhari R, Sedano MJ, Harrison AL, Long noncoding RNAs in cancer: From discovery to therapeutic targets: Adv Clin Chem, 2020; 95; 105-47

10. Salmena L, Poliseno L, Tay Y, A ceRNA hypothesis: The Rosetta Stone of a hidden RNA language?: Cell, 2011; 146(3); 353-58

11. Yang XZ, Cheng TT, He QJ, LINC01133 as ceRNA inhibits gastric cancer progression by sponging miR-106a-3p to regulate APC expression and the Wnt/β-catenin pathway: Mol Cancer, 2018; 17(1); 126

12. Miao L, Liu HY, Zhou C, LINC00612 enhances the proliferation and invasion ability of bladder cancer cells as ceRNA by sponging miR-590 to elevate expression of PHF14: J Exp Clin Cancer Res, 2019; 38(1); 143

13. Shi D, Zhang Y, Lu R, The long non-coding RNA MALAT1 interacted with miR-218 modulates choriocarcinoma growth by targeting Fbxw8: Biomed Pharmacother, 2018; 97; 543-50

14. Wang YN, Liu SY, Wang L, Long noncoding RNA PCA3 contributes to the progression of choriocarcinoma by acting as a ceRNA against miR-106b: Int J Clin Exp Pathol, 2019; 12(5); 1609-17

15. Zhang L, Wan Q, Zhou H, Targeted-regulating of miR-515-5p by LncRNA LOXL1-AS1 on the proliferation and migration of trophoblast cells: Exp Mol Pathol, 2021; 118; 104588

16. Qi X, Zhang DH, Wu N, ceRNA in cancer: possible functions and clinical implications: J Med Genet, 2015; 52(10); 710-18

17. Nishino K, Yamamoto E, Niimi K, N-acetylglucosaminyltransferase IVa promotes invasion of choriocarcinoma: Oncol Rep, 2017; 38(1); 440-48

18. Guo Y, Li S, Qu J, Let-7c inhibits metastatic ability of mouse hepatocarcinoma cells via targeting mannoside acetylglucosaminyltransferase 4 isoenzyme A: Int J Biochem Cell Biol, 2014; 53; 1-8

19. Fan J, Wang S, Yu S, N-acetylglucosaminyltransferase IVa regulates metastatic potential of mouse hepatocarcinoma cells through glycosylation of CD147: Glycoconj J, 2012; 29(5–6); 323-34

20. Gao L, Xu J, He G, CCR7 high expression leads to cetuximab resistance by cross-talking with EGFR pathway in PI3K/AKT signals in colorectal cancer: Ame J Cancer Res, 2019; 9(11); 2531-43

21. Rizeq B, Malki MI, The Role of CCL21/CCR7 chemokine axis in breast cancer progression: Cancers, 2020; 12(4); 1036

22. Sonderegger S, Pollheimer J, Knöfler M, Wnt signalling in implantation, decidualisation and placental differentiation – review: Placenta, 2010; 31(10); 839-47

23. Rena V, Angeletti S, Panzetta-Dutari G, Activation of beta-catenin signalling increases StarD7 gene expression in JEG-3 cells: Placenta, 2009; 30(10); 876-83

24. Guo L, Zhang Q, Ma X, miRNA and mRNA expression analysis reveals potential sex-biased miRNA expression: Sci Rep, 2017; 7; 39812

25. Xiao Y, Najeeb RM, Ma D, Upregulation of CENPM promotes hepatocarcinogenesis through mutiple mechanisms: J Exp Clin Cancer Res, 2019; 38(1); 458

26. Yi T, Zhou X, Sang K, MicroRNA-1270 modulates papillary thyroid cancer cell development by regulating SCAI: Biomed Pharmacother, 2019; 109; 2357-64

27. Xu N, Yang W, Liu Y, MicroRNA-411 promoted the osteosarcoma progression by suppressing MTSS1 expression: Environ Sci Pollut Res Int, 2018; 25(12); 12064-71

28. Song H, Song J, Lu L, SNHG8 is upregulated in esophageal squamous cell carcinoma and directly sponges microRNA-411 to increase oncogenicity by upregulating KPNA2: Onco Targets Ther, 2019; 12; 6991-7004

29. Liu Y, Liu T, Jin H, MiR-411 suppresses the development of bladder cancer by regulating ZnT1: Onco Targets Ther, 2018; 11; 8695-704

30. Shan D, Shang Y, Hu T, MicroRNA-411 inhibits cervical cancer progression by directly targeting STAT3: Oncol Res, 2019; 27(3); 349-58

31. Bai TL, Liu YB, Li BH, MiR-411 inhibits gastric cancer proliferation and migration through targeting SETD6: Eur Rev Med Pharmacol Sci, 2019; 23(8); 3344-50

32. Ji P, Diederichs S, Wang W, MALAT-1, a novel noncoding RNA, and thymosin beta4 predict metastasis and survival in early-stage non-small cell lung cancer: Oncogene, 2003; 22(39); 8031-41

33. Gutschner T, Hämmerle M, Diederichs S, MALAT1 – a paradigm for long noncoding RNA function in cancer: J Mol Med, 2013; 91(7); 791-801

34. Tao F, Tian X, Ruan S, miR-211 sponges lncRNA MALAT1 to suppress tumor growth and progression through inhibiting PHF19 in ovarian carcinoma: FASEB J, 2018 [Online ahead of print]

35. Shen F, Zheng H, Zhou L, Overexpression of MALAT1 contributes to cervical cancer progression by acting as a sponge of miR-429: J Cell Physiol, 2019; 234(7); 11219-26

36. Wu S, Sun H, Wang Y, MALAT1 rs664589 polymorphism inhibits binding to miR-194-5p, contributing to colorectal cancer risk, growth, and metastasis: Cancer Res, 2019; 79(20); 5432-41

37. Kim J, Piao HL, Kim BJ, Long noncoding RNA MALAT1 suppresses breast cancer metastasis: Nat Genet, 2018; 50(12); 1705-15

38. Shahabi S, Kumaran V, Castillo J, LINC00261 is an epigenetically regulated tumor suppressor essential for activation of the DNA damage response: Cancer Res, 2019; 79(12); 3050-62

39. Zhang HF, Li W, Han YD, LINC00261 suppresses cell proliferation, invasion and Notch signaling pathway in hepatocellular carcinoma: Cancer Biomark, 2018; 21(3); 575-82

40. Yan D, Liu W, Liu Y, LINC00261 suppresses human colon cancer progression via sponging miR-324-3p and inactivating the Wnt/β-catenin pathway: J Cell Physiol, 2019; 234(12); 22648-56

41. Wang Y, Xue K, Guan Y, Long noncoding RNA LINC00261 suppresses cell proliferation and invasion and promotes cell apoptosis in human choriocarcinoma: Oncol Res, 2017; 25(5); 733-42

42. Dai Y, Diao Z, Sun H, MicroRNA-155 is involved in the remodelling of human-trophoblast-derived HTR-8/SVneo cells induced by lipopolysaccharides: Hum Reprod, 2011; 26(7); 1882-91

43. Seabrook JL, Cantlon JD, Cooney AJ, Role of LIN28A in mouse and human trophoblast cell differentiation: Biol Reprod, 2013; 89(4); 95

Figures

Figure 1. (A) A1–A4 represents Pearson’s correlation analysis of GSE20510, GSE65654, GSE32346, and GSE130489 datasets, respectively. The color reflects the intensity of the correlation. The greater the absolute value of a number, the higher the correlation. (B1–B4) PCA of samples from the GSE20510, GSE65654, GSE32346, and GSE130489 datasets. In the figure, the x-axis is principal component 1 (PC1) and the y-axis is principal component 2 (PC2). (C1–C4) The volcano plot of the GSE20510, GSE65654, GSE32346, and GSE130489 datasets. (D1) Circos of GSE20510 and GSE65654 was performed using the online tool Omistudio. The outer circle represents gene lists, the inner circle represents the gene expression level, and the red curve represents the overlap between differentially expressed gene lists. (D2) Circos of GSE32346 and GSE130489. (E1, E2) 177 sharing genes between the GSE20510 and GSE65654 datasets are shown in Venn diagram, and 189 differentially expressed miRNAs were contained in the GSE32346 and GSE130489 simultaneously.Figure 2. (A) The top 14 changes in DEmRNAs’ biological processes (BP) in trophoblast and choriocarcinoma samples through GO enrichment analyses. (B) Cellular components (CC) analysis of DEmRNAs. (C) Molecular functions (MF) analysis of DEmRNAs. (D) The KEGG pathway analysis of DEmRNAs. (E) Using Cytoscape plug-in ClueGO and CluePedia to display the network of enriched terms colored by cluster identity, where nodes that share the same cluster identity are typically close to each other. DEmRNAs – differentially expressed mRNAs; GO – Gene Ontology; KEGG – Kyoto Encyclopedia of Genes and Genomes.Figure 3. (A) The protein–protein interaction (PPI) analysis in choriocarcinoma. The blue and red nodes represent downregulated and upregulated genes, respectively. The degree value is displayed by the size of the nodes. (B) The most significant modules in the PPI network.Figure 4. (A) Cytohubba a plug-in of Cytoscape identified the hub genes from the PPI network. (B) Hierarchical cluster analysis suggested that the hub genes could effectively differentiate the choriocarcinoma samples from the trophoblast samples in the GSE20510 and GSE65654 datasets. Blue represents downregulated genes, red represents upregulated genes.Figure 5. The lncRNA-miRNA-mRNA ceRNA network constructed by Cytoscape. The diamond, square, and circle represent lncRNA, miRNA, and mRNA, respectively. Blue stands for upregulated RNAs, red stands for downregulated RNAs. lncRNAs – long noncoding RNAs; miRNAs – microRNAs; mRNA – messenger RNAs; ceRNA – competing endogenous RNA.

In Press

21 Feb 2024 : Clinical Research  

Potential Value of HSP90α in Prognosis of Triple-Negative Breast Cancer

Med Sci Monit In Press; DOI: 10.12659/MSM.943049  

22 Feb 2024 : Review article  

Differentiation of Native Vertebral Osteomyelitis: A Comprehensive Review of Imaging Techniques and Future ...

Med Sci Monit In Press; DOI: 10.12659/MSM.943168  

23 Feb 2024 : Clinical Research  

A Study of 60 Patients with Low Back Pain to Compare Outcomes Following Magnetotherapy, Ultrasound, Laser, ...

Med Sci Monit In Press; DOI: 10.12659/MSM.943732  

26 Feb 2024 : Clinical Research  

Predictive Value of Combined HbA1c and Neutrophil-to-Lymphocyte Ratio for Diabetic Peripheral Neuropathy in...

Med Sci Monit In Press; DOI: 10.12659/MSM.942509  

Most Viewed Current Articles

17 Jan 2024 : Review article  

Vaccination Guidelines for Pregnant Women: Addressing COVID-19 and the Omicron Variant

DOI :10.12659/MSM.942799

Med Sci Monit 2024; 30:e942799

0:00

16 May 2023 : Clinical Research  

Electrophysiological Testing for an Auditory Processing Disorder and Reading Performance in 54 School Stude...

DOI :10.12659/MSM.940387

Med Sci Monit 2023; 29:e940387

0:00

14 Dec 2022 : Clinical Research  

Prevalence and Variability of Allergen-Specific Immunoglobulin E in Patients with Elevated Tryptase Levels

DOI :10.12659/MSM.937990

Med Sci Monit 2022; 28:e937990

0:00

01 Jan 2022 : Editorial  

Editorial: Current Status of Oral Antiviral Drug Treatments for SARS-CoV-2 Infection in Non-Hospitalized Pa...

DOI :10.12659/MSM.935952

Med Sci Monit 2022; 28:e935952

0:00

Your Privacy

We use cookies to ensure the functionality of our website, to personalize content and advertising, to provide social media features, and to analyze our traffic. If you allow us to do so, we also inform our social media, advertising and analysis partners about your use of our website, You can decise for yourself which categories you you want to deny or allow. Please note that based on your settings not all functionalities of the site are available. View our privacy policy.

Medical Science Monitor eISSN: 1643-3750
Medical Science Monitor eISSN: 1643-3750