Abstract
Preeclampsia (PE) is a pregnancy-related disease defined as onset of hypertension and proteinuria after the 20th week of pregnancy, which causes most maternal and perinatal morbidity and mortality. Although placental dysfunction is considered as the main cause of PE, the exact pathogenesis of PE is not yet fully understood. Long non-coding RNAs (lncRNAs) are implicated in a broad range of physiological and pathological processes, including the occurrence of PE. In this study, we investigated the expression and functions of HIF-1α pathway–related lncRNA-HEIPP (high expression in PE placenta) in the pathogenesis of PE. The expression of lncRNA-HEIPP in the placenta from women who underwent PE was screened by lncRNA microarray and then verified using real-time polymerase chain reaction. Then, the methylation profile of the lncRNA-HEIPP promoter and the enrichment of H3K4me3 binding were assessed by bisulfite pyrosequencing and chromatin immunoprecipitation (ChIP)–quantitative polymerase chain reaction (qPCR) assay, respectively. It was found that the level of lncRNA-HEIPP in the PE placenta was significantly higher than that in normal placenta and was increased in HTR-8/SVneo human trophoblast cells upon hypoxia treatment. Moreover, we reported that H3K4me3 manifested significantly higher promoter occupancy on lncRNA-HEIPP promoter in HTR-8/SVneo cells upon hypoxia treatment and found that the downregulation of lncRNA-HEIPP promoted trophoblast invasion. Our findings suggested that the hypoxia-induced expression of lncRNA-HEIPP mediated by H3K4me3 modification in trophoblast may contribute to the pathogenesis of PE.
Introduction
Preeclampsia (PE) is a common complication of human pregnancy defined as the occurrence of hypertension and significant proteinuria after the 20th week of gestation, occurring in approximately 3 to 5% of all pregnancies and resulting in about 50,000 maternal deaths annually worldwide (). Preeclampsia is also the leading cause of fetal morbidity and mortality (). It is characterized by maternal arterial pressure >140/90 mmHg and proteinuria >0.3 g/24 h after 20 weeks’ gestation (). Current studies suggested that the pathogenesis of PE is closely related to immune system and genetic dysfunctions (; Yong et al., 2018). However, the dysfunction of placenta plays a central role in the pathogenesis of PE (; ).
Recent studies have shown that abnormal trophoblast differentiation, migration, and apoptosis often cause decreased invasion of trophoblast cells and remodeling of uterine spiral artery, resulting in reduced uteroplacental perfusion (; ; ). This further leads to hypoxia and inflammatory changes in the placenta, which will eventually lead to PE (; ). Recent studies have shown that the important phenotypes of PE, such as hypertension and proteinuria, are closely related with hypoxia inducible factor-1α (HIF-1α) (Wang et al., 2017; ). HIF-1α is a transcription factor that controls the expression of hypoxia-induced factors such as vascular endothelial growth factor (VEGF), platelet-derived growth factor, and so on (; Zhang et al., 2019), which are considered to play important roles during the pathogenesis of PE. Recently, it has been found that long non-coding RNA (lncRNA) could regulate the invasion and apoptosis of trophoblast cell lines (; Xu et al., 2018), suggesting their potential roles in PE occurrence. However, whether hypoxia and HIF-1α regulate the expression of specific lncRNAs during PE has not been investigated yet. So, the objective of the current study is to investigate the differentially expressed lncRNAs in human placenta between normal pregnancy and PE and dissect the mechanisms underlying the regulation of lncRNAs by hypoxia as well as their functions in trophoblast cells migration.
Materials and Methods
Patient Samples
All patients provided written informed consent, and this study was approved by the Human Research Ethics Committee of the Shanghai Changzheng Hospital at the Second Military Medical University. Thirty PE patients (mean age of 33.20 ± 0.46 years) and 30 control patients (mean age of 32.60 ± 0.47 years) who received a cesarean section were included in this study. Preeclampsia was clinically diagnosed before cesarean. The main criteria for PE diagnosis were as follows: systolic blood pressure ≥140 mm Hg or diastolic blood pressure ≥90 mm Hg on two occasions at least 4 h apar, and proteinuria >0.3 g/24 h. The control pregnancy was defined as no medical complications and proteinuria, and maternal blood pressure <140/90 mm Hg. Clinical characteristics and detailed information of all patients enrolled in this study are summarized in Table 1. The original demographic data of enrolled patients can be found in Supplementary Table 1. The placenta specimen was resected from the middle of villous lobule, avoiding any visible blood clot or calcification. To minimize blood contamination, each piece of tissue was intensively washed with ice-cold phosphate-buffered saline and then immediately snap-frozen in liquid nitrogen after resection. Of the 30 pairs of placenta samples, four pairs were randomly selected for lncRNA and mRNA microarray analysis, and the other samples were used for the following verification.
TABLE 1
| Normal (n = 30) | Preeclampsia (n = 30) | P-value | |
| Maternal age (y) | 32.60 ± 0.47 | 33.20 ± 0.46 | 0.368 |
| Gestational age (wk) | 38.35 ± 0.34 | 34.44 ± 0.15b | <0.001 |
| BP (mm Hg) | |||
| Systolic | 108.13 ± 0.81 | 163.87 ± 1.30b | <0.001 |
| Diastolic | 70.17 ± 0.85 | 104.43 ± 1.23b | <0.001 |
| Proteinuria (g/24 h) | N/A | 2.62 ± 0.20 | N/A |
| Fetal birth weight (g) | 3,278.00 ± 58.01 | 2,577.33 ± 97.73b | <0.001 |
Clinical characteristics of the pregnant women enrolled in this studya.
BP, blood pressure. aData are shown as mean ± SEM. bCompared with normal pregnancy: P < 0.001.
RNA Extraction
Total RNA was extracted from snap-frozen placenta tissues obtained from PE placenta patients and control patients using TRIzol Reagent (Invitrogen, Carlsbad, CA, United States) according to the manufacturer’s instruction. The RNA integrity and concentration were evaluated using the NanoDrop ND-1000 spectrophotometer and 2100 RNA Nano 6000 Assay Kit (Agilent Technologies, Santa Clara, CA, United States).
RNA Microarray and Computational Analysis
RNA with RNA integrity number greater than 6.5 purified from total RNA following the removal of rRNA was amplified and transcribed into fluorescent cDNA along the entire length of the transcripts without 3′ bias utilizing random priming method, and cDNA was labeled and hybridized to the Human LncRNA Array V2.0 (8 × 60, Arraystar). In addition, 30,215 coding transcripts were detected using the microarray as well. The microarray was performed by KangCheng Bio-tech (Shanghai, China). The arrays were scanned using the Agilent Scanner G2505B (Agilent Technologies), and the acquired array images were analyzed using Agilent Feature Extraction software (version 11.5.1; Agilent Technologies). Quantile normalization and subsequent data processing were performed using GeneSpring GX v11.5.1 software package (Agilent Technologies).
Gene Ontology Pathway Analysis
To investigate the roles and associated pathways of differentially expressed lncRNAs and mRNAs, Gene Ontology (GO) pathway analysis was performed to annotate the transcripts with terms under the biological process, cellular component, and molecular function ontologies. GO annotations of microarray genes were downloaded from NCBI1, UniProt2, and the GO3. A Fisher exact test was performed in order to locate the significant enrichment pathway. The resulting P-values were adjusted using the Benjamini–Hochberg false discovery rate (BH FDR) algorithm. Pathway categories with a FDR < 0.05 were reported.
Quantitative Reverse Transcription Real-Time Polymerase Chain Reaction
Total RNA (2 μg) was reverse transcribed with random hexamers using the Reverse Transcription System Kit (Promega Corporation, Madison, WI). LncRNA expression in placenta tissues was quantified using a standard quantitative reverse transcription real-time PCR (qRT-PCR) on the StepOne Plus system (Applied Biosystems, Foster City, CA, United States), in a total reaction volume of 20 μL, including 10 μL SYBR Premix Ex Taq (2x), 1 μL of PCR Forward Primer (10 uM), 1 μL of PCR Reverse Primer (10 uM), 2 μL of cDNA, and 6 μL of double-distilled water. Primer sets are listed in Supplementary Table 2. The qRT-PCR was performed with an initial denaturation step of 10 min at 95°C; 95°C (15 s), 60°C (30 s), 72°C (30 s) for a total 40 cycles; and a final extension step at 72°C for 5 min. All experiments were performed in triplicate. All samples were normalized to GAPDH. The geometric mean in each triplicate was used to calculate the relative lncRNAs concentrations (ΔCt = Ct of lncRNAs − Ct of GAPDH). The fold change of expression was calculated using the 2–ΔΔCt method.
Cell Line and Hypoxia Treatment
The immortalized EVT cell line HTR-8/SVneo was used. The cells were maintained in RPMI-1640 supplemented with 5% fetal calf serum, 1% L-glutamine 200 mM, and 1% penicillin–streptomycin (all from Invitrogen) under an atmosphere of 5% CO2 at 37°C. To establish the hypoxia model, the cells were exposed to 94%N2/1% O2/5% CO2 for 2, 24, 48, and 72 h. The control group was cultured under normal air condition for the corresponding period.
HIF-1α Overexpression in HTR-8/SVneo Cells
The coding sequence of human HIF-1α was amplified by PCR from human genomic DNA libraries (GENECHEM Company, Shanghai, China) and was subcloned into the pEZ-MO2 plasmid via homologous recombination according to manufacturer’s protocol (Hieff Clone® Plus One Step Cloning Kit, YEASEN, Shanghai, China) to construct pEZ-MO2–HIF-1α overexpression. The sequence of the construct was confirmed by DNA sequencing at Sangon Biotech Company (Shanghai, China). The HTR-8/SVneo cells at 70% confluence were transiently transfected with pEZ-MO2-HIF-1α plasmid for 48 h using Lipofectamine 3000 (Invitrogen) according to manufacturer’s instructions. The empty pEZ-MO2 plasmid was also transfected as negative control (NC).
Small Interfering RNA Transfection
lncRNA-HEIPP in HTR-8/SVneo cells were knocked down by transfection of preannealed stealth small interfering RNA (siRNA) duplexes designed online by Thermo Fisher’s BLOCK-iTTM RNAi Designer4, using Lipofectamine RNAiMAX transfection reagent (Invitrogen). Silencer® Negative Control #1 siRNA (si-NC, AM4611, Ambion) was transfected as an NC. The efficiency of knockdown was determined by RT-PCR.
Cell Counting Kit-8 Cell Proliferation Assay
Cell viability was assessed using the Cell Counting Kit-8 (Dojindo Laboratories, Kumamoto, Japan) according to the manufacturer’s protocol. Cell proliferation curves were plotted using the absorbance at each time point (0, 12, 24, and 48 h). All experiments were performed in triplicate.
Matrigel Invasion Assay
The invasion of HTR-8/SVneo cells across Matrigel was objectively evaluated in an invasion chamber according to the manufacturer’s protocol (Millipore). The cells transfected with si-NC or si-HEIPP (10 and 30 μL) were plated in the upper chambers of 24-well plates at a density of 1.0 × 104 cells/well in RPMI-1640 medium containing 5% fetal bovine serum; the membrane pore of the Transwell chamber was 8 μm in diameter. The membranes were coated with Matrigel (BD, United States) to form matrix barriers at the concentration of 200 μg/mL. A 700-μL aliquot of RPMI-1640 medium containing 10% FBS was immediately placed in the lower well of the chamber as a chemoattractant. After incubation for 48 h at 37°C, cells were completely removed from the upper surface of the filter using a sterile cotton swab, and cells that had migrated to the lower surface were fixed and stained with 0.1% crystal violet. Finally, the cell migration ability was determined by counting the number of stained cells on the membranes in 10 randomly selected, non-overlapping fields at 20× magnification under microscope. The migration index was calculated as the ratio of the percentage of cell migration in the various treatments to that of the vehicle. Each experiment was performed in triplicate and repeated three times.
Bisulfite Pyrosequencing
Bisulfite pyrosequencing was used to quantity DNA methylation at two sites in the DMR of the lncRNA-HEIPP gene. Two micrograms of DNA was subjected to bisulfate conversion using the EpiTect Bisulfite Kit (59104) (Qiagen Ltd., GmbH, Hilden, Germany). Pyrosequencing primers were designed using PyroMark Assay Design 2.0 (Qiagen Ltd., GmbH, Hilden, Germany) and listed in Supplementary Table 2. One hundred nanograms of bisulfite-converted DNA was PCR amplified in a 50-μL reaction volume containing 5 × PCR buffer (KAPA Ltd., United States), 10 mM of each dNTP (KAPA Ltd., United States), 50 mM forward and reverse primers (BGI Inc., Shanghai, China), and 1 U DNA polymerase (HotStart Taq, KAPA Ltd). The PCR cycling parameters were 95°C for 3 min; 40 cycles of 95°C for 30 s, a variable annealing temperature (Ta) for 30 s, and 72°C for 1 min; and extension at 72°C for 7 min. A 25-μL aliquot of each PCR product was subjected to pyrosequencing using the PyroMark Q96 ID Pyrosequencer (Qiagen Ltd.) following the manufacturer’s recommended protocols. The degree of methylation at each CpG site was determined using Pyro Q-CpG Quantification software. All samples were analyzed in triplicate.
ChIP Assay
The HTR-8 trophoblast cell line was used for ChIP-based analysis for enrichment of H3K4me3 on lncRNA-HEIPP promoter. ChIP assays were conducted according to the manufacturer’s instruction (Millipore, 17-10086). The chromatin extracted from cells was sheared, subjected to immunoprecipitation with H3K4me3 (Abcam) or immunoglobulin G (IgG) (Millipore) antibodies, reverse cross linked, and subjected to quantitative ChIP- PCR (qChIP). The qChIP was performed on sheared DNA with following primers: SP1: F: 5′ ATTTTGCTTCC TATCCCT 3′, R:5′ ACCGACAATCTCCTCAGT 3′, SP2: F: 5′ GA CTGAGGAGATTGTCGGT 3, R: 5′ CTTATGATGGTCTGGGT GA 3′, SP3: F: 5′ TGGAGGAGGGAGATGACA 3′, R: 5′ TGGA GGAGGGAGATGACA 3′, SP4: F: 5′ CATGCCCTCTCAGC CTAAC 3′, R: 5′ AGTCCCCCCAAGAAAAACA 3′. The four pairs of primers amplified −324 to −1,288 bp region upstream from the transcriptional starting site of lncRNA-HEIPP. The results were normalized to input and expressed as a percentage of the fold difference. IgG was used as an NC. Data were obtained from at least three independent experiments.
Statistical Analysis
All data are expressed as mean ± SD (standard deviation), and all statistical analyses were performed using the SPSS statistical software package (SPSS, Inc., Chicago, IL, United States). All the data were tested for homogeneity of variance by Bartlett test before analyzing the significance. Comparisons were made using Student t-test and analysis of variance. P < 0.05 was considered statistically significant.
Results
Differential Expression of lncRNAs in Placenta Between PE and Normal Patients
The demographic data of PE women and normal pregnant women are summarized in Table 1. There was no significant difference in maternal age between normal pregnancies and pregnancies with PE (P > 0.05). The systolic and diastolic blood pressure was significantly higher in pregnant women with PE (P < 0.001). The average proteinuria is 2.62 ± 0.20 g/24 h, which means that the cases enrolled this study were severe PE patients (proteinuria >2.0 g/24 h). Women from the PE group delivered earlier than normal pregnant women mostly due to the medical termination of the pregnancy, which is the major cause of the lower infant birth weight in the PE group, as there was no intrauterine growth restriction patient in the PE group enrolled in this study.
As the transcriptome of the placenta is dynamically changed during advancing gestation (), to ensure the identical or similar baselines between the two disease groups, we selected four pairs of PE patients and controls with matched gestational age within 2 weeks (Supplementary Table 1). We aligned lncRNA array data into RefSeq_NR, UCSC, and Ensembl database and compared the lncRNA expression levels between four human PE placenta (P) and four normal samples (N). The expression profiles of lncRNAs were shown by calculating the log-fold change (P/N). Among these lncRNAs, 405 were consistently upregulated, and 344 lncRNAs were consistently downregulated (Figure 1A). The clustering analysis also showed the differential expression patterns of mRNAs between PE placentas and normal samples (Figure 1B). The complete dataset regarding the significantly upregulated and downregulated lncRNA or mRNA was deposit on the public repository5.
FIGURE 1
GO Pathway Analysis
GO pathway analysis was performed to determine the lncRNA and mRNA enrichment in biological processes, cellular components, and molecular functions. HIF-1α pathway is known to play an important role in the pathogenesis of PE and identified as one of the most significant enrichment signaling pathways under our experimental settings (Figure 2A and Supplementary Table 3). We found that a number of lncRNAs and molecules are closely related with HIF-1α pathway (Figure 2B), and the top 10 lncRNAs (AF064860.7, AC027612.3, RP11-244N9.4, RP11-574M7.1, RP11-80I15.4, RP11-939C17.2, XLOC_013276, AC023085.1, AC009236.1, and RP11-184D12.1) showing mostly dramatic different expression between PE and normal placentas were selected for further validation (Supplementary Table 4).
FIGURE 2
Validation of Differentially Expressed lncRNAs
To further validate the preliminary data from microarray, we performed qRT-PCR to determine the expression of the ten candidate lncRNAs related with HIF-1α pathway in 30 PE placenta tissues and 30 normal placenta tissues. We found that 6 of the 10 lncRNAs manifested significant upregulation in PE samples compared with normal samples (Figure 3A), which are consistent with the array results. Among the six lncRNAs, lncRNA-AF064860.7 showed the most significant differences between PE samples and normal samples, so that we named it lncRNA-HEIPP (high expression in PE placenta).
FIGURE 3
In the HTR-8/SVneo human trophoblast cell lines, we found that hypoxia could induce the expression of lncRNA-HEIPP, which reached the highest level at 24 h upon hypoxia treatment (Figure 3B). Although the other five candidate lncRNAs also showed time-dependent changes of expression upon hypoxia treatment, the increasing trend and fold change were less significant than that in lncRNA-HEIPP. Moreover, the overexpression of HIF-1α increased the expression of lncRNA-HEIPP, further suggesting the interaction between HIF-1α pathway and lncRNA-HEIPP (Figure 3C).
Hypoxia Did Not Change the Methylation Profile of lncRNA-HEIPP Promoter in Placental Trophoblast Cells
As the epigenetic regulatory elements such as DNA methylation and H3K4me3 occupancy were critical for the gene expression regulation in placental development and trophoblast invasion (; ), we analyzed the epigenetic regulatory elements of lncRNA-HEIPP in the placenta in PE versus controls. Specifically, compared to the normal oxygen condition, the DNA methylation levels of lncRNA-HEIPP promoter in HTR-8/SVneo human trophoblast cells increased along with the extension of hypoxia time, which reached the highest level at 24 h upon hypoxia treatment (Figure 4A). However, the increase of methylation level was not significant at each time point (Figure 4B).
FIGURE 4
Binding of H3K4me3 in lncRNA-HEIPP Promoter Was Increased by Hypoxia Treatment
Enrichment of histone trimethylation at lysine 4 (H3K4me3) at lncRNA-HEIPP DMRs/promoters in HTR-8/SVneo human trophoblast cells that underwent hypoxia or normoxia was assessed using ChIP-qPCR. It was shown that the binding of H3K4me3 on lncRNA-HEIPP promoter was specifically enriched compared with the binding of IgG (Figure 5A). As shown in Figure 5B, the protein expression of H3K4me3 was not affected by hypoxia treatment. However, H3K4me3 manifested significantly higher promoter occupancy on lncRNA-HEIPP promoter in HTR-8/SVneo cells that underwent hypoxia for 24 h compared to those that underwent normoxia (Figure 5C), suggesting the histone modification, instead of the DNA methylation might be the epigenetic pathway that causes the increased expression of lncRNA-HEIPP in the placenta during PE undergone hypoxia.
FIGURE 5
Downregulation of lncRNA-HEIPP Promoted Trophoblast Invasion
Next, we assessed the biological functions of lncRNA-HEIPP in trophoblast cells. We knocked down the expression of lncRNA-HEIPP in HTR-8/SVneo human trophoblast cell line using HEIPP siRNA (si-HEIPP) and found that the proliferation of HTR8/SVneo cells was not affected (data not shown). However, the invasion of HTR8/SVneo cells was significantly promoted, compared to the cells treated with NC siRNA (Figure 6).
FIGURE 6
Discussion
Preeclampsia is a unique disease during pregnancy. It is a placental syndrome characterized by high blood pressure, edema, and proteinuria after 20 weeks of pregnancy (; ). Preeclampsia is the leading cause of perinatal and maternal death. At present, the PE incidence rate is 5–8% throughout the world and is 9.4% in China (). Current studies strongly suggested that placenta play central roles in the pathogenesis of PE (). The proliferation and invasion of placental trophoblast cells are the key to successful pregnancy and embryonic development (). The earlier study found that abnormal trophoblast differentiation, migration, and apoptosis often caused the decrease in the invasion of trophoblast cells and the failure of uterine spiral artery recasting, resulting in reduced uteroplacental perfusion and subsequently leading to PE (). However, the factors that affect placental trophoblast cells and the pathogenesis of PE have not been fully elucidated.
In recent years, epigenetic regulation has been revealed to play important roles in the embryonic development and the pathogenesis in many diseases. LncRNAs are a class of transcripts whose lengths exceed 200 nt. Initially, lncRNAs were thought to be transcribed noise. However, an increasing number of studies have reported that these lncRNAs have a series of important functions and participate in the development of many diseases, such as cancer (), cardiovascular illnesses (; Wang et al., 2016), and neurological illnesses (). Recently, some lncRNAs have been shown to be associated with PE (; ; Zhang et al., 2015). However, their related pathways and regulation were not fully revealed yet. In our current study, we performed lncRNA microarray and found that 1,051 lncRNAs were upregulated, and 344 lncRNAs were downregulated in PE placenta compared with normal placenta, which is comparable to the microarray results from previous study that 738 differently expressed lncRNAs in the placentas of PE in comparison with normal pregnancies (), but is much less than another study in which 15,646 upregulated and 13,178 downregulated lncRNAs were identified in the placenta tissues of early-onset PE patients (). However, in the latter study, the tissues were collected from early-onset PE patients, and the controls were preterm birth, and the biased factor, i.e., preterm labor, may partially explain the huge differences in the numbers of differentially expressed lncRNAs between studies.
The decreased trophoblast proliferation, migration, invasion, and stimulated apoptosis constitute the pivotal reasons leading to PE (). Previous studies revealed different lncRNAs that affect cellular functions of trophoblast cells (; ; ; ). However, the hypoxia and HIF-1α–induced pathways play central roles in the pathogenesis of PE (; ; ), whereas no HIF-1α–related lncRNAs have yet been revealed in the previous studies. In the current study, we first showed that lncRNAs related with HIF-1α pathway were significantly increased in PE placenta compared to normal placenta, among which lncRNA-HEIPP manifested most significant and time-dependent upregulation upon hypoxia treatment. In tumors, it was found that microRNAs and lncRNAs are important in enabling the key hypoxia-regulated processes (; ). In placenta, some microRNAs, such as microRNA-218 and microRNA-18b, were shown to be induced by hypoxia and HIF-1a in PE patients (; Wang et al., 2017). To our knowledge, lncRNA-HEIPP is the first lncRNA that was found to be upregulated by HIF-1a in placental tissues of PE patients. Moreover, Wang et al. (2018) also found that some lncRNAs, such as NONHSAT116812 and NONHSAT145880, were also confirmed in plasma specimens, making them to be potentially utilized as indicators for PE. Whether lncRNA-HEIPP can also be detected in plasma and used as novel indicators for PE warranted further confirmation.
In recent years, epigenetic regulation of placental gene expression has been considered to play an important role in human trophoblast differentiation () and contribute to the onset of PE (; ). The studies showed higher expression levels of death domain-associated protein 6 (DAXX), MMP9, and BCL2 in placenta from PE tissue, and this increased expression was well correlated to promoter demethylation (Wang et al., 2010; ; ). However, a significant number of the observed methylation changes were not associated with corresponding changes in gene expression, and vice versa (; ), indicating that alternate methods of epigenetic regulation will need to be explored in PE. Hypoxia was shown to alter the epigenetic profile in cultured human placental trophoblasts (Yuen et al., 2013). In this study, the DNA methylation levels of lncRNA-HEIPP promoter manifested increase upon hypoxia treatment, but the increases showed no significance at each time point after hypoxia. However, it was shown that the H3K4me3 occupancy of lncRNA-HEIPP promoter was significantly enriched in the PE placenta compared to the normal placenta. H3K4me3 is a histone modification marker that usually activates the gene transcription upon its binding to the promoter regions (; Zhang et al., 2018; ). These results suggested that the histone modification, instead of the DNA methylation, might be the epigenetic pathway that causes the increased expression of lncRNA-HEIPP in the placenta during PE.
Furthermore, we demonstrated that lncRNA-HEIPP knockdown promoted the invasion of trophoblasts. LncRNAs seem to contribute to PE via different pathways. LncRNA PRNCR1 promoted the progression of eclampsia by regulating the MAPK signal pathway (), while LncRNA CCAT1 promotes the progression of PE by regulating CDK4 (). As hypoxia usually induced the activation of JAK/STAT3 signaling pathway to promote trophoblast cell viability and angiogenesis in PE (Xu et al., 2017), whether the cellular functions of lncRNA-HEIPP were also mediated by JAK/STAT3 signaling pathway needs to be further investigated.
Conclusion
Our findings suggest that HIF-1a–related lncRNA-HEIPP is upregulated in the placenta from women who underwent PE mediated by the enriched binding of H3K4me3 in lncRNA-HEIPP promoter. LncRNA-HEIPP hampers the invasion of trophoblast cells and contributes to the pathogenesis of PE. Our results demonstrated that HIF-1a pathway–related lncRNA-HEIPP expression profile is a potentially new biomarker, and the manipulation of lncRNA-HEIPP may provide theoretical basis for prevention and treatment of PE.
Statements
Data availability statement
The microarray data has been uploaded to the GEO dataset and the accession number is GSE160888.
Ethics statement
The studies involving human participants were reviewed and approved by Human Research Ethics Committee of the Shanghai Changzheng Hospital at the Second Military Medical University. The patients/participants provided their written informed consent to participate in this study.
Author contributions
NS, HC, and YM: validation, formal analysis, and investigation. YM: writing – original draft. WP: methodology and investigation. XW and QZ: validation. LG: conceptualization, writing – review and editing, project administration, and funding acquisition. WL: conceptualization, supervision, and funding acquisition. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the National Key Research and Development Project Nos. 2018YFC1002802 (WL) and 2017YFC1001404 (LG and WL), National Natural Science Foundation of China Nos. 81901482 (NS), 81622020 and 81771608 (LG), and 81873821 (WL), Program of Shanghai Academic/Technology Research Leader No. 18XD1405100 (LG), “Dawn” Program of Shanghai Education Commission No. 17SG36 (LG), Program of Shanghai Subject Chief Scientist No. 2017BR017 (LG), and the Program for Professor of Special Appointment (Eastern Scholar) at Shanghai Institutions of Higher Learning (LG).
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2020.559478/full#supplementary-material
Footnotes
1.^http://www.ncbi.nlm.nih.gov/
3.^http://www.geneontology.org/
4.^http://rnaidesigner.thermofisher.com/rnaiexpress/rnaiDesign.jsp
References
1
Alici DavutogluE.Akkaya FiratA.OzelA.YilmazN.UzunI.Temel YukselI.et al (2018). Evaluation of maternal serum hypoxia inducible factor-1alpha, progranulin and syndecan-1 levels in pregnancies with early- and late-onset preeclampsia.J. Maternal-fetal Neonatal Med.311976–1982. 10.1080/14767058.2017.1333098
2
AouacheR.BiquardL.VaimanD.MirallesF. (2018). Oxidative stress in preeclampsia and placental diseases.Int. J. Mol. Sci.19:1496. 10.3390/ijms19051496
3
ApicellaC.RuanoC. S. M.MehatsC.MirallesF.VaimanD. (2019). The role of epigenetics in placental development and the etiology of preeclampsia.Int. J. Mol. Sci.20:2837. 10.3390/ijms20112837
4
ArmalyZ.JadaonJ. E.JabbourA.AbassiZ. A. (2018). Preeclampsia: novel mechanisms and potential therapeutic approaches.Front. Physiol.9:973. 10.3389/fphys.2018.00973
5
BoriboonhirunsarnD.PradyachaipimolA.ViriyapakB. (2017). Incidence of superimposed preeclampsia among pregnant Asian women with chronic hypertension.Hypertens. Pregnancy36226–231. 10.1080/10641955.2017.1311340
6
BrosensI.PuttemansP.BenagianoG. (2019). Placental bed research: i. the placental bed: from spiral arteries remodeling to the great obstetrical syndromes.Am. J. Obstetr. Gynecol.221437–456. 10.1016/j.ajog.2019.05.044
7
CaoC.LiJ.LiJ.LiuL.ChengX.JiaR. (2017). Long non-coding RNA Uc.187 is upregulated in preeclampsia and modulates proliferation, apoptosis, and invasion of HTR-8/SVneo trophoblast cells.J. Cell. Biochem.1181462–1470. 10.1002/jcb.25805
8
ChelbiS. T.VaimanD. (2008). Genetic and epigenetic factors contribute to the onset of preeclampsia.Mol. Cell. Endocrinol.282120–129. 10.1016/j.mce.2007.11.022
9
ChenH.MengT.LiuX.SunM.TongC.LiuJ.et al (2015). Long non-coding RNA MALAT-1 is downregulated in preeclampsia and regulates proliferation, apoptosis, migration and invasion of JEG-3 trophoblast cells.Int. J. Clin. Exp. Pathol.812718–12727.
10
ChenJ.KhalilR. A. (2017). Matrix metalloproteinases in normal pregnancy and preeclampsia.Prog. Mol. Biol. Transl. Sci.14887–165. 10.1016/bs.pmbts.2017.04.001
11
ChoudhryH.HarrisA. L.McIntyreA. (2016). The tumour hypoxia induced non-coding transcriptome.Mol. Aspects Med.435–53. 10.1016/j.mam.2016.01.003
12
CotechiniT.KomisarenkoM.SperouA.Macdonald-GoodfellowS.AdamsM. A.GrahamC. H. (2014). Inflammation in rat pregnancy inhibits spiral artery remodeling leading to fetal growth restriction and features of preeclampsia.J. Exp. Med.211165–179. 10.1084/jem.20130295
13
DeyssenrothM. A.PengS.HaoK.LambertiniL.MarsitC. J.ChenJ. (2017). Whole-transcriptome analysis delineates the human placenta gene network and its associations with fetal growth.BMC Genom.18:520. 10.1186/s12864-017-3878-0
14
FangM.DuH.HanB.XiaG.ShiX.ZhangF.et al (2017). Hypoxia-inducible microRNA-218 inhibits trophoblast invasion by targeting LASP1: implications for preeclampsia development.Int. J. Biochem. Cell Biol.8795–103. 10.1016/j.biocel.2017.04.005
15
FanucchiS.FokE. T.DallaE.ShibayamaY.BornerK.ChangE. Y.et al (2019). Immune genes are primed for robust transcription by proximal long noncoding RNAs located in nuclear compartments.Nat. Genet.51138–150. 10.1038/s41588-018-0298-2
16
FerdinJ.NishidaN.WuX.NicolosoM. S.ShahM. Y.DevlinC.et al (2013). HINCUTs in cancer: hypoxia-induced noncoding ultraconserved transcripts.Cell Death. Differ.201675–1687. 10.1038/cdd.2013.119
17
FuG.YeG.NadeemL.JiL.ManchandaT.WangY.et al (2013). MicroRNA-376c impairs transforming growth factor-beta and nodal signaling to promote trophoblast cell proliferation and invasion.Hypertension61864–872. 10.1161/hypertensionaha.111.203489
18
GamageT.SchierdingW.HurleyD.TsaiP.LudgateJ. L.BhoothpurC.et al (2018). The role of DNA methylation in human trophoblast differentiation.Epigenetics131154–1173. 10.1080/15592294.2018.1549462
19
GormleyM.OnaK.KapidzicM.Garrido-GomezT.ZdravkovicT.FisherS. J. (2017). Preeclampsia: novel insights from global RNA profiling of trophoblast subpopulations.Am. J. Obstet. Gynecol.217200.e1–200.e17.
20
GutierrezJ. A.GomezI.ChiarelloD. I.SalsosoR.KleinA. D.Guzman-GutierrezE.et al (2019). Role of proteases in dysfunctional placental vascular remodelling in preeclampsia.Biochimica et Biophysica Acta. Mol. Basis Dis.1866:165448. 10.1016/j.bbadis.2019.04.004
21
Harati-SadeghM.KohanL.TeimooriB.MehrabaniM.SalimiS. (2018). The association of the placental Hypoxia-inducible factor1-alpha polymorphisms and HIF1-alpha mRNA expression with preeclampsia.Placenta6731–37. 10.1016/j.placenta.2018.05.005
22
HeX.HeY.XiB.ZhengJ.ZengX.CaiQ.et al (2013). LncRNAs expression in preeclampsia placenta reveals the potential role of LncRNAs contributing to preeclampsia pathogenesis.PLoS One8:e81437. 10.1371/journal.pone.0081437
23
HosseiniA.DolatiS.HashemiV.Abdollahpour-AlitappehM.YousefiM. (2018). Regulatory T and T helper 17 cells: their roles in preeclampsia.J. Cell. Physiol.2336561–6573. 10.1002/jcp.26604
24
JiangQ.ShanK.Qun-WangX.ZhouR. M.YangH.LiuC.et al (2016). Long non-coding RNA-MIAT promotes neurovascular remodeling in the eye and brain.Oncotarget749688–49698. 10.18632/oncotarget.10434
25
JiaoS.WangS. Y.HuangY. (2018). LncRNA PRNCR1 promoted the progression of eclampsia by regulating the MAPK signal pathway.Eur. Rev. Med. Pharmacol. Sci.223635–3642.
26
KadyrovM.KingdomJ. C.HuppertzB. (2006). Divergent trophoblast invasion and apoptosis in placental bed spiral arteries from pregnancies complicated by maternal anemia and early-onset preeclampsia/intrauterine growth restriction.Am. J. Obstet. Gynecol.194557–563. 10.1016/j.ajog.2005.07.035
27
KawasakiK.KondohE.ChigusaY.KawamuraY.MogamiH.TakedaS.et al (2019). Metabolomic profiles of placenta in preeclampsia.Hypertension73671–679. 10.1161/hypertensionaha.118.12389
28
KnoflerM.PollheimerJ. (2013). Human placental trophoblast invasion and differentiation: a particular focus on Wnt signaling.Front. Genet.4:190. 10.3389/fgene.2013.00190
29
KorkesH. A.De OliveiraL.SassN.SalahuddinS.KarumanchiS. A.RajakumarA. (2017). Relationship between hypoxia and downstream pathogenic pathways in preeclampsia.Hypertens. Pregnancy36145–150. 10.1080/10641955.2016.1259627
30
KuzmanovA.WielockxB.RezaeiM.KettelhakeA.BreierG. (2012). Overexpression of factor inhibiting HIF-1 enhances vessel maturation and tumor growth via platelet-derived growth factor-C.Int. J. Cancer131E603–E613.
31
KwakY. T.MuralimanoharanS.GogateA. A.MendelsonC. R. (2019). Human trophoblast differentiation is associated with profound gene regulatory and epigenetic changes.Endocrinology1602189–2203. 10.1210/en.2019-00144
32
LeaveyK.WilsonS. L.BainbridgeS. A.RobinsonW. P.CoxB. J. (2018). Epigenetic regulation of placental gene expression in transcriptional subtypes of preeclampsia.Clin. Epigenet.10:28.
33
LeungA.TracC.JinW.LantingL.AkbanyA.SaetromP.et al (2013). Novel long noncoding RNAs are regulated by angiotensin II in vascular smooth muscle cells.Circ. Res.113266–278. 10.1161/circresaha.112.300849
34
LiJ. L.LiR.GaoY.GuoW. C.ShiP. X.LiM. (2018). LncRNA CCAT1 promotes the progression of preeclampsia by regulating CDK4.Eur. Rev. Med. Pharmacol. Sci.221216–1223.
35
LiuD.YuX.WangS.DaiE.JiangL.WangJ.et al (2016). The gain and loss of long noncoding RNA associated-competing endogenous RNAs in prostate cancer.Oncotarget757228–57238. 10.18632/oncotarget.11128
36
LongW.RuiC.SongX.DaiX.XueX.LuY.et al (2016). Distinct expression profiles of lncRNAs between early-onset preeclampsia and preterm controls.Clin. Chimica Acta; Int. J. Clin. Chem.463193–199. 10.1016/j.cca.2016.10.036
37
MayrinkJ.CostaM. L.CecattiJ. G. (2018). Preeclampsia in 2018: revisiting concepts, physiopathology, and prediction.Sci. World J.2018:6268276.
38
Mohammadpour-GharehbaghA.JahantighD.EskandariM.SadeghM. H.NematollahiM. H.RezaeiM.et al (2019). Genetic and epigenetic analysis of the BAX and BCL2 in the placenta of pregnant women complicated by preeclampsia.Apoptosis : Int. J. Program. Cell Death24301–311. 10.1007/s10495-018-1501-8
39
NovakovicB.Evain-BrionD.MurthiP.FournierT.SafferyR. (2017). Variable DAXX gene methylation is a common feature of placental trophoblast differentiation, preeclampsia, and response to hypoxia.FASEB J.312380–2392. 10.1096/fj.201601189rr
40
SatiS.GhoshS.JainV.ScariaV.SenguptaS. (2012). Genome-wide analysis reveals distinct patterns of epigenetic features in long non-coding RNA loci.Nucleic Acids Res.4010018–10031. 10.1093/nar/gks776
41
SongX.RuiC.MengL.ZhangR.ShenR.DingH.et al (2017). Long non-coding RNA RPAIN regulates the invasion and apoptosis of trophoblast cell lines via complement protein C1q.Oncotarget87637–7646. 10.18632/oncotarget.13826
42
TongW.GiussaniD. A. (2019). Preeclampsia link to gestational hypoxia.J. Dev. Origins Health Dis.10322–333. 10.1017/s204017441900014x
43
WangS.WangX.WengZ.ZhangS.NingH.LiB. (2017). Expression and role of microRNA 18b and hypoxia inducible factor-1alpha in placental tissues of preeclampsia patients.Exp. Therapeutic Med.144554–4560.
44
WangX.ChenY.DuL.LiX.LiX.ChenD. (2018). Evaluation of circulating placenta-related long noncoding RNAs as potential biomarkers for preeclampsia.Exp. Therapeutic Med.154309–4317.
45
WangY. N.ShanK.YaoM. D.YaoJ.WangJ. J.LiX.et al (2016). Long Noncoding RNA-GAS5: a novel regulator of hypertension-induced vascular remodeling.Hypertension68736–748. 10.1161/hypertensionaha.116.07259
46
WangZ.LuS.LiuC.ZhaoB.PeiK.TianL.et al (2010). Expressional and epigenetic alterations of placental matrix metalloproteinase 9 in preeclampsia.Gynecol. Endocrinol.2696–102. 10.3109/09513590903184100
47
XuC.LiX.GuoP.WangJ. (2017). Hypoxia-Induced activation of JAK/STAT3 signaling pathway promotes trophoblast cell viability and angiogenesis in preeclampsia.Med. Sci. Monitor : Int. Med. J. Exp. Clin. Res.234909–4917. 10.12659/msm.905418
48
XuY.LianY.ZhangY.HuangS.ZuoQ.YangN.et al (2018). The long non-coding RNA PVT1 represses ANGPTL4 transcription through binding with EZH2 in trophoblast cell.J. Cell Mol. Med.221272–1282.
49
YongH. E. J.MurthiP.BrenneckeS. P.MosesE. K. (2018). Genetic approaches in preeclampsia.Methods Mol. Biol.171053–72. 10.1007/978-1-4939-7498-6_5
50
YuenR. K.ChenB.BlairJ. D.RobinsonW. P.NelsonD. M. (2013). Hypoxia alters the epigenetic profile in cultured human placental trophoblasts.Epigenetics8192–202. 10.4161/epi.23400
51
ZhangE.HeX.ZhangC.SuJ.LuX.SiX.et al (2018). A novel long noncoding RNA HOXC-AS3 mediates tumorigenesis of gastric cancer by binding to YBX1.Genome Biol.19:154.
52
ZhangY.ZhaoH. J.XiaX. R.DiaoF. Y.MaX.WangJ.et al (2019). Hypoxia-induced and HIF1alpha-VEGF-mediated tight junction dysfunction in choriocarcinoma cells: implications for preeclampsia.Clinica Chimica Acta; Int. J. Clin. Chem.489203–211. 10.1016/j.cca.2017.12.010
53
ZhangY.ZouY.WangW.ZuoQ.JiangZ.SunM.et al (2015). Down-regulated long non-coding RNA MEG3 and its effect on promoting apoptosis and suppressing migration of trophoblast cells.J. Cell. Biochem.116542–550. 10.1002/jcb.25004
Summary
Keywords
lncRNA, placenta, preeclampsia, hypoxia, H3K4me3
Citation
Sun N, Chen H, Ma Y, Pang W, Wang X, Zhang Q, Gao L and Li W (2020) H3K4me3-Mediated Upregulation of LncRNA-HEIPP in Preeclampsia Placenta Affects Invasion of Trophoblast Cells. Front. Genet. 11:559478. doi: 10.3389/fgene.2020.559478
Received
06 May 2020
Accepted
16 November 2020
Published
11 December 2020
Volume
11 - 2020
Edited by
Geng Chen, East China Normal University, China
Reviewed by
Guodong Fu, University of Toronto, Canada; Tian Fuju, Shanghai Jiao Tong University School of Medicine, China
Updates
Copyright
© 2020 Sun, Chen, Ma, Pang, Wang, Zhang, Gao and Li.
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.
*Correspondence: Lu Gao, lu.gao@smmu.edu.cnWen Li, liwen@smmu.edu.cn
†These authors have contributed equally to this work
This article was submitted to Computational Genomics, a section of the journal Frontiers in Genetics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.