ORIGINAL RESEARCH article
Identification of mRNA-, circRNA- and lncRNA- Associated ceRNA Networks and Potential Biomarkers for Preeclampsia From Umbilical Vein Endothelial Cells
- 1International Peace Maternity and Child Health Hospital, School of Medicine, Shanghai Jiao Tong University, Shanghai, China
- 2Shanghai Key Laboratory of Embryo Original Disease, Shanghai, China
Objective: The etiology and pathogenesis of preeclampsia (PE) remain unclear, and ideal biomarkers for the early detection of PE are scarce. The involvement of the competing endogenous RNA (ceRNA) hypothesis in PE is only partially understood. The present study aimed to delineate a regulatory network in PE comprised of messenger RNAs (mRNAs), circular RNAs (circRNAs), long non-coding RNAs (lncRNAs), and microRNAs (miRNAs) via ceRNA profiles from human umbilical vein endothelial cells (HUVECs) to further reveal the pathogenesis of PE and potential biomarkers.
Methods: Differentially expressed mRNAs, circRNAs, and lncRNAs were detected in HUVECs from early onset preeclampsia (EOPE) cases (n = 4) and normal pregnancies (n = 4) by microarray analysis. Bioinformatics analysis was performed to systematically analyze the data, and a relevant ceRNA network was constructed. RNAs (ANGPT2, LIPG, hsa_circ_0025992, hsa_circ_0090396, hsa_circ_0066955, hsa_circ_0041203, hsa_circ_0018116, lnc-C17orf64-1:1, lnc-SLC27A2-2:1, and lnc-UEVLD-5:1) were validated by quantitative real-time PCR (qRT-PCR) in 10 pairs of HUVECs and placental tissues from PE patients and normal pregnancies. Furthermore, expression of hsa_circ_0025992 was detected in maternal peripheral blood samples from PE patients (n = 24) and normal pregnancies (n = 30) to confirm its potential as a novel biomarker. The receiver operating characteristic (ROC) curve was applied to analyze its diagnostic value.
Results: Compared with HUVECs from normal pregnancies, HUVECs from EOPE cases had 33 differentially expressed mRNAs (DEmRNAs), 272 DEcircRNAs, and 207 DElncRNAs. GO and KEGG analyses of the DERNAs revealed the biological processes and pathways involved in PE. Based on the microarray data and the predicted miRNAs, a ceRNA network was constructed with four mRNAs, 34 circRNAs, nine lncRNAs, and 99 miRNAs. GO and KEGG analyses of the network reinforced the crucial roles of metabolic disorders, the p53 and JAK/STAT signaling pathways in PE. In addition, ROC analysis indicated that hsa_circ_0025992 could be used as a novel biomarker for PE.
Conclusion: A novel ceRNA network was revealed in PE, and the potential of hsa_circ_0025992 to serve as a new biomarker was confirmed.
Preeclampsia (PE), a systemic vascular disorder specific to pregnancy, is traditionally characterized by hypertension and proteinuria after 20 weeks of gestation (ACOG, 2020). PE affects approximately 3–5% of pregnancies, leading to maternal and fetal morbidity and mortality (Mol et al., 2016). The etiology of PE remains largely unknown; however, it is widely acknowledged that inadequate trophoblast invasion, spiral uterine artery remodeling, endothelial dysfunction, and inappropriate angiogenesis could participate in PE (Ives et al., 2020). Due to the lack of effective treatments, it is essential to identify biomarkers for early diagnosis and therapeutic targets for the precise treatment of preeclampsia (Maric-Bilkan et al., 2019). Endothelial dysfunction is a crucial pathological feature of preeclampsia, which could lead to proteinuria and coagulation dysfunction in patients (Ives et al., 2020). In addition, preeclampsia is regarded as a risk factor for long-term postpartum cardiovascular disease (CVD) (Benschop and Duvekot, 2019; Lane-Cordova et al., 2019). Growing evidence suggests that vascular endothelial dysfunction is critical for the pathogenesis of preeclampsia (Cubro et al., 2020; Mauro et al., 2020; Vicen et al., 2020). Immunological changes (Aneman et al., 2020), metabolic syndrome (Scioscia et al., 2015; Vicen et al., 2020), and endothelial dysfunction are intertwined in preeclampsia. Nevertheless, the exact mechanism of endothelial dysfunction in preeclampsia still needs to be fully elucidated.
The competing endogenous RNA (ceRNA) hypothesis surmises an interaction network between transcripts that share the same microRNA response elements (MREs) (Salmena et al., 2011). Long non-coding RNAs (lncRNAs) are transcripts longer than 200 nt that can regulate neighboring and overlapping coding gene expression in various manners, such as cis- and trans-regulation (Chen, 2016). Circular RNAs (circRNAs) are novel single-stranded and covalently closed molecules that negatively affect miRNAs by acting as miRNA sponges (Li et al., 2018). The functions of lncRNAs and circRNAs might be embodied in the corresponding mRNA genes and host genes. Additionally, lncRNAs and circRNAs can function as ceRNAs and regulate gene expression (Chen, 2016; Li et al., 2018).
Increasing evidence has revealed that lncRNAs and circRNAs are significantly involved in the pathogenesis of PE through the placental ceRNA networks. Linc00511 was reported to exert its inhibitory effect on trophoblast invasion through the miR-29b-3p/Cyr61 axis (Quan et al., 2019). LncRNA AK002210 was found to participate in the progression of PE by affecting the miR-590-3p/NAIP and ERK/MMP signaling pathways (Zhao L. et al., 2020). The higher expression of lncRNA PSG10P contributed to PE pathogenesis by regulating miR-19a-3p and IL1RAP (Wang et al., 2019). CircTRNC18 and circZDHHC20 sponge miR-762 and miR-144, respectively, affect the expression of GRHL2 and repress migration of trophoblast cells (Shen et al., 2019; Zhou et al., 2020). Due to their unique circular structure, circRNAs are resistant to ribonucleases and are more stable in body fluids or exosomes than linear RNAs. Moreover, circRNAs are abundant and specifically expressed in tissues and developmental stages (Li et al., 2018). Herein, circRNAs exhibit ideal potential for use as biomarkers for diagnosis. Hsa_circ_0036877 was identified as a ceRNA and proved to be a potential blood biomarker for PE after receiver operating characteristic (ROC) curve analysis (Hu et al., 2018). To our knowledge, differentially expressed non-coding RNAs were often identified from placental tissues, maternal blood plasma samples, and umbilical cord blood plasma samples, and could affect the angiogenesis ability of the HUVECs in vitro (Xu et al., 2017; Lip et al., 2020; Zhao X. et al., 2020). A study revealed HUVEC-originated miRNAs in PE and suggested that the PE-downregulated miR-29a/c-3p may inhibit endothelial cell migration through the FGF2-activated PI3K-AKT signaling pathway (Zhou et al., 2017). However, there was no research on HUVEC-originated ceRNA networks in PE.
To offer evidence of the relationship between endothelial dysfunction and preeclampsia, our study applied ceRNA microarray analysis in primary HUVECs extracted from PE patients and normal pregnancies. Then, we tried to illuminate the lncRNA/circRNA–miRNA–mRNA network with bioinformatics analysis. Furthermore, a genetic biomarker of PE in the blood needs to be identified.
Materials and Methods
Study Design and Sample Collection
This research was divided into two parts. First, the microarray assay was conducted with four pairs of HUVECs extracted from patients with EOPE and normal pregnancies, and 10 pairs of HUVECs and placental tissues were collected for validation. Second, a preliminary delineation of the ceRNA profiling and a prospective study on hsa_circ_0025992 in maternal whole peripheral blood from PE patients (n = 24) and normal pregnant women (n = 30) were performed to predict PE at 11–12 weeks of gestation. HUVECs and placental tissues were collected from participants who had cesarean section. Maternal whole peripheral blood samples were collected in the first trimester and selected after the diagnosis of PE. Placental tissues and blood samples were stored at −80°C until use. All samples were collected with participants’ informed consent at the International Peace Maternity and Child Health Hospital, Shanghai, China. The research was authorized and approved by the hospital ethics committee. The diagnosis of preeclampsia was based on the guidelines from the American College of Obstetricians and Gynecologists (ACOG) (ACOG, 2020), and participants with multiple gestations, chronic hypertension, gestational diabetes mellitus, thyroid dysfunctions, and kidney diseases were eliminated. The clinical characteristics are described in Tables 1, 2.
Primary HUVEC Isolation and Identification
Primary HUVECs were isolated from umbilical cords within 2 h after cesarean section under aseptic conditions. After washing the umbilical vein with PBS (HyClone, United States), 30 ml of 0.1% type I collagenase solution (Gibco, United States) was injected into the umbilical vein. Both extremities of the vein were clamped, and the cord was incubated at 37°C and 5% CO2 for 15 min. Subsequently, the cord was squeezed, the contents were transferred to a sterile 50 ml tube, and the complete endothelial cell growth medium (ECM; ScienCell, United States) was added. The tubes were centrifuged at 1000 rpm for 10 min. Finally, the cells were suspended and transferred to a culture dish precoated with 2% gelatin (Sigma-Aldrich, United States). When 50–60% confluence was achieved, the cells were identified by an immunofluorescence assay with anti-CD31 (CST, United States), anti-factor VIII (Bioss, United States), and DAPI (Yeasen, China) staining. CD31 and factor VIII expression was positive in HUVECs. The identification results were shown in Supplementary Figure 1.
RNA Purification and Labeling
Total RNA from primary HUVECs was extracted and purified utilizing mirVanaTM miRNA Isolation Kit (Ambion, United States) according to the manufacturer’s instructions. The concentration of RNA was detected by a NanoDrop ND-2000 spectrophotometer (Thermo Fisher Scientific, United States), and RNA integrity was evaluated by an Agilent Bioanalyzer 2100 (Agilent, United States) with the RNA integrity numbers (RINs) ≥7.0 and 28 s/18 s ≥0.7. Qualified total RNA was further purified by a RNeasy Mini Kit (Qiagen, Germany) and RNase-Free DNase Set (Qiagen, Germany). Total RNA was amplified and labeled using a Low Input Quick Amp WT Labeling Kit (Agilent, United States). In addition, a RNeasy Mini Kit was used to purify the labeled cRNAs.
Microarray Assay and Data Acquisition
The SBC Human (4 × 180K) ceRNA array V1.0 (SHBIO, China) assay was used to detect differentially expressed mRNA, circRNA, and lncRNA profiles in HUVECs by comparing normal pregnancies with PE patients (n = 4). According to the manufacturer’s instructions, each slide was hybridized with 1.65 μg of Cy3-labeled cRNA for 17 h using a Gene Expression Hybridization Kit (Agilent, United States) in a hybridization oven (Agilent, United States). Then, the slides were washed in staining dishes (Thermo Fisher Scientific, United States) by Gene Expression Wash Buffer Kit (Agilent, United States) and scanned exploiting an Agilent Microarray Scanner (Agilent, United States) with the default settings. Raw data were obtained by Feature Extraction software 12.0 (Agilent, United States) and normalized with the quantile algorithm of the limma package in R1. The data were analyzed with Student’s t-test, | fold change| (| FC|) ≥ 2 and P-value < 0.05 were set as the cutoff criteria.
Hierarchical Clustering Analysis
Hierarchical clustering was conducted by the Cluster and Tree View program in R1 software to identify and visualize the overview of the DERNA profiles based on the microarray dataset.
LncRNA/circRNA Targeting Genes and MREs Prediction
Long non-coding RNAs and their potential target genes were visualized based on UCSC2. Cis-regulated genes were defined as protein-coding genes located on the same chromosome within a 10 kb region upstream or downstream of the corresponding lncRNAs. Trans-regulated genes were selected by two steps: (1) Complementary or similar sequences were first selected by BLAST3 from the same species database; (2) The complementary energy between the two sequences was calculated by using RNAplex, and the genes with E ≤ −30 were selected as the trans-regulated genes of the lncRNAs. The host genes of the circRNAs obtained from the CircBase4 were regarded as their target genes. And the putative MREs of differentially expressed lncRNAs and circRNAs were identified by bioinformatics analysis with Miranda5.
GO and KEGG Pathway Analysis
Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed with Fisher’s exact-test by the clusterProfiler package in R1 software. The standard of selection was the number of genes that fell on a GO term or pathway ≥2, with P-value < 0.05. Significantly enriched GO terms and KEGG pathways are presented.
CeRNA Network Construction
Based on the ceRNA hypothesis, combined with the microarray data (DElncRNAs, DEcircRNAs, DEmRNAs, and lncRNA/circRNA targeting genes) and the intersecting predicted miRNAs, the construction of the ceRNA network was visualized and analyzed via Cytoscape_V2_8_3.
Quantitative Real-Time PCR Analysis
RNA (800 ng) was reverse transcribed by cDNA synthesis kit (TRANS, China) following the protocol. Specific gene primers in human (Table 3) were designed by NCBI Primer-BLAST3. The subsequent qRT-PCR (10 μl volume) was accomplished by a QuantStudio 7 Flex instrument (Thermo Fisher Scientific, United States) with a SYBR Premix Ex Taq Kit (Takara, Japan). Relative RNA expression was calculated by the 2–ΔΔCt method and normalized to 18S rRNA expression.
Continuous data were presented as means ± standard deviations (SD) unless specifically otherwise indicated. Statistically significant differences were estimated by Student’s t-test or the Mann–Whitney U test. The receiver operating characteristic (ROC) curve analysis and Spearman correlation analysis were used to evaluate the predictive value of the hsa_circ_0025992 and the correlation of the circRNA with blood pressure of the patients. All statistical analyses were established by GraphPad Prism 7.0 and SPSS version 20.0. P-value < 0.05 was considered as statistically significant.
Identification of DERNAs
Microarray analysis was performed to explore the pathogenesis of preeclampsia from the perspective of endothelial dysfunction. A total DERNA profile was described after stringent filtering with the following standards: absolute fold change ≥2.0 and P-value < 0.05. The analysis demonstrated that 33 DEmRNAs (10 upregulated and 23 downregulated mRNAs, Figure 1A), 272 DEcircRNAs (38 upregulated and 234 downregulated circRNAs, Figure 1B), and 207 DElncRNAs (51 upregulated and 156 downregulated lncRNAs, Figure 1C) were identified. The total microarray-detected RNAs were widely distributed on all human chromosomes, and most of the significantly dysregulated RNAs were downregulated (Figure 2).
Figure 1. Hierarchical clustering analysis of DERNAs in HUVECs between normal pregnancies and PE patients. The heatmap reveals 33 DEmRNAs (A), 272 DEcircRNAs (B), and 207 DElncRNAs (C). Each column and row correspond to a sample and a transcript, respectively. Red represents upregulation, while blue represents downregulation. CON, control; PE, preeclampsia, n = 4.
Figure 2. The distribution of DERNAs on chromosomes. Circos plot of genome-wide mRNAs (A), circRNAs (B), and lncRNAs (C). Histograms depict the distribution of DEmRNAs (D), DEcircRNAs (E), and DElncRNAs (F) mapped on each chromosome. The outermost circle represents the human chromosomes, while the inside eight circles represent eight samples from control and PE patients. The Y-axis denotes the count of transcripts mapped to the chromosome. CON, control; PE, preeclampsia, n = 4.
The variability between PE samples and controls in the microarray was shown in the volcano plot (Supplementary Figure 2). These aforementioned results illustrated that RNA expression in the PE group was different from that in the control group.
GO Enrichment Analysis of DERNAs
To gain further insight into the pathogenesis that was potentially mediated by these DERNAs in PE, GO analysis was performed. Focusing on the biological process, it was revealed that the DEmRNAs mainly affect the regulation of growth and cellular processes through “negative regulation of cell growth,” “positive/negative regulation of growth,” and “negative regulation of locomotion” in Figure 3A; it was suggested that circRNAs might participate in the “establishment of the endothelial barrier,” “de novo posttranslational protein folding,” “antigen processing and presentation,” and “mitotic process” in Figure 3B; the cis-regulated genes of lncRNAs played a pivotal role in “endothelial and epithelial cell migration” in Figure 3C, and the trans-regulated genes of lncRNAs could alert “lymphocyte apoptotic process,” “neurotransmitter uptake,” and “T cell-mediated cytotoxicity” in Figure 3D. Most of the aforementioned biological processes indicated that immune injury, intercellular signaling, and endothelial dysfunction were associated with PE.
Figure 3. Gene Ontology (GO) enrichment of DERNAs. GO enrichment of DEmRNAs (A), host genes of DEcircRNAs (B), cis-regulated genes of DElncRNAs (C), and trans-regulated genes of DElncRNAs (D).
KEGG Pathway Analysis of DERNAs
Kyoto encyclopedia of genes and genomes pathway enrichment analysis for significant DERNAs was executed to comprehend pathways and molecular interactions relevant to genes. As shown in Figure 4A, the top three pathways for DEmRNAs were “drug metabolism-other enzymes,” “malaria,” and “pathogenic Escherichia coli infection”. In Figure 4B, the top three pathways for the host genes of DEcircRNAs were “lysine biosynthesis,” “phenylalanine, tyrosine, and tryptophan biosynthesis,” and “thyroid cancer”. In Figure 4C, the top three pathways for cis-regulated genes of DElncRNAs were “RNA polymerase,” “primary immunodeficiency,” and “one-carbon pool by folate”; and in Figure 4D, the top three pathways for trans-regulated genes of DElncRNAs were “sulfur metabolism,” “ubiquinone and other terpenoid-quinone biosynthesis,” and “tryptophan metabolism”. Major nutrient metabolism disorders were also revealed in Figures 4A,C,D. Collectively, insights into the correlation of immunodeficiency, metabolic disorders, and bacterial dysbiosis with endothelial dysfunction in preeclampsia were proposed.
Figure 4. KEGG pathway analysis of DERNAs. KEGG pathway analyses of DEmRNAs (A), host genes of DEcircRNAs (B), cis-regulated genes of DElncRNAs (C), and trans-regulated genes of DElncRNAs (D).
Validation of mRNAs, circRNAs, and lncRNAs Expression
To validate the microarray profiling results, two mRNAs (ANGPT2 and LIPG), five circRNAs (hsa_circ_0025992, hsa_circ_0090396, hsa_circ_0066955, hsa_circ_0041203, and hsa_circ_0018116), and three lncRNAs (lnc-C17orf64-1:1, lnc-SLC27A2-2:1, and lnc-UEVLD-5:1) were selected for verification by quantitative real-time PCR (qRT-PCR) in the same HUVEC samples. The results in Figure 5A showed high consistency between the microarray and qRT-PCR analyses. The expression of the selected DERNAs in HUVECs and corresponding placental tissues from normal participants (n = 10) and PE patients (n = 10) was determined. The qRT-PCR results in Figures 5B,C showed that the selected genes were decreased in both HUVECs and placental tissues from the PE group; Lnc-C17orf64-1:1 expression in HUVECs and placenta of the PE group was significantly lower than that in the normal group (Figures 5B,C), and the expression levels of hsa_circ_0025992, hsa_circ_0090396, hsa_circ_0066955, hsa_circ_0041203, lnc-C17orf64-1:1, and lnc-UEVLD-5:1 were significantly decreased in placental tissues (Figure 5C). These findings indicated that some DERNAs could be implicated in the pathogenesis of PE.
Figure 5. Validation of DERNA expression by qRT-PCR. Differentially expressed mRNAs, circRNAs, and lncRNAs in PE samples by microarray analysis were validated by qRT-PCR (A), n = 4. Validation of differentially expressed mRNAs, circRNAs, and lncRNAs in 10 pairs of HUVECs (B) and corresponding placental tissues (C) of participants. | FC| = | Fold Change|. CON, control; PE, preeclampsia. Data were shown in mean ± SEM. *P < 0.05; **P < 0.01.
Global ceRNA Network Integration in PE
According to the ceRNA hypothesis, lncRNAs could compete with circRNAs for the same miRNAs and further impact gene expression. Combined with our microarray data and predicted miRNAs, a portion of the ceRNA network was constructed by Cytoscape_V2_8_3 and visualized in Figure 6A. Four key genes (CACNG8, LIPG, ACAA2, and OSMR) were identified in the network. To further uncover insights into the ceRNA network, GO, and KEGG enrichment analyses were performed as shown in Figures 6B,C. PE-associated pathways, such as the “p53 signaling pathway,” “JAK-STAT signaling pathway,” “fatty acid degradation,” and “autophagy-animal” were depicted (Gao et al., 2016; Zhao and Jiang, 2018; Nakashima et al., 2020; Khaire et al., 2021).
Figure 6. Global ceRNA network integration in PE. A portion of the theoretical ceRNA network in PE was predicted according to the microarray analysis (A). GO analysis of the molecules involved in the ceRNA network (B). KEGG pathway analysis of the molecules involved in the ceRNA network (C).
Hsa_circ_0025992 Functions as a Novel Potential Blood Biomarker
According to the microarray results and bioinformatics analysis, hsa_circ_0025992 showed a significant change, and it was predicted to sponge a reported angiogenic microRNA, hsa-miR-20a-3p (Supplementary Table 1). To confirm the potential of hsa_circ_0025992 as a biomarker, its expression in maternal whole peripheral blood was detected by qRT-PCR in normal pregnancies (n = 30) and PE patients (n = 24) at 11–12 weeks of gestation. Compared to normal patients, PE patients who were diagnosed based on ACOG at the end of gestation had a significantly higher ΔCt value of hsa_circ_0025992, indicating its lower expression level in blood samples (Figure 7A). Furthermore, receiver operating characteristic (ROC) curve analysis was performed, and the area under the ROC curve (AUC) of hsa_circ_0025992 was 0.8065 [95% confidence interval (CI) 0.6902–0.9210], P < 0.001. The cutoff value determined by the maximum Youden Index (YI) was 16.36 (ΔCt value). The sensitivity and specificity were 54.17% and 93.33%, respectively. The measurement of a ΔCt value < 16.36 demonstrated a 71.8% negative predictive value (NPV) for the diagnosis of preeclampsia. A ΔCt value above 16.36 had a relatively high positive predictive value (PPV) of 86.7% (Figure 7B). Spearman correlation analysis depicted that the ΔCt value of hsa_circ_0025992 positively correlated with the systolic blood pressure and diastolic blood pressure of the corresponding participants (Figures 7C,D). These results revealed that hsa_circ_0025992 could serve as a novel blood biomarker for the prediction of early PE.
Figure 7. Hsa_circ_0025992 serving as a potential novel blood biomarker for the prediction of early PE. Expression of hsa_circ_0025992 in blood samples from PE patients (n = 24) and normal pregnant women (n = 30) (A). Data were shown as the mean ± SEM. The ROC curve of hsa_circ_0025992 as a predictive biomarker for early PE. AUC = 0.8065 (95% CI 0.6902–0.9210, P = 0.0001); sensitivity and specificity, 54.17 and 93.33%, respectively. The cutoff value was 16.36 (ΔCt value) [the optimal cutoff point was determined at the maximum Youden Index (YI)]. The NPV of a ΔCt value ≤ 16.36 for a diagnosis of preeclampsia was 71.8%. A ΔCt value above 16.36 had a PPV of 86.7% (B). Spearman correlation analysis for hsa_circ_0025992 (ΔCt value) in blood samples with the systolic blood pressure and diastolic blood pressure of the corresponding participants (C,D). *P < 0.05, **P < 0.01, ***P < 0.001.
With the rapid development of high-throughput technology and bioinformatics analysis, increasing efforts have been made to uncover the pathology of PE, and the novel ceRNA hypothesis is becoming more noteworthy (Hu et al., 2018; Quan et al., 2019). To the best of our knowledge, there has been scarce research aimed at the ceRNA network in HUVECs to illustrate the mechanism of endothelial dysfunction in PE.
In our study, for the first time, ceRNA profiling of EOPE HUVECs was revealed, and further bioinformatics analysis was plotted to clarify the deregulated biological processes or pathways affected by ceRNAs. A set of DERNAs were identified, validated, and found to be distributed across all chromosomes, denoting that non-coding RNAs are widely implicated in gene regulation. Subsequently, GO enrichment of DERNAs suggested that immune injury substantially linked endothelial dysfunction to PE and thus might affect intercellular signaling; KEGG enrichment revealed that the microbiome and metabolic disorders had a profound influence on endothelial dysfunction and PE. These findings were consistent with those of previous studies on the mechanism of endothelial dysfunction in PE (Ahmadian et al., 2020; Aneman et al., 2020; Aplin et al., 2020; Vicen et al., 2020). In addition, it was remarkable that several reported PE-associated pathways, including the “Wnt signaling pathway” (Schlosser et al., 2020), “Toll-like receptor signaling pathway” (He et al., 2018), “mTOR signaling pathway” (He et al., 2019), “p53 signaling pathway” (Gao et al., 2016), and “JAK-STAT signaling pathway” (Zhao and Jiang, 2018) were also depicted.
Based on the microarray data and predicted miRNAs, a portion of the ceRNA network was constructed, and GO and KEGG analyses of network-related molecules were performed. The enrichment results emphasized metabolic disorders in PE, in agreement with recent cohort research on metabolomic alterations in postpartum PE plasma (Hong X. et al., 2020); therefore, more research could be done to facilitate the mechanism of PE and women’s long-term cardiovascular disease. In addition, the importance of the p53 and JAK/STAT signaling pathways was reconfirmed. There is no conflict with the former study that proved that p53 upregulation accelerated G1 arrest and apoptosis of HUVECs from preeclampsia patients (Gao et al., 2016), and a placenta microarray analysis also revealed the overrepresentation of the JAK/STAT pathway in PE (Medina-Bastidas et al., 2020), but little direct evidence explained the dysfunction of HUVECs from PE through the JAK/STAT pathway.
Moreover, the following hub genes: CACNG8, LIPG, ACAA2, and OSMR, were shown to be involved in the ceRNA network. CACNG8, which was confirmed to encode calcium channels implicated in cardiac contraction, was identified as a calcium flux-related gene that depresses ventricular function (Ortega et al., 2015) and might thus account for long-term CVD in preeclampsia patients. LIPG can be secreted by endothelial cells and is involved in lipoprotein metabolism. Cytokines greatly affect LIPG expression in endothelial cells and are involved in the inflammatory response in various diseases (Hong C. et al., 2020). ACAA2 plays an essential role in fatty acid metabolism as a key enzyme of fatty acid oxidation (Zhang et al., 2018). OSMR is a member of the interleukin 6 receptor family and mainly signals through the JAK/STAT pathway (Hermanns, 2015). Nevertheless, few studies have directly connected the mechanism of PE with these genes.
Finally, we proposed hsa_circ_0025992 as a novel potential biomarker for the early prediction of PE. According to the MREs prediction, the most significantly changed circRNA hsa_circ_0025992 might sponge hsa-miR-20a-3p, a reported angiogenic microRNA targeting vascular endothelial growth factor (Platania et al., 2019; Maisto et al., 2020), and could thus be involved in the endothelial dysfunction of PE. The area under the receiver operating characteristic (AUC) curve of hsa_circ_0025992 was 0.8065 [95% confidence interval (CI) 0.6902–0.9210], P < 0.001. Compared with hsa_circ_0036877 (Hu et al., 2018), hsa_circ_0025992 had a higher specificity of 93.33% and a ΔCt value above 16.36 with a relatively high PPV of 86.7%, whereas the sensitivity was relatively lower. A higher ΔCt value represents a lower expression level; thus, the hsa_circ_0025992 level in the blood was negatively correlated with systolic blood pressure and diastolic blood pressure. Taken together, hsa_circ_0025992 can be regarded as a potential biomarker, and the collaborative prediction of multiple factors may facilitate its predictive power.
The main shortcoming of this study was the limited cohort size due to time constraints and strict patient screening criteria. The validation of hsa_circ_0025992 as a predictive biomarker in more patients from multiple centers is necessary to explore its clinical value. Second, the relationship between hsa_circ_0025992 and hsa-miR-20a-3p needs further research; thus, we will conduct more experiments in the future to illustrate ceRNA function in vitro and in vivo.
Data Availability Statement
The datasets presented in this study can be found in the GEO repository (https://www.ncbi.nlm.nih.gov/geo/) with the accession number GSE165324.
The studies involving human participants were reviewed and approved by the International Peace Maternity and Child Health Hospital Ethics Committee. The patients/participants provided their written informed consent to participate in this study.
DC and BH analyzed microarray data and performed the experiments. SW, PZ, XZ, and JL provided the material support. DC and XY wrote the manuscript. WC and XY conceived the study and contributed valuable advice to the manuscript. All authors approved the manuscript.
This study was sponsored by the National Natural Science Foundation of China (81370727 and 81501256), the project support the guidance of medicine (Traditional Chinese and Western medicine) in Shanghai (18411963500), the project about important disease of Xuhui District Health Planning Commission (XHLHGG201805), the interdisciplinary program of Shanghai Jiao Tong University (ZH2018ZDA31), and Shanghai Sailing Program (20YF 1453800).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2021.652250/full#supplementary-material
- ^ http://cran.r-project.org/
- ^ http://genome.ucsc.edu/
- ^ https://blast.ncbi.nlm.nih.gov/
- ^ http://www.circbase.org/
- ^ http://www.miranda.org/
Ahmadian, E., Rahbar Saadat, Y., Hosseiniyan Khatibi, S. M., Nariman-Saleh-Fam, Z., Bastami, M., Zununi Vahed, F., et al. (2020). Pre-Eclampsia: microbiota possibly playing a role. Pharmacol. Res. 155:104692. doi: 10.1016/j.phrs.2020.104692
Aneman, I., Pienaar, D., Suvakov, S., Simic, T. P., Garovic, V. D., and McClements, L. (2020). Mechanisms of Key Innate Immune Cells in Early- and Late-Onset Preeclampsia. Front. Immunol. 11:1864. doi: 10.3389/fimmu.2020.01864
Benschop, L., and Duvekot, J. J. (2019). Roeters van Lennep JE. Future risk of cardiovascular disease risk factors and events in women after a hypertensive disorder of pregnancy. Heart 105, 1273–1278. doi: 10.1136/heartjnl-2018-313453
Cubro, H., Nath, K. A., Suvakov, S., Garcia-Valencia, O., Parashuram, S., White, W. M., et al. (2020). Mechanisms of vascular dysfunction in the interleukin-10-deficient murine model of preeclampsia indicate nitric oxide dysregulation. Kidney Int. 99, 646–656. doi: 10.1016/j.kint.2020.09.034
Gao, Q., Zhu, X., Chen, J., Mao, C., Zhang, L., and Xu, Z. (2016). Upregulation of P53 promoted G1 arrest and apoptosis in human umbilical cord vein endothelial cells from preeclampsia. J. Hypertens. 34, 1380–1388. doi: 10.1097/hjh.0000000000000944
He, B., Yang, X., Li, Y., Huang, D., Xu, X., Yang, W., et al. (2018). TLR9 (Toll-Like Receptor 9) Agonist Suppresses Angiogenesis by Differentially Regulating VEGFA (Vascular Endothelial Growth Factor A) and sFLT1 (Soluble Vascular Endothelial Growth Factor Receptor 1) in Preeclampsia. Hypertension 71, 671–680. doi: 10.1161/hypertensionaha.117.10510
He, C., Shan, N., Xu, P., Ge, H., Yuan, Y., Liu, Y., et al. (2019). Hypoxia-induced Downregulation of SRC-3 Suppresses Trophoblastic Invasion and Migration Through Inhibition of the AKT/mTOR Pathway: implications for the Pathogenesis of Preeclampsia. Sci. Rep. 9:10349.
Hong, X., Zhang, B., Liang, L., Zhang, Y., Ji, Y., Wang, G., et al. (2020). Postpartum plasma metabolomic profile among women with preeclampsia and preterm delivery: implications for long-term health. BMC Med. 18:277. doi: 10.1186/s12916-020-01741-4
Hu, X., Ao, J., Li, X., Zhang, H., Wu, J., and Cheng, W. (2018). Competing endogenous RNA expression profiling in pre-eclampsia identifies hsa_circ_0036877 as a potential novel blood biomarker for early pre-eclampsia. Clin. Epigenetics 10:48.
Ives, C. W., Sinkey, R., Rajapreyar, I., Tita, A. T. N., and Oparil, S. (2020). Preeclampsia-Pathophysiology and Clinical Presentations: JACC State-of-the-Art Review. J. Am. Coll. Cardiol. 76, 1690–1702. doi: 10.1016/j.jacc.2020.08.014
Lane-Cordova, A. D., Khan, S. S., Grobman, W. A., Greenland, P., and Shah, S. J. (2019). Long-Term Cardiovascular Risks Associated With Adverse Pregnancy Outcomes: JACC Review Topic of the Week. J. Am. Coll. Cardiol. 73, 2106–2116. doi: 10.1016/j.jacc.2018.12.092
Lip, S. V., Boekschoten, M. V., Hooiveld, G. J., van Pampus, M. G., Scherjon, S. A., Plösch, T., et al. (2020). Early-onset preeclampsia, plasma microRNAs, and endothelial cell function. Am. J. Obstet. Gynecol. 222, 497.e1–497.e12.
Maisto, R., Trotta, M. C., Petrillo, F., Izzo, S., Cuomo, G., Alfano, R., et al. (2020). Resolvin D1 Modulates the Intracellular VEGF-Related miRNAs of Retinal Photoreceptors Challenged With High Glucose. Front. Pharmacol. 11:235. doi: 10.3389/fphar.2020.00235
Maric-Bilkan, C., Abrahams, V. M., Arteaga, S. S., Bourjeily, G., Conrad, K. P., Catov, J. M., et al. (2019). Research Recommendations From the National Institutes of Health Workshop on Predicting, Preventing, and Treating Preeclampsia. Hypertension 73, 757–766. doi: 10.1161/hypertensionaha.118.11644
Mauro, A. K., Khurshid, N., Berdahl, D. M., Ampey, A. C., Adu, D., Shah, D. M., et al. (2020). Cytokine concentrations direct endothelial function in pregnancy and preeclampsia. J. Endocrinol. 248, 107–117. doi: 10.1530/joe-20-0397
Medina-Bastidas, D., Guzmán-Huerta, M., Borboa-Olivares, H., Ruiz-Cruz, C., Parra-Hernández, S., Flores-Pliego, A., et al. (2020). Placental Microarray Profiling Reveals Common mRNA and lncRNA Expression Patterns in Preeclampsia and Intrauterine Growth Restriction. Int. J. Mol. Sci. 21:3597. doi: 10.3390/ijms21103597
Nakashima, A., Cheng, S.-B., Ikawa, M., Yoshimori, T., Huber, W. J., Menon, R., et al. (2020). Evidence for lysosomal biogenesis proteome defect and impaired autophagy in preeclampsia. Autophagy 16, 1771–1785. doi: 10.1080/15548627.2019.1707494
Ortega, A., Tarazón, E., Roselló-Lletí, E., Gil-Cayuela, C., Lago, F., González-Juanatey, J.-R., et al. (2015). Patients with Dilated Cardiomyopathy and Sustained Monomorphic Ventricular Tachycardia Show Up-Regulation of KCNN3 and KCNJ2 Genes and CACNG8-Linked Left Ventricular Dysfunction. PLoS One 10:e0145518. doi: 10.1371/journal.pone.0145518
Platania, C. B. M., Maisto, R., Trotta, M. C., D’Amico, M., Rossi, S., Gesualdo, C., et al. (2019). Retinal and circulating miRNA expression patterns in diabetic retinopathy: an in silico and in vivo approach. Br. J. Pharmacol. 176, 2179–2194.
Quan, X., Zhao, M., Yang, X., Zhu, Y., and Tian, X. (2019). AP2γ mediated downregulation of lncRNA LINC00511 as a ceRNA suppresses trophoblast invasion by regulating miR-29b-3p/Cyr61 axis. Biomed. Pharmacother. 120:109269. doi: 10.1016/j.biopha.2019.109269
Schlosser, K., Kaur, A., Dayan, N., Stewart, D. J., Pilote, L., and Delles, C. (2020). Circulating miR-206 and Wnt-signaling are associated with cardiovascular complications and a history of preeclampsia in women. Clin. Sci. 134, 87–101. doi: 10.1042/cs20190920
Scioscia, M., Karumanchi, S. A., Goldman-Wohl, D., and Robillard, P.-Y. (2015). Endothelial dysfunction and metabolic syndrome in preeclampsia: an alternative viewpoint. J. Reprod. Immunol. 108, 42–47. doi: 10.1016/j.jri.2015.01.009
Shen, X. Y., Zheng, L. L., Huang, J., Kong, H. F., Chang, Y. J., Wang, F., et al. (2019). CircTRNC18 inhibits trophoblast cell migration and epithelial-mesenchymal transition by regulating miR-762/Grhl2 pathway in pre-eclampsia. RNA Biol. 16, 1565–1573. doi: 10.1080/15476286.2019.1644591
Vicen, M., Igreja, S. áI. C., Tripská, K., Vitverová, B., Najmanová, I., Eissazadeh, S., et al. (2020). Membrane and soluble endoglin role in cardiovascular and metabolic disorders related to metabolic syndrome. Cell. Mol. Life Sci. doi: 10.1007/s00018-020-03701-w [Epub online ahead of print]Online ahead of print.]
Wang, N., Li, R., and Xue, M. (2019). Potential regulatory network in the PSG10P/miR-19a-3p/IL1RAP pathway is possibly involved in preeclampsia pathogenesis. J. Cell. Mol. Med. 23, 852–864. doi: 10.1111/jcmm.13985
Xu, Y., Ge, Z., Zhang, E., Zuo, Q., Huang, S., Yang, N., et al. (2017). The lncRNA TUG1 modulates proliferation in trophoblast cells via epigenetic suppression of RND3. Cell Death Dis. 8:e3104. doi: 10.1038/cddis.2017.503
Zhang, Y., Wang, Y., Wang, X., Ji, Y., Cheng, S., Wang, M., et al. (2018). Acetyl-coenzyme A acyltransferase 2 promote the differentiation of sheep precursor adipocytes into adipocytes. J. Cell. Biochem. 120, 8021–8031. doi: 10.1002/jcb.28080
Zhao, L., Liang, X., Ma, Y., Li, J., Liao, S., Chen, J., et al. (2020). AK002210 promotes the proliferation, migration and invasion of trophoblast cell through regulating miR-590/NAIP signal axis. Arch. Biochem. Biophys. 688:108366. doi: 10.1016/j.abb.2020.108366
Zhao, X., Li, Y., Chen, S., Chen, D., Shen, H., Yang, X., et al. (2020). Exosomal encapsulation of miR-125a-5p inhibited trophoblast cell migration and proliferation by regulating the expression of VEGFA in preeclampsia. Biochem. Biophys. Res. Commun. 525, 646–653. doi: 10.1016/j.bbrc.2020.02.137
Zhou, B., Zhang, X., Li, T., Xie, R., Zhou, J., Luo, Y., et al. (2020). CircZDHHC20 represses the proliferation, migration and invasion in trophoblast cells by miR-144/GRHL2 axis. Cancer Cell Int. 20:19.
Zhou, C., Zou, Q. Y., Li, H., Wang, R. F., Liu, A. X., Magness, R. R., et al. (2017). Preeclampsia Downregulates MicroRNAs in Fetal Endothelial Cells: roles of miR-29a/c-3p in Endothelial Function. J. Clin. Endocrinol. Metab. 102, 3470–3479. doi: 10.1210/jc.2017-00849
Keywords: preeclampsia, HUVECs, circRNAs, ceRNA network, biomarker
Citation: Chen D, He BW, Zheng PC, Wang SY, Zhao XY, Liu JY, Yang XY and Cheng WW (2021) Identification of mRNA-, circRNA- and lncRNA- Associated ceRNA Networks and Potential Biomarkers for Preeclampsia From Umbilical Vein Endothelial Cells. Front. Mol. Biosci. 8:652250. doi: 10.3389/fmolb.2021.652250
Received: 12 January 2021; Accepted: 23 March 2021;
Published: 20 April 2021.
Edited by:Mahendra Pratap Kashyap, University of Alabama at Birmingham, United States
Reviewed by:Xiaotian Li, Fudan University, China
Nanbert Zhong, New York State Institute for Basic Research in Developmental Disabilities, United States
Sanjay Rathod, University of Pittsburgh, United States
Copyright © 2021 Chen, He, Zheng, Wang, Zhao, Liu, Yang and Cheng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
†These authors have contributed equally to this work