Exploration of the molecular mechanisms of cervical cancer based on mRNA expression profiles and predicted microRNA interactions

mRNA-lncRNA gene expression signature in HPV-associated neoplasia and cervical cancer

Cervical cancer is one of the most frequent cancers in women and is associated with human papillomavirus (HPV) in 70 % of cases. Cervical cancer occurs because of progression of low-differentiated cervical intraepithelial neoplasia through grade 2 and 3 lesions. Along with the protein-coding genes, long noncoding RNAs (lncRNAs) play an important role in the development of malignant cell transformation. Although human papillomavirus is widespread, there is currently no well-characterized transcriptomic signature to predict whether this tumor will develop in the presence of HPV-associated neoplastic changes in the cervical epithelium. Changes in gene activity in tumors reflect the biological diversity of cellular phenotype and physiological functions and can be an important diagnostic marker. We performed comparative transcriptome analysis using open RNA sequencing data to assess differentially expressed genes between normal tissue, neoplastic epithelium, and cervical cancer. Raw data were preprocessed using the Galaxy platform. Batch effect correction, identification of differentially expressed genes, and gene set enrichment analysis (GSEA) were performed using R programming language packages. Subcellular localization of lncRNA was analyzed using Locate-R and iLoc-LncRNA 2.0 web services. 1,572 differentially expressed genes (DEGs) were recorded in the "cancer vs. control" comparison, and 1,260 DEGs were recorded in the "cancer vs. neoplasia" comparison. Only two genes were observed to be differentially expressed in the "neoplasia vs. control" comparison. The search for common genes among the most strongly differentially expressed genes among all comparison groups resulted in the identification of an expression signature consisting of the CCL20, CDKN2A, CTCFL, piR-55219, TRH, SLC27A6 and EPHA5 genes. The transcription level of the CCL20 and CDKN2A genes becomes increased at the stage of neoplastic epithelial changes and stays so in cervical cancer. Validation on an independent microarray dataset showed that the differential expression patterns of the CDKN2A and SLC27A6 genes were conserved in the respective gene expression comparisons between groups.

Anti-apoptotic Effect of MiR-223-3p Suppressing PIK3C2A in Cardiomyocytes from Myocardial Infarction Rat Through Regulating PI3K/Akt Signaling Pathway

We aimed to explore the regulatory mechanism of the axis of miR-223-3p-PIK3C2A-PI3K/Akt on cardiomyocyte apoptosis in rats with myocardial infarction. Thirty 8-week-old healthy male SD rats were used for establishing the sham group and the model group, with HE staining, TUNEL staining, and TTC staining performed. After the identification of the targeting relationship between PIK3C2A and miR-223-3p, experimental rats were randomly divided into seven groups by plasmid transfection, including the Blank group, negative control (NC) group, miR-223-3p mimic group, miR-223-3p inhibitor group, siRNA-PIK3C2A group, oe-PIK3C2A group, and miR-223-3p inhibitor + oe-PIK3C2A group. Four weeks after transfection, the expression levels of miR-223-3p and PIK3C2A in tissues as well as PI3K, Akt, Bax, and bcl-2 mRNA in cells were detected by qRT-PCR and western blot, in combination with the detection of apoptosis rate by flow cytometry. Compared with the sham group, the model group showed typical myocardial injury and abnormal staining, higher apoptotic index, and larger myocardial infarction area (all P < 0.05). PIK3C2A was the target gene of miR-223-3p. The expression level of miR-223-3p in model group was significantly lower than that in sham group, while the mRNA and protein expression levels of PIK3C2A increased significantly (all P < 0.05). In cell tests, the expression level of miR-223-3p increased significantly in miR-223-3p mimic group (P < 0.05), which, however, showed no significant change in siRNA-PIK3C2A group (P > 0.05). MiR-223-3p inhibitor group and siRNA-PIK3C2A group had obviously increased PI3K, Akt, mTOR and Bcl-2 mRNA, and protein expression, while decreased mRNA and protein expression of PIK3C2A and Bax (all P < 0.05); miR-223-3p mimic groups had the opposite trends (all P < 0.05). siRNA-PIK3C2A + miR-223-3p mimic showed no obvious change relative to the control groups (all P > 0.05). Low expression of miR-223-3p may downregulate PIK3C2A expression, resulting in the inhibition of myocardial cell apoptosis in rats with myocardial infarction via the activation of PI3K/Akt signaling pathway.

Bioinformatics analysis of differentially expressed genes and pathways in the development of cervical cancer

Background

This study aimed to explore and identify key genes and signaling pathways that contribute to the progression of cervical cancer to improve prognosis.

Methods

Three gene expression profiles (GSE63514, GSE64217 and GSE138080) were screened and downloaded from the Gene Expression Omnibus database (GEO). Differentially expressed genes (DEGs) were screened using the GEO2R and Venn diagram tools. Then, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed. Gene set enrichment analysis (GSEA) was performed to analyze the three gene expression profiles. Moreover, a protein–protein interaction (PPI) network of the DEGs was constructed, and functional enrichment analysis was performed. On this basis, hub genes from critical PPI subnetworks were explored with Cytoscape software. The expression of these genes in tumors was verified, and survival analysis of potential prognostic genes from critical subnetworks was conducted. Functional annotation, multiple gene comparison and dimensionality reduction in candidate genes indicated the clinical significance of potential targets.

Results

A total of 476 DEGs were screened: 253 upregulated genes and 223 downregulated genes. DEGs were enriched in 22 biological processes, 16 cellular components and 9 molecular functions in precancerous lesions and cervical cancer. DEGs were mainly enriched in 10 KEGG pathways. Through intersection analysis and data mining, 3 key KEGG pathways and related core genes were revealed by GSEA. Moreover, a PPI network of 476 DEGs was constructed, hub genes from 12 critical subnetworks were explored, and a total of 14 potential molecular targets were obtained.

Conclusions

These findings promote the understanding of the molecular mechanism of and clinically related molecular targets for cervical cancer.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-08412-4.

Comprehensive analysis of aberrantly expressed long non‑coding RNAs, microRNAs, and mRNAs associated with the competitive endogenous RNA network in cervical cancer

Cervical cancer is a common malignant disease that poses a serious health threat to women worldwide. Growing research efforts have focused on protein‑coding and non‑coding RNAs involved in the tumorigenesis and prognosis of various types of cancer. The potential molecular mechanisms and the interaction among long non‑coding RNAs (lncRNAs), microRNAs (miRNAs), and mRNAs require further investigation in cervical cancer. In the present study, lncRNA, miRNA, and mRNA expression profiles of 304 primary tumor tissues from patients with cervical cancer and 3 solid normal tissues from The Cancer Genome Atlas (TCGA) dataset were studied via RNA sequencing (RNA‑seq). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed using R package clusterProfiler to annotate the principal functions of differentially expressed (DE) mRNAs. Kaplan‑Meier analysis was also conducted to investigate the effects of DElncRNAs, DEmiRNAs, and DEmRNAs on overall survival. A total of 2,255 mRNAs, 133 miRNAs, and 150 lncRNAs that were differentially expressed were identified with a threshold of P<0.05 and |fold change (FC)|>2. Functional enrichment analysis indicated that DEmRNAs were enriched in cancer‑associated KEGG pathways. Furthermore, 255 mRNAs, 15 miRNAs, and 12 lncRNAs that were significantly associated with overall survival in cervical carcinoma were also identified. Importantly, an miRNA‑mediated competitive endogenous RNA (ceRNA) network was successfully constructed based on the expression profiles of DElncRNAs and DEmRNAs. More importantly, it was found that the lncRNA EPB41L4A‑AS1 may function as a pivotal regulator in carcinoma of the uterine cervix. Taken together, the present study has provided novel insights into investigating the potential mechanisms underlying tumorigenesis, development, and prognosis of cervical cancer, and presented new potential avenues for cancer therapeutics.

The New Biomarker for Cervical Squamous Cell Carcinoma and Endocervical Adenocarcinoma (CESC) Based on Public Database Mining

To reconstruct the ceRNA biological network of cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC) and to select an appropriate mRNA as a biomarker that could be used for CESC early diagnosis and prognosis evaluation. We downloaded CESC data from the TCGA public database, and statistical analysis was conducted with the R software to find out differential expressed genes encoding for lncRNAs, miRNAs, and mRNAs. The differentially expressed mRNAs (DEmRNAs) screened in the ceRNA network were analyzed for survival to find the mRNAs with significantly linked to the survival prognosis. These mRNAs were searched in the Pathological Atlas to identify the final appropriate mRNAs. Differential expression analysis revealed 773 lncRNAs, 94 miRNAs, and 2466 mRNAs. Survival analysis of DEmRNAs in the ceRNA network indicated that ADGRF4, ANXA8L1, HCAR3, IRF6, and PDE2A (P < 0.05) were negatively correlated with survival time. Verification of these six DEmRNAs in the Pathology Atlas indicated that PDE2A was a possible biomarker for CESC patients. PDE2A might be a biomarker for early diagnosis and prognosis evaluation of CESC patients, but due to the lack of available data, further studies may be needed for confirmation.

Aldehyde dehydrogenase 1, a target of miR-222, is expressed at elevated levels in cervical cancer

The aim of the present study was to investigate the expression of microRNA-222 (miR-222) and aldehyde dehydrogenase 1 (ALDH1) in tissues and peripheral blood of cervical cancer patients, and to elucidate their underlying mechanisms of action. Tumor tissues and tumor-adjacent tissues were obtained from 33 cervical cancer patients and peripheral blood was obtained from these patients and 28 healthy subjects. The expression of miR-222 and ALDH1 mRNA was evaluated by reverse transcription-quantitative polymerase chain reaction (RT-qPCR). To examine the levels of ALDH1 protein in tissues and blood, western blotting and ELISA were used. To confirm a direct interaction between miR-222 and ALDH1 mRNA, a dual luciferase reporter assay was performed. HeLA cells were transfected with agomiR-222 and expression of ALDH1 in the cells was measured by RT-qPCR and western blotting. MTT assay was preform to investigate the proliferation of HeLA cells. Expression of ALDH1 mRNA and protein was elevated in cervical cancer tissues and peripheral blood from patients compared with tumor-adjacent tissues and healthy controls, while the expression of miR-222 was reduced. Upregulation of miR-222 inhibited HeLA cell proliferation possibly due to a reduction in the expression of ALDH1. A dual luciferase reporter assay showed that miR-222 can bind with the 3′-untranslated seed region of ALDH1 mRNA to regulate its expression. miR-222 regulation of ALDH1 expression may play a role in the prevention of cervical cancer.

The expression of special AT-rich binding protein 1 in cervical cancer and its clinical significance

Background

The oncogenic potential of special AT-rich binding protein 1 (SATB1) has been reported in various types of cancer, but its function in cervical cancer remains not fully investigated. This study aimed to investigate the effect of SATB1 mRNA expression on tumor progression and outcomes in the cervical cancer patients.

Methods

A total of 33 cervical cancer patients treated in our hospital from September 2012 to December 2015 were included. The mRNA expression level of STAB1 in cervical cancer tissue was determined by real-time PCR, and the patients were divided into dichotomous groups based on their SATB1 expression level. Clinical characteristics, recurrence, and survival outcomes were compared between groups.

Results

Compared with the SATB1-low group, the SATB1-high group had significantly advanced International Federation of Gynecology and Obstetrics (FIGO) stages (P=0.037) and histologic grade (P=0.036). Kaplan–Meier analysis showed that SATB1-high group had a worse overall survival (P=0.078, marginal significant). In the subgroup analysis of pathological types, adenocarcinomas group (n=8) had a significantly higher SATB1 expression level as compared with the squamous cell carcinomas (n=18) and adenosquamous carcinomas (n=7) groups (both P<0.05). Cervical squamous cell carcinomas patients with a high-expression SATB1 (n=8) had more advanced FIGO stages (P=0.015) and histologic grades (P=0.060, marginal significant) as well as a higher (P=0.069, marginal significant) incidence of lymphatic metastasis than those with a low expression of SATB1 (n=10).

Conclusion

These results showed that expression of SATB1 may have an effect on the disease progression and survival outcome of cervical cancer.

MicroRNA-218-5p inhibits cell growth and metastasis in cervical cancer via LYN/NF-κB signaling pathway

Background

We are committed to investigate miR-218-5 effects on the progression of cervical cancer (CC) cell and find out the molecular mechanism.

Methods

GSE9750 was obtained from GEO database and R Limma package was applied to filter out dysregulated genes. The pathways were enriched by GSEA software, ClusterProfiler and enrichplot packages to predict the function of DEGs. The binding sites of LYN were detected by miRanda and TargetScan. The miR2Disease database was used to find miRNAs related with CC. The expression of miR-218-5p and LYN were quantified by qRT-PCR and that of LYN protein was measured by western blot. The targeted relationships between miR-218-5p and LYN were verified by dual-luciferase reporter assay. Colony formation assays, wound healing, transwell invasion assay and flow cytometer analysis were performed to investigate the roles that miR-218-5p and LYN played in migration, invasion and death of cervical carcinoma. Xenografts established in nude mice were used to assess tumor growth in vivo.

Results

The highly expressed mRNA LYN was selected by microarray analysis in GSE9750. NF-κB signaling pathway was enriched base on GSEA results. The expression of miR-218-5p was lower but LYN was higher in CC primary tumors compared with normal control. In addition, miR-218-5p could regulate the expression of LYN in HeLa cells negatively. Overexpression of LYN could promote cell migration and invasion, but inhibit cell death in vitro, and also promote tumor formation in vivo via activating NF-κB signaling pathway which could be reversed by miR-218-5p.

Conclusions

MiR-218-5p suppressed the progression of CC via LYN/NF-κB signaling pathway.

Integrated Analysis of a Competing Endogenous RNA Network Revealing a Prognostic Signature for Cervical Cancer

Given the high morbidity and the trend of younger individuals being affected observed in cervical cancer, it is important to identify sensitive and effective biomarkers for predicting the survival outcome of patients. Based on data from 307 cervical cancer cases acquired from The Cancer Genome Atlas portal, 1920 differentially expressed mRNAs, 70 microRNAs(miRNAs), and 493 long non-coding(lncRNAs) were screened by comparing cervical cancer tissues with paracancerous tissues. A competing endogenous (ceRNA) network containing 50 lncRNAs, 16 miRNAs, and 81 mRNAs was constructed. Eighteen RNAs, comprising 13 mRNAs, 2 miRNAs, and 3 lncRNAs, were identified as significant prognostic factors by univariate Cox proportional hazards regression. ETS-related gene and fatty acid synthase signatures were discovered using a multivariate Cox regression model built to identify independent prognostic factors in cervical cancer patients. Receiver operating characteristic (ROC) analysis was used to determine the optimal cut-off value for distinguishing the risk level of cervical cancer patients. High-risk patients exhibited a poorer prognosis than low-risk patients did. This study focused on ceRNA networks to provide a novel perspective and insight into cervical cancer and suggested that the identified signature can serve as an independent prognostic biomarker in cervical cancer.