- Open Access
H3K27 acetylation and gene expression analysis reveals differences in placental chromatin activity in fetal growth restriction
© The Author(s). 2018
Received: 17 October 2017
Accepted: 29 May 2018
Published: 26 June 2018
Posttranslational modification of histone tails such as histone 3 lysine 27 acetylation (H3K27ac) is tightly coupled to epigenetic regulation of gene expression. To explore whether this is involved in placenta pathology, we probed genome-wide H3K27ac occupancy by chromatin immunoprecipitation sequencing (ChIP-seq) in healthy placentas and placentas from pathological pregnancies with fetal growth restriction (FGR). Furthermore, we related specific acetylation profiles of FGR placentas to gene expression changes.
Analysis of H3K27ac occupancy in FGR compared to healthy placentas showed 970 differentially acetylated regions distributed throughout the genome. Principal component analysis and hierarchical clustering revealed complete segregation of the FGR and control group. Next, we identified 569 upregulated genes and 521 downregulated genes in FGR placentas by RNA sequencing. Differential gene transcription largely corresponded to expected direction based on H3K27ac status. Pathway analysis on upregulated transcripts originating from hyperacetylated sites revealed genes related to the HIF-1-alpha transcription factor network and several other genes with known involvement in placental pathology (LEP, FLT1, HK2, ENG, FOS). Downregulated transcripts in the vicinity of hypoacetylated sites were related to the immune system and growth hormone receptor signaling. Additionally, we found enrichment of 141 transcription factor binding motifs within differentially acetylated regions. Of the corresponding transcription factors, four were upregulated, SP1, ARNT2, HEY2, and VDR, and two downregulated, FOSL and NR4A1.
We demonstrate a key role for genome-wide alterations in H3K27ac in FGR placentas corresponding with changes in transcription profiles of regions relevant to placental function. Future studies on the role of H3K27ac in FGR and placental-fetal development may help to identify novel targets for therapy of this currently incurable disease.
The dynamics of histone 3 lysine 27 acetylation (H3K27ac) in DNA regulatory regions is one of the components playing a fundamental role in the precise timing and level of gene transcription [1, 2]. Consequently, aberrant H3K27ac has been suggested to be involved in disease pathology by eliciting pathological gene expression programs [3, 4]. H3K27ac marks both active promoters and distal enhancers, the most important and best understood regulatory domains. To study involvement of this regulatory level/layer in placental pathology, we mapped H3K27ac occupancy in healthy placentas and placentas from pregnancies with fetal growth restriction (FGR).
FGR, a condition in which the fetus is unable to achieve its full growth potential through inadequate supply of nutrients and growth factors, occurs in approximately 5% of pregnancies . FGR imposes a major risk of perinatal morbidity and mortality  and programs the health of the fetus throughout life, by being associated with a future risk of type 2 diabetes and cardiovascular and renal disease [7, 8]. At the histopathological level, FGR placenta exhibit signs of disrupted placental development characterized by increased infarction area, increased syncytial knotting, inflammation, and impaired trophoblast invasion into the spiral arteries of the uterus due to inappropriate maternal-fetal immune interaction [9–11].
Previous studies of placentas of FGR pregnancies have reported differences in gene transcription across many regions across the genome. Pathways associated with altered placental gene expression in FGR include angiogenesis, immune modulation, energy production, and growth signaling [12–15]. Most adaptive responses of the placenta, e.g., to support restricted fetal growth, are thought to result from changes in epigenetic regulation [16–18]. Recently, DNA methylation has been mapped in the human placenta  and a number of studies have suggested differences in genome-wide methylation profiles in FGR placentas [20–22]. We assume that disruptions in other epigenetic systems (e.g., histone modifications and other posttranslational chromatin modifications) regulating placenta gene expression are also likely to be involved [23, 24].
In this study, we mapped differential H3K27ac profiles in DNA regulatory regions in relation to disrupted development of the placenta, by exploring the presence of H3K27ac using chromatin immunoprecipitation sequencing (ChIP-seq) in FGR and control placentas. Next, we performed RNA-seq to examine whether differences in H3K27ac also reflects gene expression levels. With this approach, we identified previously unstudied alterations in promoter and enhancer activity related to placenta pathology.
Detection of regions with differential H3K27ac occupancy
Genes annotated to within 20 kb of differentially acetylated regions
Pathway analysis of differential acetylated regions in FGR placentas
Cross-validation of functional effects by combining H3K27ac and mRNA profiles
To cross-validate functional meaning of the differentially acetylated peaks and to examine whether hyper- or hypoacetylated state of H3K27 was also accompanied by differential gene expression, we performed RNA-seq. We identified 569 upregulated genes and 521 downregulated genes in FGR vs. control (full list supplied in Additional file 1: Table S4, heatmap and MA and V plots in Additional file 3). These gene expression profiles largely overlapped with findings from earlier gene expression studies in FGR placentas using microarrays [12–15]. This is also reflected in the biological processes and pathways associated with the up- and downregulated genes supplied in Additional file 1: Table S5. These findings confirm disruption of gene expression in placenta from FGR pregnancies involved in important processes of placental development in angiogenesis and immune modulation.
Identification of pathways by combining H3K27ac profiles and gene expression levels
Pathway with TOPPfun on upregulated genes near hyperacetylated sites revealed genes to be related to the HIF-1-alpha transcription factor network (Fig. 4c, Additional file 1: Table S6). In addition, the downregulated transcripts and hypoacetylated regions could be grouped as genes encoding secreted soluble factors, extracellular matrix proteins, and molecules involved in growth hormone receptor signaling (Fig. 4c, Additional file 1: Table S6).
Transcription factor binding motif analysis
In this study, we used ChIP-seq and RNA-seq to perform an in-depth analysis of DNA regions with differential chromatin activity in healthy human placentas and pathological placentas of pregnancies affected by FGR. With this combined approach, we were able to identify H3K27 acetylation as a key additional layer involved in epigenetic regulation of gene expression in placental function which in turn might be related to impaired fetal growth.
In summary, our findings confirm that, in FGR, the placenta exhibits substantial genome-wide alterations in H3K27 acetylation, with corresponding changes in transcription profiles in several regions pertinent to placental development and presumably function. The identified acetylation patterns show a clear distinction between FGR and healthy placentas, and the genes annotated to differentially acetylated regions are involved in pathways known to be affected in FGR [27–29]. We confirmed functionality of the differentially acetylated regions by showing a clear correlation between the acetylation profiles with differential gene expression. In-depth analyses of the differences in FGR placenta that revealed candidate genes and pathways that fit previously reported histopathological and protein data. For example, we found hyperacetylation and higher mRNA of the sFlt-1 region in FGR placentas, consistent with studies reporting upregulation of sFlt-1 protein in FGR placentas [36, 37]. Furthermore, sFlt-1 is known to be enriched in syncytial knots, which are also observed in FGR placenta . Similarly, the leptin protein was showed to be upregulated in FGR placentas [39, 40] and GH and CSH proteins are often downregulated in FGR placentas [41, 42]. Additionally, we unmasked an interplay between differentially transcribed transcription factors with enriched TFBM within the differentially acetylated regions. Together, our data suggest an important role of H3K27ac as an additional layer/level in the regulation of placental gene expression. Our findings confirm that a combined ChIP-seq and RNA-seq approach may provide a useful approach to probe the pathophysiology of placental disease and discover novel targets to improve placental health and thus support fetal growth and development.
To date, we are not aware of any previous reports on histone modification in FGR placentas. Other groups have focused attention on epigenetic regulation of placentas in FGR involving DNA methylation. Genome-wide methylation studies in placenta of FGR suggest involvement of pathways associated with lipid metabolism, transcription, and cadherin and Wnt signaling [20–22]. Others have reported on DNA methylation of imprinted loci, since unbalanced placental expression of imprinted genes has been reported in FGR placenta [12, 43–45]. However, methylation actually fluctuates very little at these loci . While DNA methylation may contribute to placental gene regulation, it can only in part explain differences in gene expression [23, 24]. Together with our results, this emphasizes the need to include other epigenetic processes to fully appreciate the complexity of mechanisms underlying placental gene expression and (primary or secondary adaptive) gene responses during (disrupted) early human development.
Mechanisms involved in differential acetylation of histones within the human placenta have not been studied. In view of the “fetal origins of adult disease” hypothesis (also known as the Barker hypothesis), i.e., the concept that exposure to intrauterine conditions may have long-lasting effects on fetal development, these epigenetic alterations might well be the result of changes in the intrauterine environment associated with FGR . One of the intrauterine factors that might play a role in changes in epigenetic profiles in the placenta of FGR is the timing and degree of hypoxia. For example, it was showed previously that the proliferation of trophoblasts is highly dependent on the crosstalk between HIFs and histone deacetylases (HDACs) in response to hypoxia . In addition, differences in histone acetylation may be the in part driven by variation in the genetic code itself, e.g., single-nucleotide polymorphisms, that shape chromatin architecture and thereby disrupt normal placental development [47, 48]. Considering potential intervention, it is of great interest to study the contribution of both environmental and intrinsic factors in more detail.
Our study has several strengths. First, our data provide new and comprehensive information on epigenetic regulation of placental gene expression in FGR. Our data suggest that H3K27ac profiles have a role in genome-wide regulation of gene expression in FGR placenta and point towards disruptions in important pathways involved in placental development. Histone modifications represent mid-long-term effects  and are likely more stable and independent of sample moment compared to RNA, and by combining CHIP-seq data with RNA-seq data, we were able to relate epigenetic marks to functional gene expression patterns, confirming relevance of identified pathways. Another strength is our careful collection and detailed phenotyping of FGR placentas during Cesarean sections, excluding those women who have labored. Thus, we avoided confounding effects of temporal changes associated with parturition .
There are some limitations in our approach that also need to be addressed. First, we could not avoid the limitation of the use of nearest gene approach for functional annotation. Here, we used a 20-kb window to annotate genes which, although quite wide, still allowed identification of specific pathways and genes that corresponded with previous literature. With regard to the collected material, there were unavoidable difference in gestational ages in material from FGR and control due to the severity of the selected FGR cases. It is possible that gestational age in itself has an effect on H3K27ac and gene expression as shown in early mouse placenta  and for other epigenetic marks such as methylation as shown in human placenta . While gestational age might be a confounder, the strongest determinant in our cohort is likely to be the disease, given that all placentas were from third trimester pregnancies while majority of gene expression changes occur earlier in development  and the confirmed disease-specific pathways. In addition, largest differences in epigenetics in the third trimester are mainly induced by parturition , which we ruled out by collecting samples from C-sections only.
With this study, we show that genome-wide H3K27ac profiles might be useful to better understand pathophysiology of FGR. Therefore, we would recommend future studies to confirm the identified acetylation targets in a larger population as this might be helpful to address association of the differentially acetylated regions with different subtypes of FGR, gestational-age specific effects on acetylation within the third trimester of healthy and FGR pregnancy, and, most importantly, to explore cause-effect relationships. The latter might be very valuable as the plasticity of histone marks form an attractive strategy for intervention, especially since no therapy is currently available to improve placental function. In addition, it will be valuable to study the combination of other epigenetic marks, such as DNA methylation and other histone marks, to better discriminate promoters from enhancers (H3K4me3, H3K4me1) or repressed chromatin (H3K27me3) . Moreover, since placental samples represent an average state across different cellular compartments with trophoblast being the predominant cell type , newer techniques looking at subpopulations of cells might help to unravel cell-specific deregulated pathways .
We demonstrate involvement of H3K27ac in key regions related to the placental function and placental pathology associating with impaired fetal growth. Our findings underscore that an approach using combined CHIP-seq and RNA-seq analyses facilitates discovery of novel genes involved in FGR, and identification of molecular pathways and processes associated with the placental function and early (disrupted) human growth and development. Future studies on unraveling the role of H3K27ac in FGR pregnancies and placental-fetal development may help to identify novel targets for therapy of this currently incurable disease.
Study design and sample collection
Placenta biopsies were collected immediately after birth from women with a FGR pregnancy and women with uncomplicated pregnancy. Each placenta was sampled at four random locations. For this study, we only included women undergoing a primary Cesarean section. We defined FGR as an estimated fetal weight <p3 and only included cases with pulsatility index of the uterine artery >p97.5 since we were particularly interested in cases with FGR based on placenta insufficiency. Controls were included when planned for Cesarean section because of either breech presentation or history of Caesarian section. Additional file 4 provides an overview of characteristics from the pregnancies of which the placenta were derived. The material was snap-frozen directly after sampling and stored at − 80 °C. Before further processing, all four biopsies per placenta were pooled and ground into powder using mortar and pestle cooled with liquid nitrogen.
ChIP-sequencing for H3K27ac occupancy
The powdered samples, two scoops diluted in 700 μl sPBS, were used for chromatin isolation using the MAGnify™ Chromatin Immunoprecipitation System kit (Life Technologies, Thermo Fisher Scientific, Carlsbad, CA) according to manufacturer’s instructions. In brief, the tissue was crosslinked with 1% formaldehyde and the crosslinking was stopped by adding 1.25 M glycine. Cells were lysed using the kit-provided lysis buffer, and nuclei were sonicated using Covaris microTUBE (duty cycle 5%, intensity 2, 200 cycles per burst, 60 s cycle time, 10 cycles). We aimed for DNA fragments of 200–2000 bp long. We continued with one fifth of the volume, and the sheared chromatin was diluted and then incubated with 1 μl anti-H3K27ac antibody (ab4729, Abcam) pre-coupled to magnetic beads for 2 h at 4 °C. Beads were extensively washed, and crosslinking was reversed by the kit-provided reverse crosslinking buffer with proteinase K. DNA was purified using ChIP DNA Clean & Concentrator kit (Zymo Research). Libraries were prepared using the NEXTflex™ Rapid DNA Sequencing Kit (Bioo Scientific). Samples were PCR-amplified and checked for the proper size range and for the absence of adaptor dimers on a 2% agarose gel, and barcoded libraries were sequenced 75 bp single-end on Illumina NextSeq500 sequencer (Utrecht Sequencing Facility).
Mapping of ChIP-sequencing reads
Sequencing reads were mapped against the reference genome (hg19 assembly, GRCh37) using BWA package (mem –t 7 –c 100 –M –R) . Multiple reads mapping to the same location and strand have been collapsed to single reads, and only uniquely placed reads were used for peak/region calling. Regions were called using Cisgenome 2.0 ( –e 150 -maxgap 200 –minlen 200) . Subsequently, to obtain a common reference, region coordinates from all FGR and control samples were stretched to at least 2000 bp and collapsed into a single common list. Overlapping regions were merged based on their outmost coordinates. Only the autosomal regions supported by at least two independent datasets were further analyzed. Sequencing reads from each ChIP-seq library were overlapped with the common region list, to set the H3K27ac occupancy for every region-sample pair.
ChIP data analysis
Regions with differential H3K27ac occupancy between FGR samples and controls were identified using DESeq2 (p < 0.05 by Wald test)  and are referred to as “differentially acetylated regions.” Hierarchical clustering based on differentially acetylated regions was performed with quantile-normalized, log2-transformed, and median-centered read counts per common region. To avoid log2 transformation of zero values, one read was added to each region. A Manhattan plot was created representing the distribution of regions detected to be differentially acetylated in FRG vs. controls across autosomal chromosomes. The unsupervised PCA analysis was performed for the top 500 most variable regions using DESeq2. Based on the initial analysis, we excluded one control sample because of being an extreme outlier in the PCA analysis and continued the analysis with n = 5 FGR and n = 4 control samples.
To assign biological meaning to differentially acetylated peaks, we used three approaches. First, we annotated genes to the differentially acetylated sites in silico using a window of ± 20 kb from transcription start site (TSS). Second, GREAT software (Stanford) was used to assign biological meaning to a set of non-coding genomic regions by analyzing the annotations of genes flanking differentially acetylated regions . GREAT incorporates annotations from 20 ontologies and accounts for the length of gene regulatory domains. ToppFun was used for gene list enrichment analysis and candidate gene prioritization based on functional annotations and protein interaction networks (accessed in March 2017). In addition, we identified whether the differentially acetylated regions contain enrichment of specific transcription factor binding domains by overlapping the differentially hypo- and hyperacetylated regions in the FGR vs. control group to placenta DNAse hypersensitivity site (DHS) datasets obtained from the ENCODE database (ENCFF203HVV, ENCFF249GZW, ENCFF919NRH). The genomic sequence of overlapping DHS was repeat masked, and the enrichment of TFBM was calculated against a random set of non-overlapping DHS sequences using the Analysis Motif of Enrichment (AME tool) of the MEME Suite with default settings using human (HOCOMOCO v9) motif database.
Total RNA was extracted from placental powder using RNeasy® (Qiagen, Hilden, Germany) according to the manufacturer’s instructions by which all RNA molecules longer than 200 nucleotides are purified. Polyadenylated mRNA fraction was isolated using Poly(A) Beads (NEXTflex, San Jose, CA), and sequencing libraries were constructed using the Rapid Directional RNA-seq kit (NEXTflex, San Jose, CA). Libraries were sequenced on the Nextseq500 platform (Illumina, San Diego, CA), producing single-end reads of 75 bp (Utrecht Sequencing Facility). Reads were aligned to the human reference genome GRCh37 using STAR version 2.4.2a. Picard’s AddOrReplaceReadGroups (v1.98) was used to add read groups to the BAM files, which were sorted with Sambamba v0.4.5, and transcript abundances were quantified with HTSeq-count version 0.6.1p1 using the union mode. Subsequently, reads per kilobase million reads sequenced (RPKMs) were calculated with edgeR’s RPKM function.
RNA-seq data analysis
Differentially expressed genes were identified using the DESeq2 package with standard settings. Genes with padj < 0.05 were considered as differentially expressed. Again, ToppFun was used for gene list enrichment analysis and candidate gene prioritization (accessed in March 207). The list of up- and downregulated genes were tested separately using probability density function p value calculation, FDR B&H correction, p value cutoff 0.05.
ChIP-seq and RNA-seq overlap analyses
To validate functional consequences of differentially acetylated regions, we overlapped the genes annotated to differentially hyperacetylated regions with upregulated genes derived from the mRNA-seq and genes annotated to differentially hypoacetylated regions with downregulated genes derived from the mRNA-seq. Functional pathway analysis with TOPPFUN were performed on both sets of overlapping genes. In addition, we searched the up- and downregulated genes for transcription factors (TF) with enriched TFBM within the differentially acetylated regions.
Data are shown as mean ± SD or median (range) and analyzed with Student’s t test or χ2 test.
The authors acknowledge N. van Dungen and N. Lansu for their technical support.
This study was funded by Fonds Gezond Geboren (personal grant awarded to BBR) and Dutch Kidney Foundation (15O141).
Availability of data and materials
The datasets generated and/or analyzed are partly available within the article as supplemental files. Other raw data is available from the corresponding author on reasonable request.
BBR, AF, and ATL contributed to the research idea. BBR and ATL contributed to the study design and sample collection. NDP contributed to the sample processing. MM and NDP contributed to the analysis of data. MM contributed to the statistics. NDP contributed to the drafting of the manuscript. BBR, PN, ATL, JAJ, MM, and AF contributed to the editing of the manuscript. All authors read and approved the final manuscript.
Ethics approval and consent to participate
This study was approved by the Medical Ethics Committee of the University Medical Center Utrecht, Utrecht, the Netherlands (protocol number NL57021.041.16, METC 16-198), and written informed consent was obtained.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Ernst J, Kheradpour P, Mikkelsen TS, Shoresh N, Ward LD, Epstein CB, et al. Mapping and analysis of chromatin state dynamics in nine human cell types. Nature. 2011;473:43–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Creyghton MP, Cheng AW, Welstead GG, Kooistra T, Carey BW, Steine EJ, et al. Histone H3K27ac separates active from poised enhancers and predicts developmental state. Proc Natl Acad Sci U S A. 2010;107:21931–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Mirabella AC, Foster BM, Bartke T. Chromatin deregulation in disease. Chromosoma. 2016;125:75–93.View ArticlePubMedGoogle Scholar
- Hnisz D, Abraham BJ, Lee TI, Lau A, Saint-Andre V, Sigova AA, et al. Super-enhancers in the control of cell identity and disease. Cell. 2013;155:934–47.View ArticlePubMedGoogle Scholar
- Romo A, Carceller R, Tobajas J. Intrauterine growth retardation (IUGR): epidemiology and etiology. Pediatr Endocrinol Rev. 2009;6(3):332–6.PubMedGoogle Scholar
- Bernstein IM, Horbar JD, Badger GJ, Ohlsson A, Golan A. Morbidity and mortality among very-low-birth-weight neonates with intrauterine growth restriction. The Vermont Oxford Network. Am J Obstet Gynecol. 2000;182(1 1):198–206.View ArticlePubMedGoogle Scholar
- Barker DJ, Osmond C, Golding J, Kuh D, Wadsworth ME. Growth in utero, blood pressure in childhood and adult life, and mortality from cardiovascular disease. BMJ. 1989;298:564–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Kermack AJ, Van Rijn BB, Houghton FD, Calder PC, Cameron IT, Macklon NS. The “developmental origins” hypothesis: relevance to the obstetrician and gynecologist. J Dev Orig Health Dis. 2015;6:415–24.View ArticlePubMedGoogle Scholar
- Chaddha V, Viero S, Huppertz B, Kingdom J. Developmental biology of the placenta and the origins of placental insufficiency. Semin Fetal Neonatal Med. 2004;9:357–69.View ArticlePubMedGoogle Scholar
- Veerbeek JHW, Brouwers L, Koster MPH, Koenen SV, van Vliet EOG, Nikkels PGJ, et al. Spiral artery remodeling and maternal cardiovascular risk: the spiral artery remodeling (SPAR) study. J Hypertens. 2016;34:1570–7.View ArticlePubMedGoogle Scholar
- Veerbeek JHW, Nikkels PGJ, Torrance HL, Gravesteijn J, Post Uiterweer ED, Derks JB, et al. Placental pathology in early intrauterine growth restriction associated with maternal hypertension. Placenta. 2014;35:696–701.View ArticlePubMedGoogle Scholar
- McMinn J, Wei M, Schupf N, Cusmai J, Johnson EB, Smith AC, et al. Unbalanced placental expression of imprinted genes in human intrauterine growth restriction. Placenta. 2006;27:540–9.View ArticlePubMedGoogle Scholar
- McCarthy C, Cotter FE, McElwaine S, Twomey A, Mooney EE, Ryan F, et al. Altered gene expression patterns in intrauterine growth restriction: potential role of hypoxia. Am J Obstet Gynecol. 2007;196:70.e1–6.View ArticleGoogle Scholar
- Struwe E, Berzl G, Schild R, Blessing H, Drexel L, Hauck B, et al. Microarray analysis of placental tissue in intrauterine growth restriction. Clin Endocrinol. 2010;72:241–7.View ArticleGoogle Scholar
- Madeleneau D, Buffat C, Mondon F, Grimault H, Rigourd V, Tsatsaris V, et al. Transcriptomic analysis of human placenta in intrauterine growth restriction. Pediatr Res. 2015;77:799–807.View ArticlePubMedGoogle Scholar
- Nelissen ECM, van Montfoort APA, Dumoulin JCM, Evers JLH. Epigenetics and the placenta. Hum Reprod Update. 2011;17:397–417.View ArticlePubMedGoogle Scholar
- Bianco-Miotto T, Mayne BT, Buckberry S, Breen J, Rodriguez Lopez CM, Roberts CT. Recent progress towards understanding the role of DNA methylation in human placental development. Reproduction. 2016;152:R23–30.View ArticlePubMedPubMed CentralGoogle Scholar
- Lewis RM, Cleal JK, Hanson MA. Review: placenta, evolution and lifelong health. Placenta. 2012;33:S28–32. https://doi.org/10.1016/j.placenta.2011.12.003.View ArticlePubMedGoogle Scholar
- Schroeder DI, Blair JD, Lott P, Yu HOK, Hong D, Crary F, et al. The human placenta methylome. Proc Natl Acad Sci U S A. 2013;110:6037–42.View ArticlePubMedPubMed CentralGoogle Scholar
- Roifman M, Choufani S, Turinsky AL, Drewlo S, Keating S, Brudno M, et al. Genome-wide placental DNA methylation analysis of severely growth-discordant monochorionic twins reveals novel epigenetic targets for intrauterine growth restriction. Clin Epigenetics. 2016;8:70. https://doi.org/10.1186/s13148-016-0238-x.View ArticlePubMedPubMed CentralGoogle Scholar
- Hillman SL, Finer S, Smart MC, Mathews C, Lowe R, Rakyan VK, et al. Novel DNA methylation profiles associated with key gene regulation and transcription pathways in blood and placenta of growth-restricted neonates. Epigenetics. 2015;10:50–61. https://doi.org/10.4161/15592294.2014.989741.View ArticlePubMedPubMed CentralGoogle Scholar
- Lambertini L, Lee T-L, Chan W-Y, Lee M-J, Diplas A, Wetmur J, et al. Differential methylation of imprinted genes in growth-restricted placentas. Reprod Sci. 2011;18:1111–7.View ArticlePubMedGoogle Scholar
- Joo JE, Hiden U, Lassance L, Gordon L, Martino DJ, Desoye G, et al. Variable promoter methylation contributes to differential expression of key genes in human placenta-derived venous and arterial endothelial cells. BMC Genomics. 2013;14:475.View ArticlePubMedPubMed CentralGoogle Scholar
- Lopez-Abad M, Iglesias-Platas I, Monk D. Epigenetic characterization of CDKN1C in placenta samples from non-syndromic intrauterine growth restriction. Front Genet. 2016;7:62.View ArticlePubMedPubMed CentralGoogle Scholar
- Kaartokallio T, Cervera A, Kyllonen A, Laivuori K, Kere J, Laivuori H. Gene expression profiling of pre-eclamptic placentae by RNA sequencing. Sci Rep. 2015;5:14107.View ArticlePubMedPubMed CentralGoogle Scholar
- Kleinrouweler CE, van Uitert M, Moerland PD, Ris-Stalpers C, van der Post JAM, Afink GB. Differentially expressed genes in the pre-eclamptic placenta: a systematic review and meta-analysis. PLoS One. 2013;8:e68991.View ArticlePubMedPubMed CentralGoogle Scholar
- Gourvas V, Dalpa E, Konstantinidou A, Vrachnis N, Spandidos DA, Sifakis S. Angiogenic factors in placentas from pregnancies complicated by fetal growth restriction (review). Mol Med Rep. 2012;6:23–7.PubMedGoogle Scholar
- Kimura C, Watanabe K, Iwasaki A, Mori T, Matsushita H, Shinohara K, et al. The severity of hypoxic changes and oxidative DNA damage in the placenta of early-onset preeclamptic women and fetal growth restriction. J Matern Fetal Neonatal Med. 2013;26:491–6.View ArticlePubMedGoogle Scholar
- Prins JR, Faas MM, Melgert BN, Huitema S, Timmer a HMN, et al. Altered expression of immune-associated genes in first-trimester human decidua of pregnancies later complicated with hypertension or foetal growth restriction. Placenta. 2012;33:453–5. https://doi.org/10.1016/j.placenta.2012.02.010.View ArticlePubMedGoogle Scholar
- McLeay RC, Bailey TL. Motif enrichment analysis: a unified framework and an evaluation on ChIP data. BMC Bioinformatics. 2010;11:165.View ArticlePubMedPubMed CentralGoogle Scholar
- Mandl M, Depping R. Hypoxia-inducible aryl hydrocarbon receptor nuclear translocator (ARNT) (HIF-1beta): is it a rare exception? Mol Med. 2014;20:215–20.View ArticlePubMedPubMed CentralGoogle Scholar
- Diez H, Fischer A, Winkler A, Hu C-J, Hatzopoulos AK, Breier G, et al. Hypoxia-mediated activation of Dll4-Notch-Hey2 signaling in endothelial progenitor cells and adoption of arterial cell fate. Exp Cell Res. 2007;313:1–9.View ArticlePubMedGoogle Scholar
- Ma R, Gu Y, Zhao S, Sun J, Groome LJ, Wang Y. Expressions of vitamin D metabolic components VDBP, CYP2R1, CYP27B1, CYP24A1, and VDR in placentas from normal and preeclamptic pregnancies. Am J Physiol Endocrinol Metab. 2012;303:E928–35.View ArticlePubMedPubMed CentralGoogle Scholar
- Kent LN, Rumi MAK, Kubota K, Lee D-S, Soares MJ. FOSL1 is integral to establishing the maternal-fetal interface. Mol Cell Biol. 2011;31:4801–13.View ArticlePubMedPubMed CentralGoogle Scholar
- Soares MJ, Chakraborty D, Renaud SJ, Kubota K, Bu P, Konno T, et al. Regulatory pathways controlling the endovascular invasive trophoblast cell lineage. J Reprod Dev. 2012;58:283–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Nevo O, Many A, Xu J, Kingdom J, Piccoli E, Zamudio S, et al. Placental expression of soluble FMS-like tyrosine kinase 1 is increased in singletons and twin pregnancies with intrauterine growth restriction. J Clin Endocrinol Metab. 2008;93:285–92.View ArticlePubMedGoogle Scholar
- Hoeller A, Ehrlich L, Golic M, Herse F, Perschel FH, Siwetz M, et al. Placental expression of sFlt-1 and PlGF in early preeclampsia vs. early IUGR vs. age-matched healthy pregnancies. Hypertens pregnancy. 2017;36:151–60.View ArticlePubMedGoogle Scholar
- Rajakumar A, Cerdeira AS, Rana S, Zsengeller Z, Edmunds L, Jeyabalan A, et al. Transcriptionally active syncytial aggregates in the maternal circulation may contribute to circulating soluble FMS-like tyrosine kinase 1 in preeclampsia. Hypertens (Dallas, Tex 1979). 2012;59:256–64.View ArticleGoogle Scholar
- Schrey S, Kingdom J, Baczyk D, Fitzgerald B, Keating S, Ryan G, et al. Leptin is differentially expressed and epigenetically regulated across monochorionic twin placenta with discordant fetal growth. Mol Hum Reprod. 2013;19:764–72.View ArticlePubMedGoogle Scholar
- Li RHW, Poon SCS, Yu MY, Wong YF. Expression of placental leptin and leptin receptors in preeclampsia. Int J Gynecol Pathol. 2004;23:378–85.View ArticlePubMedGoogle Scholar
- Velegrakis A, Sfakiotaki M, Sifakis S. Human placental growth hormone in normal and abnormal fetal growth. Biomed reports. 2017;7:115–22.View ArticleGoogle Scholar
- Mannik J, Vaas P, Rull K, Teesalu P, Laan M. Differential placental expression profile of human growth hormone/chorionic somatomammotropin genes in pregnancies with pre-eclampsia and gestational diabetes mellitus. Mol Cell Endocrinol. 2012;355:180–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Diplas AI, Lambertini L, Lee M-J, Sperling R, Lee YL, Wetmur J, et al. Differential expression of imprinted genes in normal and IUGR human placentas. Epigenetics. 2009;4:235–40.View ArticlePubMedGoogle Scholar
- Iglesias-Platas I, Martin-Trujillo A, Petazzi P, Guillaumet-Adkins A, Esteller M, Monk D. Altered expression of the imprinted transcription factor PLAGL1 deregulates a network of genes in the human IUGR placenta. Hum Mol Genet. 2014;23:6275–85.View ArticlePubMedPubMed CentralGoogle Scholar
- Camprubi C, Iglesias-Platas I, Martin-Trujillo A, Salvador-Alarcon C, Rodriguez MA, Barredo DR, et al. Stability of genomic imprinting and gestational-age dynamic methylation in complicated pregnancies conceived following assisted reproductive technologies. Biol Reprod. 2013;89:50.View ArticlePubMedGoogle Scholar
- Maltepe E, Krampitz GW, Okazaki KM, Red-Horse K, Mak W, Simon MC, et al. Hypoxia-inducible factor-dependent histone deacetylase activity determines stem cell fate in the placenta. Development. 2005;132:3393–403.View ArticlePubMedGoogle Scholar
- Handy DE, Castro R, Loscalzo J. Epigenetic modifications: basic mechanisms and role in cardiovascular disease. Circulation. 2011;123:2145–56.View ArticlePubMedPubMed CentralGoogle Scholar
- Miguel-Escalada I, Pasquali L, Ferrer J. Transcriptional enhancers: functional insights and role in human disease. Curr Opin Genet Dev. 2015;33:71–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Turner BM. Defining an epigenetic code. Nat Cell Biol. 2007;9:2–6.View ArticlePubMedGoogle Scholar
- Lee KJ, Shim SH, Kang KM, Kang JH, Park DY, Kim SH, et al. Global gene expression changes induced in the human placenta during labor. Placenta. 2010;31:698–704.View ArticlePubMedGoogle Scholar
- Tuteja G, Chung T, Bejerano G. Changes in the enhancer landscape during early placental development uncover a trophoblast invasion gene-enhancer network. Placenta. 2016;37:45–55.View ArticlePubMedGoogle Scholar
- Uusküla L, Männik J, Rull K, Minajeva A, Kõks S, Vaas P, et al. Mid-gestational gene expression profile in placenta and link to pregnancy complications. PLoS One. 2012;7:e49248.View ArticlePubMedPubMed CentralGoogle Scholar
- Karlic R, Chung H-R, Lasserre J, Vlahovicek K, Vingron M. Histone modification levels are predictive for gene expression. Proc Natl Acad Sci U S A. 2010;107:2926–31.View ArticlePubMedPubMed CentralGoogle Scholar
- Bonet B, Brunzell JD, Gown AM, Knopp RH. Metabolism of very-low-density lipoprotein triglyceride by human placental cells: the role of lipoprotein lipase. Metabolism. 1992;41:596–603.View ArticlePubMedGoogle Scholar
- Gormley M, Ona K, Kapidzic M, Garrido-Gomez T, Zdravkovic T, Fisher SJ. Preeclampsia: novel insights from global RNA profiling of trophoblast subpopulations. Am J Obstet Gynecol. 2017;217(2):200.e1–200.e17.View ArticleGoogle Scholar
- Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25:1754–60.View ArticlePubMedPubMed CentralGoogle Scholar
- Ji H, Jiang H, Ma W, Johnson DS, Myers RM, Wong WH. An integrated software system for analyzing ChIP-chip and ChIP-seq data. Nat Biotechnol. 2008;26:1293–300.View ArticlePubMedPubMed CentralGoogle Scholar
- Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.View ArticlePubMedPubMed CentralGoogle Scholar
- McLean CY, Bristor D, Hiller M, Clarke SL, Schaar BT, Lowe CB, et al. GREAT improves functional interpretation of cis-regulatory regions. Nat Biotechnol. 2010;28:495–501.View ArticlePubMedPubMed CentralGoogle Scholar