Next Article in Journal
Overdose of Quetiapine—A Case Report with QT Prolongation
Next Article in Special Issue
Prenatal Exposure to Locally Emitted Air Pollutants Is Associated with Birth Weight: An Administrative Cohort Study from Southern Sweden
Previous Article in Journal
Exposure to Sub-Lethal Doses of Permethrin Is Associated with Neurotoxicity: Changes in Bioenergetics, Redox Markers, Neuroinflammation and Morphology
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Toxics
Volume 9
Issue 12
10.3390/toxics9120338
toxics-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

Early Pregnancy Exposure to Ambient Air Pollution among Late-Onset Preeclamptic Cases Is Associated with Placental DNA Hypomethylation of Specific Genes and Slower Placental Maturation

by
Karin Engström
1,*,†,
Yumjirmaa Mandakh
2,†,
Lana Garmire
3,
Zahra Masoumi
4,
Christina Isaxon
5,
Ebba Malmqvist
2,
Lena Erlandsson
4,‡ and
Stefan R. Hansson
4,6,‡
1
EPI@LUND, Division of Occupational and Environmental Medicine, Department of Laboratory Medicine, Lund University, SE-22362 Lund, Sweden
2
Environment Society and Health, Division of Occupational and Environmental Medicine, Department of Laboratory Medicine, Lund University, SE-22362 Lund, Sweden
3
Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, MI 48104, USA
4
Division of Obstetrics and Gynecology, Department of Clinical Sciences Lund, Lund University, SE-22185 Lund, Sweden
5
Ergonomics and Aerosol Technology, Department of Design Sciences, Lund University, SE-22362 Lund, Sweden
6
Obstetrics and Gynecology, Skåne University Hospital, SE-22184 Lund, Sweden
*
Author to whom correspondence should be addressed.
†
These authors have contributed equally to this work and share first authorship.
‡
These authors have contributed equally to this work and share last authorship.
Toxics 2021, 9(12), 338; https://doi.org/10.3390/toxics9120338
Submission received: 12 November 2021 / Revised: 25 November 2021 / Accepted: 29 November 2021 / Published: 6 December 2021
(This article belongs to the Special Issue Update Reproductive Toxicology Research Associated with Air Pollutions)
Browse Figures
Versions Notes

Abstract

:
Exposure to ambient air pollution during pregnancy has been associated with an increased risk of preeclampsia (PE). Some suggested mechanisms behind this association are changes in placental DNA methylation and gene expression. The objective of this study was to identify how early pregnancy exposure to ambient nitrogen oxides (NOx) among PE cases and normotensive controls influence DNA methylation (EPIC array) and gene expression (RNA-seq). The study included placentas from 111 women (29 PE cases/82 controls) in Scania, Sweden. First-trimester NOx exposure was assessed at the participants’ residence using a dispersion model and categorized via median split into high or low NOx. Placental gestational epigenetic age was derived from the DNA methylation data. We identified six differentially methylated positions (DMPs, q < 0.05) comparing controls with low NOx vs. cases with high NOx and 14 DMPs comparing cases and controls with high NOx. Placentas with female fetuses showed more DMPs (N = 309) than male-derived placentas (N = 1). Placentas from PE cases with high NOx demonstrated gestational age deceleration compared to controls with low NOx (p = 0.034). No differentially expressed genes (DEGs, q < 0.05) were found. In conclusion, early pregnancy exposure to NOx affected placental DNA methylation in PE, resulting in placental immaturity and showing sexual dimorphism.

1. Introduction

Preeclampsia (PE) affects, in average, 0.2%–9% of pregnancies worldwide [1,2] and has been associated with increased risk of maternal mortality [3,4] and maternal cardiovascular disease (CVD) in later life. In addition, offspring born to PE mothers are at high risk of low birth weight at birth and CVD later in life [5,6]. Preeclampsia presents with new onset hypertension after 20 weeks of gestation accompanied by proteinuria and/or dysfunction of one or more maternal organs [7]. The pathogenesis of PE is not fully understood, but placental dysfunction due to syncytiotrophoblast stress explains the underlying mechanism of the clinical presentation of PE [8]. Preeclampsia is commonly divided into early-onset PE (manifesting before week 34) and late-onset PE (at or after 34 weeks) [9]. Early-onset PE is characterized by poor placentation and higher maternal and fetal morbidity and mortality. Late-onset PE makes up 80% of the PE cases and is associated with maternal risk factors such as metabolic syndrome, obesity, and impaired glucose tolerance. Studies have implicated that a female fetus is more common in early-onset PE [10] and a male fetus in late-onset PE [11].
PE is a “multifactorial syndrome” in which ambient air pollution may play a role [12,13,14,15,16]. Ambient particulate air pollution is the fourth leading risk factor for attributable deaths and the seventh leading risk factor for attributable burden of disease [17]. Systematic reviews and meta-analyses suggest that exposure to ambient air pollution during pregnancy is associated with increased risk of PE [13,18], especially for the first trimester [18]. However, the underlying mechanism of the association between ambient air pollution during pregnancy and increased risk of PE is not clearly elucidated.
Placental DNA methylation may be one of the underlying mechanisms linking PE to ambient air pollution. DNA methylation occurs predominantly on cytosines (C) that are followed by guanine (G) residues, referred to as CpG sites, that, along with other epigenetic modifications, regulate gene expression [19]. DNA methylation plays a critical role in cell differentiation and epigenetic reprogramming during embryonic development in the first trimester [20]. Environmental exposures occurring at this critical period have a potential to induce additional epigenetic changes and alter gene expression, resulting in increased susceptibility to disease [21,22]. In previous studies, altered DNA methylation in PE placentas has been seen, which was presumed to be associated with placental development and function [23]. In addition, a differential DNA methylation placental profile has been found in relation to exposure to ambient air pollution during pregnancy [24].
Nitrogen oxide (NOx) is a generic term for gaseous oxides containing nitrogen including nitric oxide (NO) and nitrogen dioxide (NO2) gas. The main source of these gases is internal combustion engines that burn fuels such as oil, diesel, gas, and organic matter at high temperatures. The health effects of NO2 have recently been reviewed by WHO with strong evidence of mortality effects at levels above 10 µg/m3 [25]. Prenatal exposure to NOx during the first and second trimesters has been associated with hypomethylation of two CpGs located at the Adenosine receptor A2b (ADORA2B) gene, suggested to be involved in PE pathogenesis [26]. However, no studies have yet evaluated the association between exposure to NOx during pregnancy and DNA methylation among PE cases and normotensive pregnancies.
The placenta grows rapidly throughout gestation and matures within the limited time in a physiologic phenomenon known as placental ageing [27]. Timely ageing of the placenta is needed for optimal fetal growth and development. The biological age of a tissue can be estimated with high accuracy using epigenetic data from whole-genome DNA methylation arrays [28], and there are specific epigenetic clocks for the placenta that calculates gestational aging [29].
The present study aimed to identify how early pregnancy exposure to ambient nitrogen oxides (NOx) among PE cases and normotensive controls influence DNA methylation and gene expression.

2. Materials and Methods

2.1. Study Population

This study included placentas from pregnant women enrolled in a biobank specifically created for the research on PE at the Department of Obstetrics and Gynecology, Lund University, Sweden. The biobank consists of a biomedical database with biological research samples and detailed clinical data linked to Swedish personal identification numbers for each participant. The placental biopsies used in this study were collected from both Caesarean and vaginal deliveries in 2008–2015 at Skåne University Hospital in Lund and Malmö, Sweden (Figure 1). After excluding those who had multiple pregnancies, smoked during pregnancy, and lived outside the territory of our dispersion model, a total of 111 women were included in the study. Of these, 29 were preeclamptic (one was early-onset and 28 late-onset) and 82 were normotensive. Complete data on NOx exposure during the first trimester, PE status, and potential confounding variables are presented in Table 1. In the descriptive statistics (Table 1), the χ2 test for categorical variables and t-test for continuous variables were used to assess the differences between PE cases and normotensive controls. A p-value of less than 0.05 was considered statistically significant. The study protocol was approved by Lund University Ethical Review Board, LU803-2 (2015/14), Lund, Sweden. Written informed consent was obtained from all participants.

2.2. Air Pollution Exposure Assessment

After selecting study participants with available placental samples in the biobank, we obtained the geocoded residential addresses in Scania, Sweden from the regional population registry using the personal identification numbers. Hourly concentrations of NOx during pregnancy were modeled at geocoded residence level using a modified Gaussian-plume dispersion model with local emission data in EnviMan software [30]. We were not able to model NO and NO2 separately, because only NOx emission data was available at the time of exposure assessment. Detailed exposure assessment using the dispersion model has been described previously [31]. The model provides hourly concentrations on 100 × 100 m grids. To account for long-range transboundary air pollution, a background concentration of 2.5 µg/m3 was added to the modeled NOx concentrations. The mean exposure to NOx during the first trimester for each person was calculated as the average of the modeled hourly NOx concentrations. Mean NOx exposure was then used as a dichotomized variable via the median cut-off (12.1 μg/m3; Interquartile Range: 7.9, 19.1) in high and low exposure groups.

2.3. Study Design

The study participants were divided into groups based on a combination of pregnancy status (PE or normotensive control) and exposure to ambient NOx (low or high NOx exposure based on median split). The study participants were thus divided into four categories: controls with low NOx exposure (N = 45), controls with high NOx exposure (N = 37), PE cases with low NOx exposure (N = 10), and PE cases with high NOx exposure (N = 19). DNA methylation was compared between these groups.
In the gene expression analyses, we only compared PE cases to controls, since the number of individuals (N = 17) included were too low for more in-depth evaluation of differential gene expression.

2.4. Association between NOx Exposure and PE

Associations between NOx exposure during the first trimester (low or high exposure) and PE status was estimated with crude and adjusted logistic regression. The final model was adjusted for clinical data, such as maternal age, pregestational body mass index (BMI), gestational age (linear), and fetal sex [31,32,33,34].

2.5. Biological Sample Collection

The placentas were collected within 4 hours after delivery and stored at 4 °C until biopsies were taken. Briefly, a placental biopsy of 1 × 1 × 1 cm size was excised 7 cm from the umbilical cord insertion or from a central part of the villi region if cord insertion was asymmetrical. Vertically, the biopsy was taken centrally at the thickest point to not include chorionic plate or maternal tissue. Due to the heterogeneity of placental tissue, we also adjusted for the estimated proportion of placental cell types in the statistical analyses (described in Section 2.10). After rinsing in phosphate-buffered saline solution, the placental biopsies were patted dry with a napkin and stored on dry ice for 30 min prior to transfer to −80 °C until further analysis.

2.6. DNA and RNA Extraction

Total DNA and RNA from placental biopsies were extracted using Qiagen AllPrep DNA/RNA/Protein Mini Kit (Qiagen, Hilden, Germany). The DNA and RNA quality (A260/280 and A260/230) were evaluated using a NanoDrop Spectrophotometer ND-1000 (NanoDrop technologies, Wilmington, DE, USA). The RNA integrity was evaluated with LabChip (LabChip GX/GXII systems, PerkinElmer, MA, USA) resulting in 17 samples with an RNA integrity number (RIN) above 7. Only samples with RIN >7 were included in the RNA-seq analyses.

2.7. Placental DNA Methylation Laboratory Analyses

A total amount of 500 ng DNA was bisulfite-treated using the EZ DNA Methylation kit v 1.1 (Zymo Research, Irvine, MA, USA). DNA samples were randomized for distribution in two 96-well analysis plates prior to analysis. Genome-wide DNA methylation was determined at the Center for Translational Genomics (CTG), Lund University, Sweden, using the Infinium MethylationEPIC BeadChip (lllumina, San Diego, CA, USA), analyzing approximately 850,000 specific markers of DNA methylation. All beadchips were from the same batch.

2.8. Placental DNA Methylation Preprocessing

The statistical software R, version 4.0.208 (R Project for Statistical Computing), was used for thte preprocessing and analysis of the data. The R packages ChAMP [35,36] and minfi [37,38] were used for image processing, quality control, filtering, and normalization. Normalization was performed using the beta-mixture quantile (BMIQ) normalization [35,39]. Non-CpG probes and multi-hit probes were filtered. The detection p-value was below 0.01 for at least >98% of CpGs for all samples. The CpGs with detection p-value above 0.01 in more than 20% of the samples were removed.
Analyses of differentially methylated positions (DMPs) were first performed including all samples irrespective of fetal sex. Group comparisons that had significant DMPs were then analysed and stratified for fetal sex. For the analysis in which all individual samples were included, probes in the sex chromosomes were removed. Overall, 805185 CpGs were retained in the analyses including all individual samples, and 823536 CpGs were analyzed in the sex-stratified analyses.
The placental DNA methylation at specific CpG sites was used to estimate the proportion of different placental cell types (Trophoblasts, Stromal, Hofbauer, Endothelial, nucleated RBC, Syncytiotrophoblast) using the R package Planet [40]. Here, reference data from the first trimester were used to infer cell composition using the constrained projection approach implemented in minfi. These estimated cell type proportions were later evaluated as potential covariates in the statistical models.

2.9. Predicted Gestational Age Acceleration/Deceleration and Associations with PE and NOx Exposure

The predicted gestational age (GA) was estimated using GA clocks for placental DNA methylation data [29] using the package Planet [40]. By regressing the predicted GA to the chronological GA, the residual gives each placenta a value for GA acceleration/deceleration [41]. A positive value indicates that the placenta ages faster than expected, i.e., an accelerating placental epigenetic clock or GA acceleration. A negative value indicates that the placenta ages slower than expected, which is GA deceleration. The so-called “control placental clock” (CPC) was used [40].
The correlation between predicted GA and chronological GA was evaluated using Spearman correlation for all individuals as well as stratified for PE status. The differences in GA acceleration/deceleration between PE cases and controls, as well as combined PE status and NOx exposure groups (categorical variable), were analysed using linear regression.

2.10. DNA Methylation Preprocessing and Data Analysis

Singular value decomposition (SVD) was performed in ChAMP to identify the technical and biological variables associated with DNA methylation (data not shown). These variables were then included in ComBat adjustments [42,43] using ChAMP (for technical variables) or considered as covariates in linear models. The SVD analyses showed that the technical variables “sentrix position” and “analysis plate” accounted for a substantial fraction of variation (p < 0.05) in principal components (PCs) one and two, respectively. ComBat adjustments were performed to remove the technical variation. The SVD plots also showed that all estimated cell type proportions had a p < 0.05 in the first and/or second PC.
We used estimated fractions of syncytiotrophoblast and Hofbauer cells in the final model, considering that these two cell types form a majority of the placental cell population. Villous syncytiotrophoblasts make up the multinucleated syncytia and are the major cell type in the placenta [44]. Hofbauer cells are located in the placental villous core and are suggested to be involved in immune tolerance during pregnancy [45].
Pregestational BMI, gestational age, fetal sex, and DNA concentration also showed a p < 0.05 in the first or second PCs and were thus included in the model. For the sex-stratified analyses, we performed a model with fewer covariates due to a lower number of individual samples included. Sex-specific models were adjusted for estimated fractions of syncytiotrophoblast, Hofbauer cells, DNA concentration, and gestational age.
The DMPs were evaluated by fitting a linear regression model to each CpG using the R package limma [46], with adjustments as described above. M-values (the log2 ratio of the intensities of methylated probe versus unmethylated probe) were employed. Pair-wise comparisons between the combined PE status and NOx exposure groups were performed using a contrast matrix (for a summary of the different comparisons made, see Supplemental Material, Figure S1). Empirical Bayes smoothing was applied to the standard errors with a robust selection against outlier sample variances. Further, p-values were adjusted for multiple comparisons for all CpGs by the Benjamini–Hochberg false discovery rate (FDR) method to obtain q-values. A q-value of 0.05 or lower was considered statistically significant. Analyses were performed including samples from: (1) both fetal sexes in the same analysis and (2) stratified for fetal sex for group comparisons that had significant DMPs.
Imprinted genes are expressed in a parent-of-origin-specific manner that was established in the parental germline and is caused by epigenetic processes involving DNA methylation and histone methylation. Many imprinted genes are involved in embryogenesis and placental development. All imprinted genes were removed from the results, and we used the following dataset to confirm imprinted genes (https://www.geneimprint.com/site/genes-by-species, accessed on 5 November 2021).

2.11. In Silico Analyses: Interaction Network Analyses and Gene Ontology Analyses for DNA Methylation

We used ConsensusPathDB (CPDB) (http://cpdb.molgen.mpg.de/, accessed on 5 November 2021) to search for potential gene regulatory interactions (induced network molecule analysis). In the analyses, we included genes that had CpGs with q < 0.1 in the DMP analyses. When comparing cases and controls with high NOx exposure for female placentas, we had to use a cut-off value of q value < 0.05 due to the high number of DMPs.
Gene ontology (GO) analyses were run for genes that had CpGs with q < 0.1 in the DMP analyses. This was performed using the gometh (DMPs) function in the R package missMethyl [47]. A q-value of 0.05 or lower for the term in the GO analyses was considered statistically significant.

2.12. RNA Sequencing Analysis

Total RNA from placental biopsies, prepared as described above, was used to prepare libraries using the TruSeq® Stranded mRNA Library Prep (20020594, Illumina, San Diego, CA, USA). Library preparation quality control (QC) was performed using LabChip DNA High Sensitivity Reagent kit (CLS138948, Perkin Elmer, Waltham, MA, USA) and DNA 1K/12K/Hi Sensitivity Assay LabChip (760517, Perkin Elmer). The library pool(s) were quantified using the QuantIT® dsDNA HS Assay Kit (Q33120, Thermo Fischer, Waltham, MA, USA). The library pool was then sequenced as paired-end, 75-bp reads on a NovaSeq 6000 (Illumina).
Raw data were demultiplexed, and FASTQ files for each sample were generated using the bcl2fastq software (Illumina, San Diego, CA, USA). The FASTQ data were checked using the FastQC tool [48]. Alignment of reads were made to the reference genome sequence from the Ensembl database, the Human GRCh38. Assembly of the alignments into full transcripts and quantification of the expression levels of each gene and transcript were performed using the StringTie software [49]. There were 57,492 input genes, of which 1243 had zero expression, and 37,725 (65%) showed low expression and were removed using the filterByExpr function in DSEq2 [50]. Thus, 18,524 genes were included in the final analysis. Counts were normalized to Trimmed Mean of M-values (TMM) using limma. Differentially expressed genes (DEG) analysis was performed using Limma voom [46]. Since there were few individuals (N = 17), we only compared PE cases (N = 10) and controls (N = 7). Models were adjusted for fetal sex. A q-value of 0.05 or lower was considered statistically significant. Since no genes reached q < 0.1, we did not perform any GO analyses.

3. Results

3.1. Characteristics of the Study Participants

First, we summarized the descriptive data of the participants and evaluated differences in characteristics between PE cases and controls using the χ2 test for categorical variables and a t-test for continuous variables (Table 1). There was a statistically significant lower mean gestational age among PE cases compared to controls (38.5 vs. 39.5 weeks, p = 0.014), and PE cases were also more likely to have had PE previously (p = 0.046).

3.2. Association between NOx Exposure and PE Status

Associations between NOx exposure during the first trimester (low or high exposure, low exposure is reference) and PE status were estimated with crude and adjusted logistic regression. The crude odds ratio (OR) for the association between exposure to ambient NOx during the first trimester and PE was 2.3 (95% confidence interval (CI): 0.96, 5.6, p = 0.062). After controlling for maternal age, pregestational BMI, gestational age (linear), and fetal sex, the OR was 3.1 (95% CI: 1.1, 8.6), p = 0.028.

3.3. Association between Exposure to Ambient NOx and PE Status with Gestational Age Acceleration/Deceleration

We investigated whether the combined PE status and NOx exposure group was associated with GA acceleration/deceleration. By regressing the predicted GA to the chronological GA, the residual gives each placenta a value for GA acceleration/deceleration [41]. A positive value indicates that the placenta ages faster than expected, i.e., an accelerating placental epigenetic clock or GA acceleration, while a negative value indicates a GA deceleration.
Spearman correlations (Rs) between predicted GA and chronological GA showed somewhat stronger correlations for PE cases (Rs = 0.64, p < 0.001) than for controls (Rs = 0.56, p < 0.001).
The PE cases demonstrated GA deceleration compared to controls, although this was not statistically significant (beta = −0.43, p = 0.057, controls are reference) (Figure 2A). Placentas from PE cases with high NOx exposure showed significant GA deceleration compared to placentas from controls with low NOx exposure (beta = −0.61, p = 0.035, model R2 = 0.03) (Figure 2B). However, the variable of combined PE status and NOx group (as a categorical variable) was not statistically significant in the model (p = 0.20, model R2 = 0.04) (Supplemental Material, Table S1).

3.4. Differentially Methylated Positions Associated with PE Status and NOx Exposure

We investigated whether the combined PE status and NOx exposure group was associated with differential methylation of specific positions in the genome.
In the analyses including placentas of both fetal sexes, we identified six DMPs in the comparison between controls with low NOx exposure vs. PE cases with high NOx exposure, one DMP in the comparison between PE cases with low NOx exposure vs. controls with high NOx exposure and fourteen DMPs in the comparison between cases vs. controls, both with high NOx exposure (q < 0.05, Table 2), making up a total of 19 CpGs showing differential methylation.
Two CpGs were statistically significant in both group comparisons: cg15999356 in the Yes1-associated transcriptional regulator gene (YAP1) and cg26672098 (located at chromosome 2p21, not annotated to any gene) showed lower methylation in the PE cases with high NOx exposure group compared to the reference groups in both comparisons. All statistically significant CpG sites displayed lower methylation in the PE cases with high NOx exposure group compared to the reference groups, except for one CpG site, which instead showed a higher methylation: cg02404739 in FAM111 trypsin-like peptidase B (FAM111B). To note is that two of the top DMPs, cg17283620 in 3-hydroxyanthranilate 3,4-dioxygenase (HAAO), and cg02171814 in serpin family F member 2 (SERPINF2), had nearby single-nucleotide polymorphisms (SNPs) with a frequency above 10% in the CEU reference population (Utah Residents with Northern and Western European Ancestry) [51,52]. These SNPs were rs12617051 (HAAO, minor allele frequency [MAF] 43%) and rs77305322 (SERPINF2, MAF 27%). HAAO rs12617051 was situated in the CpG site, while the SERPINF2 rs77305322 was situated 35 base pairs away from the CpG site.
When stratifying for fetal sex, one DMP was seen for placentas with male fetuses in the comparison between PE cases with low NOx exposure vs. controls with high NOx exposure. This DMP was situated in the signal transducer and activator of transcription 3 (STAT3). For placentas with female fetuses, we identified 309 DMPs (q < 0.05) when comparing cases and controls with high NOx exposure (Table 2). A total of 280 out of 309 (90%) statistically significant CpG sites (q < 0.05) displayed lower methylation in the PE cases with high NOx exposure group compared to controls with low NOx exposure. The top 50 statistically significant DMPs in the comparison between PE cases with low NOx exposure vs. controls with high NOx exposure for placentas with female fetuses are shown in Table 2 (top 13, based on q-values) and the Supplemental Material, Table S2 (ranked 14-50, based on q-values). Notably, cg09234983 in large tumor suppressor kinase 2 (LATS2) had a nearby SNP with a MAF above 10% in the CEU population. This SNP was rs9506597 and had a MAF of 28%. The distance between rs9506597 and the query base of the probe was 27 base pairs. There was also one DMP (cg16162930, not annotated to any gene) in the comparison between PE cases with low NOx exposure vs. controls with high NOx exposure for placentas with female fetuses. However, this CpG had a nearby SNP (rs2105225) with a MAF of 44% in the CEU population, situated at the CpG site.

3.5. In Silico Analyses: Interaction Network Analyses and Gene Ontology Analyses for DNA Methylation

The results from the analysis of gene regulatory interactions pointed to specific transcription factors (TFs). A central TF when comparing controls with low NOx exposure vs. PE cases with high NOx exposure was hepatocyte nuclear factor 4 alpha (HNF4A) (Supplemental Material, Figure S2A). When comparing cases and controls, both with high NOx exposure, E2F transcription factor 4 (E2F4) and TATA-box binding protein associated factor 1 (TAF1) were central TFs (Supplemental Material, Figure S2B). Lastly, TAF1 was a central TF when comparing low-exposed controls vs high-exposed PE in placentas from female fetuses (Supplemental Material, Figure S2C).
We performed GO analyses, including CpGs with q < 0.1 in the DMP analyses. The number of CpGs annotated to genes and thus included in the GO analyses were 123 in the comparison between controls with low NOx exposure vs. PE cases with high NOx exposure, and 264 in the comparison between cases and controls, both with high NOx exposure. In the comparison between PE cases with low NOx exposure vs. controls with high NOx exposure, there was just one DMP with q < 0.1, and, thus, we did not perform any GO analysis. We found no significant hits in the GO analyses after the correction for multiple testing in any of these two group comparisons (FDR < 0.05) (Supplemental Material, Table S3).
For the sex-stratified analyses, for placentas from female fetuses, there were no DMPs with q < 0.1 in the DMP analyses comparing controls with high NOx exposure vs. PE cases with high exposure, while the number of genes included in the GO analyses in the comparison between cases and controls with high NOx exposure were 1465. For placentas from male fetuses, there was just one DMP with q < 0.1 in the comparison between PE cases with low NOx exposure vs. controls with high NOx exposure, and there were no DMPs with q < 0.1 in any of the other two group comparisons. Thus, no GO analyses were performed for placentas from male fetuses.
When stratifying for fetal sex, for placentas from female fetuses, one enriched pathway in the GO analyses, structural constituents of postsynapse, had an FDR < 0.1 when comparing cases and controls with high NOx exposure (Supplemental Material, Table S3).

3.6. RNA Sequencing: Analysis of Differentially Expressed Genes and Gene Ontology

We investigated whether PE status was associated with differential gene expression.
Only samples with RIN > 7 were included in the RNA sequencing; thus, the sample size for the RNA-seq was only 17 individuals, and, due to the small sample size, no gender-stratified analyses were performed. Descriptive data for study participants included in RNA-seq are shown in Supplemental Material, Table S4. No genes reached statistical significance in the DEG analyses (q < 0.05). Top genes in the DEG analyses, ranked on unadjusted p-values, are shown in Table 3. None of the genes with DMPs (for both sexes together) showed a p < 0.05 in the DEG analyses.

4. Discussion

Our study has identified several differentially methylated placental genes that were associated with PE in relation to exposure to ambient NOx during the first trimester. We also found sex differences, where placentas with female fetuses showed more differentially methylated placental genes than placentas with male fetuses. Placentas from PE cases with high NOx exposure showed a significant gestational age deceleration compared to placentas from normotensive controls with low NOx exposure, suggesting a state of placental villous immaturity [53].
The DNA methylation analysis indicated DNA hypomethylation in locations associated with several genes, of which many have been investigated in tumor growth and metastasis. Their specific role in placentation may remain unclear, but reactivation of embryonic developmental processes such as those involved in placental growth and invasion have been suggested to explain the molecular mechanisms in tumor growth and metastasis, depicting many similarities between the two phenomena [54,55,56]. Accordingly, it is interesting to find hypomethylation of two members of the SH3 domain protein family (SH3D21 and YAP1) as well as RPRD1B in placentas from PE cases with high NOx exposure compared to controls with low NOx exposure. These three genes have been shown to play significant roles in regulation of cell proliferation, migration, and even differentiation in cancer. In trophoblasts, YAP1 has been shown to regulate villous cytotrophoblasts proliferation and stemness, thus suppressing trophoblast fusion (to create syncytium) and even limiting extra villous trophoblast differentiation [57]. Inhibition or degradation of YAP1, on the other hand, can lead to trophoblast apoptosis [58] or decreased trophoblast invasion [59], indicating the significance of balance in YAP1 activity in different trophoblast cell types for effective placentation. Recent studies found significantly lower levels of YAP1 mRNA expression [60] and protein levels [61] in severe PE compared to uncomplicated pregnancies. Taken together, hypomethylation of YAP1 may play a role in trophoblast fusion and proliferation in PE, thereby potentially also contributing to defective villous maturation and impaired placental function [62].
Another gene correlated with a hypomethylated position was CACNA2D3, which regulates voltage-dependent calcium channels. By introducing mitochondria-mediated apoptosis, CACNA2D3 also plays a critical role in tumor suppression [63]. While calcium signaling and transfer in the placenta are complicated and significant for fetal development [64], altered expression of L-type voltage-dependent calcium channels in pregnancy disorders such as PE has been associated with depolarization in arterial smooth muscles and hypoxic vasoconstriction in the placenta [64,65]. Furthermore, intracellular alterations in calcium homeostasis have been shown in PE, where reports show reduced levels of the calcium-binding messenger calmodulin [66] and the downstream second messenger calcium/calmodulin-dependent kinase IV (CAMK4) [67]. Reduced levels of CAMK4 suppresses the proliferation and migratory capacity of trophoblasts, along with increased apoptosis [67]. Lastly, the hypomethylated DMP cg26672098 is not linked to any gene but located on chromosome 2p21. It was significant in the comparison between PE cases and controls, both with high NOx exposure. This CpG has previously been found to be hypomethylated in association with early-onset PE [68].
When comparing PE cases and controls, both with high NOx exposure, many of the genes were hypomethylated and previously associated with changes in actin cytoskeleton and cell proliferation, as well as migration and metastasis in cancer. For instance, studying PEBP1 in cancer showed that it controls the cell cytoskeleton by interacting with actin, as well as regulating various signaling pathways and GTPases upstream of actin cytoskeleton regulation [69]. These changes in the cytoskeleton can be associated with membrane changes in processes from cell migration to vesicle formation and endocytosis [69]. In the placenta, PEBP1 has been suggested to regulate trophoblast migration and is expressed in villous cytotrophoblasts while its expression is translocated to the syncytium in PE placentas [70]. Due to syncytiotrophoblast stress in PE, there is an increased production and shedding of vesicles compared to normal pregnancies [71]. A key depolymerizer of actin filaments also found in a hypomethylated region in PE is MICAL1, which regulates actin cytoskeleton dynamics in a variety of processes from cytokinesis [72] to invasion and mobility in cancer [73]. Additionally, downregulation of MICAL1 in PE placentas has been observed along with increased trophoblast proliferation [74]. In this comparison, we also saw YAP1, which was discussed earlier, and MB21D2, whose overexpression facilitates cell proliferation and invasion in cancer, while its suppression leads to apoptosis [75]. Another gene in a hypomethylated position is CD163, which is mostly known for being expressed on macrophages and endocytoses hemoglobin/haptoglobin complexes [76]. In the placenta, CD163 is expressed on the anti-inflammatory placental Hofbauer cells [77], and reduced numbers of these cells as well as down-regulation of CD163 expression have previously been reported in PE [78,79]. This would result in a reduction of macrophage functions, such as promoting anti-inflammatory and pro-angiogenic responses in the placenta.
Our analysis also indicated three specific hypomethylated positions associated with HAAO, SUPT3H, and SERPINF2, which were not directly relevant to the above-mentioned mechanisms but still intriguing. For instance, in pregnancy, tryptophan catabolism is required to regulate the maternal immune reaction towards an implanting fetus [80]. In PE, tryptophan catabolism is deregulated, and HAAO mRNA levels are decreased [81,82], which may play some role in maternal–fetal tolerance. Tryptophan levels have also been shown to be decreased in PE, along with kynurenine pathway enzymes that break down tryptophan [82]. Tryptophan is also a precursor for the neurotransmitter serotonin, which is a critical player in fetal brain development [83]. Another gene, SUPT3H, considered to be important in the regulation of transcription [84], can, if mutated, cause placental defect and fetal mortality in mice [85]. The other gene with a unique role was SERPINF2, which is a plasmin inhibitor [86], but it also has a role in vascular remodeling [87]. Interestingly, it has been reported to be increased in maternal plasma in PE and co-localized with fibrin deposits in the PE placenta [88]. The last one, PLCXD2, is a phospholipase located in the nucleus with no known connection to the placenta or pregnancy, but a novel intronic SNP has recently been associated with risk for several cancers [89]. The only hypermethylated position in this comparison was associated with FAM111B, a nuclear serine protease that can promote proliferation, migration, and invasion in lung cancer, while its suppression induces cell apoptosis in cancer [90,91].
There were also a few DMPs in the comparison between PE cases with low NOx exposure and controls with high NOx exposure. Since this comparison looks at cases and controls with different exposure levels, there is no rationale or relevant scientific question for any further evaluation. However, the gene found was sodium channel and clathrin linker 1 (SCLT1), which is an adaptor protein that is involved in linking clathrin to sodium channels and therefore necessary for receptor endocytosis [92]. It has not previously been associated to PE or the placenta, but placenta-derived extracellular vesicles have been shown to be taken up by endothelial cells through clathrin-dependent endocytosis [93].
Overall, DNA methylation analysis of PE placentas in relation to exposure to high ambient NOx indicate a potential deregulation in cell proliferation, differentiation, and migration associated with altered cytoskeletal dynamics and cell signaling. These changes may underlie the lower invasive capacity of the trophoblasts seen during implantation at Stage 1 of PE development, resulting in impaired placental perfusion [94]. Later, these changes can contribute to defective morphology and altered metabolism and immune response that are hallmarks of Stage 2 in PE [8]. In addition, our DNA methylation analysis indicated DNA hypomethylation, and this is in line with a recent systematic review demonstrating an association between exposure to ambient air pollution during pregnancy and global loss of methylation in the umbilical cord blood and placenta [95]. In addition, this report indicated that early pregnancy is the most vulnerable period when it comes to exposure.
For the female placentas comparing PE cases and controls, both with high NOx exposure, our analysis mainly indicated hypomethylated positions. We identified several genes about which very little is known in the literature regarding their function and expression, and many of them had no known connection to pregnancy, the placenta, or PE. Information regarding expression in different tissues was instead available to us at the online databases The Human Protein Atlas (www.proteinatlas.org [96]) and the GeneCards (www.genecards.org, accessed on 5 November 2021 [97]). Several of the genes were connected to cell and/or vesicle migration and to the cytoskeleton. Again, making a comparison to tumor cells, the process of metastasis and tumor mobility has similarities to placental growth and invasion. One of the key elements in tumor mobility is the cytoskeleton that undergoes active change during this process. This results in changes in plasma membrane and in cellular balance between adhesion and mobility. Members of the tetraspanin family have been shown to play important roles in cytoskeleton dynamics during tumor metastasis and in the regulation of tumor microenvironment [98]. Interestingly, we found hypomethylation of TSPAN4, a member of the tetraspanin family, that has recently been shown to play a vital role in the formation of migrasomes, which are cellular organelles that form in migrating cells [99]. Additionally, TSPAN4 has been shown to be expressed in placenta in different types of trophoblasts, as well as Hofbauer cells [96]. Hofbauer cells are of placental origin and have recently been shown to have microbicidal capacity as well as a possible role in placental angiogenesis and remodeling [100], which are functions that would require mobility. Two other hypomethylated genes, ACTR3C and KRT24, are both involved with the cytoskeleton. Protein ARP11, produced by ACTR3C, is an intracellular protein expressed in the placenta [96]. It is involved in regulation of the actin cytoskeleton [101] and in early formation of endosomes by being part of the dynein–dynactin motor complex [102]. Also involved with motor complexes is TRAK1, a mitochondrial adaptor protein that is essential for mitochondrial trafficking. When mitochondria move along microtubules, they are attached to TRAK proteins that link them to the two motor complexes kinesin-1 and dynein–dynactin driving the transport [103,104]. Another hypomethylated gene, KRT24, is highly expressed in trophoblast cells in the placenta [105] and forms the intermediate filaments that is part of the cytoskeleton. Yet another hypomethylated gene found to be expressed in the placenta is RAB5C [96]. It is an intracellular protein that is a member of the RAS oncogene family, located in endosomes and thought to be involved in docking and/or fusion of vesicles [106]. Not directly associated with the above-mentioned mechanisms, the peroxisomal membrane protein SLC25A17, known as a coenzyme A transporter and vital for functional peroxisomes [107], was also hypomethylated and belongs to mitochondrial solute carriers. It is found to have cytoplasmic and membrane expression in all tissues, including trophoblasts in placenta [96]. Finally, our analysis showed hypomethylation at a position associated to LATS2, which is a primary upstream negative regulator of YAP1 [108], again pointing at the significance of balance in YAP1 activity for normal placenta development.
Hypermethylation was found at positions associated with genes RTN4RL1 and PRAC. The first one, RTN4RL1, is a cell surface protein that is expressed at low levels in the placenta and predominantly in the brain [96]. Signaling via RTN4RL1 leads to activation of Rho and downstream restriction of actin polymerization, thereby regulating actin cytoskeleton during neural development to control the number of synapses formed [109]. The second one, PRAC, is a small nuclear protein and not reported to be expressed in the placenta. Instead, it is previously identified as prostate cancer susceptibility genes PRAC1 and PRAC2, specifically over-expressed in human prostate and colon cancer [110].
For pregnancies with male fetuses, one DMP was detected in the comparison between PE cases with low NOx exposure and controls with high NOx exposure. Although we see no rationale or relevant scientific question for evaluating this group comparison, the gene STAT3 is still of interest for PE. STAT3 is a transcription factor and mediates expression of a variety of genes in response to cell stimuli, and thus plays a key role in many cellular processes such as cell growth and apoptosis. Decreased STAT3 expression and activation have been suggested to play an important role in the pathogenesis of PE [111,112,113].
Our data further supports sexual dimorphism in the placenta and an association between PE and fetal sex. Female placentas are suggested to have a growth trajectory and invasion that is generally slower compared to male placentas [114,115], which is connected to cytoskeleton and migration. It is also suggested that disturbances, or in-utero insults, during early placenta development may affect growth of the female fetus more than the male, which is due to differences in placental growth [115]. Additionally, placental function during gestation is different between female and male placentas when it comes to strategies for metabolism and transport of nutrients, where male placentas have higher nutrient transfer capacity [114,116].
A gestational age deceleration among PE placentas exposed to high NOx resembles reports about placental delayed villous maturation or villous immaturity associated with placental insufficiency, fetal hypoxia, term fetal death, and fetal growth restriction [117]. Normal maturation of the placenta involves formation of the new terminal villi and syncytium in response to increasing demands from the fetus. These are all activities that require changes in the cellular cytoskeleton. Villous immaturity is characterized by reduced formation of terminal villi [118], resulting in impaired placental function. In PE, there is evidence of altered syncytium homeostasis with reduced trophoblast proliferation, a reduced level of trophoblast fusion to create new syncytium, increased syncytiotrophoblast apoptosis, and reduced production of placental growth factor (PlGF) [119,120]. This is in accordance with our DNA methylation results, indicating that cytoskeletal dynamics as well as trophoblast function/apoptosis seem to be altered in PE cases exposed to high levels of NOx and might be contributing to reduced maturation of the placenta. Placental age deceleration has previously been associated with maternal cardio-metabolic factors such as pre-pregnancy obesity and high blood pressure in women carrying male fetuses [121], further supporting sexual dimorphism. Interestingly, these are also known maternal risk factors for developing PE [122]. Additionally, exposure to a negative pregnancy environment associated with insulin-treated gestational diabetes or maternal Sjögren’s syndrome have been linked to epigenetic age deceleration in umbilical cord blood [123]. Early-onset, but not late-onset, PE has previously been associated with placental age acceleration [19,124]. However, in this study, very few PE cases belong to early-onset PE (1 out of 29 for all cases, and 1 out 19 for high NOx exposure). Taken together, our data suggest a scenario where early-onset PE with poor placentation, hypoxia, and oxidative stress result in compensatory mechanisms such as accelerated villous maturation and epigenetic age acceleration, while late-onset PE, on the other hand, affected by maternal risk factors and air pollution exposure, results in placental immaturity and epigenetic age deceleration.
The CPDB analyses highlighted three specific transcription factors. The first one, HNF4A, is involved in several forms of cancer [125,126] and implicated in the formation of progenitor cells from stem cells. The next one, E2F4, is a transcriptional repressor involved in many different functions in the cell, such as cell cycle regulation and regenerative processes in stem cells and in cancer [127]. Both HNF4A and E2F4 are implicated in signaling pathways related to tumorigenesis such as proliferation, invasion, self-renewal, and apoptosis [128]. The last one, TAF1, is an X-linked gene and the largest component of the basal transcription factor complex and therefore important for transcription in general. It has a binding site in the progesterone receptor gene promoter [129]. Variants of the gene or disrupted expression are implicated in neurodevelopment [130,131] but mutated variants of TAF1 have also been implicated in cancer [132]. Additionally, here, we suggest that the impact these transcriptions factors have on molecular mechanisms in tumor growth and metastasis can be extrapolated to similar processes during placenta formation.
Several important limitations of this study need to be considered. Firstly, due to the low quality of RNA for most of the placenta samples, the power of the DEG analyses was very low, and it was not possible to expand our analysis of DNA methylation by correlating it to gene expression on a genome-wide scale. Secondly, although our findings of differential DNA methylation associated with ambient NOx exposure during the first trimester and PE were assessed on mostly term placentas, it has been postulated and confirmed that DNA methylation levels increase during the first trimester to reach a consistent level during second and third trimesters [133,134]. Thirdly, the current investigation was limited by modeling only NOx exposure and no other pollutants such as PM2.5 was considered due to practical constraints in modeling other ambient air pollutants during exposure assessment. Moreover, the modeled concentration of NOx pertains to ambient air pollution at residence level but not accounting for total exposure, including indoor, occupational, and commuting exposures, which is a common and inherent limitation of epidemiologic studies of air pollution and PE [13]. Fourthly, we did not validate the findings in a validation cohort. However, it was not possible to find a relevant validation cohort consisting of placentas from PE cases and controls with data on NOx exposure. Lastly, it was not possible to do sub-analysis by differentiating between early- and late-onset PE [8] due to a small sample size for the early-onset group. Additionally, four of the DMPs contained SNPs in and around the CpG site, which could bias the results due to impaired base extension since the probe cannot hybridize properly. This influence is stronger the closer the SNP is located to the CpG site. The SNPs in the HAAO cg17283620 and in the non-annotated cg16162930 were situated by the CpG sites, and this warrants caution. Thus, the results for these CpG sites are likely to be influenced by some genetic variation, but we estimate it to be like the CEU population. However, the two other SNPs (in SERPINF2 and LATS2) did not include any SNPs nearby the CpG site.
The key strengths of this study are its access to a unique biobank designed to study PE and the individual exposure assessment by a local-scale high-resolution dispersion modeling. Additionally, since maternal smoking during pregnancy alters DNA methylation [135,136], we only included mothers who never smoked by their self-report. In addition, placenta is regarded as an appropriate tissue for DNA methylation analysis [137].

5. Conclusions

This study has enhanced the scientific understanding of how early pregnancy exposure to high concentrations of ambient NOx affect placental DNA methylation in PE. These changes resulted in placental immaturity and showed that sexual dimorphism may be an important factor in the etiology.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/toxics9120338/s1, Figure S1: Summary of groups based on a combination of pregnancy status (PE or normotensive control) and exposure to ambient NOx (low or high NOx exposure based on median split) (Table A) and the pair-wise comparisons performed in the DMP analyses, where comparisons that had DMPs (q < 0.05) are marked in bold (Table B), Figure S2: Potential gene regulatory interactions (induced network molecule analysis) including genes with a q < 0.1 in the DMP analyses. Gene interaction networks are shown in purple, while protein interaction networks are shown in yellow. Figure A shows a comparison of controls with low NOx exposure vs. PE cases with high NOx exposure, figure B shows a comparison of cases and controls, both with high NOx exposure, and figure C shows a comparison of controls with low NOx exposure vs. PE cases with high NOx exposure in placentas from female fetuses. Table S1: Differences in gestational age between groups based on combinations of PE status and NOx exposure, analysed using linear regression, Table S2: Differentially methylated positions (DMPs) when comparing controls with high NOx exposure, (reference) vs. PE cases with high NOx exposure in pregnancies with female fetuses. The twelve DMPs with lowest q-value are already presented in Table 2 in the main text, thus the DMPs presented here are the DMPs ordered 14–50 ranked by q-value. Table S3: Top Gene Ontology (GO) terms for analyses including genes with CpGs with q < 0.1 in the DMP analyses for different combined PE status and NOx group comparisons, Table S4: Comparisons between the study populations included in the placental DNA methylation analyses and RNA sequencing analyses.

Author Contributions

Conceptualization, C.I., E.M., L.E. and S.R.H.; data curation, E.M., L.E. and S.R.H.; formal analysis, K.E., Y.M., L.G. and Z.M.; funding acquisition, C.I., E.M., L.E. and S.R.H.; investigation, Y.M.; methodology, K.E., Y.M., L.G., E.M., L.E. and S.R.H.; project administration, E.M. and S.R.H.; resources, Z.M., E.M. and S.R.H.; supervision, L.E. and S.R.H.; visualization, K.E.; writing—original draft, K.E., Y.M., Z.M. and L.E.; writing—review and editing, K.E., Y.M., L.G., Z.M., C.I., E.M., L.E. and S.R.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Swedish Research Council (DNR 2016-01003, 2018-04840, 2018-03156), FORTE (DNR 2016-00572) ALF, Region Skåne, and SUS foundations and with support from the Erasmus + Programme of the European Union (Framework Agreement number: 2013-0040). This publication reflects only the views of the authors, and the Commission cannot be held responsible for any use that may be made of the information contained therein.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by Lund University Ethical Review Board, LU803-2 (2015/14), Lund, Sweden.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data underlying this article cannot be shared publicly due to the privacy of individuals that participated in the study. The datasets related to the clinical parameters and the DNA analysis are available upon request to the corresponding author. However, the datasets related to individual exposure levels of air pollution are not readily available because of limitations set in the ethical permission stating that only the researchers involved in the project are allowed to access data. Requests to access the dataset should be directed to data-holders, Statistic Sweden, and The Swedish National Board of Health and Welfare for birth registry data.

Acknowledgments

We are grateful to Eva Hansson for her support in DNA/RNA extractions, Ralf Rittner for his support in exposure assessment, and Kristoffer Mattisson for his support in creating an exposure map. We thank tbe Center for Translational Genomics, Lund University, and Clinical Genomics Lund, SciLifeLab, for providing sequencing and DNA methylation services.

Conflicts of Interest

Stefan R. Hansson holds a patent related to diagnosis and treatment of PE and is a co-founder of Guard Therapeutics International AB (formerly named A1M Pharma AB). The rest of the authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Khan, N.; Andrade, W.; De Castro, H.; Wright, A.; Wright, D.; Nicolaides, K.H. Impact of New Definitions of Pre-Eclampsia on Incidence and Performance of First-Trimester Screening. Ultrasound Obstet. Gynecol. 2020, 55, 50–57. [Google Scholar] [CrossRef] [PubMed]
  2. Umesawa, M.; Kobashi, G. Epidemiology of Hypertensive Disorders in Pregnancy: Prevalence, Risk Factors, Predictors and Prognosis. Hypertens. Res. 2017, 40, 213–220. [Google Scholar] [CrossRef] [PubMed]
  3. Thornton, C.; Dahlen, H.; Korda, A.; Hennessy, A. The Incidence of Preeclampsia and Eclampsia and Associated Maternal Mortality in Australia from Population-Linked Datasets: 2000–2008. Am. J. Obstet. Gynecol. 2013, 208, 476.e1–476.e5. [Google Scholar] [CrossRef]
  4. MacKay, A.P.; Berg, C.J.; Atrash, H.K. Pregnancy-Related Mortality from Preeclampsia and Eclampsia. Obstet. Gynecol. 2001, 97, 533–538. [Google Scholar] [CrossRef] [PubMed]
  5. Burton, G.J.; Fowden, A.L.; Thornburg, K.L. Placental Origins of Chronic Disease. Physiol. Rev. 2016, 96, 1509–1565. [Google Scholar] [CrossRef]
  6. Davis, E.F.; Lazdam, M.; Lewandowski, A.J.; Worton, S.A.; Kelly, B.; Kenworthy, Y.; Adwani, S.; Wilkinson, A.R.; McCormick, K.; Sargent, I.; et al. Cardiovascular Risk Factors in Children and Young Adults Born to Preeclamptic Pregnancies: A Systematic Review. Pediatrics 2012, 129, e1552–e1561. [Google Scholar] [CrossRef] [PubMed]
  7. 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. Hypertension 2018, 72, 24–43. [Google Scholar] [CrossRef] [Green Version]
  8. Redman, C.W.G.; Staff, A.C.; Roberts, J.M. Syncytiotrophoblast Stress in Preeclampsia: The Convergence Point for Multiple Pathways. Am. J. Obstet. Gynecol. 2021, S0002–9378, 31115-7. [Google Scholar] [CrossRef] [PubMed]
  9. Raymond, D.; Peterson, E. A Critical Review of Early-Onset and Late-Onset Preeclampsia. Obstet. Gynecol. Surv. 2011, 66, 497–506. [Google Scholar] [CrossRef] [PubMed]
  10. Taylor, B.D.; Ness, R.B.; Klebanoff, M.A.; Tang, G.; Roberts, J.M.; Hougaard, D.M.; Skogstrand, K.; Haggerty, C.L. The Impact of Female Fetal Sex on Preeclampsia and the Maternal Immune Milieu. Pregnancy Hypertens. 2018, 12, 53–57. [Google Scholar] [CrossRef] [PubMed]
  11. Elsmén, E.; Källén, K.; Maršál, K.; Hellström-Westas, L. Fetal Gender and Gestational-Age-Related Incidence of Pre-Eclampsia. Acta Obstet. Gynecol. Scand. 2006, 85, 1285–1291. [Google Scholar] [CrossRef] [PubMed]
  12. Gyselaers, W. Preeclampsia Is a Syndrome with a Cascade of Pathophysiologic Events. J. Clin. Med. 2020, 9, 2245. [Google Scholar] [CrossRef]
  13. Pedersen, M.; Stayner, L.; Slama, R.; Sorensen, M.; Figueras, F.; Nieuwenhuijsen, M.J.; Raaschou-Nielsen, O.; Dadvand, P. Ambient Air Pollution and Pregnancy-Induced Hypertensive Disorders: A Systematic Review and Meta-Analysis. Hypertension 2014, 64, 494–500. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Malmqvist, E.; Jakobsson, K.; Tinnerberg, H.; Rignell-Hydbom, A.; Rylander, L. Gestational Diabetes and Preeclampsia in Association with Air Pollution at Levels below Current Air Quality Guidelines. Environ. Health Perspect. 2013, 121, 488–493. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Nääv, Å.; Erlandsson, L.; Isaxon, C.; Åsander Frostner, E.; Ehinger, J.; Sporre, M.K.; Krais, A.M.; Strandberg, B.; Lundh, T.; Elmér, E.; et al. Urban PM2.5 Induces Cellular Toxicity, Hormone Dysregulation, Oxidative Damage, Inflammation, and Mitochondrial Interference in the HRT8 Trophoblast Cell Line. Front. Endocrinol. 2020, 11, 75. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Familari, M.; Nääv, Å.; Erlandsson, L.; de Iongh, R.U.; Isaxon, C.; Strandberg, B.; Lundh, T.; Hansson, S.R.; Malmqvist, E. Exposure of Trophoblast Cells to Fine Particulate Matter Air Pollution Leads to Growth Inhibition, Inflammation and ER Stress. PLoS ONE 2019, 14, e0218799. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. World Health Organization. Ambient Air Pollution: A Global Assessment of Exposure and Burden of Disease. Clean Air J. 2016, 26, 6. [Google Scholar] [CrossRef]
  18. Bai, W.; Li, Y.; Niu, Y.; Ding, Y.; Yu, X.; Zhu, B.; Duan, R.; Duan, H.; Kou, C.; Li, Y.; et al. Association between Ambient Air Pollution and Pregnancy Complications: A Systematic Review and Meta-Analysis of Cohort Studies. Environ. Res. 2020, 185, 109471. [Google Scholar] [CrossRef] [PubMed]
  19. Bianco-Miotto, T.; Mayne, B.T.; Buckberry, S.; Breen, J.; Rodriguez Lopez, C.M.; Roberts, C.T. Recent Progress towards Understanding the Role of DNA Methylation in Human Placental Development. Reproduction 2016, 152, R23–R30. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Collas, P.; Noer, A.; Timoskainen, S. Programming the Genome in Embryonic and Somatic Stem Cells. J. Cell. Mol. Med. 2007, 11, 602–620. [Google Scholar] [CrossRef] [Green Version]
  21. Jirtle, R.L.; Skinner, M.K. Environmental Epigenomics and Disease Susceptibility. Nat. Rev. Genetics 2007, 8, 253–262. [Google Scholar] [CrossRef]
  22. Saenen, N.D.; Martens, D.S.; Neven, K.Y.; Alfano, R.; Bové, H.; Janssen, B.G.; Roels, H.A.; Plusquin, M.; Vrijens, K.; Nawrot, T.S. Air Pollution-Induced Placental Alterations: An Interplay of Oxidative Stress, Epigenetics, and the Aging Phenotype? Clin. Epigenet. 2019, 11, 124. [Google Scholar] [CrossRef] [PubMed]
  23. Cirkovic, A.; Garovic, V.; Milin Lazovic, J.; Milicevic, O.; Savic, M.; Rajovic, N.; Aleksic, N.; Weissgerber, T.; Stefanovic, A.; Stanisavljevic, D.; et al. Systematic Review Supports the Role of DNA Methylation in the Pathophysiology of Preeclampsia: A Call for Analytical and Methodological Standardization. Biol. Sex Differ. 2020, 11, 36. [Google Scholar] [CrossRef]
  24. Neven, K.Y.; Saenen, N.D.; Tarantini, L.; Janssen, B.G.; Lefebvre, W.; Vanpoucke, C.; Bollati, V.; Nawrot, T.S. Placental Promoter Methylation of DNA Repair Genes and Prenatal Exposure to Particulate Air Pollution: An ENVIR ON AGE Cohort Study. Lancet Planet. Health 2018, 2, e174–e183. [Google Scholar] [CrossRef] [Green Version]
  25. World Health Organization. WHO Global Air Quality Guidelines: Particulate Matter (PM2.5 and PM10), Ozone, Nitrogen Dioxide, Sulfur Dioxide and Carbon Monoxideitle; World Health Organization: Geneva, Switzerland, 2021; ISBN 9789240034228. [Google Scholar]
  26. Abraham, E.; Rousseaux, S.; Agier, L.; Giorgis-Allemand, L.; Tost, J.; Galineau, J.; Hulin, A.; Siroux, V.; Vaiman, D.; Charles, M.-A.; et al. Pregnancy Exposure to Atmospheric Pollution and Meteorological Conditions and Placental DNA Methylation. Environ. Int. 2018, 118, 334–347. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Cox, L.S.; Redman, C. The Role of Cellular Senescence in Ageing of the Placenta. Placenta 2017, 52, 139–145. [Google Scholar] [CrossRef] [Green Version]
  28. Horvath, S. DNA Methylation Age of Human Tissues and Cell Types. Genome Biol. 2013, 14, R115. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Lee, Y.; Choufani, S.; Weksberg, R.; Wilson, S.L.; Yuan, V.; Burt, A.; Marsit, C.; Lu, A.T.; Ritz, B.; Bohlin, J.; et al. Placental Epigenetic Clocks: Estimating Gestational Age Using Placental DNA Methylation Levels. Aging 2019, 11, 4238–4253. [Google Scholar] [CrossRef]
  30. Opsis. 2006. Available online: http://Www.Opsis.Se/ (accessed on 5 November 2021).
  31. Mandakh, Y.; Rittner, R.; Flanagan, E.; Oudin, A.; Isaxon, C.; Familari, M.; Hansson, S.R.; Malmqvist, E. Maternal Exposure to Ambient Air Pollution and Risk of Preeclampsia: A Population-Based Cohort Study in Scania, Sweden. Int. J. Environ. Res. Public Health 2020, 17, 1744. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Iodice, S.; Hoxha, M.; Ferrari, L.; Carbone, I.F.; Anceschi, C.; Miragoli, M.; Pesatori, A.C.; Persico, N.; Bollati, V. Particulate Air Pollution, Blood Mitochondrial DNA Copy Number, and Telomere Length in Mothers in the First Trimester of Pregnancy: Effects on Fetal Growth. Oxidative Med. Cell. Longev. 2018, 2018, 5162905. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Martens, D.S.; Cox, B.; Janssen, B.G.; Clemente, D.B.P.; Gasparrini, A.; Vanpoucke, C.; Lefebvre, W.; Roels, H.A.; Plusquin, M.; Nawrot, T.S. Prenatal Air Pollution and Newborns’ Predisposition to Accelerated Biological Aging. JAMA Pediatrics 2017, 171, 1160. [Google Scholar] [CrossRef] [PubMed]
  34. Melody, S.M.; Wills, K.; Knibbs, L.D.; Ford, J.; Venn, A.; Johnston, F. Maternal Exposure to Ambient Air Pollution and Pregnancy Complications in Victoria, Australia. Int. J. Environ. Res. Public Health 2020, 17, 2572. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Morris, T.J.; Butcher, L.M.; Feber, A.; Teschendorff, A.E.; Chakravarthy, A.R.; Wojdacz, T.K.; Beck, S. ChAMP: 450 k Chip Analysis Methylation Pipeline. Bioinformatics 2014, 30, 428–430. [Google Scholar] [CrossRef] [PubMed]
  36. Tian, Y.; Morris, T.J.; Webster, A.P.; Yang, Z.; Beck, S.; Feber, A.; Teschendorff, A.E. ChAMP: Updated Methylation Analysis Pipeline for Illumina BeadChips. Bioinformatics 2017, 33, 3982–3984. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Fortin, J.-P.; Triche, T.J.; Hansen, K.D. Preprocessing, Normalization and Integration of the Illumina HumanMethylationEPIC Array. bioRxiv 2017, 15, 558–560. [Google Scholar] [CrossRef] [PubMed]
  38. Aryee, M.J.; Jaffe, A.E.; Corrada-Bravo, H.; Ladd-Acosta, C.; Feinberg, A.P.; Hansen, K.D.; Irizarry, R.A. Minfi: A Flexible and Comprehensive Bioconductor Package for the Analysis of Infinium DNA Methylation Microarrays. Bioinformatics 2014, 30, 1363–1369. [Google Scholar] [CrossRef] [Green Version]
  39. Teschendorff, A.E.; Marabita, F.; Lechner, M.; Bartlett, T.; Tegner, J.; Gomez-Cabrero, D.; Beck, S. A Beta-Mixture Quantile Normalization Method for Correcting Probe Design Bias in Illumina Infinium 450 k DNA Methylation Data. Bioinformatics 2013, 29, 189–196. [Google Scholar] [CrossRef] [Green Version]
  40. Yuan, V.; Hui, D.; Yin, Y.; Peñaherrera, M.S.; Beristain, A.G.; Robinson, W.P. Cell-Specific Characterization of the Placental Methylome. BMC Genom. 2021, 22, 6. [Google Scholar] [CrossRef]
  41. Khouja, J.N.; Simpkin, A.J.; O’Keeffe, L.M.; Wade, K.H.; Houtepen, L.C.; Relton, C.L.; Suderman, M.; Howe, L.D. Epigenetic Gestational Age Acceleration: A Prospective Cohort Study Investigating Associations with Familial, Sociodemographic and Birth Characteristics. Clin. Epigenet. 2018, 10, 86. [Google Scholar] [CrossRef]
  42. Johnson, W.E.; Li, C.; Rabinovic, A. Adjusting Batch Effects in Microarray Expression Data Using Empirical Bayes Methods. Biostatistics 2007, 8, 118–127. [Google Scholar] [CrossRef]
  43. Leek, J.T.; Johnson, W.E.; Parker, H.S.; Fertig, E.J.; Jaffe, A.E.; Storey, J.D.; Zhang, Y.; Torres, L.C. Sva: Surrogate Variable Analysis. R Package Version 3.24.4. 2017. Available online: https://bioconductor.org/packages/release/bioc/html/sva.html (accessed on 18 November 2021).
  44. Gude, N.M.; Roberts, C.T.; Kalionis, B.; King, R.G. Growth and Function of the Normal Human Placenta. Thromb. Res. 2004, 114, 397–407. [Google Scholar] [CrossRef]
  45. Reyes, L.; Golos, T.G. Hofbauer Cells: Their Role in Healthy and Complicated Pregnancy. Front. Immunol. 2018, 9, 2628. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Ritchie, M.E.; Phipson, B.; Wu, D.; Hu, Y.; Law, C.W.; Shi, W.; Smyth, G.K. Limma Powers Differential Expression Analyses for RNA-Sequencing and Microarray Studies. Nucleic Acids Res. 2015, 43, e47. [Google Scholar] [CrossRef] [PubMed]
  47. Phipson, B.; Maksimovic, J.; Oshlack, A. MissMethyl: An R Package for Analyzing Data from Illumina’s HumanMethylation450 Platform. Bioinformatics 2016, 32, 286–288. [Google Scholar] [CrossRef]
  48. Andrews, S. FastQC: A Quality Control Tool for High Throughput Sequence Data. 2010. Available online: http://www.bioinformatics.babraham.ac.uk/projects/fastqc (accessed on 19 November 2021).
  49. Pertea, M.; Pertea, G.M.; Antonescu, C.M.; Chang, T.-C.; Mendell, J.T.; Salzberg, S.L. StringTie Enables Improved Reconstruction of a Transcriptome from RNA-Seq Reads. Nat. Biotechnol. 2015, 33, 290–295. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Love, M.I.; Huber, W.; Anders, S. Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Frazer, K.A.; Ballinger, D.G.; Cox, D.R.; Hinds, D.A.; Stuve, L.L.; Gibbs, R.A.; Belmont, J.W.; Boudreau, A.; Hardenbol, P.; Leal, S.M.; et al. A Second Generation Human Haplotype Map of over 3.1 Million SNPs. Nature 2007, 449, 851–861. [Google Scholar] [CrossRef] [PubMed]
  52. Altshuler, D.M.; Gibbs, R.A.; Peltonen, L.; Schaffner, S.F.; Yu, F.; Dermitzakis, E.; Bonnen, P.E.; De Bakker, P.I.W.; Deloukas, P.; Gabriel, S.B.; et al. Integrating Common and Rare Genetic Variation in Diverse Human Populations. Nature 2010, 467, 52–58. [Google Scholar] [CrossRef]
  53. Seidmann, L.; Suhan, T.; Unger, R.; Gerein, V.; Kirkpatrick, C.J. Transient CD15-Positive Endothelial Phenotype in the Human Placenta Correlates with Physiological and Pathological Fetoplacental Immaturity. Eur. J. Obstet. Gynecol. Reprod. Biol. 2014, 180, 172–179. [Google Scholar] [CrossRef] [PubMed]
  54. Costanzo, V.; Bardelli, A.; Siena, S.; Abrignani, S. Exploring the Links between Cancer and Placenta Development. Open Biol. 2018, 8, 180081. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Holtan, S.G.; Creedon, D.J.; Haluska, P.; Markovic, S.N. Cancer and Pregnancy: Parallels in Growth, Invasion, and Immune Modulation and Implications for Cancer Therapeutic Agents. Mayo Clin. Proc. 2009, 84, 985–1000. [Google Scholar] [CrossRef] [Green Version]
  56. Oliveira, M.M.S.; Salgado, C.M.; Viana, L.R.; Gomes-Marcondes, M.C.C. Pregnancy and Cancer: Cellular Biology and Mechanisms Affecting the Placenta. Cancers 2021, 13, 1667. [Google Scholar] [CrossRef] [PubMed]
  57. Meinhardt, G.; Haider, S.; Kunihs, V.; Saleh, L.; Pollheimer, J.; Fiala, C.; Hetey, S.; Feher, Z.; Szilagyi, A.; Than, N.G.; et al. Pivotal Role of the Transcriptional Co-Activator YAP in Trophoblast Stemness of the Developing Human Placenta. Proc. Natl. Acad. Sci. USA 2020, 117, 13562–13570. [Google Scholar] [CrossRef] [PubMed]
  58. Zha, W.; Guan, S.; Liu, N.; Li, Y.; Tian, Y.; Chen, Y.; Wang, Y.; Wu, F. Let-7a Inhibits Bcl-Xl and YAP1 Expression to Induce Apoptosis of Trophoblast Cells in Early-Onset Severe Preeclampsia. Sci. Total Environ. 2020, 745, 139919. [Google Scholar] [CrossRef] [PubMed]
  59. Wang, H.; Xu, P.; Luo, X.; Hu, M.; Liu, Y.; Yang, Y.; Peng, W.; Bai, Y.; Chen, X.; Tan, B.; et al. Phosphorylation of Yes-Associated Protein Impairs Trophoblast Invasion and Migration: Implications for the Pathogenesis of Fetal Growth Restriction†. Biol. Reprod. 2020, 103, 866–879. [Google Scholar] [CrossRef] [PubMed]
  60. Sun, M.; Na, Q.; Huang, L.; Song, G.; Jin, F.; Li, Y.; Hou, Y.; Kang, D.; Qiao, C. YAP Is Decreased in Preeclampsia and Regulates Invasion and Apoptosis of HTR-8/SVneo. Reprod. Sci. 2018, 25, 1382–1393. [Google Scholar] [CrossRef]
  61. Liu, R.; Wei, C.; Ma, Q.; Wang, W. Hippo-YAP1 Signaling Pathway and Severe Preeclampsia (SPE) in the Chinese Population. Pregnancy Hypertens. 2020, 19, 1–10. [Google Scholar] [CrossRef]
  62. Turowski, G.; Vogel, M. Re-View and View on Maturation Disorders in the Placenta. Apmis 2018, 126, 602–612. [Google Scholar] [CrossRef] [PubMed]
  63. Jin, Y.; Cui, D.; Ren, J.; Wang, K.; Zeng, T.; Gao, L. CACNA2D3 Is Downregulated in Gliomas and Functions as a Tumor Suppressor. Mol. Carcinog. 2017, 56, 945–959. [Google Scholar] [CrossRef]
  64. Baczyk, D.; Kingdom, J.C.P.; Uhlén, P. Calcium Signaling in Placenta. Cell Calcium 2011, 49, 350–356. [Google Scholar] [CrossRef] [PubMed]
  65. Jakoubek, V.; Bíbová, J.; Hampl, V. Voltage-Gated Calcium Channels Mediate Hypoxic Vasoconstriction in the Human Placenta. Placenta 2006, 27, 1030–1033. [Google Scholar] [CrossRef]
  66. Hansson, S.R.; Chen, Y.; Brodszki, J.; Chen, M.; Hernandez-Andrade, E.; Inman, J.M.; Kozhich, O.A.; Larsson, I.; Marsal, K.; Medstrand, P.; et al. Gene Expression Profiling of Human Placentas from Preeclamptic and Normotensive Pregnancies. Mol. Hum. Reprod. 2006, 12, 169–179. [Google Scholar] [CrossRef] [PubMed]
  67. Li, H.Q.; Fan, J.J.; Li, X.H.; Bao, D. MiR-507 Inhibits the Growth and Invasion of Trophoblasts by Targeting CAMK4. Eur. Rev. Med. Pharmacol. Sci. 2020, 24, 5856–5862. [Google Scholar] [CrossRef] [PubMed]
  68. Wang, T.; Xiang, Y.; Zhou, X.; Zheng, X.; Zhang, H.; Zhang, X.; Zhang, J.; He, L.; Zhao, X. Epigenome-Wide Association Data Implicate Fetal/Maternal Adaptations Contributing to Clinical Outcomes in Preeclampsia. Epigenomics 2019, 11, 1003–1019. [Google Scholar] [CrossRef] [Green Version]
  69. Schoentgen, F.; Jonic, S. PEBP1/RKIP Behavior: A Mirror of Actin-Membrane Organization. Cell. Mol. Life Sci. 2020, 77, 859–874. [Google Scholar] [CrossRef] [PubMed]
  70. Ciarmela, P.; Marzioni, D.; Islam, M.S.; Gray, P.C.; Terracciano, L.; Lorenzi, T.; Todros, T.; Petraglia, F.; Castellucci, M. Possible Role of RKIP in Cytotrophoblast Migration: Immunohistochemical and in Vitro Studies. J. Cell. Physiol. 2012, 227, 1821–1828. [Google Scholar] [CrossRef] [PubMed]
  71. Goswami, D.; Tannetta, D.S.; Magee, L.A.; Fuchisawa, A.; Redman, C.W.; Sargent, I.L.; von Dadelszen, P. Excess Syncytiotrophoblast Microparticle Shedding Is a Feature of Early-Onset Pre-Eclampsia, but Not Normotensive Intrauterine Growth Restriction. Placenta 2006, 27, 56–61. [Google Scholar] [CrossRef] [PubMed]
  72. Frémont, S.; Hammich, H.; Bai, J.; Wioland, H.; Klinkert, K.; Rocancourt, M.; Kikuti, C.; Stroebel, D.; Romet-Lemonne, G.; Pylypenko, O.; et al. Oxidation of F-Actin Controls the Terminal Steps of Cytokinesis. Nat. Commun. 2017, 8, 14528. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Deng, W.; Wang, Y.; Gu, L.; Duan, B.; Cui, J.; Zhang, Y.; Chen, Y.; Sun, S.; Dong, J.; Du, J. MICAL1 Controls Cell Invasive Phenotype via Regulating Oxidative Stress in Breast Cancer Cells. BMC Cancer 2016, 16, 489. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Li, W.; Li, Y. The Pseudogene URAHP Promotes Proliferation and Regulates the Pathogenesis of Preeclampsia. Am. J. Transl. Res. 2020, 12, 4715–4727. [Google Scholar] [PubMed]
  75. Gracilla, D.E.; Korla, P.K.; Lai, M.; Chiang, A.; Liou, W.; Sheu, J.J. Overexpression of Wild Type or a Q311E Mutant MB21D2 Promotes a Pro-oncogenic Phenotype in HNSCC. Mol. Oncol. 2020, 14, 3065–3082. [Google Scholar] [CrossRef] [PubMed]
  76. Fabriek, B.O.; Dijkstra, C.D.; van den Berg, T.K. The Macrophage Scavenger Receptor CD163. Immunobiology 2005, 210, 153–160. [Google Scholar] [CrossRef] [PubMed]
  77. Joerink, M.; Rindsjo, E.; van Riel, B.; Alm, J.; Papadogiannakis, N. Placental Macrophage (Hofbauer Cell) Polarization Is Independent of Maternal Allergen-Sensitization and Presence of Chorioamnionitis. Placenta 2011, 32, 380–385. [Google Scholar] [CrossRef] [PubMed]
  78. Tang, Z.; Buhimschi, I.A.; Buhimschi, C.S.; Tadesse, S.; Norwitz, E.; Niven-Fairchild, T.; Huang, S.T.; Guller, S. Decreased Levels of Folate Receptor-Beta and Reduced Numbers of Fetal Macrophages (Hofbauer Cells) in Placentas from Pregnancies with Severe Pre-Eclampsia. Am. J. Reprod. Immunol. 2013, 70, 104–115. [Google Scholar] [CrossRef] [Green Version]
  79. Yang, S.W.; Cho, E.H.; Choi, S.Y.; Lee, Y.K.; Park, J.H.; Kim, M.K.; Park, J.Y.; Choi, H.J.; Lee, J.I.; Ko, H.M.; et al. DC-SIGN Expression in Hofbauer Cells May Play an Important Role in Immune Tolerance in Fetal Chorionic Villi during the Development of Preeclampsia. J. Reprod. Immunol. 2017, 124, 30–37. [Google Scholar] [CrossRef]
  80. Munn, D.H.; Zhou, M.; Attwood, J.T.; Bondarev, I.; Conway, S.J.; Marshall, B.; Brown, C.; Mellor, A.L. Prevention of Allogeneic Fetal Rejection by Tryptophan Catabolism. Science 1998, 281, 1191–1193. [Google Scholar] [CrossRef]
  81. Broekhuizen, M.; Klein, T.; Hitzerd, E.; de Rijke, Y.B.; Schoenmakers, S.; Sedlmayr, P.; Danser, A.H.J.; Merkus, D.; Reiss, I.K.M. L-Tryptophan–Induced Vasodilation Is Enhanced in Preeclampsia. Hypertension 2020, 76, 184–194. [Google Scholar] [CrossRef]
  82. Keaton, S.A.; Heilman, P.; Bryleva, E.Y.; Madaj, Z.; Krzyzanowski, S.; Grit, J.; Miller, E.S.; Jalmby, M.; Kalapotharakos, G.; Racicot, K.; et al. Altered Tryptophan Catabolism in Placentas From Women With Pre-Eclampsia. Int. J. Tryptophan Res. 2019, 12, 1178646919840321. [Google Scholar] [CrossRef]
  83. Hanswijk, S.I.; Spoelder, M.; Shan, L.; Verheij, M.M.M.; Muilwijk, O.G.; Li, W.; Liu, C.; Kolk, S.M.; Homberg, J.R. Gestational Factors throughout Fetal Neurodevelopment: The Serotonin Link. Int. J. Mol. Sci. 2020, 21, 5850. [Google Scholar] [CrossRef]
  84. Yu, J.; Madison, J.M.; Mundlos, S.; Winston, F.; Olsen, B.R. Characterization of a Human Homologue of TheSaccharomyces CerevisiaeTranscription Factor Spt3 (SUPT3H). Genomics 1998, 53, 90–96. [Google Scholar] [CrossRef] [Green Version]
  85. Perez-Garcia, V.; Fineberg, E.; Wilson, R.; Murray, A.; Mazzeo, C.I.; Tudor, C.; Sienerth, A.; White, J.K.; Tuck, E.; Ryder, E.J.; et al. Placentation Defects Are Highly Prevalent in Embryonic Lethal Mouse Mutants. Nature 2018, 555, 463–468. [Google Scholar] [CrossRef] [PubMed]
  86. Carpenter, S.L.; Mathew, P. Alpha2-Antiplasmin and Its Deficiency: Fibrinolysis out of Balance. Haemophilia 2008, 14, 1250–1254. [Google Scholar] [CrossRef] [PubMed]
  87. Hou, Y.; Okada, K.; Okamoto, C.; Ueshima, S.; Matsuo, O. Alpha2-Antiplasmin Is a Critical Regulator of Angiotensin II-Mediated Vascular Remodeling. Arterioscler. Thromb. Vasc. Biol. 2008, 28, 1257–1262. [Google Scholar] [CrossRef] [Green Version]
  88. Shao, X.; Wang, Y.; Liu, Y.; Guo, X.; Li, D.; Huo, R.; Jia, W.; Cao, G.; Li, Y.-X.; Liu, M.; et al. Association of Imbalanced Sex Hormone Production with Excessive Procoagulation Factor SerpinF2 in Preeclampsia. J. Hypertens. 2019, 37, 197–205. [Google Scholar] [CrossRef]
  89. Jin, G.; Ma, H.; Wu, C.; Dai, J.; Zhang, R.; Shi, Y.; Lu, J.; Miao, X.; Wang, M.; Zhou, Y.; et al. Genetic Variants at 6p21.1 and 7p15.3 Are Associated with Risk of Multiple Cancers in Han Chinese. Am. J. Hum. Genet. 2012, 91, 928–934. [Google Scholar] [CrossRef] [Green Version]
  90. Sun, H.; Liu, K.; Huang, J.; Sun, Q.; Shao, C.; Luo, J.; Xu, L.; Shen, Y.; Ren, B. FAM111B, a Direct Target of P53, Promotes the Malignant Process of Lung Adenocarcinoma. OncoTargets Ther. 2019, 12, 2829–2842. [Google Scholar] [CrossRef] [Green Version]
  91. Kawasaki, K.; Nojima, S.; Hijiki, S.; Tahara, S.; Ohshima, K.; Matsui, T.; Hori, Y.; Kurashige, M.; Umeda, D.; Kiyokawa, H.; et al. FAM111B Enhances Proliferation of KRAS-Driven Lung Adenocarcinoma by Degrading P16. Cancer Sci. 2020, 111, 2635–2646. [Google Scholar] [CrossRef]
  92. Liu, C.; Cummins, T.R.; Tyrrell, L.; Black, J.A.; Waxman, S.G.; Dib-Hajj, S.D. CAP-1A Is a Novel Linker That Binds Clathrin and the Voltage-Gated Sodium Channel Nav1.8. Mol. Cell. Neurosci. 2005, 28, 636–649. [Google Scholar] [CrossRef] [PubMed]
  93. Cronqvist, T.; Erlandsson, L.; Tannetta, D.; Hansson, S.R. Placental Syncytiotrophoblast Extracellular Vesicles Enter Primary Endothelial Cells through Clathrin-Mediated Endocytosis. Placenta 2020, 100, 133–141. [Google Scholar] [CrossRef]
  94. Fisher, S.J. Why Is Placentation Abnormal in Preeclampsia? Am. J. Obstet. Gynecol. 2015, 213, S115–S122. [Google Scholar] [CrossRef] [Green Version]
  95. Isaevska, E.; Moccia, C.; Asta, F.; Cibella, F.; Gagliardi, L.; Ronfani, L.; Rusconi, F.; Stazi, M.A.; Richiardi, L. Exposure to Ambient Air Pollution in the First 1000 Days of Life and Alterations in the DNA Methylome and Telomere Length in Children: A Systematic Review. Environ. Res. 2021, 193, 110504. [Google Scholar] [CrossRef]
  96. Uhlen, M.; Fagerberg, L.; Hallstrom, B.M.; Lindskog, C.; Oksvold, P.; Mardinoglu, A.; Sivertsson, A.; Kampf, C.; Sjostedt, E.; Asplund, A.; et al. Proteomics. Tissue-Based Map of the Human Proteome. Science 2015, 347, 1260419. [Google Scholar] [CrossRef]
  97. Stelzer, G.; Rosen, N.; Plaschkes, I.; Zimmerman, S.; Twik, M.; Fishilevich, S.; Stein, T.I.; Nudel, R.; Lieder, I.; Mazor, Y.; et al. The GeneCards Suite: From Gene Data Mining to Disease Genome Sequence Analyses. Curr. Protoc. Bioinform. 2016, 54, 1.30.1–1.30.33. [Google Scholar] [CrossRef]
  98. Detchokul, S.; Williams, E.D.; Parker, M.W.; Frauman, A.G. Tetraspanins as Regulators of the Tumour Microenvironment: Implications for Metastasis and Therapeutic Strategies. Br. J. Pharmacol. 2014, 171, 5462–5490. [Google Scholar] [CrossRef] [Green Version]
  99. Huang, Y.; Zucker, B.; Zhang, S.; Elias, S.; Zhu, Y.; Chen, H.; Ding, T.; Li, Y.; Sun, Y.; Lou, J.; et al. Migrasome Formation Is Mediated by Assembly of Micron-Scale Tetraspanin Macrodomains. Nat. Cell Biol. 2019, 21, 991–1002. [Google Scholar] [CrossRef]
  100. Thomas, J.R.; Appios, A.; Zhao, X.; Dutkiewicz, R.; Donde, M.; Lee, C.Y.C.; Naidu, P.; Lee, C.; Cerveira, J.; Liu, B.; et al. Phenotypic and Functional Characterization of First-Trimester Human Placental Macrophages, Hofbauer Cells. J. Exp. Med. 2021, 218, e20200891. [Google Scholar] [CrossRef]
  101. Shindo-Okada, N.; Iigo, M. Expression of the Arp11 Gene Suppresses the Tumorigenicity of PC-14 Human Lung Adenocarcinoma Cells. Biochem. Biophys. Res. Commun. 2003, 312, 889–896. [Google Scholar] [CrossRef] [PubMed]
  102. Zhang, J.; Yao, X.; Fischer, L.; Abenza, J.F.; Penalva, M.A.; Xiang, X. The P25 Subunit of the Dynactin Complex Is Required for Dynein-Early Endosome Interaction. J. Cell Biol. 2011, 193, 1245–1255. [Google Scholar] [CrossRef] [Green Version]
  103. Henrichs, V.; Grycova, L.; Barinka, C.; Nahacka, Z.; Neuzil, J.; Diez, S.; Rohlena, J.; Braun, M.; Lansky, Z. Mitochondria-Adaptor TRAK1 Promotes Kinesin-1 Driven Transport in Crowded Environments. Nat. Commun. 2020, 11, 3123. [Google Scholar] [CrossRef] [PubMed]
  104. Fenton, A.R.; Jongens, T.A.; Holzbaur, E.L.F. Mitochondrial Adaptor TRAK2 Activates and Functionally Links Opposing Kinesin and Dynein Motors. Nat. Commun. 2021, 12, 4578. [Google Scholar] [CrossRef] [PubMed]
  105. Sprecher, E.; Itin, P.; Whittock, N.V.; McGrath, J.A.; Meyer, R.; DiGiovanna, J.J.; Bale, S.J.; Uitto, J.; Richard, G. Refined Mapping of Naegeli-Franceschetti-Jadassohn Syndrome to a 6 CM Interval on Chromosome 17q11.2-Q21 and Investigation of Candidate Genes. J. Investig. Dermatol. 2002, 119, 692–698. [Google Scholar] [CrossRef]
  106. Han, H.J.; Sudo, K.; Inazawa, J.; Nakamura, Y. Isolation and Mapping of a Human Gene (RABL) Encoding a Small GTP-Binding Protein Homologous to the Ras-Related RAB Gene. Cytogenet. Cell Genet. 1996, 73, 137–139. [Google Scholar] [CrossRef] [PubMed]
  107. Kim, Y.I.; Nam, I.K.; Lee, D.K.; Bhandari, S.; Charton, L.; Kwak, S.; Lim, J.Y.; Hong, K.; Kim, S.J.; Lee, J.N.; et al. Slc25a17 Acts as a Peroxisomal Coenzyme A Transporter and Regulates Multiorgan Development in Zebrafish. J. Cell. Physiol. 2020, 235, 151–165. [Google Scholar] [CrossRef] [PubMed]
  108. He, C.; Lv, X.; Huang, C.; Hua, G.; Ma, B.; Chen, X.; Angeletti, P.C.; Dong, J.; Zhou, J.; Wang, Z.; et al. YAP1-LATS2 Feedback Loop Dictates Senescent or Malignant Cell Fate to Maintain Tissue Homeostasis. EMBO Rep. 2019, 20, e44948. [Google Scholar] [CrossRef]
  109. Wills, Z.P.; Mandel-Brehm, C.; Mardinly, A.R.; McCord, A.E.; Giger, R.J.; Greenberg, M.E. The Nogo Receptor Family Restricts Synapse Number in the Developing Hippocampus. Neuron 2012, 73, 466–481. [Google Scholar] [CrossRef] [Green Version]
  110. Liu, X.F.; Olsson, P.; Wolfgang, C.D.; Bera, T.K.; Duray, P.; Lee, B.; Pastan, I. PRAC: A Novel Small Nuclear Protein That Is Specifically Expressed in Human Prostate and Colon. Prostate 2001, 47, 125–131. [Google Scholar] [CrossRef] [PubMed]
  111. Zhang, Z.; Yang, X.; Zhang, L.; Duan, Z.; Jia, L.; Wang, P.; Shi, Y.; Li, Y.; Gao, J. Decreased Expression and Activation of Stat3 in Severe Preeclampsia. J. Mol. Histol. 2015, 46, 205–219. [Google Scholar] [CrossRef]
  112. Zhang, Z.; Wang, X.; Wang, J.; Zhang, L. The Decreased Expression of Stat3 and P-Stat3 in Preeclampsia-like Rat Placenta. J. Mol. Histol. 2018, 49, 175–183. [Google Scholar] [CrossRef]
  113. Zhang, X.L.; Pan, J.; Wang, G.; Li, H.; Wang, J.L.; Jiang, X.M.; Zhao, X.Y. Effect of STAT3 Mediated Epigenetic Regulation in Pre-Eclampsia: An Analysis of Trial Data. Eur. Rev. Med. Pharmacol. Sci. 2020, 24, 5843–5855. [Google Scholar] [CrossRef]
  114. Kalisch-Smith, J.I.; Simmons, D.G.; Dickinson, H.; Moritz, K.M. Review: Sexual Dimorphism in the Formation, Function and Adaptation of the Placenta. Placenta 2017, 54, 10–16. [Google Scholar] [CrossRef] [PubMed]
  115. Kalisch-Smith, J.I.; Simmons, D.G.; Pantaleon, M.; Moritz, K.M. Sex Differences in Rat Placental Development: From Pre-Implantation to Late Gestation. Biol. Sex Differ. 2017, 8, 17. [Google Scholar] [CrossRef] [Green Version]
  116. Rosenfeld, C.S. Sex-Specific Placental Responses in Fetal Development. Endocrinology 2015, 156, 3422–3434. [Google Scholar] [CrossRef] [Green Version]
  117. Seidmann, L.; Kamyshanskiy, Y.; Martin, S.Z.; Fruth, A.; Roth, W. Immaturity for Gestational Age of Microvasculature and Placental Barrier in Term Placentas with High Weight. Eur. J. Obstet. Gynecol. Reprod. Biol. 2017, 215, 134–140. [Google Scholar] [CrossRef]
  118. Benirschke, K.; Kaufmann, P. Three-Dimensional Aspects of Villous Maldevelopment. In Pathology of the Human Placenta; Springer: New York, NY, USA, 1995; pp. 167–181. [Google Scholar]
  119. Fujii, T.; Nagamatsu, T.; Morita, K.; Schust, D.J.; Iriyama, T.; Komatsu, A.; Osuga, Y.; Fujii, T. Enhanced HIF2α Expression during Human Trophoblast Differentiation into Syncytiotrophoblast Suppresses Transcription of Placental Growth Factor. Sci. Rep. 2017, 7, 12455. [Google Scholar] [CrossRef] [Green Version]
  120. Roland, C.S.; Hu, J.; Ren, C.-E.; Chen, H.; Li, J.; Varvoutis, M.S.; Leaphart, L.W.; Byck, D.B.; Zhu, X.; Jiang, S.-W. Morphological Changes of Placental Syncytium and Their Implications for the Pathogenesis of Preeclampsia. Cell. Mol. Life Sci. 2016, 73, 365–376. [Google Scholar] [CrossRef] [Green Version]
  121. Workalemahu, T.; Ouidir, M.; Shrestha, D.; Wu, J.; Grantz, K.L.; Tekola-Ayele, F. Differential DNA Methylation in Placenta Associated With Maternal Blood Pressure During Pregnancy. Hypertension 2020, 75, 1117–1124. [Google Scholar] [CrossRef] [PubMed]
  122. Fox, R.; Kitt, J.; Leeson, P.; Aye, C.Y.L.; Lewandowski, A.J. Preeclampsia: Risk Factors, Diagnosis, Management, and the Cardiovascular Impact on the Offspring. J. Clin. Med. 2019, 8, 1625. [Google Scholar] [CrossRef] [Green Version]
  123. Girchenko, P.; Lahti, J.; Czamara, D.; Knight, A.K.; Jones, M.J.; Suarez, A.; Hämäläinen, E.; Kajantie, E.; Laivuori, H.; Villa, P.M.; et al. Associations between Maternal Risk Factors of Adverse Pregnancy and Birth Outcomes and the Offspring Epigenetic Clock of Gestational Age at Birth. Clin. Epigenet. 2017, 9, 49. [Google Scholar] [CrossRef] [Green Version]
  124. Mayne, B.T.; Leemaqz, S.Y.; Smith, A.K.; Breen, J.; Roberts, C.T.; Bianco-Miotto, T. Accelerated Placental Aging in Early Onset Preeclampsia Pregnancies Identified by DNA Methylation. Epigenomics 2016, 9, 279–289. [Google Scholar] [CrossRef] [Green Version]
  125. Zhang, X.; Du, L.; Qiao, Y.; Zhang, X.; Zheng, W.; Wu, Q.; Chen, Y.; Zhu, G.; Liu, Y.; Bian, Z.; et al. Ferroptosis Is Governed by Differential Regulation of Transcription in Liver Cancer. Redox Biol. 2019, 24, 101211. [Google Scholar] [CrossRef]
  126. Brunton, H.; Caligiuri, G.; Cunningham, R.; Upstill-Goddard, R.; Bailey, U.-M.; Garner, I.M.; Nourse, C.; Dreyer, S.; Jones, M.; Moran-Jones, K.; et al. HNF4A and GATA6 Loss Reveals Therapeutically Actionable Subtypes in Pancreatic Cancer. Cell Rep. 2020, 31, 107625. [Google Scholar] [CrossRef]
  127. Hsu, J.; Sage, J. Novel Functions for the Transcription Factor E2F4 in Development and Disease. Cell Cycle 2016, 15, 3183–3190. [Google Scholar] [CrossRef] [PubMed]
  128. Bhawe, K.; Roy, D. Interplay between NRF1, E2F4 and MYC Transcription Factors Regulating Common Target Genes Contributes to Cancer Development and Progression. Cell. Oncol. 2018, 41, 465–484. [Google Scholar] [CrossRef]
  129. Tseng, L.; Tang, M.; Wang, Z.; Mazella, J. Progesterone Receptor (HPR) Upregulates the Fibronectin Promoter Activity in Human Decidual Fibroblasts. DNA Cell Biol. 2003, 22, 633–640. [Google Scholar] [CrossRef]
  130. O’Rawe, J.A.; Wu, Y.; Dorfel, M.J.; Rope, A.F.; Au, P.Y.; Parboosingh, J.S.; Moon, S.; Kousi, M.; Kosma, K.; Smith, C.S.; et al. TAF1 Variants Are Associated with Dysmorphic Features, Intellectual Disability, and Neurological Manifestations. Am. J. Hum. Genet. 2015, 97, 922–932. [Google Scholar] [CrossRef] [Green Version]
  131. Cheng, H.; Capponi, S.; Wakeling, E.; Marchi, E.; Li, Q.; Zhao, M.; Weng, C.; Piatek, S.G.; Ahlfors, H.; Kleyner, R.; et al. Missense Variants in TAF1 and Developmental Phenotypes: Challenges of Determining Pathogenicity. Hum. Mutat. 2020, 41, 449–464. [Google Scholar] [CrossRef]
  132. Oh, H.R.; An, C.H.; Yoo, N.J.; Lee, S.H. Frameshift Mutations in the Mononucleotide Repeats of TAF1 and TAF1L Genes in Gastric and Colorectal Cancers with Regional Heterogeneity. Pathol. Oncol. Res. 2017, 23, 125–130. [Google Scholar] [CrossRef]
  133. Fuke, C.; Shimabukuro, M.; Petronis, A.; Sugimoto, J.; Oda, T.; Miura, K.; Miyazaki, T.; Ogura, C.; Okazaki, Y.; Jinno, Y. Age Related Changes in 5-Methylcytosine Content in Human Peripheral Leukocytes and Placentas: An HPLC-Based Study. Ann. Hum. Genet. 2004, 68, 196–204. [Google Scholar] [CrossRef]
  134. Novakovic, B.; Yuen, R.K.; Gordon, L.; Penaherrera, M.S.; Sharkey, A.; Moffett, A.; Craig, J.M.; Robinson, W.P.; Saffery, R. Evidence for Widespread Changes in Promoter Methylation Profile in Human Placenta in Response to Increasing Gestational Age and Environmental/Stochastic Factors. BMC Genom. 2011, 12, 529. [Google Scholar] [CrossRef] [Green Version]
  135. Maccani, J.; Maccani, M. Altered Placental DNA Methylation Patterns Associated with Maternal Smoking: Current Perspectives. Adv. Genom. Genet. 2015, 2015, 205. [Google Scholar] [CrossRef] [Green Version]
  136. Suter, M.; Ma, J.; Harris, A.; Patterson, L.; Brown, K.A.; Shope, C.; Showalter, L.; Abramovici, A.; Aagaard-Tillery, K.M. Maternal Tobacco Use Modestly Alters Correlated Epigenome-Wide Placental DNA Methylation and Gene Expression. Epigenetics 2011, 6, 1284–1294. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  137. Mikeska, T.; Craig, J.M. DNA Methylation Biomarkers: Cancer and Beyond. Genes 2014, 5, 821–864. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Toxics 09 00338 g001 550
Figure 1. Study setting in Scania, Sweden, showing the two recruitment hospitals.
Figure 1. Study setting in Scania, Sweden, showing the two recruitment hospitals.
Toxics 09 00338 g001
Toxics 09 00338 g002 550
Figure 2. Gestational age (GA) acceleration/deceleration in (A) cases and controls and in (B) the different combinations of PE status and NOx exposure group.
Figure 2. Gestational age (GA) acceleration/deceleration in (A) cases and controls and in (B) the different combinations of PE status and NOx exposure group.
Toxics 09 00338 g002
Table
Table 1. Description of study participants a and test of differences in characteristics between preeclampsia (PE) cases and controls b.
Table 1. Description of study participants a and test of differences in characteristics between preeclampsia (PE) cases and controls b.
CharacteristicsPE Cases, N = 29 (%)Controls, N = 82 (%)p-Value
Maternal age (years)31.2 ± 5.630.2 ± 3.90.30
<3520 (69.0)70 (85.4)0.053
35≤9 (31.0)12 (14.6)
Maternal pregestational body mass index (BMI) (kg/m2)
<18.50 (0)1 (1.2)0.062
18.5–24.911 (37.9)44 (53.7)
25–29.97 (24.1)25 (30.5)
30≤11 (37.9)12 (14.6)
Maternal parity
Nullipara17 (58.6)60 (73.2)0.26
Primipara10 (34.5)16 (19.5)
Multipara2 (6.9)6 (7.3)
Previous gestational hypertension0 (0)2 (2.4)0.40
Previous preeclampsia6 (20.7)6 (7.3)0.046
Mode of delivery
Vaginal22 (75.9)67 (81.7)0.50
Caesarean7 (24.1)15 (18.3)
Fetal sex
Male16 (55.2)43 (52.4)0.80
Female13 (44.8)39 (47.6)
Preterm birth (<week 37)4 (13.8)4 (4.9)0.11
Fetal weight (g)3485 ± 5623282 ± 753
Year of birth
20086 (20.7)15 (18.3)0.11
20098 (27.6)21 (25.6)
20102 (6.9)0
20112 (6.9)19 (23.2)
20145 (17.2)13 (15.9)
20156 (20.7)14 (17.1)
Gestational age by week38.5 ± 2.139.5 ± 1.70.014
a Mean and standard deviation (SD) are presented for continuous variables, and numbers (percentages) are presented for categorical variables. b Differences between PE cases and controls were analysed with an χ2 test for categorical variables and with a t-test for continuous variables.
Table
Table 2. Differentially methylated positions (ranked by q-value) between comparisons of combinations of preeclampsia (PE) status and ambient nitrogen oxides (NOx) exposure group.
Table 2. Differentially methylated positions (ranked by q-value) between comparisons of combinations of preeclampsia (PE) status and ambient nitrogen oxides (NOx) exposure group.
CpGChrGeneGene Name2log FC
(95% CI) a		Beta bq-Value
Controls with low NOx exposure (reference) vs. PE cases with high NOx exposure
cg275967795NAc 0.78 (0.72, 0.86)0.640.028
cg266720982NA 0.75 (0.68, 0.83)0.340.028
cg078063613CACNA2D3Calcium voltage-gated channel auxiliary subunit alpha 2 delta 30.62 (0.53, 0.73)0.720.031
cg2393522020RPRD1BRegulation of nuclear pre-mRNA domain containing 1B0.77 (0.70, 0.85)0.180.034
cg056943311SH3D21SH3 domain-containing protein 210.78 (0.71, 0.85)0.680.034
cg1599935611YAP1Yes1-associated transcriptional regulator0.77 (0.70, 0,85)0.660.036
Controls with high NOx exposure (reference) vs. PE cases with high NOx exposure
cg172836202HAAO3-hydroxyanthranilate 3,4-dioxygenase0.1 (0.05, 0.21)0.140.005
cg266720982NA 0.74 (0.67, 0.81)0.340.006
cg2483245715NA 0.73 (0.66, 0.81)0.560.010
cg1553446112PEBP1phosphatidylethanolamine binding protein 10.70 (0.62, 0.8)0.690.021
cg1853135115NA 0.77 (0.7, 0.84)0.610.021
cg0240473911FAM111BFAM111 trypsin-like peptidase B2.58 (1.85, 3.63)0.160.021
cg077153793PLCXD2phosphatidylinositol specific phospholipase C X domain containing 20.76 (0.68, 0.84)0.50.042
cg1823527412CD163CD163 molecule, also known as M1300.68 (0.59, 0.78)0.830.042
cg235022956MICAL1microtubule associated monooxygenase, calponin, and LIM domain containing 10.77 (0.70, 0.85)0.830.042
cg1599935611YAP1Yes1 associated transcriptional regulator0.77 (0.74, 0.87)0.660.042
cg072495173MB21D2Mab-21 domain containing 2 also known as:C3orf590.80 (0.74, 0.87)0.560.047
cg084935906SUPT3HSPT3 homolog, SAGA, and STAGA complex component0.74 (0.66, 0.83)0.660.047
cg0217181417SERPINF2serpin family F member 20.74 (0.66, 0.82)0.610.047
Controls with high NOx exposure (reference) vs. PE cases with low NOx exposure
cg189047844SCLT1sodium channel and clathrin linker 10.54 (0.45, 0.67)0.020.031
Female pregnancies only—Controls with high NOx exposure (reference) vs. PE cases with high NOx exposure
cg1139840017RTN4RL1reticulon 4 receptor like 12.71 (2, 3.66)0.950.011
cg0923498313LATS2large tumor suppressor kinase 20.59 (0.5, 0.69)0.390.011
cg223434761NA 0.46 (0.36, 0.59)0.240.011
cg0856855011C11orf16chromosome 11 open reading frame 160.65 (0.56, 0.74)0.570.011
cg244744097ACTR3Cactin related protein 3C0.47 (0.37, 0.59)0.740.011
cg177649463TRAK1trafficking kinesin protein 10.46 (0.36, 0.59)0.880.011
cg0421954417KRT24keratin 240.57 (0.48, 0.68)0.30.011
cg2455145917RAB5Cmember RAS oncogene family0.46 (0.35, 0.59)0.740.011
cg146202349NA 0.67 (0.59, 0.76)0.310.015
cg1117324611TSPAN4tetraspanin 40.44 (0.33, 0.57)0.30.015
cg0815719422SLC25A17Solute Carrier Family 25 Member 170.6 (0.51, 0.71)0.550.016
cg0457403417PRACPRAC small nuclear protein2.88 (2.03, 4.08)0.80.016
cg248740908NA 0.62 (0.53, 0.73)0.640.019
Female pregnancies only—Controls with high NOx exposure (reference) vs. PE cases with low NOx exposure
cg161629301NA 0.10 (0.05, 0.20)0.570.015
Male pregnancies only—Controls with high NOx exposure (reference) vs. PE cases with low NOx exposure
cg0980443917STAT3Signal transducer and activator of transcription 30.15 (0.09, 0.26)0.020.004
Abbreviations: Chr, chromosome; CI, confidence interval; q-value, false discovery rate (FDR)—adjusted p-value using the Benjamini–Hochberg method; FC, fold change. a 2logFC, binary logarithmic fold change. LogFC denotes β1 from the following robust regression model: M-value = β1 × PE/NOx group + β2 × DNA concentration + β3 × pregestational BMI + β4 × gestational age + β5 × fetal sex + β6 × estimated fraction Hofbauer cells + β7 × estimated fraction Syncytiotrophoblast cells. For female placentas, the models look as follows: M-value = β1 × PE/NOx group + β2 × DNA concentration + β3 × gestational age + β4 × estimated fraction Hofbauer cells + β5 × estimated fraction Syncytiotrophoblast cells. b Average methylation state, expressed as Beta-value, for all study participants, ranging from 0 to 1 (1 means fully methylated). c NA, not annotated, i.e., the CpG is not present in any known gene.
Table
Table 3. Differentially expressed genes (ranked by q-value) between comparisons of combinations of PE status and NOx group.
Table 3. Differentially expressed genes (ranked by q-value) between comparisons of combinations of PE status and NOx group.
GeneGene Name2logFC
(95% CI) a		p-Value
(Unadjusted)		q-Value
ESCO2establishment of sister chromatid cohesion N-acetyltransferase 22.4 (1.7, 3.3)0.000020.33
MUC20P1mucin 20, cell surface associated pseudogene 119.0 (6.1, 59.4)0.000030.33
LINC01003long intergenic non-protein coding RNA 10030.58 (0.47, 0.73)0.00010.64
KCNA4potassium voltage-gated channel subfamily A member 43.6 (1.9, 6.7)0.00040.94
TTC3tetratricopeptide repeat domain 31.4 (1.2, 1.6) 0.00050.94
ABCC13ATP binding cassette subfamily C member 13 (pseudogene)0.36 (0.21, 0.60)0.00060.94
GNAI1G protein subunit alpha i10.72 (0.61, 0.86)0.00090.94
SRRM3serine/arginine repetitive matrix 30.44 (0.27, 0.69)0.00130.94
AL451007.2 0.29 (0.15, 0.58)0.00140.94
PCDHB18Pprotocadherin beta 18 pseudogene 3.3 (1.7, 6.7)0.00160.94
a 2logFC, binary logarithmic fold change. LogFC denotes β1 from the following robust regression model: TMM = β1 × PE status + β2 × Sex.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Engström, K.; Mandakh, Y.; Garmire, L.; Masoumi, Z.; Isaxon, C.; Malmqvist, E.; Erlandsson, L.; Hansson, S.R. Early Pregnancy Exposure to Ambient Air Pollution among Late-Onset Preeclamptic Cases Is Associated with Placental DNA Hypomethylation of Specific Genes and Slower Placental Maturation. Toxics 2021, 9, 338. https://doi.org/10.3390/toxics9120338

AMA Style

Engström K, Mandakh Y, Garmire L, Masoumi Z, Isaxon C, Malmqvist E, Erlandsson L, Hansson SR. Early Pregnancy Exposure to Ambient Air Pollution among Late-Onset Preeclamptic Cases Is Associated with Placental DNA Hypomethylation of Specific Genes and Slower Placental Maturation. Toxics. 2021; 9(12):338. https://doi.org/10.3390/toxics9120338

Chicago/Turabian Style

Engström, Karin, Yumjirmaa Mandakh, Lana Garmire, Zahra Masoumi, Christina Isaxon, Ebba Malmqvist, Lena Erlandsson, and Stefan R. Hansson. 2021. "Early Pregnancy Exposure to Ambient Air Pollution among Late-Onset Preeclamptic Cases Is Associated with Placental DNA Hypomethylation of Specific Genes and Slower Placental Maturation" Toxics 9, no. 12: 338. https://doi.org/10.3390/toxics9120338

APA Style

Engström, K., Mandakh, Y., Garmire, L., Masoumi, Z., Isaxon, C., Malmqvist, E., Erlandsson, L., & Hansson, S. R. (2021). Early Pregnancy Exposure to Ambient Air Pollution among Late-Onset Preeclamptic Cases Is Associated with Placental DNA Hypomethylation of Specific Genes and Slower Placental Maturation. Toxics, 9(12), 338. https://doi.org/10.3390/toxics9120338

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