Next Article in Journal
Ocular Manifestations in a Chinese Pedigree of Familial Amyloidotic Polyneuropathy Carrying the Transthyretin Mutation c.401A>G (p.Tyr134Cys)
Previous Article in Journal
Mobilome of Apicomplexa Parasites
Previous Article in Special Issue
Loss of Heterozygosity in the Tumor DNA of De Novo Diagnosed Patients Is Associated with Poor Outcome for B-ALL but Not for T-ALL
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 13
Issue 5
10.3390/genes13050888
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Preliminary Interpretations of Epigenetic Profiling of Cord Blood in Preeclampsia

by
Junrui Ma
1,2,†,
Zhongqun Zhan
3,†,
Ning Li
4,†,
Yanli Huang
5,
Yan Li
5,
Lu Liu
5,
Qi Shen
3,
Qiao Chu
1,
Xiaonan Wang
1,
Benqing Wu
3,6,* and
Hui Zhang
3,*
1
School of Public Health, Shanghai Jiao Tong University School of Medicine, Shanghai 200025, China
2
Faculty of Medical Laboratory Science, College of Health Science and Technology, Shanghai Jiao Tong University School of Medicine, Shanghai 200025, China
3
Institute of Translational Medicine, University of Chinese Academy of Sciences Shenzhen Hospital, Shenzhen 518106, China
4
Cytotherapy Laboratory, The First Affiliated Hospital (Shenzhen People’s Hospital) Southern University of Science and Technology, Shenzhen 518020, China
5
Department of Obstetrics, University of Chinese Academy of Sciences Shenzhen Hospital, Shenzhen 518106, China
6
Department of Neonatology, Shenzhen Guangming Maternity & Child Healthcare Hospital, Shenzhen 518107, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Genes 2022, 13(5), 888; https://0-doi-org.brum.beds.ac.uk/10.3390/genes13050888
Submission received: 12 April 2022 / Revised: 1 May 2022 / Accepted: 14 May 2022 / Published: 16 May 2022
(This article belongs to the Special Issue Epigenomics, Epitranscriptomics and Novel Therapy Approaches in Hematology)
Browse Figures
Versions Notes

Abstract

:
Preeclampsia (PE) is characterized by new-onset hypertension after 20 weeks of pregnancy and results in high maternal and fetal mortality worldwide. It has been reported that PE is associated with abnormalities in the umbilical cord and cord blood. However, previous studies were focused primarily on the transcriptomics level, while the underlying gene regulatory landscapes are still unclear. Thus, we performed the Assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC-seq) using the umbilical cord blood samples collected from a patient with superimposed PE and three healthy donors to uncover the chromatin accessibility changes attributed to PE. We have identified genes associated with immunomodulation and hypoxia response that have higher chromatin accessibility close to their transcription start sites. Motif analysis indicated that the GATA family transcription factor binding was enriched in PE and may play an essential regulatory role in the disease progression. Overall, our findings provide an overview of gene regulatory programs and the corresponding downstream pathways associated with PE that may influence the placenta function and fetal growth.

1. Introduction

Preeclampsia (PE), one of the leading causes of maternal and fetal mortality, is a pregnancy-associated complication characterized by new-onset hypertension and is often associated with proteinuria [1]. The International Society for the Study of Hypertension in Pregnancy (ISSHP) defines PE as de-novo hypertension (systolic blood pressure (SBP) ≥ 140 mmHg or diastolic blood pressure (DBP) ≥ 90 mmHg) present after 20 weeks of gestation accompanied by ≥1 of the following new-onset conditions: proteinuria (≥1+, 30 mg/dL; urine Protein/Creatinine Ratio (PCR) ≥ 30 mg/mmol (0.3 mg/mg)), material organ dysfunction, including acute kidney injury (creatinine ≥ 90 umol/L; 1 mg/dL), liver involvement, neurological complications, hematological complications, uteroplacental dysfunction [2,3]. Proteinuria is no longer required for a diagnosis of PE. The most life-threatening complication is HELLP (hemolysis, elevated liver enzymes, and low platelets) syndrome, leading to high mortality in patients with PE [4]. Depending on the gestational age when hypertension arises, it can be classified into either early-onset (EO, <34 weeks gestation) or late-onset (LO, ≥34 weeks gestation) subtypes [5]. Hypertension, as the main criteria in a diagnosis of PE, is strongly associated with the initiation and progression of PE. There are two types of hypertension. Gestational hypertension occurs after 20 weeks of pregnancy in the absence of other symptoms of PE. Chronic hypertension occurs before pregnancy of other causes rather than pregnancy itself and is associated with the development of superimposed PE [6].
Cord blood refers to the blood remains in the placenta and the umbilical cord after delivery. Thus, it has been considered a valuable source to study if PE has adverse effects on the placenta and if these can be delivered to the fetus. As previously reported, PE is associated with abnormalities in the cord blood in different aspects, including changes in Biochemistry, immune cell subtypes, Hematology, and Epigenetics [7,8,9,10]. In addition, fetuses exposed to PE havean increased risk of neonatal death, cardiovascular diseases [11], neurodevelopmental disorders [12], thrombocytopenia [13], and hypertension [14]. However, the underlying molecular mechanisms of PE and the corresponding regulatory programs are still largely unclear.
It has been reported that PE is associated with enhanced oxidative stress (OS) [15] and inflammatory response [16] contributed by distinct T helper cell subtypes [7]. Th2 and Treg cells arelargely diminished in PE compared to normal pregnancies. The shift in Th1/Th2 and Th17/Treg balance may further promote inflammation in the umbilical cord blood in PE [8]. RNA-sequencing/microarray was performed using placenta and plasma cell-free RNA to study pregnancies complicated by PE at the transcriptomic level [15,17]. Up-regulated genes in PE are associated with cell proliferation and differentiation, immunity and defense, cell structure, and lipid metabolism and transport, compared to healthy donors. Several genes were suggested as early markers for PE detection, including CLDN7, PAPPA2, SNORD14A, PLEKHH1, MAGEA10, TLE6, and FABP1, identified from plasma samples. Higher up in the hierarchy, gene expression during the development and progression of PE is regulated by epigenetic mechanisms, including DNA methylation, histone modification, and microRNAs (miRNAs) [18]. However, investigations on the epigenetic changes in PE are still limited. DNA methylation landscape isthe best-studied epigenetic mechanism in PE [19,20]. Genome-scale DNA methylation profile of cord blood samples from patients with early on-set PE revealed a trend of hypomethylation. Differential expression analysis identified a group of genes involved in inflammation, lipid metabolism, and proliferation, including RUNX3, LINE-1, ADORA2B, VHL, MPP-9, TBXAS1, ERVW-1, SPESP1, WNT2, AGT, DDAH1, CALCA, IL12B, FAS, PIK31, IGF1, GNAS, IGF2, and HSD11B2 [18,19,21].
Assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC-seq) is a method for mapping chromatin accessibility genome-wide [22]. Compared to DNA methylation profiling, ATAC-Seq gives us a more precise overview of open chromatin regions and the underlying gene regulatory networks by peak calling since chromatin accessibility is determined by a dynamic interplay among histones, TFs, and active chromatin remodelers [23]. Thus, we performed ATAC-Seq on cord blood samples from three healthy pregnant donors and a patient with superimposed PE with high sequencing coverage (~33 k reads/sample on average). We identified genes associated with immunomodulation and hypoxia response that have higher chromatin accessibility close to their transcription start sites. Motif analysis indicated that the GATA family transcription factor binding was enriched in PE and may play an important regulatory role. Overall, our findings provide an overview of gene regulatory programs and the corresponding downstream pathways associated with PE that may influence the placenta function and fetal growth. This study will be used as a preliminary investigation of a large-scale epigenetic profiling study of PE with different subtypes.

2. Materials and Methods

2.1. Patient Recruitment Criteria

This study has been approved by the Ethics Committee of the University of Chinese Academy of Sciences Shenzhen Hospital. The research adhered to the tenets of the Declaration of Helsinki as amended in 2013 [24]. All participants recruited have gone through the informed consent process and agreed to donate their cord blood samples. To minimize interpersonal variability, participants were carefully chosen as follows:
PE group—excluded if met any of the following criteria: abnormal pregnancy; gestational diabetes mellitus (GDM); severe organ dysfunction, such as liver, kidney, and heart; co-infection or other inflammatory diseases; history of abnormal pregnancy or delivery (including spontaneous abortion, stillbirth, malformation, etc.); a family history of psychosis or mental illness; a history of alcohol or drug abuse.
Control group: healthy pregnant women who participated in the prenatal examination on time and met the above exclusion criteria with age ≤ 35.

2.2. Sample Preparation for ATAC-seq

Cord blood (approximately 50–75 mL) was collected from each donor during delivery. Lymphoprep (Serumwerk Bernburg cat. 1858) was used to isolate mononuclear cells (MNCs) from the cord blood. Red blood cells were lysed using RBC lysing buffer (Biolegend cat. 420310). Cells were recovered from liquid nitrogen. Approximately 10 K cells were processed following the established ATAC-Seq protocol [25]. Samples were sequenced using the MGI2000 platform.

2.3. Peak Calling

Raw ATAC-seq reads were mapped to the hg38 reference genome (BioProject: PRJNA31257) using Bowtie (version: 1.3.1). Duplicates were removed with the MarkDuplicates function from Picard Tools (version: 2.26.11) with default parameters. Low-quality reads with mapping quality (mapQ) lower than 30 were removed using samtools. Reads mapped to ChrX, Y, and Random were filtered. MACS2 (version: 2.2.7.1) was used for peak calling with parameters -nomodel -shift -100 -extsize 200. The blacklisted genomic (obtained from https://github.com/Boyle-Lab/Blacklist/blob/master/lists/hg38-blacklist.v2.bed.gz, accessed on 3 March 2022) regions were removed from the peak files using bedtools. Peaks from all samples were merged using bedtools. HTseq-count (version: 1.99.2) was used for counting reads within each peak. Then R package ‘EdgeR’ was used to identify differentially expressed peaks with a false discovery rate of less than 0.1 and |log fold change (logFC)| > 2. To identify peaks with high confidence, we further strengthened the selection criteria by adding an additional cutoff of log count per million (logCPM) > 0.

2.4. Peak Annotation

The annotatePeaks.pl function of HOMER [26] (version: 4.11.1) was utilized to associate peaks with nearby genes and to annotate the genomic location of each peak.

2.5. Peak Visualization

Integrative Genomics Viewer (IGV) [27] (2.12.3) was utilized to visualize the peaks.

2.6. Pathway Enrichment Analysis

The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis was performed by KOBAS [28] (http://kobas.cbi.pku.edu.cn/, accessed on 28 March 2022). Enriched pathways were selected as FDR < 0.1(calculated by hypergeometric test/Fisher’s exact test).

3. Results

3.1. Case Presentation

The patient with PE recruited in this study was a 30-year-old pregnant woman with a history of chronic hypertension for 8 years who developed chronic hypertension associated with superimposed PE with proteinuria at 35 weeks of pregnancy. Before labor, the last blood test results showed abnormal coagulation (PT: 9.3 s, FIB: 3.88 g/L, PT-R: 121.1% and D-Dimer: 0.79 mg /L). However, proteinuria and most of the abnormal blood test results were improved after delivery. (Table 1). To understand if PE would induce any epigenetic changes that may impact the placenta function and fetal growth, we performed ATAC-seq, which allowed us to uncover gene regulatory programs using the cord blood sample from the patient with PE. Three cord blood samples collected from healthy pregnant women were used as negative controls. (Supplementary Table S1; The main test result of healthy gravidae (before labor)).

3.2. Immune, Myocardiopathy, and Hypoxia Response-Related Pathways Are Up-Regulated in Superimposed PE

ATAC-Seq identifies the gnome-wide set of cis-regulatory elements and transcription factor binding profiles from open chromatin regions using a hyperactive Tn5 transposase that cuts and ligates sequencing adaptor into the regions. Hence, these regions are enriched with sequencing fragments that are also referred to as peaks. MACS2 is a widely used peak calling tool that defines regions where transcription factors bind.
In this study, 34,304 peaks were identified from the four ATAC-Seq samples using MACS2. Although peak intensities of the four samples at TSS were similar in general (Supplementary Figure S1A,B), differential expression analysis, using edgeR with cutoffs of |logFC| > 2 and FDR < 0.1, identified 786 differential expressed peaks comparing between the patient with superimposed PE and the three healthy controls, among which 767 were up-regulated and 19 were down-regulated. (Supplementary Figure S1C). Peaks were then mapped to their closest genes based on the peak-to-TSS (Transcription Start Site, TSS) distances.
To further investigate the biological functions of genes with higher chromatin accessibility at their TSS in PE, we conducted pathway analysis using the KEGG database to annotate the up-regulated genes. Only unique genes were extracted for the pathway analysis. Results indicated that the enriched pathways fell into four broad groups: (1) immune (in purple), including Th1 and Th2 cell differentiation, T cell receptor signaling, Toll-like receptor signaling, leukocyte transendothelial migration, and inflammatory mediator regulation of TRP channels; (2) myocardiopathy (in blue), including vascular smooth muscle contraction, dilated cardiomyopathy and adrenergic signaling in cardiomyocytes); (3) hematopoiesis (in red), including Platelet activation; and (4) hypoxia response (in green), including HIF-1 signaling (Figure 1a). Full lists of pathways with adjusted p-values are available in Supplementary Table S2 (up-regulated pathway). Genes contributing to these four pathways were highlighted in the same color, respectively (Figure 1b–e).

3.3. Higher Chromatin Accessibility across TSS of Genes of Differentially Expressed Peaks in Superimposed PE

We further investigated the genomic locations of the 786 differentially expressed peaks as detected above. Since these peaks represent either cis-regulatory elements or transcription factor binding sites, they may locate in different regions relative to a gene:(1) Promoter (≤1 kb): less than 1 kilobase (kb) away from the TSS; (2) Promoter (1–2 kb): 1–2 kb away from the TSS; (3) Promoter (2–3 kb): 2–3 kb away from the TSS; (4) 5′ UTR: untranslated region (UTR) at the 5′ of a gene; (5) 3′ UTR: UTR at the 3′ of a gene; (6) 1st Exon: the first exon of a gene; (7) Other Exon: other exons apart from the first one; (8) 1st Intron: the first intron of a gene; (9) Other Intron: other introns apart from the first one; (10) Distal Intergenic: the distal intergenic regions. Thus, each peak was annotated with a unique genomic location. More than 50% of the peaks, either up-regulated or down-regulated, were at the non-promoter regions, indicating potential cis-regulatory element binding sites, which could either be enhancers or silencers that regulate genes nearby (Figure 2a,b).

3.4. GATA-Binding Motifs Are Enriched in Superimposed PE

HOMER contains a motif discovery algorithm that was designed for regulatory element analysis in genomics applications [26], which can be applied to nucleic acid sequences. Thus, it was used to identify transcription factors that may have a regulatory role in the progression of PE by finding the enriched motifs in the defined open chromatin regions that were only detected in the PE sample. The results revealed enrichment of five GATA family transcription factors—GATA 1,2,3,4,6 and a co-binding motif of both GATA and SCL. (Figure 3). We have extracted the peak regions in various types of blood cells from the available ChIP-Seq/ChIP-ChIP datasets of GATA factors, predominately GATA1 and GATA2. Associating these peaks within nearby genes, we confirmed that the GATA family binds to different genomic locations, among which the majority are not located inthe promoter regions [29,30,31,32,33,34] (Supplementary Figure S2), indicating that the GATAtranscription factors were able to bind to canonical GATA sequence motifs at cis-regulatory regions of many PE-associated genes during the disease progression.

3.5. Identification of Novel PE-Associated Genes

We further filtered the differentially expressed peaks with an additional cutoff of logCPM> 0. As expected, only a few peaks remained but were considered with higher confidence in terms of reproducibility. There were eight peaks up-regulated, among which five are close to TSS of five protein-coding genes—PTGIS, SERINC2, PRR25, TNFRSF6B, and PCMTD2, that have not been previously reported in PE (Figure 4a). PTGIS is a well-known marker gene for cardiovascular diseases and hypertension. PRR25 is involved in cardiomyopathy. SERINC2, TNFRSF6B, and PCMTD2 are associated with immune and pro-inflammatory responses (Supplementary Figure S3). In addition, we observed a similar expression pattern of well-known PE-associated genes, including PAPPA2 and GNAS, although they were not selected as high confidence (Figure 4b).

4. Discussion

PE is a common hypertensive disorder in 2–8% of pregnancies and may lead to severe outcomes, even maternal and fetal deaths. Fetuses after exposure to PE are at a higher risk of developing cardiovascular diseases [11] and various neurodevelopmental disorders, particularly autism spectrum disorder and attention-deficit/hyperactivity disorder (ADHD) [12], attributed to oxidative stress and inflammation. Thus, an in-depth understanding of the underlying molecular mechanisms, especially the transcriptional regulation, is urgently needed to prevent PE and provide further guidance in treating PE and PE-associated disorders. Cord blood collected from the umbilical cord after delivery provides an environment shared between the placenta and the fetus. Thus, it has been considered a valuable source to study if PE has adverse effects on the placenta and if these can be delivered to the fetus. Hence, ATAC-seq was administrated to reveal the epigenetic changes in cord blood from a patient with PE and to interrogate the underlying gene regulatory programs.
Previous studies have already claimed the association of hypoxia with PE. HIF-1α, as an oxygen sensor that regulates oxygen homeostasis in the human placenta, plays a key role in the pathogenesis of PE [35,36,37]. Furthermore, maternal PE can also result in fetal hypoxia [38]. An accumulation of pro-inflammatory cytokines, an increase in reactive oxygen species (ROS), and a reduced concentration of antioxidants and antioxidant enzymes in PE have been reported, indicating enhanced inflammation and oxidative stress [12]. Our ATAC-seq data showed higher chromatin accessibility of hypoxia response-related genes in the patient with PE. This points to hypoxia in PE in the cord blood, which may further pass to the fetus under the regulation of the HIF-1 pathway. In terms of cardiovascular function, the up-regulation of Dilated cardiomyopathy (DCM), vascular smooth muscle contraction, and Adrenergic signaling pathways may be owing to uncontrolled hypertension that leads to damaged and narrowed coronary arteries, heart attack, or even heart failure.
Intriguingly, our results showed an up-regulation of the platelet activation pathway in the cord blood from the epigenetic landscape, which is consistent with abnormal coagulation, as indicated in the patient’s regular blood test result. It has been reported that platelet activation is involved in the progression of PE and contributes to the thrombotic and coagulopathic complications of the disease [39,40]. Our results confirmed this statement and provided further evidence that the platelet activation itself can potentially be used as a biomarker for the prediction of PE.
Altered immunomodulation is the first pathogenic sign of PE, which refers to an abnormal immune response to the allogeneic fetus [41,42,43]. Pathway enrichment analysis discovered an enrichment of immune response pathways that involved various subtypes of T helper cells. T helper cells are CD4 positive lymphocytes that can differentiate into several subtypes, including Th1, Th2, Th17, etc. These subtypes modulate immune responses by secreting a different panel of cytokines that trigger or suppress inflammation. It has been found that the Th1/Th2 balance is altered in PE, which is characterized by a higher ratio of circulating Th1/Th2 lymphocytes [44]. Our ATAC-Seq results also provided a hint of changes in the T helper cell subtype ratios. We found that the T helper cell subtype differentiation-related pathway and the Th1 and Th2 cell differentiation pathways were up-regulated in PE, which is in line with previous findings [41,42,43]. In addition, immune response pathways, including Toll-like receptor signaling pathway, T cell receptor signaling pathway, and Inflammatory mediator regulation of TRP channels, were observed, suggesting a pro-inflammatory outcome contributed by the imbalanced subtypes of T helper cells. We also detected an up-regulation of leukocyte transendothelial migration in PE, indicating an inflammatory infiltration. All the above results reflect an inflammatory microenvironment of cord blood in PE.
Furthermore, motif analysis of the up-regulated peak regions revealed enrichment of GATA family transcription factors. GATA factors are associated with the development and differentiation of T helper cells and other types of immune cells, such as dendritic cells and macrophages [45]. The GATA family, especially GATA3, regulates Th1/Th2 cell differentiation [46,47]. In addition, the GATA family is also associated with the hypoxia response via HIF-1 pathways. GATA3 has the ability to interact with and stabilize HIF-1α protein under the hypoxia condition, resulting in prolonged hypoxia response [48]. Thus, the GATA family, including GATA1, GATA2, GATA3, and GATA4, are likely to be key regulatory elements that regulate the differentially expressed genes involved in the pathogenesis of PE.
In agreement with previous studies, our findings provide further evidence that hypoxia, oxidative stress, inflammation, and cardiovascular dysfunction are mutually linked, resulting in the pathophysiological progression of PE (Figure 5). Although the number of samples involved in this study is limited, especially for the PE condition, our ATAC-seq results provided an overview of gene regulatory programs and the corresponding downstream pathways associated with PE that may influence placenta function and fetal growth. Additionally, we discovered novel genes with high confidence that were not previously reported and could be potentially used as biomarkers for the prediction and phenotyping of PE. This case study can be used as a preliminary input and provides guidance to any future large-scale epigenetic profiling studies of PE.

5. Conclusions

Our ATAC-seq results provide an overview of gene regulatory programs and the corresponding downstream pathways associated with PE that may influence placenta function and fetal growth. In addition, we discovered novel genes with high confidence that were not previously reported and could be potentially used as biomarkers for the prediction and phenotyping of PE.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/genes13050888/s1. The main test results of healthy gravidae (before labor) are submitted as Table S1. The lists of pathways and differentially expressed peaks are submitted as Table S2 and Table S3, respectively. The heatmap of distribution of transcription factor binding loci relative to TSS, the distribution of transcription factor binding loci relative to TSS, and volcano plot of differential expressed genes are submitted as Figure S1A–C, respectively. The binding site distribution is presented in Figure S2; the ChIP-Seq data were obtained from GSM970258, GSM1278240, GSM970257, GSM651546, GSM651547, GSM1097883, GSM1816080, GSM607949, GSM607950. All bed files were merged into one bed file and then were annotated. Peaks across the promoter regions of PGTIS, SERINC2, PCMTD2, TNFRSF6B, and PRR25 are presented in Figure S3.

Author Contributions

Conceptualization, J.M., Z.Z., X.W. and H.Z.; Data curation, Y.H., Y.L., L.L., Q.S. and Q.C.; Investigation, J.M., Z.Z. and N.L.; Supervision, X.W., B.W. and H.Z.; Writing—original draft, J.M., Z.Z., X.W., B.W. and H.Z.; Writing—review and editing, J.M., Z.Z., N.L., X.W., B.W. and H.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Natural Science Foundation of Guangdong Province (2020A1515010714, 2019A1515011068); the Shenzhen Basic Research Program (JCYJ20190813160607211, JCYJ20190813132201654, JCYJ20210324122606016, JCYJ20210324114009026); the Science, Technology, and Innovation Bureau of Shenzhen Guangming District (2020R01074, 2020R01076); Shanghai Pujiang Program [grant number: 2020PJC080]; National Natural Science Foundation of China [grant number: 72004133].

Institutional Review Board Statement

This study was performed according to a protocol approved by the Ethics Committee at the University of Chinese Academy of Sciences Shenzhen Hospital. The research adhered to the tenets of the Declaration of Helsinki as amended in 2013.

Informed Consent Statement

Informed consent was obtained from all subjects involved in this study. Written informed consent has been obtained from the patients and all participants to publish this paper.

Data Availability Statement

The ATAC-seq sequence is available in GSE199479. ChIP-Seq data of GATA family transcriptional factors binding sites were obtained from GSM970258, GSM1278240, GSM970257, GSM651546, GSM651547, GSM1097883, GSM1816080, GSM607949, GSM607950.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mol, B.W.; Roberts, C.T.; Thangaratinam, S.; Magee, L.A.; De Groot, C.J.; Hofmeyr, G.J. Pre-eclampsia. Lancet 2016, 387, 999–1011. [Google Scholar] [CrossRef]
  2. Brown, M.A.; Magee, L.A.; Kenny, L.C.; Karumanchi, S.A.; McCarthy, F.P.; Saito, S.; Hall, D.R.; Warren, C.E.; Adoyi, G.; Ishaku, S. Hypertensive Disorders of Pregnancy: ISSHP Classification, Diagnosis, and Management Recommen-dations for International Practice. Hypertension 2018, 72, 24–43. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Tanner, M.S.; Davey, M.-A.; Mol, B.W.; Rolnik, D.L. The evolution of the diagnostic criteria of preeclampsia-eclampsia. Am. J. Obstet. Gynecol. 2022, 226, S835–S843. [Google Scholar] [CrossRef] [PubMed]
  4. Socha, M.; Malinowski, B.; Puk, O.; Dubiel, M.; Wiciński, M. The NLRP3 Inflammasome Role in the Pathogenesis of Pregnancy Induced Hypertension and Preeclampsia. Cells 2020, 9, 1642. [Google Scholar] [CrossRef] [PubMed]
  5. Kawasaki, K.; Kondoh, E.; Chigusa, Y.; Kawamura, Y.; Mogami, H.; Takeda, S.; Horie, A.; Baba, T.; Matsumura, N.; Mandai, M.; et al. Metabolomic Profiles of Placenta in Preeclampsia. Hypertension 2019, 73, 671–679. [Google Scholar] [CrossRef]
  6. American College of Obstetricians and Gynecologists. Hypertension in pregnancy. Report of the American College of Obstetricians and Gynecologists’ Task Force on Hypertension in Pregnancy. Obstet. Gynecol. 2013, 122, 1122–1131. [Google Scholar]
  7. Seshadri Reddy, V.; Duggina, P.; Vedhantam, M.; Manne, M.; Varma, N.; Nagaram, S. Maternal serum and fetal cord-blood ische-mia-modified albumin concentrations in normal pregnancy and preeclampsia: A systematic review and meta-analysis. J. Matern. Fetal Neonatal Med. 2018, 31, 3255–3266. [Google Scholar] [CrossRef]
  8. Vargas-Rojas, M.I.; Solleiro-Villavicencio, H.; Soto-Vega, E. Th1, Th2, Th17 and Treg levels in umbilical cord blood in preeclamp-sia. J. Matern. Fetal Neonatal Med. 2016, 29, 1642–1645. [Google Scholar] [CrossRef]
  9. Gumina, D.L.; Black, C.P.; Balasubramaniam, V.; Winn, V.D.; Baker, C.D. Umbilical Cord Blood Circulating Progenitor Cells and Endothelial Colony-Forming Cells Are Decreased in Preeclampsia. Reprod. Sci. 2017, 24, 1088–1096. [Google Scholar] [CrossRef]
  10. Sukenik-Halevy, R.; Amiel, A.; Kidron, D.; Liberman, M.; Ganor-Paz, Y.; Biron-Shental, T. Telomere homeostasis in trophoblasts and in cord blood cells from pregnancies complicated with preeclampsia. Am. J. Obstet. Gynecol. 2016, 214, 283.e1–283.e7. [Google Scholar] [CrossRef]
  11. Wojczakowski, W.; Kimber-Trojnar, Ż.; Dziwisz, F.; Słodzińska, M.; Słodziński, H.; Leszczyńska-Gorzelak, B. Preeclampsia and Cardiovascular Risk for Offspring. J. Clin. Med. 2021, 10, 3154. [Google Scholar] [CrossRef] [PubMed]
  12. Barron, A.; McCarthy, C.M.; O’Keeffe, G.W. Preeclampsia and Neurodevelopmental Outcomes: Potential Pathogenic Roles for Inflammation and Oxidative Stress? Mol. Neurobiol. 2021, 58, 2734–2756. [Google Scholar] [CrossRef] [PubMed]
  13. Holden, C.; Moses, N.; Fox, M.; Glose, G.; Vaughn, B.C.; Marshall-Taylor, S. Collaborating to Address Infant Mortality: Lessons Learned from the Brownsville Action Community for Health Equality. Prog. Community Health Partnersh. Res. Educ. Action 2011, 5, 281–288. [Google Scholar] [CrossRef] [PubMed]
  14. Teng, R.J.; Wu, T.J.; Sharma, R.; Garrison, R.D.; Hudak, M.L. Early neonatal hypotension in premature infants born to preeclamptic mothers. J. Perinatol. Off. J. Calif. Perinat. Assoc. 2006, 26, 471–475. [Google Scholar] [CrossRef]
  15. Sitras, V.; Paulssen, R.; Grønaas, H.; Leirvik, J.; Hanssen, T.; Vårtun, Å.; Acharya, G. Differential Placental Gene Expression in Severe Preeclampsia. Placenta 2009, 30, 424–433. [Google Scholar] [CrossRef]
  16. Catarino, C.; Santos-Silva, A.; Belo, L.; Rocha-Pereira, P.; Rocha, S.; Patrício, B.; Quintanilha, A.; Rebelo, I. Inflammatory Disturbances in Preeclampsia: Relationship between Maternal and Umbilical Cord Blood. J. Pregnancy 2012, 2012, 684384. [Google Scholar] [CrossRef]
  17. Rasmussen, M.; Reddy, M.; Nolan, R.; Camunas-Soler, J.; Khodursky, A.; Scheller, N.M.; Cantonwine, D.E.; Engelbrechtsen, L.; Mi, J.D.; Dutta, A.; et al. RNA profiles reveal signatures of future health and disease in pregnancy. Nature 2022, 601, 422–427. [Google Scholar] [CrossRef]
  18. Kamrani, A.; Alipourfard, I.; Ahmadi-Khiavi, H.; Yousefi, M.; Rostamzadeh, D.; Izadi, M.; Ahmadi, M. The role of epigenetic changes in preeclampsia. BioFactors 2019, 45, 712–724. [Google Scholar] [CrossRef]
  19. Apicella, C.; Ruano, C.S.M.; Méhats, C.; Miralles, F.; Vaiman, D. The Role of Epigenetics in Placental Development and the Etiology of Preeclampsia. Int. J. Mol. Sci. 2019, 20, 2837. [Google Scholar] [CrossRef] [Green Version]
  20. Espinoza, J. Abnormal fetal-maternal interactions: An evolutionary value? Obs. Gynecol. 2012, 120 Pt 1, 370–374. [Google Scholar] [CrossRef]
  21. Ching, T.; Ha, J.; Song, M.-A.; Tiirikainen, M.; Molnár, J.; Berry, M.J.; Towner, D.; Garmire, L.X. Genome-scale hypomethylation in the cord blood DNAs associated with early onset preeclampsia. Clin. Epigenet. 2015, 7, 21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Buenrostro, J.D.; Wu, B.; Chang, H.Y.; Greenleaf, W.J. ATAC-seq: A Method for Assaying Chromatin Accessibility Genome-Wide. Curr. Protoc. Mol. Biol. 2015, 109, 21.29.1–21.29.9. [Google Scholar] [CrossRef] [PubMed]
  23. Klemm, S.L.; Shipony, Z.; Greenleaf, W.J. Chromatin accessibility and the regulatory epigenome. Nat. Rev. Genet. 2019, 20, 207–220. [Google Scholar] [CrossRef]
  24. World Medical Association. World Medical Association Declaration of Helsinki: Ethical principles for medical research in-volving human subjects. JAMA 2013, 310, 2191. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Buenrostro, J.D.; Giresi, P.G.; Zaba, L.C.; Chang, H.Y.; Greenleaf, W.J. Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position. Nat. Methods 2013, 10, 1213–1218. [Google Scholar] [CrossRef] [PubMed]
  26. Heinz, S.; Benner, C.; Spann, N.; Bertolino, E.; Lin, Y.C.; Laslo, P.; Cheng, J.X.; Murre, C.; Singh, H.; Glass, C.K. Simple Combinations of Lineage-Determining Transcription Factors Prime cis-Regulatory Elements Required for Macrophage and B Cell Identities. Mol. Cell 2010, 38, 576–589. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Robinson, J.T.; Thorvaldsdóttir, H.; Turner, D.; Mesirov, J.P. igv.js: An embeddable JavaScript implementation of the Integrative Genomics Viewer (IGV). bioRxiv 2020. [Google Scholar] [CrossRef]
  28. Bu, D.; Luo, H.; Huo, P.; Wang, Z.; Zhang, S.; He, Z.; Wu, Y.; Zhao, L.; Liu, J.; Guo, J.; et al. KOBAS-i: Intelligent prioritization and exploratory visualization of biological functions for gene enrichment analysis. Nucleic Acids Res. 2021, 49, W317–W325. [Google Scholar] [CrossRef]
  29. Romano, O.; Miccio, A. GATA factor transcriptional activity: Insights from genome-wide binding profiles. IUBMB Life 2020, 72, 10–26. [Google Scholar] [CrossRef]
  30. Tijssen, M.R.; Cvejic, A.; Joshi, A.; Hannah, R.L.; Ferreira, R.; Forrai, A.; Bellissimo, D.C.; Oram, S.H.; Smethurst, P.A.; Wilson, N.K.; et al. Genome-wide Analysis of Simultaneous GATA1/2, RUNX1, FLI1, and SCL Binding in Megakaryocytes Identifies Hematopoietic Regulators. Dev. Cell 2011, 20, 597–609. [Google Scholar] [CrossRef] [Green Version]
  31. Beck, D.; Thoms, J.A.I.; Perera, D.; Schütte, J.; Unnikrishnan, A.; Knezevic, K.; Kinston, S.J.; Wilson, N.K.; O’Brien, T.A.; Göttgens, B.; et al. Genome-wide analysis of transcriptional regulators in human HSPCs reveals a densely interconnected network of coding and noncoding genes. Blood 2013, 122, e12–e22. [Google Scholar] [CrossRef] [PubMed]
  32. Kanhere, A.; Hertweck, A.; Bhatia, U.; Gökmen, M.R.; Perucha, E.; Jackson, I.; Lord, G.M.; Jenner, R.G. T-bet and GATA3 orchestrate Th1 and Th2 differentiation through lineage-specific targeting of distal regulatory elements. Nat. Commun. 2012, 3, 1268. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. He, A.; Kong, S.W.; Ma, Q.; Pu, W.T. Co-occupancy by multiple cardiac transcription factors identifies transcriptional enhancers active in heart. Proc. Natl. Acad. Sci. USA 2011, 108, 5632–5637. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Sulahian, R.; Casey, F.; Shen, J.; Qian, Z.R.; Shin, H.; Ogino, S.; A Weir, B.; Vazquez, F.; Liu, X.S.; Hahn, W.C.; et al. An integrative analysis reveals functional targets of GATA6 transcriptional regulation in gastric cancer. Oncogene 2013, 33, 5637–5648. [Google Scholar] [CrossRef] [Green Version]
  35. Tal, R. The Role of Hypoxia and Hypoxia-Inducible Factor-1Alpha in Preeclampsia Pathogenesis. Biol. Reprod. 2012, 87, 134. [Google Scholar] [CrossRef]
  36. Soleymanlou, N.; Jurisica, I.; Nevo, O.; Ietta, F.; Zhang, X.; Zamudio, S.; Post, M.; Caniggia, I. Molecular evidence of placental hypoxia in preeclampsia. J. Clin. Endocrinol. Metab. 2005, 90, 4299–4308. [Google Scholar] [CrossRef] [Green Version]
  37. Hu, X.-Q.; Zhang, L. Hypoxia and the integrated stress response promote pulmonary hypertension and preeclampsia: Implications in drug development. Drug Discov. Today 2021, 26, 2754–2773. [Google Scholar] [CrossRef]
  38. Backes, C.H.; Markham, K.; Moorehead, P.; Cordero, L.; Nankervis, C.A.; Giannone, P.J. Maternal Preeclampsia and Neonatal Outcomes. J. Pregnancy 2011, 2011, 214365. [Google Scholar] [CrossRef]
  39. Lee, A.; Chow, B.L.; Skeith, L.; A Nicholas, B.J.; Poon, M.-C.; Poole, M.A.W.; O Agbani, P.E. Platelet Membrane Procoagulation in Preeclampsia. Blood 2020, 136 (Suppl. S1), 7. [Google Scholar] [CrossRef]
  40. Macey, M.; Bevan, S.; Alam, S.; Verghese, L.; Agrawal, S.; Beski, S.; Thuraisingham, R.; MacCallum, P. Platelet activation and endogenous thrombin potential in pre-eclampsia. Thromb. Res. 2010, 125, e76–e81. [Google Scholar] [CrossRef]
  41. Rambaldi, M.P.; Weiner, E.; Mecacci, F.; Bar, J.; Petraglia, F. Immunomodulation and preeclampsia. Best Pract. Research. Clin. Obstet. Gynaecol. 2019, 60, 87–96. [Google Scholar] [CrossRef]
  42. Kanellopoulos-Langevin, C.; Caucheteux, S.M.; Verbeke, P.; Ojcius, D.M. Tolerance of the fetus by the maternal immune system: Role of inflammatory mediators at the feto-maternal interface. Reprod. Biol. Endocrinol. 2003, 1, 121. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Van Nieuwenhoven, A.V.; Heineman, M.; Faas, M. The immunology of successful pregnancy. Hum. Reprod. Updat. 2003, 9, 347–357. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Saito, S.; Sakai, M.; Sasaki, Y.; Tanebe, K.; Tsuda, H.; Michimata, T. Quantitative analysis of peripheral blood Th0, Th1, Th2 and the Th1:Th2 cell ratio during normal human pregnancy and preeclampsia. Clin. Exp. Immunol. 1999, 117, 550–555. [Google Scholar] [CrossRef] [PubMed]
  45. Tindemans, I.; Serafini, N.; Di Santo, J.; Hendriks, R.W. GATA-3 Function in Innate and Adaptive Immunity. Immunity 2014, 41, 191–206. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Ranganath, S.; Murphy, K.M. Structure and Specificity of GATA Proteins in Th2 Development. Mol. Cell. Biol. 2001, 21, 2716–2725. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Romao-Veiga, M.; Ribeiro, V.R.; Matias, M.L.; Nunes, P.R.; Romagnoli, G.G.; Peracoli, J.C.; Peracoli, M.T.S. DAMPs are able to skew CD4+ T cell subsets and increase the inflammatory profile in pregnant women with preeclampsia. J. Reprod. Immunol. 2022, 149, 103470. [Google Scholar] [CrossRef] [PubMed]
  48. Lin, M.-C.; Lin, J.-J.; Hsu, C.-L.; Juan, H.-F.; Lou, P.-J.; Huang, M.-C. GATA3 interacts with and stabilizes HIF-1α to enhance cancer cell invasiveness. Oncogene 2017, 36, 4243–4252. [Google Scholar] [CrossRef]
Genes 13 00888 g001 550
Figure 1. Immune, myocardiopathy, and hypoxia response-related pathways are up-regulated in superimposed PE (a) Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of genes with higher chromatin accessibility (b–e) The volcano plot of immune, hypoxia response, hypertension, and platelet activation related gene. HIF-1 stands for Hypoxia Inducible Factor-1.
Figure 1. Immune, myocardiopathy, and hypoxia response-related pathways are up-regulated in superimposed PE (a) Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of genes with higher chromatin accessibility (b–e) The volcano plot of immune, hypoxia response, hypertension, and platelet activation related gene. HIF-1 stands for Hypoxia Inducible Factor-1.
Genes 13 00888 g001
Genes 13 00888 g002 550
Figure 2. Genomic locations of the differentially expressed peaks. (a) Proportions of genomic locations. UTR stands for the untranslated region. (b) Proportions of peak-to-TSS distances.
Figure 2. Genomic locations of the differentially expressed peaks. (a) Proportions of genomic locations. UTR stands for the untranslated region. (b) Proportions of peak-to-TSS distances.
Genes 13 00888 g002
Genes 13 00888 g003 550
Figure 3. GATA-binding motifs are enriched in superimposed PE (a) The enriched motifs. The analysis was conducted by HOMER with parameters set to ‘find regulatory element’. (b) The enriched GATA family motifs.
Figure 3. GATA-binding motifs are enriched in superimposed PE (a) The enriched motifs. The analysis was conducted by HOMER with parameters set to ‘find regulatory element’. (b) The enriched GATA family motifs.
Genes 13 00888 g003
Genes 13 00888 g004 550
Figure 4. Identification of novel PE-associated genes (a) The M-versus-A plot (MA plot) of peaks. (b) Heatmap of possible PE-associated genes.
Figure 4. Identification of novel PE-associated genes (a) The M-versus-A plot (MA plot) of peaks. (b) Heatmap of possible PE-associated genes.
Genes 13 00888 g004
Genes 13 00888 g005 550
Figure 5. The mechanism map.
Figure 5. The mechanism map.
Genes 13 00888 g005
Table
Table 1. The main abnormal test result of the patient.
Table 1. The main abnormal test result of the patient.
Before LaborAfter LaborReference Range
Urine
Protein+ +−−
24 h urine protein (mg/24 h)876.58 ↑↑−<150
Blood *
BNP (pg/mL)909 ↑↑NA<100
PT (s)9.3 ↓10.310~15
FIB (g/L)3.88 ↑5.75 ↑1.8~3.5
PT-R (%)121.10 ↑111.7075~120
D-D (mg/L)0.79 ↑1.12 ↑<0.5
* BNP stands for B-type Natriuretic Peptide. PT stands for prothrombin time. FIB stands for fibrinogen. PT-R stands for prothrombin time activity. D-D stands for D-Dimer. ↑ means higher than normal value. ↓ means lower than normal value. + means positive, − means negative.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ma, J.; Zhan, Z.; Li, N.; Huang, Y.; Li, Y.; Liu, L.; Shen, Q.; Chu, Q.; Wang, X.; Wu, B.; et al. Preliminary Interpretations of Epigenetic Profiling of Cord Blood in Preeclampsia. Genes 2022, 13, 888. https://0-doi-org.brum.beds.ac.uk/10.3390/genes13050888

AMA Style

Ma J, Zhan Z, Li N, Huang Y, Li Y, Liu L, Shen Q, Chu Q, Wang X, Wu B, et al. Preliminary Interpretations of Epigenetic Profiling of Cord Blood in Preeclampsia. Genes. 2022; 13(5):888. https://0-doi-org.brum.beds.ac.uk/10.3390/genes13050888

Chicago/Turabian Style

Ma, Junrui, Zhongqun Zhan, Ning Li, Yanli Huang, Yan Li, Lu Liu, Qi Shen, Qiao Chu, Xiaonan Wang, Benqing Wu, and et al. 2022. "Preliminary Interpretations of Epigenetic Profiling of Cord Blood in Preeclampsia" Genes 13, no. 5: 888. https://0-doi-org.brum.beds.ac.uk/10.3390/genes13050888

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop