Prenatal fine particulate exposure associated with reduced childhood lung function and nasal epithelia GSTP1 hypermethylation: Sex-specific effects.

BACKGROUND: In utero exposure to particulate matter with an aerodynamic diameter of less than 2.5 μm (PM2.5) has been linked to child lung function. Overlapping evidence suggests that child sex and exposure timing may modify effects and associations may be mediated through glutathione S-transferase P1 (GSTP1) methylation.
METHODS: We prospectively examined associations among prenatal PM2.5 exposure and child lung function and GSTP1 methylation in an urban pregnancy cohort study. We employed a validated satellite-based spatiotemporally resolved prediction model to estimate daily prenatal PM2.5 exposure over gestation. We used Baysian distributed lag interaction models (BDLIMs) to identify sensitive windows for prenatal PM2.5 exposure on child lung function and nasal epithelia GSTP1 methylation at age 7 years, and to examine effect modification by child sex.
RESULTS: BDLIMs identified a sensitive window for prenatal PM2.5 exposure at 35-40 weeks gestation [cumulative effect estimate (CEE) = - 0.10, 95%CI = - 0.19 to - 0.01, per μg/m3 increase in PM2.5] and at 36-40 weeks (CEE = - 0.12, 95%CI = - 0.20 to - 0.01) on FEV1 and FVC, respectively, in boys. BDLIMs also identified a sensitive window of exposure at 37-40 weeks gestation between higher prenatal PM2.5 exposure and increased GSTP1 percent methylation. The association between higher GSTP1 percent methylation and decreased FEV1 was borderline significant in the sample as a whole (β = - 0.37, SE = 0.20, p = 0.06) and in boys in stratified analyses (β = - 0.56, SE = 0.29, p = 0.05).
CONCLUSIONS: Prenatal PM2.5 exposure in late pregnancy was associated with impaired early childhood lung function and hypermethylation of GSTPI in DNA isolated from nasal epithelial cells. There was a trend towards higher GSTP1 percent methylation being associated with reduced FEV1. All findings were most evident among boys.

Background

The importance of early life lung development on future respiratory health is well established [1, 2]. Understanding environmental exposures and mechanisms that lead to and maintain this early predisposition is key to identifying children at increased risk of future respiratory disease and, ultimately, reducing chronic pulmonary disease morbidity and mortality. Genetic variants and early life factors including tobacco smoke, birth weight, gestational age, and asthma are associated with reduced lung function over the life course, however these factors account for a relatively small proportion of the risk suggesting that unidentified factors remain.

Animal models have highlighted the need to understand PM effects on lung development starting prenatally [3, 4]. Recent human data corroborates the association between in utero PM2.5 exposure and asthma risk in children with an increased focus on sensitive windows for exposure effects [5]. To our knowledge, only one study has reported associations between prenatal exposure to PM2.5 and reduced child lung function, however exposure was measured at a single time point and extrapolated over pregnancy and thus could not flexibly assess sensitive windows of effect [6]. Because the well-orchestrated processes involved in programming lung growth and development over gestation do not necessarily occur within clinically defined trimesters, research that allows flexibility in identifying sensitive windows may be particularly informative [7]. Further, prior research suggests that boys may be more vulnerable to the effects of prenatal environmental exposures, including air pollution exposure [5], perhaps mediated through increased vulnerability to oxidative stress.

Epigenetic programing, including DNA methylation, starting during fetal life is one pathway by which environmental factors may influence gene expression thus programing future disease risk [8, 9]. Global and gene-specific methylation is altered in response to prenatal exposures and these changes appear stable in early childhood [10, 11]. For example, an epigenome-wide study of cord blood DNA methylation in children prenatally exposed to nitrogen dioxide (NO2) demonstrated differential methylation of genes involved in mitochondrial function [12], and a study examining exposure to PM10 averaged over trimesters in relation to loci specific methylation found that early pregnancy PM10 exposure was associated with placental DNA methylation of LINE1 and HSD11B2 [13].

Particulate matter is a strong oxidant able to generate reactive oxygen species (ROS); oxidative stress pathways are believed to be central in the association between air pollution exposure and respiratory outcomes [14, 15]. The glutathione S-transferase (GST) superfamily, specifically the P1 isoform (GSTP1), is expressed in the respiratory tract and functions in oxidant defenses, xenobiotic metabolism, and detoxification of hyperperoxides; GSTP1 variants have been found to increase susceptibility to tobacco smoke and air pollution in several studies [16-18]. GSTP1 CpG promoter hypermethylation may inactivate GSTP1 [19], thereby reducing cellular detoxification capabilities and increasing susceptibility to oxidative stress [20].

The airway epithelium is an important driver of pulmonary disease – for example, the airway epithelium regulates inflammatory responses in asthma [21]. Lower airway, or bronchial, epithelial cell sampling requires invasive procedures (e.g., bronchoscopy or lung tissue biopsy), which are not suitable for healthy children. As nasal epithelial cells (NECs) are surrogates for bronchial epithelial cells and more readily accessible, it has been proposed that examination of associations between environmental pollutants and epigenetic changes in DNA isolated from NECs may provide insight into the mechanisms that underlie associations between PM2.5 and airway diseases [22-24].

We examined the impact of prenatal PM2.5 exposure on children’s lung function in an urban, ethnically mixed longitudinal pregnancy cohort. We employed Bayesian distributed lag interaction models (BDLIMs) to estimate windows of sensitivity between prenatal PM2.5 exposure and children’s pulmonary function measured at age 7 years and to examine effect modification by child sex. We also examined associations between prenatal PM2.5 exposure and GSTP1 percent methylation in DNA from NECs isolated in children at time of spirometry. Finally we explored associations between GSTP1 percent methylation and children’s lung function measures. We hypothesized that children born to mothers with higher PM2.5 exposure would have lower spirometry outcomes and hypermethylation of NEC GSTP1 at age 7 years, relative to those born to mothers with low PM2.5 exposure. We posited that the highest levels of GSTP1 percent methylation would be associated with greater reductions in spirometry measures. We hypothesized that boys would be particularly vulnerable.

Methods

Study participants

Subjects were from the Asthma Coalition on Community, Environment and Social Stress (ACCESS) project, an urban ethnically-diverse pregnancy cohort designed to examine the effects of environmental exposures on childhood respiratory outcomes [25]. From August 2002 through January 2007, n = 500 pregnant women were recruited from two Boston hospitals and affiliated health centers at 28.4 ± 7.9 weeks gestation and 455 gave birth to a live singleton infant. Of those approached and eligible, 78% were enrolled; no significant differences in education, race/ethnicity, or income were observed in those who enrolled versus those who were not enrolled. A subset (230 of 375) of children actively followed at age 6.99 ± 0.89 years participated in a pulmonary function visit, during which time NECs were also collected. Procedures were approved by human studies committees at the Brigham and Women’s Hospital and Boston Medical Center and written consent was obtained from all mothers and assent was obtained for children age ≥7 years.

Daily prenatal PM2.5 levels

Maternal residence was geocoded at baseline and updated throughout gestation if the subject moved; geocoded addresses were used to estimate residence-specific daily prenatal PM2.5 exposure over the pregnancy, as detailed previously [26]. Briefly, a novel spatiotemporal model estimated daily high resolution PM2.5 by regressing daily surface PM2.5 measurements, taken from the U.S. Environmental Protection Agency Air Quality System and Interagency Monitoring of Protected Visual Environments Network, with daily aerosol optical depth measurements, land-use terms (elevation, distance to major roads, percent open space, point emissions, and area emissions), and meteorological variables (temperature, wind speed, visibility). The model linked moderate resolution imaging spectroradiometer satellite-derived aerosol optical depth measurements at a 1 km × 1 km spatial resolution and ground measurement covering the Northeast USA were calibrated on a daily basis and validated with robust out of sample 10-fold cross-validation. The mean cross-validated R2 for the New England sub-region that includes the greater Boston area included in this study was 0.88. To reduce noise created by the day-to-day PM2.5 variation, women’s gestational PM2.5 exposure estimates were calculated by averaging daily predictions over each week of the pregnancy for these analyses.

Pulmonary function testing

Research assistants trained by an experienced pediatric respiratory therapist and pulmonologists (WJM, RJW) on our team measured child height, weight, and lung function with over-reading performed for all spirometry tests to ensure quality control. Height was measured to the nearest 0.1 cm using a stadiometer and weight was measured to the nearest 0.1 kg using an electronic scale. Spirometry was performed in participant homes with a portable MedGraphics™ laptop supported spirometer, which displays real-time flow-volume plots to facilitate testing. Testing procedures met American Thoracic Society guidelines [27, 28] with techniques modified for children ≤8 years of age [29, 30]. Subjects without acute respiratory symptoms for ≥3 weeks were eligible for testing. Short-acting beta-agonists, atropinics and theophylline preparations were withheld for 4 h and long-acting beta-agonists for 12 h before testing. Forced vital capacity (FVC, milliliters), forced expiratory volume in 1 s (FEV1, milliliters), and forced expiratory flow between 25 and 75% of the FVC (FEF25–75, milliliters per second) were recorded from a minimum of 3 (no more than 8) maneuvers. Lung function measures, height and weight were all approximately normally distributed. Raw FEV1, FVC, FEF25–75 and FEV1/FVC values were adjusted for age, sex, height, and race/ethnicity using multivariable regression, and then converted to z-scores with a mean of 0 and a standard deviation of 1 to describe each child’s position relative to that of other individuals in the distribution, as done in previous studies [31, 32].

GSTP-1 percent methylation

Epithelial cells were collected from the anterior nares. Briefly, the child was seated in a chair with head tilted back. Each nostril was cleaned with a cotton swab wet in saline solution. A 1 mm Microinvasive Boston specimen brush was inserted into each nostril until the tip was immediately inferior to the nasal bone and vigorously rotated along the nares avoiding the nasal septum. Notably, prior work by our group has demonstrated methylation concordance (R2 = 0.93) between inferior turbinate and anterior nares cells, providing an acceptable and safer, less-invasive sampling procedure for young children [33]. The brushes were immediately submersed in 6 mL of 2% acetylcysteine-saline solution. Samples were incubated at room temperature with gentle shaking for 30 min, then the brush was removed and the solution was centrifuged at 4500 rpm for 10 min. DNA was extracted from the cell pellet using the Promega MaxwellR 16 Buccal Swab LEV DNA Kit following manufacturer instructions and stored at − 80 °C prior to bisulfite conversion.

GSTP1 methylation analyses

The nasal cell DNA was bisulfite converted using the EZ DNA Methylation-Gold Kit (Zymo Research, CA, USA) according to the manufacturer’s protocol. Polymerase chain reaction (PCR) and pyrosequencing was used to quantify the percentage methylation at each site of interest. Control samples were used to verify bisulfite conversion efficiency. The measured degree of methylation is presented as the percentage of methylated cytosines divided by the sum of methylated and unmethylated cytosines (%5mC). Pyrosequencing primers were designed not to overlap with any single-nucleotide polymorphism or repeated elements [34]. The forward PCR primer sequence was TTTGGGAAAGAGGGAAAGGT and the reverse 5′ end biotin labeled primer was AACCTTATAAAAATAATCCC. PCR cycling conditions were a 5 min hold at 95 °C followed by 95 °C for 45 s, 50 °C for 45 s, and 72 °C for 45 s for 45 cycles, followed by a 5 min hold at 72 °C. PCR products were purified and sequenced on a Q96 MD pyrosequencing system (Qiagen), as previously described [35]. The sequencing primer sequence was AGAGGGAAAGGTTTTTT and the sequence entry was CGGTTAGTTGCGCGGCGATTTCGGGGATTTTAGGGCGTTTTTTTGCGGTCGACGTTCGGGGTGTAGCGGTCGTCGGGGTTGGGGTCGGCGGGAGTTCGCGGGATT. Seventeen CpG sites were analyzed on chromosome 11 (GRCh37/hg19).

Covariates

Maternal race, age at enrollment, and education were obtained through questionnaires. Mothers who reported smoking at baseline or in the third trimester of pregnancy were classified as smokers; postnatal secondhand tobacco smoke exposure was classified based on maternal-reported smoking or report of others smoking in the home at each postpartum interview. Maternal-reported clinician-diagnosed asthma was determined at 3-month intervals for the first 24 months of life and annually thereafter until age 6 years. Mothers were asked, “Has a doctor or nurse ever said that your child had asthma?”. The majority of children received a diagnosis of asthma after 3 years of age (78.7%) [5].

Statistical analysis

Of the 230 children recruited for the spirometry visit, 211 (92%) provided acceptable spirometry. Complete data on prenatal daily PM2.5 levels was available in 192 (91%); of these, 171 were born full-term (born at ≥37 weeks gestation) and were included in the analysis examining associations between prenatal PM2.5 and spirometry outcomes. In this analytic sample, 131/171 (76.6%) also had adequate nasal cell DNA for GSTP1 assays.

In order to identify sensitive windows for the effects of prenatal PM2.5 in relation to spirometry outcomes and GSTP1 DNA methylation, as well as effect modification by sex, we implemented Bayesian distributed lag interaction models (BDLIMs) as detailed previously [7]. Using a BDLIM, we estimated the time-varying association for each participant’s weekly exposures throughout the gestational period and z-scores of spirometry outcomes. BDLIM extends the traditional constrained distributed lag model framework that identifies sensitive windows [36], and also accounts for within window effects and tests for effect modifications. Significant sensitive exposure windows were identified as weeks during pregnancy with a statistically significant association.

We first conducted BDLIM in the overall sample and then examined effect modification by child’s sex. BDLIM partitions the distributed lag function into two components: 1) the weights that identify sensitive windows of susceptibility, and 2) the coefficients that identify the magnitude of the within-window effects. Thus, BDLIM is able to detect whether boys and girls have either the same, or different, sensitive windows (weights) and magnitude of within-window effects (effects). In other words, models with 4 types of patterns (i.e., same weight and same effect, same weight and different effects, different windows and same effect, different windows and different effects) were tested to determine whether the associations between weekly PM2.5 and the outcomes were modified by sex. Subsequently, in addition to estimating the time-varying associations, BDLIM is also able to estimate the cumulative effects of PM2.5 exposure over pregnancy that account for both sensitive windows and within-window effects corresponding to each sex. Deviance information criterion (DIC) was used to determine the best fitting model with optimal number of knots and whether the weights and effects are different across sex.

For analysis on spirometry outcomes (i.e., FEV1, FVC, FEF25–75, FEV1/FVC), in addition to the covariates adjusted in z-scores including age, sex, height and race/ethnicity, we also adjusted for maternal age and education in the main analysis. Analysis on GSTP1 percent methylation were adjusted for child’s age, sex, race/ethnicity, maternal age and education.

The distribution of GSTP1 percent methylation was right-skewed (median 2.26, IQR 1.56–3.02). We hypothesized that high percent methylation (hypermethylation) would be associated with impaired lung function. Therefore, GSTP1 percent methylation was a priori dichotomized into high [fourth quartile; ≥3.02 (n = 33)] versus low [< 3.02 (n = 98)]. Multivariable linear regression models were employed to examine the association between low versus high GSTP1 methylation on PFT z-scores with additional adjustment for maternal age and education. Sex-stratified analyses were also explored. In sensitivity analysis, we also further adjusted for pre- and postnatal smoking, as well as children’s asthma which is likely on the pathway between air pollution and lung function.

BDLIM analyses were implemented in R statistical software (v3.3.1, Vienna, Austria) and descriptive analyses as well as linear regression models were implemented in SAS (v9.4, SAS Institute, Cary, NC).

Results

Table 1 summarizes participant characteristics. The distribution of covariates among those included in the spirometry analysis and those in the GSTP1 analysis were similar. Most mothers were ethnic minority (predominantly Hispanic, followed by Black), had ≤12 years of education (~ 60%), and never smoked (~ 70%). Average prenatal PM2.5 levels were similar for both groups and there were no significant differences in terms of maternal age at enrollment, race/ethnicity, education, pre- and postnatal smoking, and child’s asthma status (Table 1).

Table 1

ACCESS participant characteristics

	Included in spirometry analysis		Included in GSTP1 analysis
	(n = 171)		(n = 131)
Prenatal PM2.5 level (μg/m3; median, IQR) a	10.9	10.2─11.7	11.0	10.2─11.8
Child sex (n, %)
Girls	82	48.0	63	48.1
Boys	89	52.1	68	51.9
Race/Ethnicity (n, %)
Black	39	22.8	28	21.4
Hispanic	112	65.5	85	64.9
White/Other	20	11.7	18	13.7
Maternal education (n, %)
> 12 yrs	57	33.3	53	40.5
≤ 12 yrs	114	66.7	78	59.5
Child age at spirometry measure (yr; mean, SD)	6.9	0.8	6.9	0.8
Maternal age at enrollment (yr; mean, SD)	27.2	5.7	27.6	5.7
Pre- and postnatal smoking status (n, %) b
Never smoked	122	71.4	95	72.5
Smoked prenatally, but not postnatally	14	8.2	11	8.4
Did not smoke prenatally, but smoked postnatally	24	14.0	16	12.2
Smoked both pre- and postnatally	11	6.4	9	6.9
Child asthma status (n, %)
No	144	84.2	109	83.2
Yes	27	15.8	22	16.8
Spirometry outcomes
FEV1 raw value (L; mean, SD)	1.44	0.25	1.42	0.25
FVC raw value (L; mean, SD)	1.58	0.29	1.55	0.28
FEV1/FVC ratio (mean, SD)	0.92	0.05	0.92	0.05
FEF25–75 raw value (L/s; mean, SD)	1.87	0.46	1.88	0.46
z-score of FEV1 (mean, SD) c	0.01	0.99	0.00	1.03
z-score of FVC (mean, SD) c	−0.04	0.99	−0.06	1.02
z-score of FEV1/FVC ratio (mean, SD) c	0.04	1.00	0.06	0.98
z-score of FEF25–75 (mean, SD) c	0.08	0.96	0.07	0.98

aAveraged over entire pregnancy

bCombination of prenatal maternal smoking and postnatal household smoking status

cAdjusted for age, sex, height, race

Effects of prenatal PM2.5 exposure on child spirometry

In the sample considered as a whole, BDLIMs identified a statistically significant sensitive window of exposure (36–39 weeks gestation) during which children exposed to higher prenatal PM2.5 had an increased risk of reduced FEV1 z-score, after adjusting for maternal age and education and z-score adjustment for child age, sex, height and weight (Fig. 1a). In sex-stratified analyses, we observed a similar significant sensitive window of prenatal PM2.5 exposure at 35–40 weeks gestation in boys whereas no significant sensitive window was seen in girls (Fig. 1b). The results from BDLIMs suggested that the interaction between prenatal PM2.5 and sex was attributable to both different sensitive windows and different within-window effects (the normalized posterior density was 0.86 for the model with this assumption, which can be interpreted as a probability that this was the best fitting pattern of effect modification). In order to further assess the time-weighted associations over the entire pregnancy, we also estimated the cumulative effects accounting for identified sensitive windows and within-window associations. The estimated cumulative effect of prenatal PM2.5 per 1 μg/m3 increase in PM2.5 was significant for boys (cumulative effect estimate = − 0.10, 95%CI = − 0.19 to − 0.01), but not in girls (Table 2).

Fig. 1

Associations between weekly PM2.5 levels over gestation and FEV1 z-score. This figure demonstrates the association between PM2.5 exposure over pregnancy and FEV1 z-scores using a BDLIM assuming week-specific effects, for (a) overall sample, and (b) interaction by sex. Models were adjusted for maternal age and education, and z-scores were adjusted for child’s age at spirometry test, sex, race/ethnicity, and height. The y-axis represents the change in FEV1 z-scores corresponding to a 1 μg/m3 increase in PM2.5; the x-axis is gestational age in weeks. Solid lines show the predicted change in FEV1 z-score. Gray areas indicate 95% confidence intervals (CIs). A sensitive window is identified for the weeks where the estimated pointwise 95% CI (shaded area) does not include zero

Table 2

Estimated cumulative effects of prenatal PM2.5 exposure over gestation: accounting for interaction by sex, sensitive windows and within-window effects identified by BDLIM

Spirometry (z-scores)a	Overall			Interaction Model
				Girls			Boys
	Cumulative Effect	95% CI		Cumulative Effect	95% CI		Cumulative Effect	95% CI
FEV1	− 0.05	− 0.13	0.02	0.10	0.00	0.18	− 0.10	− 0.19	− 0.01
FVC	− 0.09	− 0.17	− 0.02	0.03	− 0.05	0.14	− 0.12	− 0.20	−0.01
FEV1/FVC	0.03	−0.05	0.09	0.02	−0.04	0.11	0.02	−0.04	0.10
FEF25–75	0.03	−0.05	0.09	0.08	−0.01	0.17	0.00	−0.07	0.07

aModels were adjusted for maternal age and education, and z-scores were adjusted for child’s age at spirometry test, sex, height, race/ethnicity; corresponding to per 1 μg/m3 increase in PM2.5

When considering the sample as a whole, the estimated cumulative effect of PM2.5 exposure across pregnancy on FVC, accounting for time-varying effects determined by BDLIMs was also significant (cumulative effect estimate = − 0.09 per 1 μg/m3 increase in PM2.5, 95%CI = − 0.17 to − 0.02; Table 2). BDLIMs did not identify a statistically significant sensitive window of exposure between prenatal PM2.5 and FVC in the overall sample, however models did identify a sensitive window when also considering infant sex (Fig. 2a). Specifically, the BDLIM accounting for sex demonstrated that boys born to mothers with increased prenatal PM2.5 exposure at 36–40 weeks gestation were at increased risk of having reduced FVC z-scores (Fig. 2b), and the estimated cumulative effect was also statistically significant in boys (cumulative effect estimate = − 0.12 per 1 μg/m3 increase in PM2.5, 95%CI = − 0.20 to − 0.01; Table 2). No significant window or cumulative effect was found among girls.

Fig. 2

Associations between weekly PM2.5 levels over gestation and FVC z-score. This figure demonstrates the association between PM2.5 exposure over pregnancy and FVC z-scores using a BDLIM assuming week-specific effects, for (a) overall sample, and (b) interaction by sex. Models were adjusted for maternal age and education, and z-scores were adjusted for child’s age at spirometry test, sex, race/ethnicity, and height. The y-axis represents the change in FVC z-scores corresponding to a 1 μg/m3 increase in PM2.5; the x-axis is gestational age in weeks. Solid lines show the predicted change in FVC z-score. Gray areas indicate 95% confidence intervals (CIs). A sensitive window is identified for the weeks where the estimated pointwise 95% CI (shaded area) does not include zero

For the FEV1/FVC ratio, BDLIMs identified a statistically significant sensitive window of exposure to prenatal PM2.5 at 34–36 weeks gestation and reduced FEV1/FVC z-score. However, the overall cumulative effect was not significant (Table 2), and no significant effect modification by child sex was found. For FEF25–75, BDLIMs did not identify a statistically significant sensitive window of exposure between prenatal PM2.5 and reduced FEF25–75 z-score; the estimated cumulative effect of PM2.5 exposure across pregnancy on FEF25–75 z-score was not significant (Table 2), nor was effect modification by child sex.

Effect of prenatal PM2.5 exposure on GSTP1 percent methylation

We found a statistically significant sensitive window of exposure between prenatal PM2.5 at 37–40 weeks gestation and increased GSTP1 percent methylation (Fig. 3), although the overall estimated cumulative effect across pregnancy was not significant (cumulative effect estimate =0.01 per 1 μg/m3 increase in PM2.5, 95%CI = − 0.78 to 0.76). No sex-specific effects were identified.

Fig. 3

Associations between weekly PM2.5 levels over gestation and GSTP1 percent methylation. This figure demonstrates the association between PM2.5 exposure over pregnancy and GSTP1 DNA methylation using a BDLIM assuming week-specific effects. The model was adjusted for maternal age, education, child’s age at spirometry test, sex, and race/ethnicity. The y-axis represents the change in GSTP1 percent methylation corresponding to a 1 μg/m3 increase in PM2.5; the x-axis is gestational age in weeks. Solid lines show the predicted change in GSTP1 percent methylation. Gray areas indicate 95% confidence intervals (CIs). A sensitive window is identified for the weeks where the estimated pointwise 95% CI (shaded area) does not include zero

Associations between GSTP1 percent methylation and child spirometry

In the sample as a whole, high percent methylation was associated with decreases in FEV1 (z-score = − 0.38, p = 0.06) and FEF25–75 (z-score = − 0.38, p = 0.05), adjusting for maternal age and education (Table 3, Model 1) although associations did not quite reach statistical significance at the p < 0..05 level. Sensitivity analysis with additional adjustment for child asthma and tobacco smoke exposure did not significantly change the point estimates (Table 3, Models 2 and 3). In sex-stratified analyses, associations between high GSTP1 percent methylation and reduced FEV1 (z-score = − 0.56, p = 0.05) and FEF25–75 (z-score = − 0.54, p = 0.06) were only seen among boys and associations were again borderline significant (Table 4). Interaction terms did not approach significance.

Table 3

Associations between GSTP1%methylationa (%m) and childhood PFT z-scoresb: Linear regression

Spirometry (z-scores)a	Univariate Model			Multivariable-adjusted Modelsc
				Model 1d			Model 2e			Model 3f
	β	s.e.	p	Β	s.e.	p	β	s.e.	p	β	s.e.	p
FEV 1
GSTP1%m low	Ref	–	–	Ref	–	–	Ref	–	–	Ref	–	–
GSTP1%m high	−0.38	0.20	0.06	−0.37	0.20	0.06	−0.40	0.20	0.05	−0.41	0.20	0.04
FVC
GSTP1%m low	Ref	–	–	Ref	–	–	Ref	–	–	Ref	–	–
GSTP1%m high	−0.20	0.20	0.32	−0.19	0.21	0.36	−0.21	0.21	0.29	−0.22	0.21	0.27
FEV 1 /FVC
GSTP1%m low	Ref	–	–	Ref	–	–	Ref	–	–	Ref	–	–
GSTP1%m high	−0.27	0.19	0.16	−0.27	0.19	0.16	−0.27	0.20	0.16	−0.27	0.20	0.16
FEF 25–75
GSTP1%m low	Ref	–	–	Ref	–	–	Ref	–	–	Ref	–	–
GSTP1%m high	−0.37	0.19	0.05	−0.38	0.19	0.05	−0.40	0.19	0.04	−0.39	0.19	0.04

aGSTP1 average methylation ≥3.02 [n = 33 (25%)] vs < 3.02; dichotomized around 4th quartile

bPFT z-scores adjusted for age, sex, height, race

cMultivariable-adjusted linear regressions (Models 1–3) predicting PFT z-scores (dependent variables)

dModel 1 additionally adjusted for maternal age and education

eModel 2 (sensitivity model) additionally adjusted for child asthma

fModel 3 (sensitivity model) additionally adjusted for child asthma and tobacco smoke exposure

Table 4

Sex-stratified associations between GSTP1%methylationa (%m) and childhood PFT z-scoresb: Linear regressionc

Spirometry (z-scores)	Girls (low n = 65; high, n = 16)			Boys (low n = 67, high n = 17)			p-interaction
	Β	s.e.	p	β	s.e.	p
FEV 1
GSTP1%m low	Ref	–	–	Ref	–	–
GSTP1%m high	−0.26	0.29	0.38	−0.56	0.29	0.05	0.37
FVC
GSTP1%m low	Ref	–	–	Ref	–	–
GSTP1%m high	−0.14	0.28	0.61	−0.31	0.31	0.32	0.53
FEV 1 /FVC
GSTP1%m low	Ref	–	–	Ref	–
GSTP1%m high	−0.15	0.23	0.51	−0.38	0.33	0.26	0.67
FEF 25–75
GSTP1%m low	Ref	–	–	Ref	–	–
GSTP1%m high	−0.26	0.27	0.34	−0.54	0.28	0.06	0.54

aGSTP1 average methylation ≥3.02 vs < 3.02; dichotomized around 4th quartile

bPFT z-scores adjusted for age, sex, height, race

cMultivariable-adjusted linear regressions additionally adjusted for maternal age and education predicting PFT z-scores (dependent variables)

Discussion

This is the first prospective study to combine advanced statistical modeling with highly resolved ambient PM2.5 exposure estimates to objectively determine susceptible windows of exposure between prenatal PM2.5 exposure and early childhood lung function. Children born to mothers exposed to higher PM2.5 in late pregnancy (≥ 35 weeks gestation) were at increased risk for impaired lung function with male children in particular being at increased risk. Higher prenatal PM2.5 exposure was also associated with increased GSTP1 methylation in children’s nasal cell DNA with models showing a similar window of susceptibility (37–40 weeks gestation). Linear regression models showed an inverse relationship between increased GSTP1 DNA methylation and reduced lung function, specifically amongst boys. These findings combined with prior literature showing that GSTP1 methylation may decrease GSTP1 expression and increase susceptibility to oxidative stress [19], suggest that epigenetic regulation of GSTP1 may be one pathway underlying the association between prenatal PM2.5 exposure and reduced lung function outcomes in early childhood.

A more refined temporal understanding of the impact of prenatal PM2.5 exposure on lung function can better delineate underlying mechanisms. Only one prior study has tried to address effects of exposure timing across pregnancy in this context. Morales and colleagues investigated associations between average trimester benzene and nitrogen dioxide (NO2) exposure and child lung function at age 4.5 years, demonstrating an inverse relationship between second trimester benzene and NO2 exposures and lung function [37]. Our study adds to the literature by investigating a different pollutant, PM2.5. Further, studies considering relatively arbitrary assignment of exposure windows (e.g., averaged over pregnancy or within clinically defined trimesters) rather than being grounded in an understanding of well-orchestrated and timed developmental processes relevant to lung development progressing over gestation may introduce bias or miss effect altogether [7]. By using data-driven methods that incorporate both current and prior gestational exposure, we were able to examine associations between prenatal PM2.5 exposure and child lung function and identify an unbiased estimate of the window of susceptibility on lung function. Results suggest that the effect of PM2.5 on future child lung function largely occurs during the saccular and alveolar phases of lung development. During these phases, the lung undergoes rapid growth and remodeling, including secondary septation of primitive saccules into alveoli, to support efficient gas exchange [38]. Alveologenesis is characterized by proliferation of alveolar type II (ATII) cells, stem cells for alveolar type I cells that line the alveolar surface and produce surfactant. Concurrent elastogenesis spatially instructs future bud formation and it, along with angiogeneisis, is essential to alveolar septation. While our study is the first prospective human study to identify a specific window of susceptibility to prenatal PM2.5 for lung function, mouse models of prenatal PM2.5 on lung function similarly found that pups exposed to in utero PM2.5 had impaired lung function and significantly altered alveolar structure and elastic properties [4].

Lung development and subsequent function is determined by a large number of biological pathways. While heritable genetics play a role, the largest genome-wide association study suggests that 26 loci and > 100 variants may explain 7.5% of FEV1/FVC variance and 3.4% of FEV1 variance [39] highlighting the need to consider gene by environment interactions. Environmental exposures can induce epigenetic modifications, including DNA methylation, that alter gene expression without changing the underlying genetic code [20]. Indeed, one point of epigenome establishment is during the prenatal period and evidence suggests that these changes may persist over the lifecourse or even through generations [10, 12, 40, 41].

Air pollution exposure generates reactive oxygen species (ROS) and the resulting oxidant imbalance is thought to be a central pathway mediating the association between air pollution exposure and impaired lung function [42, 43]. GSTP1 controls enzymes involved in the detoxification of ROS [44] - GSTP1 variants are more susceptible to the health effects of air pollution, including lung function, [45] and GSTP1 methylation is associated with oxidative stress [46]. The developing fetus is especially prone to oxidative injury, as antioxidant defenses remain immature through pregnancy; male fetuses may be at even increased risk. In vitro work exposing ATII cells to PM2.5 demonstrates an increase in ROS with subsequent change in ATII phenotype [47]. Pups of pregnant rats exposed to PM2.5 during gestation demonstrate impaired lung function with alveolar destruction and thickened alveoli septum demonstrated on histology; maternal PM2.5 exposure correlated with increasing levels of pup oxidative stress and markers of oxidative stress in pup lung tissue [48]. While our study demonstrates that prenatal PM2.5 is associated with increased GSTP1 methylation during the same window of sensitivity as reduced lung function and that elevated GSTP1 methylation is associated with reduced lung function, future studies with larger sample size are needed to corroborate these findings and to more formally examine mediation.

We note several strengths to our study. We prospectively captured prenatal environmental exposures, a number of important confounders and covariates, and lung function - an objective measure of respiratory outcomes in early childhood - an urban cohort largely composed of ethnic minorities at increased risk for ambient air pollution exposure and impaired lung function. We assessed prenatal PM2.5 exposure using a validated hybrid spatiotemporal LUR model incorporating highly-resolved satellite-derived AOD measures. We then objectively identified susceptible windows of exposure for PM2.5 on child lung function by using data-driven, advanced statistical techniques to more objectively determine sensitive windows of effect. Also, this is the first study to combine these methods with gene-specific methylation to explore mechanisms mediating the association between prenatal air pollution and child lung function. Also, while studies indicate that tobacco smoke exposure is a major contributor to indoor air pollution in inner-city homes [49, 50], variations in indoor source particles (e.g., from smoking and/or cooking sources) are largely uncorrelated with variations in outdoor source particles, which is the main exposure considered herein [51]. Thus, while particles of indoor origin are an important predictor of fetal development in and of themselves, they are unlikely to confound associations between ambient particulate matter and child lung function. Indeed, investigators have begun to conceptualize tobacco smoke exposure as an effect modifier rather than a confounder in this context [52]. Observed effects remained significant even in sensitivity analyses adjusting for prenatal maternal smoking and second-hand postnatal tobacco smoke exposure, another major determinant of childhood lung function [53].

We also acknowledge limitations. While we adjust for several confounding factors, we did not have data on dietary or other environmental factors, such as temperature, which may co-vary with PM2.5. However, as with smoking exposures, maternal diet and temperature are increasingly conceptualized as effect modifiers in the context of health effects of particulate matter rather than confounders [54, 55]. This should be explored in future studies. Our epigenetic analyses involve specimens collected at the time of lung function testing; future studies with serially collected airway epithelial samples will allow formal mediation analyses and expand our understanding of the influence of prenatal environmental exposures on, and stability of, the epigenome. Also, while previous studies show that GSTP1 methylation decreases GSTP1 expression [19, 56], we were unable to measure GSPT1 expression due to logistical limitations that made collecting nasal cell RNA unfeasible. Future studies that include more direct measurement of oxidative stress could more definitively assess whether the observed associations are mediated through oxidative stress pathways. Given that our study examines these effects in an ethnically diverse population, additional studies are warranted to generalize our results to other populations. Finally, as our statistical modeling techniques evolve, it will be important to note that co-occurring pollutants may influence different periods of lung development and interact, thus highlighting the need for future multi-pollutant study.

In summary, we use data-driven statistical techniques to demonstrate that increased PM exposure during the final phase of gestation may increase risk for future impaired lung function, possibly mediated through epigenetic regulation of GSTP1.

Increased prenatal PM2.5 exposure in later pregnancy (≥ 35 weeks gestation) was associated with increased risk for impaired lung function in 7 year olds and increased GSTP1 methylation. Boys were especially susceptible. Elucidating timing of exposure may better inform mechanistic underpinnings and identify those at heightened risk.

Conclusions

Data-drive statistical techniques demonstrate that prenatal PM2.5 exposure in late pregnancy was associated with impaired early childhood lung function and DNA hypermethylation of nasal epithelia GSTPI, especially in boys. Higher GSTP1 percent methylation was significantly associated with reduced FEV1. Future studies with larger sample size are needed to corroborate these findings and more formally examine mediation of the link between prenatal air pollution exposure and childhood lung function by GSTP1.