Abstract
BACKGROUND This study aimed to identify significantly altered circRNAs/lncRNAs/miRNAs/mRNAs pathways in preeclampsia (PE), investigate their target relationships, and determine their biological functions. MATERIAL AND METHODS Base on RNA-seq technique and the GEO database, expression profiles of circRNAs/lncRNAs/miRNAs/mRNAs related to PE were obtained. Differentially expressed RNAs were determined using the Limma package in R. Gene set enrichment analysis (GSEA) was performed using GSEA software (v. 3.0) and illustrated by ClusterProfiler and ggplot2 package in R. DAVID database (v. 6.8) was implemented to analyze functional categories and the association between genes and the corresponding Gene Ontology (GO) classification. The R visualization package GOPlot was used to get a better visualization of the relationships between genes and the selected functional categories. CeRNA networks which visualized the correlations between circRNA/lncRNA-miRNA-mRNA were constructed using Cytoscape software (v. 3.6.0). Targetscan and miRanda database were used to predict target relationships between circRNA/lncRNA-miRNA-mRNA. QRT-PCR and luciferase reporter assay were used to verify the expression and target relationship of has_circ_0088196/LINC01492/miR-100-5p/LIF (leukemia inhibitory factor). RESULTS The jak-stat signaling pathway was activated and miR-100-5p was downregulated in PE compared with normal tissues both in collected placental tissue samples and GEO database. Upregulated LIF, LINC01492, and hsa_circ_0088196 were negatively correlated with miR-100-5p expression and had a targeted relationship with miR-100-5p. CONCLUSIONS miR-100-5p may suppress PE development, while LIF, LINC01492, and hsa_circ_0088196 may promote it though inhibiting miR-100-5p. The jak-stat signaling pathway was activated and involved in PE progression.