Front. Toxicol. Frontiers in Toxicology Front. Toxicol. 2673-3080 Frontiers Media S.A. 10.3389/ftox.2021.659407 Toxicology Original Research Association of Prenatal Ambient Air Pollution Exposure With Placental Mitochondrial DNA Copy Number, Telomere Length and Preeclampsia Mandakh Yumjirmaa 1 * Oudin Anna 1 2 Erlandsson Lena 3 Isaxon Christina 4 Hansson Stefan R. 3 5 Broberg Karin 6 † Malmqvist Ebba 1 † 1Environment Society and Health, Division of Occupational and Environmental Medicine, Department of Laboratory Medicine, Lund University, Lund, Sweden 2Section of Sustainable Health, Department of Public Health and Clinical Medicine, Umeå University, Umeå, Sweden 3Division of Obstetrics and Gynecology, Department of Clinical Sciences Lund, Lund University, Lund, Sweden 4Ergonomics and Aerosol Technology, Department of Design Sciences, Lund University, Lund, Sweden 5Department of Obstetrics & Gynaecology, Skåne University Hospital, Malmö, Sweden 6Division of Occupational and Environmental Medicine, Department of Laboratory Medicine, Lund University, Lund, Sweden

Edited by: Phoebe Stapleton, Rutgers, The State University of New Jersey, United States

Reviewed by: Colette Miller, United States Environmental Protection Agency (EPA), United States; Sabbir Khan, University of Texas MD Anderson Cancer Center, United States

 *Correspondence: Yumjirmaa Mandakh yumjirmaa.mandakh@med.lu.se

This article was submitted to Developmental and Reproductive Toxicology, a section of the journal Frontiers in Toxicology

†These authors have contributed equally to this work and share last authorship

 26 05 2021 2021 3 659407 27 01 2021 19 04 2021 Copyright © 2021 Mandakh, Oudin, Erlandsson, Isaxon, Hansson, Broberg and Malmqvist. 2021 Mandakh, Oudin, Erlandsson, Isaxon, Hansson, Broberg and Malmqvist

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

Background: Studies have shown that ambient air pollution is linked to preeclampsia (PE), possibly via generation of oxidative stress in the placenta. Telomere length and mitochondrial DNA copy number (mtDNAcn) are sensitive to oxidative stress damage.

Objective: To study the association between prenatal exposure to ambient nitrogen oxides (NOx, a marker for traffic-related air pollution), and PE, as well as potential mediation effects by placental telomere length and mtDNAcn.

Methods: This is a cross-sectional study of 42 preeclamptic and 95 arbitrarily selected normotensive pregnant women with gestational ambient NOx exposure assessment in southern Scania, Sweden. Hourly concentrations of NOx were estimated at the residential addresses by a Gaussian-plume dispersion model with 100 × 100 m spatial resolutions and aggregated into trimester-specific mean concentrations. Placental relative mtDNAcn and telomere length were measured using qPCR. Linear and logistic regression models were used to investigate associations, adjusted for perinatal and seasonal characteristics.

Results: Exposure was categorized into low and high exposures by median cut-offs during first [11.9 μg/m3; interquartile range (IQR) 7.9, 17.9], second (11.6 μg/m3; IQR: 7.1, 21.1), third trimesters (11.9 μg/m3; IQR: 7.7, 19.5) and entire pregnancy (12.0 μg/m3; IQR: 7.6, 20.1). Increased risk of PE was found for high prenatal NOx exposure during the first trimester (OR 4.0; 95% CI: 1.4, 11.1; p = 0.008), and entire pregnancy (OR 3.7; 95% CI: 1.3, 10.4; p = 0.012). High exposed group during the first trimester had lower placental relative mtDNAcn compared with low exposed group (−0.20; 95% CI: −0.36, −0.04; p = 0.01). Changes in relative mtDNAcn did not mediate the association between prenatal NOx exposure and PE. No statistically significant association was found between placental relative telomere length, prenatal NOx exposure and PE.

Conclusion: In this region with relatively low levels of air pollution, ambient NOx exposure during the first trimester was associated with reduced placental relative mtDNAcn and an increased risk of PE. However, we did not find any evidence that mtDNAcn or TL mediated the association between air pollution and PE. Future research should further investigate the role of mtDNAcn for pregnancy complications in relation to exposure to ambient air pollution during pregnancy.

 preeclampsia mitochondrial DNA copy number telomere length nitrogen oxides placenta ambient air pollution Forskningsrådet om Hälsa, Arbetsliv och Välfärd10.13039/501100006636 Introduction

Ambient air pollution is a major global contributor to morbidity and mortality globally, accounting for 7.6% of all deaths worldwide in 2015 (Cohen et al., 2017). Exposure to ambient air pollution during pregnancy is particularly worrisome due to the health threat posed to both pregnant women and their fetuses. Several epidemiologic reviews have highlighted the possible association between maternal exposure to ambient air pollutants and pregnancy complications, such as gestational diabetes mellitus (Zhang et al., 2020), hypertensive disorders of pregnancy (Bai et al., 2020; Yu et al., 2020) and fetal growth restriction (FGR) (Fu et al., 2019). According to the Barker hypothesis, different environmental conditions present already at the time of conception may pose an increased risk of developing diseases later in life (Barker, 2004). Apparently, FGR associated with prenatal exposure to ambient air pollution appears to be a predisposing factor for impaired physical growth in the early years after birth (Clemente et al., 2017; Fossati et al., 2020).

The pathophysiological mechanisms linking ambient air pollution with pregnancy complications are not completely understood. One reason is the complex composition of air pollution containing a mixture of components, such as particulate matter (PM) of different sizes and compositions, ozone, nitrogen oxides (NOx), carbon monoxide and sulfur dioxide. Another reason is that pregnancy complications, especially preeclampsia (PE), tend to be caused by a dysfunctional placenta. According to the recommendation of International Society for the Study of Hypertension in Pregnancy (ISSHP), PE is defined as de-novo hypertension combined with maternal organ dysfunction and/or utero-placental dysfunction including FGR presenting after 20 weeks of gestation (Brown et al., 2018). More specifically, PE is more of a syndrome with different clinical manifestations than a well-defined disease (Gyselaers, 2020). It is classified into early (occurring before 34 gestational weeks) and late (occurring after 34 weeks) onset PE (Tranquilli et al., 2013).

The etiology and pathophysiology of early- and late-onset PE are still not fully understood, however, placental dysfunction due to syncytiotrophoblast stress may explain the underlying pathophysiologic mechanism of PE syndrome (Redman et al., 2021). Syncytiotrophoblast is the outer cell layer on placental villae in direct contact with maternal blood (Manna et al., 2019) and is formed by underlying cytotrophoblast fusion and senescence (Redman et al., 2021). The key characteristic of trophoblast senescence is a telomere shortening (Burton et al., 2017) and trophoblast fusion dependent on mitochondrial dynamics (Holland et al., 2017). Syncytiotrophoblast stress is suggested to be caused by insufficient placental perfusion and hypoxia due to incomplete placentation in early-onset PE or due to placental growth and maternal constitutional factors such as obesity, autoimmune disease and infection in late-onset PE (Redman et al., 2021).

Exposure to ambient air pollution is potentially associated with increased oxidative stress damage to nuclear and mitochondrial DNA (mtDNA), as well as lipid peroxidation, resulting in DNA mutations and cell death (Møller and Loft, 2010). Previous studies have reported a dose-dependent increase of a DNA adduct 8-hydroxy-deoxyguanosine in DNA isolated from lung, venous blood and urine from laboratory animals and humans exposed to diesel exhaust gases, such as NOx and polycyclic aromatic hydrocarbons, and/or diesel exhaust particles, both the fine and ultrafine PM (Ichinose et al., 1997; Harri et al., 2005; Vinzents et al., 2005; Lee et al., 2012). Particularly, polycyclic aromatic hydrocarbons on ambient PM has been indicated to induce upregulation of cytochrome P450, which generates reactive oxygen species (ROS) (Bonvallot et al., 2001). Furthermore, recent evidence has suggested that elevated levels of malondialdehyde, a lipid peroxidation product of syncytiotrophoblast membranes, are linked to both short-term exposure to ambient air pollution (Li et al., 2020) and PE (Ramiro-Cortijo et al., 2020). Taken together, air pollution-induced oxidative stress may cause placental dysfunction (Saenen et al., 2019). Interestingly, telomeres and mitochondria are both considered to be important targets of oxidative stress induced by environmental pollutants (Martens and Nawrot, 2016).

Telomeres are specialized structures located at the end of each chromosome. In human somatic cells, telomeres gradually shorten during cell division until they reach a critical point that triggers replicative cell senescence or apoptosis (Hayflick and Moorhead, 1961; Victorelli and Passos, 2017; Pańczyszyn et al., 2020). The high guanine content in the telomeres render them susceptible to ROS and accelerated telomere length shortening has been found to be associated with oxidative stress (Richter and Zglinicki, 2007; Pańczyszyn et al., 2020). A growing body of literature suggests that exposure to environmental pollutants during pregnancy is associated with changes in placental telomere length (Lin et al., 2013; Bijnens et al., 2015; Herlin et al., 2019). Moreover, it has been indicated that placentas from pregnancies complicated with PE and FGR have shorter telomere length as compared with uncomplicated pregnancies (Biron-Shental et al., 2010).

Mitochondria are cell organelles that are responsible for ATP production and several metabolic pathways, such as the heme synthesis. Mitochondria in placental cytotrophoblast and syncytiotrophoblast have different structure and functions (Bustamante et al., 2014). Due to its location in the mitochondrial matrix and its structure without histones, mtDNA are vulnerable to oxidative stress damage induced by both environmental pollutants and endogenous ROS production formed in the electron transport chain of the mitochondrial matrix (Burton et al., 2017; Roubicek and Souza-Pinto, 2017). Several studies have shown that maternal exposure to ambient air pollutants is associated with reduced placental mtDNA copy number (mtDNAcn) (Janssen et al., 2012, 2015; Clemente et al., 2017). In contrast, placentas from early-onset PE had increased mtDNAcn accompanied with increased mitochondrial mass when compared to normal pregnancies, whereas there was no difference between mtDNAcn in late-onset PE and control placentas (Vishnyakova et al., 2016). Moreover, placentas from pregnancies complicated by PE have shown mitochondrial swelling and broken cristae along with excess endogenous production of ROS (Muralimanoharan et al., 2012).

Here, we hypothesize that exposure to ambient air pollution during pregnancy induces placental dysfunction seen in PE, measured as shorter telomere length and lower mtDNAcn. The aim of this study was therefore to explore the association between prenatal exposure to ambient NOx and PE, as well as potential mediation effects of placental telomere length and mtDNAcn.

 Materials and Methods Study Design and Population

This cross-sectional study used biobank data of placental biopsies collected from Cesarean and vaginal births between 2008 and 2015 from normotensive and PE pregnancies at Skåne University Hospital in Lund and Malmö, Sweden (Figure 1). The study protocol of placental sample collection (Dnr 243/2005) and air pollution exposure assessment at geocoded residential locations of participants (Dnr 696/2014) were approved by Lund University Ethical Review Board, LU803-2 (2015/14) and a written informed consent was obtained from all participants.

Map of study setting in Scania, Sweden.

PE is recorded in the Swedish Medical Birth Register in accordance with the Swedish adaptation of the 10th version of the International Statistical Classification of Diseases and Related Health Problems (ICD-10) by the World Health Organization (WHO) (Socialstyrelsen, 2018). We selected all available PE cases and chose random controls arbitrarily for the present study. Using the unique personal identification numbers in Sweden, the biobank consists of a compiled database of medical information regarding prenatal visits, delivery and perinatal outcomes. As PE or no PE was the prerequisite for collection of placental samples, this data was also inserted and double-checked with data from prenatal visits.

A flow chart of this study is shown in Figure 2. A total of 169 frozen placental tissue samples were available in the biobank for retrospective NOx exposure assessment during pregnancy of the participants. It was, however, not possible to estimate the exposure for eight participants due to their residential location outside the territory of our dispersion model and they were thus excluded. In addition, we excluded participants who had twins or smoked during pregnancy. Two participants with missing covariate information were also excluded. The remaining 137 participants were included for further analysis.

Flowchart of the study.

 Air Pollution Exposure

Using the Swedish personal identification numbers and linking them to the regional population registry, we obtained geocoded residential locations of each participant in Scania, Sweden. Ambient concentrations of NOx during pregnancy were estimated at outdoor ground level at their geocoded residential locations using EnviMan software (Opsis, 2006). The Gaussian-plume dispersion model, AERMOD developed by the United States Environmental Protection Agency (USEPA, 1998), was modified and incorporated into EnviMan. The dispersion model used in this study is explained in detail elsewhere (Rittner et al., 2020). Briefly, EnviMan relies on the emission database created by the local Air Quality Board and maintained by the Environmental Department at the City of Malmö. Sources in the emission database are road traffic, shipping, aviation, rail transport, industry, power plants, small-scale heating, machinery, as well as local emissions from Zealand, Denmark.

EnviMan calculates hourly-average concentrations of NOx at a spatial resolution of 100 × 100 m from an emission rate and meteorological conditions including the wind speed, wind direction, temperature and solar radiation (Stroh et al., 2007). Our Gaussian dispersion model was only modeled for southwestern Scania. In a similar setting, Stroh et al. (2012) have found the modeled residential NO2 levels were strongly correlated (rS = 0.8, p > 0.001) with measured façade NO2 levels. According to our previous studies, a total NOx concentrations were used in the data analysis by adding a background concentration of 2.5 μg/m3 on the modeled hourly NOx concentrations to account for transboundary air pollution (Malmqvist et al., 2013; Mandakh et al., 2020). We averaged hourly NOx concentrations to monthly mean NOx concentrations for each pregnant women's pregnancy period. The mean of at least 2 monthly average concentrations of NOx were used as the corresponding trimester mean NOx concentrations. Furthermore, third trimester mean NOx concentrations were estimated based on their gestational age of delivery. In each trimester there were different number of pregnant women with complete data on exposure. In the first trimester there were 129, in the second and third trimesters there were 124 out of 137 total cases. One hundred seventeen women had complete data on exposure for all three trimesters. As our study sample size was relatively small, we used a standard method of dichotomizing a continuous variable by median to categorize statistically the trimester-specific NOx concentrations to low and high exposure groups.

This study was limited to a single-pollutant analysis of ambient NOx although our dispersion model had been used in the similar setting to model other ambient air pollutants, such as PM2.5 (particulate matter with an aerodynamic diameter smaller than 2.5 μm), PM10 and black carbon as reported in our population-based study (Mandakh et al., 2020). The participants given birth prior to 2009 in the present study were overlapping with Mandakh et al. (2020) and had data on locally generated PM2.5 and black carbon. We could not model these pollutants for those given birth after 2009 due to some missing data of PM2.5 measurement at local monitoring stations and economic constraint during the retrospective air pollution exposure assessment. As previously reported, modeled ambient NOx was highly and positively correlated with modeled local PM2.5 (r = 0.90, p < 0.01) (Mandakh et al., 2020).

 Biological Sample Collection

Following a standardized protocol as recommended by Burton et al. (2014), biopsies (1 × 1 × 1 cm) from a central part of the villi region of the placenta (within a 7 cm radius of the umbilical cord) were collected within 30 min after delivery. Samples were rinsed in a sodium chloride solution, dried, and then frozen on dry ice and stored at −80°C until further analysis.

 DNA Extraction

Total DNA was extracted from the frozen placental biopsies kept on dry ice using Qiagen AllPrep DNA/RNA/Protein Mini Kit (Qiagen, Hilden, Germany) according to manufacturer's instruction, and thereafter stored at −80°C until further use. The DNA samples were thawed only once in preparation for the present study. The purity ratios [A260/280 1.87 (IQR 1.86, 1.88) and A260/230 2.05 (IQR 1.73, 2.20)] and DNA yield (ng/μl) were determined by using a NanoDrop Spectrophotometer ND-1,000 (NanoDrop technologies, Wilmington, USA).

 Quantification of Mitochondrial DNA Copy Number and Telomere Length

To quantify relative mtDNAcn and telomere length, real-time polymerase chain reactions (qPCR) were performed using an Applied Biosystems 7900HT Fast Real-Time PCR System (Applied Biosystems, Foster City, CA, USA). The PCR products of mtDNAcn and telomere length were related to PCR product of a single copy gene [hemoglobin beta (HBB)]. MT-TL1 encoding the Mitochondrially Encoded TRNA-Leu (UUA/G) 1 was the target gene that was used to determine mtDNAcn (Venegas and Halberg, 2012). Master mixes for mtDNAcn quantification was prepared with PerfeCTa SYBR Green FastMix PCR reagent (ROX) (2x) (QuantaBio, Beverly, MA, USA), 10 μM mtDNA forward primer tRNA F3212 (5′-CAC CCA AGA ACA GGG TTT GT-3′), and 10 μM reverse primer tRNA R3319 (5′-TGG CCA TGG GTA TGT TGT TA-3′) (Custom DNA Oligos Synthesis Services, ThermoFisher Scientific). HBB master mix consisted of Fast SYBR Green Master Mix (2x) (Applied Biosystems, Waltham, USA), 10 μM HBB forward primer (5′-TGT GCT GGC CCA TCA CTT TG-3′), and 10 μM reverse primer (5′-ACC AGC CAC CAC TTT CTG ATA GG-3′) (Xu et al., 2018). Master mixes for telomere length runs were prepared with PCR buffer (10x) (Invitrogen, Carlsbad, CA, USA), 50 mM MgCl2, 10 mM dNTPs, 20 μM TEL1 primer (5'-CGG TTT GTT TGG GTT TGG GTT TGG GTT TGG GTT TGG GTT-3'), 20 μM TEL2 primer (5'-GGC TTG CCT TAC CCT TAC CCT TAC CCT TAC CCT TAC CCT-3'), 10x SYBR Green I (Invitrogen, Carlsbad, CA, USA) (50x), ROX (1x) (Invitrogen, Carlsbad, CA, USA), and 0.5 U Taq Platina DNA polymerase (Invitrogen, Carlsbad, CA, USA) (Herlin et al., 2019). Each reaction consisted of 2.5 μl DNA (4 ng/μl for mtDNA and HBB, 2 ng/μl for telomere length) and 7.5 μl master mix, resulting in total volume of 10 μl.

The thermal cycling profile of mtDNAcn PCR started with 3 min at 95°C followed by 35 cycles of 15 s of denaturation at 95°C and 1 min of annealing/extension at 60°C. HBB PCR condition was 20 s at 95°C, followed by 40 cycles of 1 s at 95 °C and 20 s at 60°C (Xu et al., 2018). PCR conditions for telomere length were 3 min at 95°C followed by 25 cycles of 15 s at 95°C and 1 min at 56°C (Herlin et al., 2019). The quantification of mtDNAcn, telomere length, HBB gene copy number and the standard curves were run in triplicates at separate runs. To generate a standard curve ranging from 0.25 to 16 ng/μl, a series of two-fold dilutions were carried out for three individual reference DNA samples. R2 of each standard curve was >0.99. Relative quantification of the concentration of mtDNAcn was determined by standard curve analysis and the ratio between mtDNAcn and HBB gene copy number (Xu et al., 2018). Similarly, the relative telomere length was the quotient of telomere length and HBB gene copy number (Herlin et al., 2019).

 Statistical Analysis

All statistical analyses were performed using the IBM SPSS Statistics software version 25 (IBM, Chicago, IL, USA). Continuous data was checked for normality and presented as means ± SD or as median with interquartile range (IQR) when data were not normally distributed. Categorical data are presented as numbers and percentage. A p-value of <0.05 was considered statistically significant. After NOx exposure assessment, 20 subjects were missing one (n = 7) or two (n = 13) trimester-specific NOx exposure data although they were term pregnancies. To account for a cumulative effect of exposure during the entire period of pregnancy, complete case analysis (n = 117) on complete exposure data at all trimesters was performed and reported.

For analysis of associations between prenatal NOx exposure and risk of PE, adjusted odds ratios (OR) with their 95% confidence intervals (95% CI) were estimated using logistic regression, after controlling for maternal age, pregestational body mass index (BMI), parity, gestational age and fetal sex. These covariates were defined based on a priori knowledge (Pedersen et al., 2014; Clemente et al., 2017; Martens et al., 2017; Iodice et al., 2018; Mandakh et al., 2020; Melody et al., 2020). The season of birth was also chosen as a covariate to control due to the seasonal variations in both PE prevalence (Wacker et al., 1998; Phillips et al., 2004; Rudra and Williams, 2005) and NOx concentrations (Pedersen et al., 2014; Choi et al., 2015; Silva et al., 2020). Linear regression was used to identify how prenatal NOx exposure was associated with placental relative mtDNAcn and telomere length, respectively, after adjusting for the same set of confounders adjusted for the logistic regression model mentioned above.

The standard Baron and Kenny (1986) method of mediation analysis was employed for investigating whether placental relative mtDNAcn or telomere length mediated the association between prenatal NOx exposure and PE. The Baron and Kenny method is a three step method where in the present study: (1) association between prenatal NOx exposure and PE was investigated; (2) association between prenatal NOx exposure and placental relative mtDNAcn or telomere length was investigated; and (3) if the association between prenatal NOx exposure and placental relative mtDNAcn or telomere length was statistically significant, mtDNAcn or telomere length were entered in to the regression model of the first step (i.e., prenatal NOx exposure in association with PE). In all models, the same set of confounding factors were included into the models. To be considered a mediating variable in the third step: (1) placental relative mtDNAcn or telomere length should be statistically significant; and (2) the coefficient for prenatal NOx exposure should be substantially diminished as compared to in the first step.

 Results Characteristics

A total of 42 PE cases consisted of three early-onset PE and 39 late-onset PE. Table 1 shows characteristics of the study participants stratified by PE status. Preeclamptic women were more likely to be older than 35 years, primiparous, obese and had PE or gestational hypertension in their previous pregnancies, as compared to normotensive pregnant women. Also, they were more likely to give birth during the autumn. Aspirin and antihypertensive drug use during pregnancy was more prevalent among preeclamptic women compared with normotensive women, since aspirin had been routinely prescribed for high-risk pregnancies to reduce the risk of PE. There were no differences in placental relative telomere length and mtDNAcn between PE and normotensive pregnant women.

Characteristics of study participants stratified by preeclampsia status.

Characteristics Total (n = 137) PE (n = 42) Normotensive (n = 95) p-value
Maternal age in years [mean (SD)] 30.1 ± 4.7 30.5 ± 5.9 29.97 ± 4.1
  <35 112 (81.8) 31 (73.8) 81 (85.3) 0.110
 ≥35 25 (18.2) 11 (26.2) 14 (14.7)
Parity [n (%)]
  Nullipara 92 (67.2) 24 (57.1) 68 (71.6) 0.228
  Primipara 34 (24.8) 13 (31.0) 21 (22.1)
  Multipara 11 (8.0) 5 (11.9) 6 (6.3)
Pregestational BMI [kg/m3] 25.7 ± 5.3 27.6 ± 6.3 24.9 ± 4.5
  <18.5 3 (2.2) 1 (2.4) 2 (2.1) 0.008*
  18.5–24.9 66 (48.2) 16 (38.1) 50 (52.6)
  25–29.9 40 (29.2) 9 (21.4) 31 (32.6)
 ≥30 28 (20.4) 16 (38.1) 12 (12.6)
Fetal sex [n (%)]
  Male 73 (53.3) 23 (54.8) 50 (52.6) 0.818
  Female 64 (46.7) 19 (45.2) 45 (47.4)
Gestational age, days [mean (SD)] 274.2 ± 15.8 269.3 ± 14.6 276.3 ±15.9
Birth weight [g] 3455.1 ± 702.2 3323.1 ± 748.9 3513.4 ± 676.4
Previous PE [n (%)] 16 (11.7) 8 (19.0) 8 (8.4) 0.074
Previous gestational HTNa [n (%)] 7 (5.1) 3 (7.1) 4 (4.2) 0.472
Medication of importance [n (%)]
  Aspirin 6 (4.4) 5 (11.9) 1 (1.1) 0.01*b
  Antihypertensive drugs 19 (13.9) 16 (38.1) 3 (3.2) <0.001*
Mode of delivery [n (%)]
  Vaginal 106 (79.1) 31 (75.6) 75 (80.6) 0.515
  Cesarean 28 (20.9) 10 (24.4) 18 (19.4)
Municipality [n (%)]
  Malmö 35 (25.9) 8 (20.0) 27 (28.4) 0.505
  Lund 46 (34.1) 15 (37.5) 31 (32.6)
  Other 54 (40.0) 17 (42.5) 37 (38.9)
Season of birth [n (%)]
  Winter 18 (13.1) 6 (14.3) 12 (12.6) 0.055
  Spring 26 (19.0) 4 (9.5) 22 (23.2)
  Summer 39 (28.5) 9 (21.4) 30 (31.6)
  Autumn 54 (39.4) 23 (54.8) 31 (32.6)
Placental mtDNAcn
  Mean ± SD 0.88 ± 0.39 0.87 ± 0.40 0.88 ± 0.39
  Median (IQR) 0.82 (0.57–1.08) 0.82 (0.49–1.14) 0.83 (0.59–1.05)
Placental telomere length
  Mean ± SD 1.39 ± 0.36 1.36 ± 0.33 1.41 ± 0.38
  Median (IQR) 1.36 (1.16–1.62) 1.35 (1.12–1.58) 1.37 (1.16–1.63)

Pearson Chi-Square p < 0.05;

HTN, Hypertension;

Fisher's Exact test.

The mean NOx concentrations varied between 14.0 and 16.2 μg/m3 over the three trimesters and entire period of pregnancy (Table 2). As shown in Figure 3, the trimester-specific and total NOx concentrations had moderate (rp = 0.65, p < 0.001) to high (rp = 0.90, p < 0.001) correlations that were statistically significant. The range in ambient NOx was wide, modeled NOx concentrations were split to low and high exposure groups by median cut-offs during first (11.9 μg/m3; IQR: 7.9, 17.9), second (11.6 μg/m3; IQR: 7.1, 21.1), third trimesters (11.9 μg/m3; IQR: 7.7, 19.5), and entire period of pregnancy (12.0 μg/m3; IQR: 7.6, 20.1) (Table 2). The characteristics of study population stratified by ambient NOx exposure group is presented in the Supplementary Table 1. Those exposed to high NOx concentrations during pregnancy were more likely to be primipara, overweight and to give birth in the spring and summer.

Descriptive statistics of prenatal NOx exposure (μg/m3).

Time window N Mean ± SD Min P25 P50 P75 Max
Trimester 1 129 14.0 ± 7.6 5.1 7.9 11.9 17.9 49.5
Trimester 2 124 16.2 ± 12.4 4.4 7.1 11.6 21.1 62.5
Trimester 3 124 14.7 ± 8.6 4.3 7.7 11.9 19.5 41.9
Entire pregnancy 137 14.8 ± 8.4 4.6 7.6 12.0 20.1 38.9

Correlation plot of NOx concentrations at all trimesters and entire period of pregnancy.

 Prenatal NO<sub>x</sub> Exposure and PE

Independent of maternal age, pregestational BMI, parity, season of birth and fetal sex, statistically significant associations were found between prenatal NOx exposure and PE when looking at exposure during first trimester (OR 4.0; 95% CI: 1.4, 11.1; p = 0.008) and entire pregnancy (OR 3.7; 95% CI: 1.3, 10.4; p = 0.012) (Table 3). Effect estimates associated with NOx exposures during the second and third trimesters were decreased, but statistically not significant.

Odds ratios (OR) and their 95% Confidence Intervals (CI) for PE in association with high NOx exposure as compared with low NOx exposure (median cut-off).

Time window Exposure group B Adjusted OR§ 95% CI for AOR p
Lower Upper
First trimester Low Reference
High 1.38 4.0 1.4 11.1 0.008
Second trimester Low Reference
High 1.02 2.8 0.97 8.0 0.056
Third trimester Low Reference
High 0.91 2.5 0.93 6.7 0.071
Entire pregnancy Low Reference
High 1.31 3.7 1.3 10.4 0.012

Adjusted for maternal age, pregestational BMI, parity, gestational age, season of birth and fetal sex.

 Prenatal NO<sub>x</sub> Exposure, Placental Relative mtDNAcn and PE

Women exposed to high concentrations of NOx had lower placental mtDNAcn (−0.20, 95% CI: −0.36, −0.04; p = 0.01) as compared with women exposed to low NOx concentrations during the first trimester (Table 4). The association between prenatal NOx exposure and placental relative mtDNAcn during other time windows were in the same direction, but the effect estimates were slightly lower and not statistically significant. Placental mtDNAcn was not associated with PE neither in unadjusted (OR 0.8; 95% CI: 0.3, 2.1; p = 0.59) nor in adjusted analyses (OR 0.6; 95% CI: 0.2, 2.1; p = 0.46).

Linear regression analysis of placental relative mtDNAcn in women with high NOx exposure as compared with women with low NOx exposure (median cut-off).

Time window Exposure group Unadjusted change in mtDNAcn, (95% CI) p Adjusted§ change in mtDNAcn, (95% CI) p
First trimester Low Reference Reference
High −0.17 (−0.31, −0.03) 0.02 −0.20 (−0.36, −0.04) 0.01
Second trimester Low Reference Reference
High −0.15 (−0.30, −0.01) 0.04 −0.16 (−0.33, 0.01) 0.06
Third trimester Low Reference Reference
High −0.16 (−0.30, −0.01) 0.03 −0.15 (−0.31, 0.01) 0.06
Entire pregnancy Low Reference Reference
High −0.13 (−0.27, 0.02) 0.09 −0.14 (−0.31, −0.02) 0.08

Adjusted for maternal age, pregestational BMI, parity, gestational age, season of birth and fetal sex.

 Mediation Analysis

There was an association between prenatal NOx exposure during first trimester and PE, as well as an association between prenatal NOx exposure during first trimester and placental relative mtDNAcn. To investigate if placental mtDNAcn was mediating this association, we performed a mediation analysis.

In Table 5, Y is the dependent variable, i.e., PE; X the independent variable, prenatal NOx exposure; and Me the proposed mediator variable, mtDNAcn. In order for a variable to be considered a mediator, β32 should be statistically significant, however this was not the case (p = 0.954). Furthermore, β31 should be either substantially smaller than β11, which it was slightly smaller (β11 = 1.382 > β31 = 1.375), or statistically not significant, which it was not (p = 0.010). Therefore, mtDNAcn was not a mediator for the association between prenatal NOx exposure and PE (Table 5).

Mediation analysis to determine the causal mediated effect of placental relative mtDNAcn on the association between prenatal NOx exposure and PE.

Model βi p-value (s)
Step 1 Y = β10 + β11X + ε1 β11 = 1.382 0.008
Step 2 Me = β20 + β21X + ε2 β21 = −0.201 0.014
Step 3 Y = β30 + β31X + β32Me + ε3 β31 = 1.375 β32 = −0.037 0.010 0.954

βi, b-coefficients;

εi, standard errors.

 Prenatal NO<sub>x</sub> Exposure, Placental Relative Telomere Length and PE

There was no statistically significant difference in placental relative telomere length between exposure groups (Table 6). There was no significant association between placental relative telomere length and PE neither in unadjusted (OR 0.6; 95% CI: 0.2, 2.0; p = 0.43) nor in the adjusted analyses (OR 0.4; 95% CI: 0.1, 1.4; p = 0.14). Since the assumption of mediation analysis did not meet, we did not further investigate any mediation effect of placental relative telomere length on the association between prenatal NOx exposure and PE.

Linear regression analysis of placental telomere length in women with high NOx exposure as compared with women with low NOx exposure (median cut-off).

Time window Exposure group Unadjusted change in placental relative telomere length (95% CI) p Adjusted§ change in placental relative telomere length (95% CI) p
First trimester Low Reference Reference
High 0.01 (−0.11, 0.14) 0.82 0.03 (−0.11, 0.17) 0.67
Second trimester Low Reference Reference
High 0.04 (−0.08, 0.17) 0.51 0.07 (−0.08, 0.21) 0.37
Third trimester Low Reference Reference
High 0.04 (−0.09, 0.16) 0.57 0.07 (−0.06, 0.21) 0.30
Entire pregnancy Low Reference Reference
High 0.06 (−0.07, 0.18) 0.37 0.08 (−0.06, 0.21) 0.27

Adjusted for maternal age, pregestational BMI, parity, gestational age, season of birth and fetal sex.

 Discussion

We found that prenatal exposure to ambient NOx during the first trimester was associated with reduced placental relative mtDNAcn and increased risk of PE, which indicate that early pregnancy is the most vulnerable period for damage of the placenta by ambient air pollution. In contrary to our hypothesis, we found no evidence that placental relative mtDNAcn mediate the association between ambient NOx and PE. Placental relative telomere length was not associated with prenatal NOx exposure or with PE.

Early pregnancy exposure to ambient NOx coincides with epigenetic reprogramming of germ cells and preimplantation embryos by DNA methylation occurring in the first trimester (Reik et al., 2001). A plausible mechanism linking prenatal NOx exposure and PE may be abnormal trophoblast fusion, which is a defining characteristic of PE. The maternal-fetal interface of a placenta consists of syncytiotrophoblasts that are formed through trophoblast fusion (Roland et al., 2016). Normal trophoblast syncytialization requires fully functional mitochondria (Poidatz et al., 2015). However, reduced placental mtDNAcn due to prenatal NOx exposure reflects mitochondrial dysfunction which may interfere with trophoblast fusion (Walker et al., 2020). Although the exact mechanisms remain to be clarified, the reduced placental mtDNAcn in association with higher exposure to NOx in the present study could give an indication that air pollution may be a contributing factor for dysfunctional mitochondrial biogenesis during pregnancy. As previously shown, exposure to ambient PM during pregnancy has been found to be associated with increased levels of mitochondrial 8-hydroxy 2′-deoxyguanine, which is considered to be a marker of oxidative DNA damage (Grevendonk et al., 2016) which was found to be correlated with loss of mtDNA (Qian and Van Houten, 2010). Taken together, these findings suggest that ambient NOx exposure during first trimester may play a role in placental dysfunction seen in PE.

However, mediation analysis of placental mtDNAcn undertaken here, did not appear to influence the relationship between prenatal NOx exposure and PE. Interestingly, the only previous study analyzing mtDNAcn as a mediator found that 10% of the association between maternal exposure to nitrogen dioxide (NO2) during pregnancy and birth weight, was mediated by a decrease in placental mtDNAcn (Clemente et al., 2016). The sample size was larger (n = 926) in that study and the mean NO2 concentrations ranged between 21.1 and 25.5 μg/m3 which is higher than the mean NOx concentrations (14–16.2 μg/m3) in our study, which could possibly explain the effect of mediation. Future research in a larger cohort is therefore needed to assess the value of this biomarker for assessing mitochondrial dysfunction in placental pathology.

Although direct deleterious effects of ambient NOx on placenta have not previously been studied, the detrimental effect of fine and ultrafine particles on placenta have been reported. Ambient black carbon particles have been found to be translocated to the fetal side of the placenta collected from full term pregnancies (Bové et al., 2019). A recent in-vitro study of first trimester trophoblast cells exposed to PM2.5 indicated morphological changes, accompanied by altered protease and H2O2 production in mitochondria (Nääv et al., 2020). In another recent in-vitro study, exogenous NOx has been found in the placental homogenates from healthy term pregnancies in the form of iron nitrosyl complexes (Mukosera et al., 2020). Furthermore, placental villous tissue collected from women with pregnancies complicated by PE-FGR had increased levels of NOx as compared to normotensive pregnancies with FGR (Mukosera et al., 2020). Whether prenatal exposure to NOx leads to placental mitochondrial dysfunction, irrespective of the simultaneous effect of PM, remains to be determined in future studies.

In the present study, the effect of prenatal exposure to PM on placental mtDNAcn or PE was not evaluated. Both NOx and PM mainly, however, originate from various combustion processes, such as a combustion engine. In our previous population-based study in southern Scania, there was a statistically significant, high positive correlation between mean NOx and local PM2.5 concentrations (rp = 0.90, p < 0.01) and mean NOx and total PM2.5 concentrations (rp = 0.58, p < 0.01) during the entire pregnancy (Mandakh et al., 2020). Therefore, it can be hypothesized that exposure to NOx during pregnancy is representative also for exposure to locally emitted PM2.5, however, given the difference in biological fate of PM2.5 and NOx in human body, it is not possible to account the change in placental mtDNAcn observed in the NOx exposure model to PM2.5. Nevertheless, this study do not have the data to fully support this hypothesis since we could only model NOx exposure for the majority of participants.

In contrast to earlier findings, however, we found no association between prenatal NOx exposure and placental telomere length, nor between placental telomere length and PE. Previous studies reported associations between telomere shortening in placenta and nearest proximity to major roads (Bijnens et al., 2015) and with maternal exposure to PM2.5 during pregnancy (Martens et al., 2017). A 5 μg/m3 increment in PM2.5 exposure during the entire pregnancy was associated with 13.2% decrease in placental telomere length (Martens et al., 2017). As compared to our study, Martens et al. (2017) had a larger sample size (n = 641) and a larger variation (4.3–32.5 μg/m3) in weekly PM2.5 exposure. Furthermore, shorter telomeres were found in placental trophoblasts from pregnancies complicated by PE, FGR and PE with FGR as compared to uncomplicated pregnancies (Biron-Shental et al., 2010). However, the finding that telomere shortening was not found in cord blood cells from pregnancies complicated with PE suggests that air pollution-related shorter telomere length might be restricted to placental trophoblasts (Sukenik-Halevy et al., 2016).

To our knowledge, this is the first study to report an association between prenatal exposure to ambient NOx, placental dysfunction measured with biomarkers of toxicity, and PE. Another key strength of this study is the access to a unique biobank specially created for studying preeclampsia with extensive biological samples linked with detailed clinical data. In addition, the individual exposure was assessed at residence level by high resolution dispersion model after obtaining the geocoded residential locations of each participant by linking the unique personal identification number with the well-documented regional population registry. Further studies, which take these biomarkers of toxicity into account, should also measure oxidative stress biomarkers in the placenta to better explain the association between prenatal ambient air pollution exposure and the risk of PE.

The generalizability of these results is subject to at least four limitations. Firstly, our study might be predisposed to exposure misclassification due to a residual confounder related to unmeasured socioeconomic status of the participants and not accounting for exposure during commuting, work, or indoor exposure. We were not able to obtain data on building characteristics and socioeconomic status of the participants. Difficulties to assess the total exposure is, however, a common limitation in air pollution epidemiologic studies (Pedersen et al., 2014). Secondly, it is beyond the scope of this study to further examine associations with sub-clinical diagnosis of PE such as early- and late-onset PE, or PE with or without FGR, due to the limited sample size. In addition, we were not able to match PE cases with normotensive controls based on covariates such as gestational age due to the fact that they were not equally distributed among the term and preterm births after the exposure assessment. Moreover, we could not exclude the fact that prematurity can be a result of other pathologic pregnancies and therefore a preterm normotensive pregnant women were not selected as a preterm control to match preterm PE cases. Lastly, the current study has only examined NOx exposure due to the technical error and economic constraint to model other ambient air pollutants, such as PM2.5 during exposure assessment.

Our data show that exposure to ambient air pollution (NOx) during the first trimester was associated with reduced placental relative mtDNAcn and an increased risk of PE. However, we did not find any evidence that mtDNAcn or TL mediated the association between air pollution and PE. Future research should further investigate the role of mtDNAcn for pregnancy complications in relation to exposure to ambient air pollution during pregnancy.

 Data Availability Statement

The datasets related to the clinical parameters and the DNA-analysis are all available upon request. However the dataset 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.

 Ethics Statement

The studies involving human participants were reviewed and approved by Lund University Ethical Review Board, LU803-2 (2015/14). The patients/participants provided their written informed consent to participate in this study.

 Author Contributions

EM, YM, SH, KB, AO, CI, and LE: conceptualization, methodology, and writing—review and editing. YM, KB, and AO: formal analysis and investigation. KB, SH, and LE: resources. EM, SH, and LE: data curation. YM: writing—original draft preparation. EM, AO, KB, and SH: supervision. EM: project administration. EM, CI, and SH: funding acquisition. All authors have read and approved the submitted version of the manuscript.

 Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

We thank Eva Hansson for her support on DNA and RNA extractions, Ralf Rittner for his support on exposure assessment, Daniela Pineda and Eleni Korre Nte Paoula for their support with genetic analyses, and Kristoffer Mattisson for his support on creating a map.

 Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/ftox.2021.659407/full#supplementary-material

 References Bai W. Li Y. Niu Y. Ding Y. Yu X. Zhu B. . (2020). Association between ambient air pollution and pregnancy complications: a systematic review and meta-analysis of cohort studies. Environ. Res. 185:109471. 10.1016/j.envres.2020.10947132276169 Barker D. J. P. (2004). The developmental origins of adult disease. J. Am. College Nutr. 23, 588S−595S. 10.1080/07315724.2004.1071942827374623 Baron R. M. Kenny D. A. (1986). The moderator-mediator variable distinction in social psychological research: conceptual, strategic, and statistical considerations. J. Pers. Soc. Psychol. 51, 1173–1182. 10.1037/0022-3514.51.6.11733806354 Bijnens E. Zeegers M. P. Gielen M. Kicinski M. Hageman G. J. Pachen D. . (2015). Lower placental telomere length may be attributed to maternal residential traffic exposure; a twin study. Environ. Int. 79, 1–7. 10.1016/j.envint.2015.02.00825756235 Biron-Shental T. Sukenik-Halevy R. Sharon Y. Goldberg-Bittman L. Kidron D. Fejgin M. D. . (2010). Short telomeres may play a role in placental dysfunction in preeclampsia and intrauterine growth restriction. Am. J. Obstet. Gynecol. 202, 381.e381–387. 10.1016/j.ajog.2010.01.03620350645 Bonvallot V. Baeza-Squiban A. Baulig A. Brulant S. Boland S. Muzeau F. . (2001). Organic compounds from diesel exhaust particles elicit a proinflammatory response in human airway epithelial cells and induce cytochrome p450 1A1 expression. Am. J. Resp. Cell Mol. Biol. 25, 515–521. 10.1165/ajrcmb.25.4.451511694458 Bové H. Bongaerts E. Slenders E. Bijnens E. M. Saenen N. D. Gyselaers W. . (2019). Ambient black carbon particles reach the fetal side of human placenta. Nat. Commun. 10:3866. 10.1038/s41467-019-11654-331530803 Brown M. A. Magee L. A. Kenny L. C. Karumanchi S. A. McCarthy F. P. Saito S. . (2018). The hypertensive disorders of pregnancy: ISSHP classification, diagnosis & management recommendations for international practice. Pregnancy Hypertension 13, 291–310. 10.1016/j.preghy.2018.05.004 Burton G. J. Sebire N. J. Myatt L. Tannetta D. Wang Y. L. Sadovsky Y. . (2014). Optimising sample collection for placental research. Placenta 35, 9–22. 10.1016/j.placenta.2013.11.00524290528 Burton G. J. Yung H. W. Murray A. J. (2017). Mitochondrial – Endoplasmic reticulum interactions in the trophoblast: stress and senescence. Placenta 52, 146–155. 10.1016/j.placenta.2016.04.00127091649 Bustamante J. Ramírez-Vélez R. Czerniczyniec A. Cicerchia D. Aguilar de Plata A. C. Lores-Arnaiz S. (2014). Oxygen metabolism in human placenta mitochondria. J. Bioenerg. Biomembr. 46, 459–469. 10.1007/s10863-014-9572-x25234730 Choi H. Melly S. Spengler J. (2015). Intraurban and longitudinal variability of classical pollutants in Kraków, Poland 2000–2010. Int. J. Environ. Res. Public Health 12, 4967–4991. 10.3390/ijerph12050496725955530 Clemente D. Casas M. Janssen B. Lertxundi A. Santa Marina L. Iñiguez C. . (2017). Prenatal ambient air pollution exposure, infant growth and placental mitochondrial DNA content in the INMA birth cohort. Environ. Res. 157, 96–102. 10.1016/j.envres.2017.05.01828535425 Clemente D. B. Casas M. Vilahur N. Begiristain H. Bustamante M. Carsin A. E. . (2016). Prenatal ambient air pollution, placental mitochondrial DNA content, and birth weight in the INMA (Spain) and ENVIRONAGE (Belgium) birth cohorts. Environ. Health Perspect. 124, 659–665. 10.1289/ehp.140898126317635 Cohen A. J. Brauer M. Burnett R. Anderson H. R. Frostad J. Estep K. . (2017). Estimates and 25-year trends of the global burden of disease attributable to ambient air pollution: an analysis of data from the Global Burden of Diseases Study 2015. Lancet 389, 1907–1918. 10.1016/S0140-6736(17)30505-628408086 Fossati S. Valvi D. Martinez D. Cirach M. Estarlich M. Fernández-Somoano A. . (2020). Prenatal air pollution exposure and growth and cardio-metabolic risk in preschoolers. Environ. Int. 138:105619. 10.1016/j.envint.2020.10561932193046 Fu L. Chen Y. Yang X. Yang Z. Liu S. Pei L. . (2019). The associations of air pollution exposure during pregnancy with fetal growth and anthropometric measurements at birth: a systematic review and meta-analysis. Environ. Sci. Pollut. Res. Int. 26, 20137–20147. 10.1007/s11356-019-05338-031111384 Grevendonk L. Janssen B. G. Vanpoucke C. Lefebvre W. Hoxha M. Bollati V. . (2016). Mitochondrial oxidative DNA damage and exposure to particulate air pollution in mother-newborn pairs. Environ. Health 15:10. 10.1186/s12940-016-0095-226792633 Gyselaers W. (2020). Preeclampsia is a syndrome with a cascade of pathophysiologic events. J. Clin. Med. 9:2245. 10.3390/jcm907224532679789 Harri M. Svoboda P. Mori T. Mutanen P. Kasai H. Savela K. (2005). Analysis of 8-hydroxydeoxyguanosine among workers exposed to diesel particulate exhaust: comparison with urinary metabolites and PAH air monitoring. Free Radic. Res. 39, 963–972. 10.1080/1071576050019011516087477 Hayflick L. Moorhead P. S. (1961). The serial cultivation of human diploid cell strains. Exp. Cell Res. 25, 585–621. 10.1016/0014-4827(61)90192-613905658 Herlin M. Broberg K. Igra A. M. Li H. Harari F. Vahter M. (2019). Exploring telomere length in mother-newborn pairs in relation to exposure to multiple toxic metals and potential modifying effects by nutritional factors. BMC Med. 17:77. 10.1186/s12916-019-1309-630971237 Holland O. Dekker Nitert M. Gallo L. A. Vejzovic M. Fisher J. J. Perkins A. V. (2017). Review: Placental mitochondrial function and structure in gestational disorders. Placenta 54, 2–9. 10.1016/j.placenta.2016.12.01228024805 Ichinose T. Yajima Y. Nagashima M. Takenoshita S. Nagamachi Y. Sagai M. (1997). Lung carcinogenesis and formation of 8-hydroxy-deoxyguanosine in mice by diesel exhaust particles. Carcinogenesis 18, 185–192. 10.1093/carcin/18.1.1859054605 Iodice S. Hoxha M. Ferrari L. Carbone I. F. Anceschi C. Miragoli M. . (2018). Particulate air pollution, blood mitochondrial DNA copy number, and telomere length in mothers in the first trimester of pregnancy: effects on fetal growth. Oxid. Med. Cell Longev. 2018:5162905. 10.1155/2018/516290530524658 Janssen B. G. Byun H.-M. Gyselaers W. Lefebvre W. Baccarelli A. A. Nawrot T. S. (2015). Placental mitochondrial methylation and exposure to airborne particulate matter in the early life environment: an ENVIRONAGE birth cohort study. Epigenetics 10, 536–544. 10.1080/15592294.2015.104841225996590 Janssen B. G. Munters E. Pieters N. Smeets K. Cox B. Cuypers A. . (2012). Placental mitochondrial DNA content and particulate air pollution during in utero life. Environ. Health Perspect. 120, 1346–1352. 10.1289/ehp.110445822626541 Lee M.-W. Chen M.-L. Lung S.-C. C. Tsai C.-J. Lai C.-F. S. Yang S.-C. . (2012). Increase of urinary concentrations of 8-hydroxy-2′-deoxyguanosine in diesel exhaust emission inspector exposed to polycyclic aromatic hydrocarbons. Int. Arch. Occupational Environ. Health 85, 273–282. 10.1007/s00420-011-0663-221674251 Li Z. Liu Q. Xu Z. Guo X. Wu S. (2020). Association between short-term exposure to ambient particulate air pollution and biomarkers of oxidative stress: a meta-analysis. Environ. Res. 191:110105. 10.1016/j.envres.2020.11010532835677 Lin S. Huo X. Zhang Q. Fan X. Du L. Xu X. . (2013). Short placental telomere was associated with cadmium pollution in an electronic waste recycling town in China. PLoS ONE 8:e60815. 10.1371/journal.pone.006081523565277 Malmqvist E. Jakobsson K. Tinnerberg H. Rignell-Hydbom A. Rylander L. (2013). Gestational diabetes and preeclampsia in association with air pollution at levels below current air quality guidelines. Environ. Health Perspect. 121, 488–493. 10.1289/ehp.120573623563048 Mandakh Y. Rittner R. Flanagan E. Oudin A. Isaxon C. Familari M. . (2020). Maternal exposure to ambient air pollution and risk of preeclampsia: a population-based cohort study in Scania, Sweden. Int. J. Environ. Res. Public Health 17:1744. 10.3390/ijerph1705174432155988 Manna S. McCarthy C. McCarthy F. P. (2019). Placental ageing in adverse pregnancy outcomes: telomere shortening, cell senescence, and mitochondrial dysfunction. Oxid. Med. Cell. Longev. 2019:3095383. 10.1155/2019/309538331249642 Martens D. S. Cox B. Janssen B. G. Clemente D. B. P. Gasparrini A. Vanpoucke C. . (2017). Prenatal air pollution and newborns' predisposition to accelerated biological aging. JAMA Pediatr. 171, 1160–1167. 10.1001/jamapediatrics.2017.302429049509 Martens D. S. Nawrot T. S. (2016). Air pollution stress and the aging phenotype: the telomere connection. Curr. Environ. Health Rep. 3, 258–269. 10.1007/s40572-016-0098-827357566 Melody S. M. Wills K. Knibbs L. D. Ford J. Venn A. Johnston F. (2020). Maternal exposure to ambient air pollution and pregnancy complications in Victoria, Australia. Int. J. Environ. Res. Public Health 17:2572. 10.3390/ijerph1707257232283665 Møller P. Loft S. (2010). Oxidative damage to DNA and lipids as biomarkers of exposure to air pollution. Environ. Health Perspect. 118, 1126–1136. 10.1289/ehp.090172520423813 Mukosera G. T. Clark T. C. Ngo L. Liu T. Schroeder H. Power G. G. . (2020). Nitric oxide metabolism in the human placenta during aberrant maternal inflammation. J. Physiol. 598, 2223–2241. 10.1113/JP27905732118291 Muralimanoharan S. Maloyan A. Mele J. Guo C. Myatt L. G. Myatt L. (2012). MIR-210 modulates mitochondrial respiration in placenta with preeclampsia. Placenta 33, 816–823. 10.1016/j.placenta.2012.07.00222840297 Nääv Å. Erlandsson L. Isaxon C. Åsander Frostner E. Ehinger J. Sporre M. K. . (2020). Urban PM2.5 induces cellular toxicity, hormone dysregulation, oxidative damage, inflammation, and mitochondrial interference in the HRT8 trophoblast cell line. Front. Endocrinol. 11:75. 10.3389/fendo.2020.0007532226408 Opsis (2006). Available online at: http://www.opsis.se/ (accessed August 09, 2020). Pańczyszyn A. Boniewska-Bernacka E. Goc A. (2020). The role of telomeres and telomerase in the senescence of postmitotic cells. DNA Repair. 95:102956. 10.1016/j.dnarep.2020.10295632937289 Pedersen M. Stayner L. Slama R. Sorensen M. Figueras F. Nieuwenhuijsen M. J. . (2014). Ambient air pollution and pregnancy-induced hypertensive disorders: a systematic review and meta-analysis. Hypertension 64, 494–500. 10.1161/HYPERTENSIONAHA.114.0354524935943 Phillips J. K. Bernstein I. M. Mongeon J. A. Badger G. J. (2004). Seasonal variation in preeclampsia based on timing of conception. Obstetrics Gynecol. 104, 1015–1020. 10.1097/01.AOG.0000143306.88438.cf15516394 Poidatz D. Dos Santos E. Gronier H. Vialard F. Maury B. De Mazancourt P. . (2015). Trophoblast syncytialisation necessitates mitochondrial function through estrogen-related receptor-γ activation. Mol. Hum. Rep. 21, 206–216. 10.1093/molehr/gau10225376642 Qian W. Van Houten B. (2010). Alterations in bioenergetics due to changes in mitochondrial DNA copy number. Methods 51, 452–457. 10.1016/j.ymeth.2010.03.00620347038 Ramiro-Cortijo D. de la Calle M. Rodriguez-Rodriguez P. Phuthong S. López de Pablo Á. L. Martín-Cabrejas M. A. . (2020). First trimester elevations of hematocrit, lipid peroxidation and nitrates in women with twin pregnancies who develop preeclampsia. Pregnancy Hypertension 22, 132–135. 10.1016/j.preghy.2020.09.00532971355 Redman C. W. G. Staff A. C. Roberts J. M. (2021). Syncytiotrophoblast stress in preeclampsia: the convergence point for multiple pathways. Am. J. Obstetrics Gynecol. Available online at: https://www.ajog.org/article/S0002-9378(20)31115-7/fulltext. 10.1016/j.ajog.2020.09.04733546842 Reik W. Dean W. Walter J. (2001). Epigenetic reprogramming in mammalian development. Science 293, 1089–1093. 10.1126/science.1063443 Richter T. Zglinicki T. V. (2007). A continuous correlation between oxidative stress and telomere shortening in fibroblasts. Exp. Gerontol. 42, 1039–1042. 10.1016/j.exger.2007.08.00517869047 Rittner R. Gustafsson S. Spanne M. Malmqvist E. (2020). Particle concentrations, dispersion modelling and evaluation in southern Sweden. SN Appl. Sci. 2:1013. 10.1007/s42452-020-2769-1 Roland C. S. Hu J. Ren C. E. Chen H. Li J. Varvoutis M. S. . (2016). Morphological changes of placental syncytium and their implications for the pathogenesis of preeclampsia. Cell Mol. Life Sci. 73, 365–376. 10.1007/s00018-015-2069-x26496726 Roubicek D. A. Souza-Pinto N. C. (2017). Mitochondria and mitochondrial DNA as relevant targets for environmental contaminants. Toxicology 391, 100–108. 10.1016/j.tox.2017.06.01228655544 Rudra C. B. Williams M. A. (2005). Monthly variation in preeclampsia prevalence: Washington state, 1987–2001. J. Maternal Fetal Neonatal Med. 18, 319–324. 10.1080/1476705050027583816390791 Saenen N. D. Martens D. S. Neven K. Y. Alfano R. Bové H. Janssen B. G. . (2019). Air pollution-induced placental alterations: an interplay of oxidative stress, epigenetics, and the aging phenotype? Clin. Epigenetics 11:124. 10.1186/s13148-019-0688-z31530287 Silva A. V. Oliveira C. M. Canha N. Miranda A. I. Almeida S. M. (2020). Long-term assessment of air quality and identification of aerosol sources at Setúbal, Portugal. Int. J. Environ. Res. Public Health 17:5447. 10.3390/ijerph1715544732731647 Socialstyrelsen (2018). ICD-10-SE. Available Online at: http://www.socialstyrelsen.se/publikationer2018/2018-4-13 (accessed August 20, 2018). Stroh E. Harrie L. Gustafsson S. (2007). A study of spatial resolution in pollution exposure modelling. Int. J. Health Geogr. 6:19. 10.1186/1476-072X-6-1917547740 Stroh E. Rittner R. Oudin A. Ardö J. Jakobsson K. Björk J. . (2012). Measured and modeled personal and environmental NO2 exposure. Population Health Metrics 10:10. 10.1186/1478-7954-10-1022681784 Sukenik-Halevy R. Amiel A. Kidron D. Liberman M. Ganor-Paz Y. Biron-Shental T. (2016). Telomere homeostasis in trophoblasts and in cord blood cells from pregnancies complicated with preeclampsia. Am. J. Obstetrics Gynecol. 214, 283.e281–283.e287. 10.1016/j.ajog.2015.08.05026321036 Tranquilli A. L. Brown M. A. Zeeman G. G. Dekker G. Sibai B. M. (2013). The definition of severe and early-onsetpreeclampsia. Statements from the International Society for the Study of Hypertension in Pregnancy (ISSHP). Pregnancy Hypertens. 3, 44–47. 10.1016/j.preghy.2012.11.00126105740 USEPA (1998). User's Guide for the AMS/EPA Regulatory Model-AERMOD. Available Online at: https://www3.epa.gov/scram001/7thconf/aermod/aermodug.pdf (accessed August 09, 2020). Venegas V. Halberg M. C. (2012). “Measurement of mitochondrial DNA copy number,” in Mitochondrial Disorders: Biochemical and Molecular Analysis, ed. Wong L.-J.C. (Houston, TX: Springer), 327–335. 10.1007/978-1-61779-504-6_2233220500 Victorelli S. Passos J. F. (2017). Telomeres and cell senescence - size matters not. EBioMedicine 21, 14–20. 10.1016/j.ebiom.2017.03.02728347656 Vinzents P. S. Møller P. Sørensen M. Knudsen L. E. Hertel O. Jensen F. P. . (2005). Personal exposure to ultrafine particles and oxidative DNA damage. Environ. Health Perspect. 113, 1485–1490. 10.1289/ehp.756216263500 Vishnyakova P. A. Volodina M. A. Tarasova N. V. Marey M. V. Tsvirkun D. V. Vavina O. V. . (2016). Mitochondrial role in adaptive response to stress conditions in preeclampsia. Sci. Rep. 6:32410. 10.1038/srep3241027573305 Wacker J. Schulz M. Frühauf J. Chiwora F. Solomayer E. Bastert G. (1998). Seasonal change in the incidence of preeclampsia in Zimbabwe. Acta Obstet. Gynecol. Scand. 77, 712–716. 10.1080/j.1600-0412.1998.770703.x9740517 Walker O. L. S. Ragos R. Wong M. K. Adam M. Cheung A. Raha S. (2020). Reactive oxygen species from mitochondria impacts trophoblast fusion and the production of endocrine hormones by syncytiotrophoblasts. PLoS ONE 15:e0229332. 10.1371/journal.pone.022933232092105 Xu Y. Lindh C. H. Jönsson B. A. G. Broberg K. Albin M. (2018). Occupational exposure to asphalt mixture during road paving is related to increased mitochondria DNA copy number: a cross-sectional study. Environ. Health 17:29. 10.1186/s12940-018-0375-029587765 Yu H. Yin Y. Zhang J. Zhou R. (2020). The impact of particulate matter 2.5 on the risk of preeclampsia: an updated systematic review and meta-analysis. Environ. Sci. Pollut. Res. 27, 37527–37539. 10.1007/s11356-020-10112-832740838 Zhang H. Wang Q. He S. Wu K. Ren M. Dong H. . (2020). Ambient air pollution and gestational diabetes mellitus: a review of evidence from biological mechanisms to population epidemiology. Sci. Total Environ. 719:137349. 10.1016/j.scitotenv.2020.13734932114225

Funding. This project has been funded by Swedish Research Council FORTE 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.