Abstract:Background Mounting evidence suggests that noncoding RNAs (lncRNAs) were involved in various human cancers. However, the role of these lncRNAs in HPV-driven cervical cancer (CC) has not been extensively studied. Considering that HR-HPV infections contribute to cervical carcinogenesis by regulating the expression of lncRNAs, miRNAs and mRNAs, we aim to systematically analyze lncRNAs and mRNAs expression profile to identify novel lncRNAs-mRNAs co-expression networks and explore their potential impact on tumorigenesis in HPV-driven CC. Methods LncRNA/mRNA microarray technology was utilized to identify the differentially expressed lncRNAs (DElncRNAs) and mRNAs (DEmRNAs) in HPV-16 and HPV-18 cervical carcinogenesis compared to normal cervical tissues. Venn diagram and weighted gene co-expression network analysis (WGCNA) were used to identify the hub DElncRNAs/DEmRNAs that were both significantly correlated with HPV-16 and HPV-18 CC patients . LncRNA-mRNA correlation analysis and functional enrichment pathway analysis were performed on these key DElncRNAs/DEmRNAs in HPV-16 and HPV-18 CC patients to explore their mutual mechanism in HPV-driven CC. A lncRNA-mRNA co-expression score (CES) model was established and validated by using the Cox regression method. Afterward, the clinicopathological characteristics were analyzed between CES-high and CES-low groups. In vitro, functional experiments were performed to evaluate the role of LINC00511 and PGK1 in cell proliferation, migration and invasion in CC cells. To understand whether LINC00511 play as an oncogenic role partially via modulating the expression of PGK1, rescue assays were used. Results We identified 81 lncRNAs and 211 mRNAs that were commonly differentially expressed in HPV-16 and HPV-18 CC tissues compared to normal tissues. The results of lncRNA-mRNA correlation analysis and functional enrichment pathway analysis showed that the LINC00511-PGK1 co-expression network may make an important contribution to HPV-mediated tumorigenesis and be closely associated with metabolism-related mechanisms. Combined with clinical survival data, the prognostic lncRNA-mRNA co-expression score (CES) model based on LINC00511 and PGK1 could precisely predict patients' overall survival (OS). CES-high patients had a worse prognosis than CES-low patients and the enriched pathways and potential targets of applicable drugs were explored in CES-high patients. In vitro experiments confirmed the oncogenic functions of LINC00511 and PGK1 in the progression of CC, and revealed that LINC00511 functions in an oncogenic role in CC cells partially via modulating the expression of PGK1. Conclusions Together, these data identify co-expression modules that provide valuable information to understand the pathogenesis of HPV-mediated tumorigenesis, which highlights the pivotal function of the LINC00511-PGK1 co-expression network in cervical carcinogenesis. Furthermore, our CES model has a reliable predicting ability that could stratify CC patients into low- and high-risk groups of poor survival. This study provides a bioinformatics method to screen prognostic biomarkers which leads to lncRNA-mRNA co-expression network identification and construction for patients' survival prediction and potential drug applications in other cancers.

Co-expression Network Analysis of the Lncrnas and Mrnas Associated with Cervical Cancer Progression

Integrated Bioinformatical Analysis of Lncrna-Mrna Co-Expression Profiles of Cervical Cancer

Identification of Lncrnas by Microarray Analysis Reveals the Potential Role of Lncrnas in Cervical Cancer Pathogenesis

Comprehensive Analysis of Lncrnas Microarray Profile and Mrna-Lncrna Co-Expression in Oncogenic HPV-positive Cervical Cancer Cell Lines.

ceRNA Network and WGCNA Analyses of Differentially Expressed Genes in Cervical Cancer Tissues for Association with Survival of Patients

Integrated Analysis of ceRNA Regulatory Network Associated With Tumor Stage in Cervical Cancer

Identification of the Coregulated Mrna and Lncrna Functional Pathways in Ovarian Cancer Using LncRNAs2Pathways by Lncrnas Based on an Integrated Network Propagation

Development Immune-Related Prognostic Model and LncRNA-miRNA-mRNA ceRNA Network for Cervical Cancer

Association of Biobehavioral Factors with Non-Coding RNAs in Cervical Cancer.

Integrated Analysis of Long Non‑coding RNA Competing Interactions Revealed Potential Biomarkers in Cervical Cancer: Based on a Public Database.

Comprehensive Analysis of lncRNA and miRNA Regulatory Network Reveals Potential Prognostic Non-coding RNA Involved in Breast Cancer Progression

Comprehensive Analysis of Differentially Expressed Profiles of Lncrnas and Construction of Mir-133B Mediated Cerna Network in Colorectal Cancer

Genome-Wide Analysis of the Expression of Circular RNA Full-Length Transcripts and Construction of the circRNA-miRNA-mRNA Network in Cervical Cancer

Long Non-Coding RNA in Cervical Cancer: from Biology to Therapeutic Opportunity

Comprehensive analysis of lncRNA-mRNA co-expression networks in HPV-driven cervical cancer reveals the pivotal function of LINC00511-PGK1 in tumorigenesis

How Long Non-Coding RNAs and MicroRNAs Mediate the Endogenous RNA Network of Head and Neck Squamous Cell Carcinoma: a Comprehensive Analysis

Construction and Integrated Analysis of a Lncrna-Associated Competing Endogenous RNA Network Reveal Functional Lncrnas in Pancreatic Cancer

Comprehensive Analysis of Lncrna–mirna– Mrna Network Ascertains Prognostic Factors in Patients with Colon Cancer

A Cerna Network of BBOX1-AS1-hsa-miR-125b-5p/hsa-miR-125a-5p-CDKN2A Shows Prognostic Value in Cervical Cancer

Comprehensive Analysis Of Long Non-Coding Rna-Messenger Rna-Microrna Co-Expression Network Identifies Cell Cycle-Related Lncrna In Hepatocellular Carcinoma

The Construction and Analysis of Cerna Networks in Invasive Breast Cancer: a Study Based on the Cancer Genome Atlas.