Air pollution exposure and gestational diabetes mellitus among pregnant women in Massachusetts: a cohort study

Background Rodent and human studies suggest an association between air pollution exposure and type 2 diabetes mellitus, but the extent to which air pollution is associated with gestational diabetes mellitus (GDM) is less clear. Methods We used the Massachusetts Registry of Vital Records to study primiparous women pregnant from 2003-2008 without pre-existing diabetes. We used satellite-based spatiotemporal models to estimate first and second trimester residential particulate (PM2.5) exposure and geographic information systems to estimate neighborhood traffic density. We obtained GDM status from birth records. We performed logistic regression analyses adjusted for sociodemographics on the full cohort and after stratification by maternal age and smoking habits. Results Of 159,373 women, 5,381 (3.4 %) developed GDM. Residential PM2.5 exposure ranged 1.3–19.3 μg/m3 over the second trimester. None of the exposures were associated with GDM in the full cohort [e.g. OR 0.99 (95 % CI: 0.95, 1.03) for each interquartile range (IQR) increment in second trimester PM2.5]. There were also no consistent associations after stratification by smoking habits. When the cohort was stratified by maternal age, women less than 20 years had 1.36 higher odds of GDM (95 % CI: 1.08, 1.70) for each IQR increment in second trimester PM2.5 exposure. Conclusions Although we found no evidence of an association between air pollution exposure and GDM among all women in our study, greater exposure to PM2.5 during the second trimester was associated with GDM in the youngest age stratum. Electronic supplementary material The online version of this article (doi:10.1186/s12940-016-0121-4) contains supplementary material, which is available to authorized users.


Background
Gestational diabetes mellitus (GDM) complicates 2-6 % of pregnancies worldwide and as many as 10-20 % of pregnancies in high-risk populations [1]. GDM increases risk of adverse perinatal outcomes such as fetal hypoglycemia and birth trauma to mother and infant. GDM additionally primes infants and mothers for higher risk of cardiometabolic disease later in life [2]. Maternal characteristics such as obesity, older age, and family history of type 2 diabetes mellitus are known to increase risk of GDM, but up to half of women with GDM do not have these classic determinants [3], suggesting a role for environmental factors. Identification of environmental triggers is critical to target at-risk women and provide opportunities for prevention.
Air pollution is one remediable environmental trigger that may predispose susceptible pregnant women to GDM. Automobiles and power plants emit both gasses and particulate air pollutants. The smallest of these particules, with an aerodynamic diameter less than 2.5 μm (PM 2.5 ), are readily inhaled. Experimental PM 2.5 exposure increases insulin resistance in rodents [4,5] through endothelial dysfunction, inflammation, and oxidative stress [6]. Consistent with the rodent findings, several epidemiologic studies have demonstrated an association between higher air pollution exposure and increased risk of type 2 diabetes mellitus [reviewed in [7]].
However, epidemiologic analyses of prenatal air pollution exposure and abnormal glucose tolerance in pregnancy are conflicting with most [8][9][10][11] but not all [12] showing an association. Two of the prior studies were limited by use of birth cohorts with relatively small sample sizes and few cases of GDM [8,12] and another did not include individual-level socioeconomic status covariate data [9].
In the present analysis, our objective was to use the Massachusetts Registry of Vital Records to evaluate the extent to which first and second trimester residential exposure to PM 2.5 and neighborhood traffic density were associated with GDM in a large cohort of pregnant women. We hypothesized that prenatal air pollution exposure would be associated with GDM.

Study population and design
We obtained data on registered live births in Massachusetts from January 1, 2003 through December 31, 2008 and latitude and longitude of each residential address at the time of delivery from the Massachusetts Registry of Vital Records and Statistics (http://www.mass.gov/ eohhs/gov/departments/dph/programs/admin/dmoa/ vitals/). We obtained data only for births associated with a Massachusetts residential address. Daily PM 2.5 exposure estimates were available as of January 1, 2003, so we included mothers in this analysis whose last menstrual period (LMP) occurred on or after January 1, 2003, which enabled us to create the early pregnancy PM 2.5 exposure estimates detailed below. Of the 362,148 women who met these inclusion criteria, we restricted the dataset to 161,144 nulliparous women to ascertain only independent observations in our analysis, as the Massachusetts Registry of Vital Records provides a unique ID for each birth rather than each woman. We further restricted our analysis to mothers who delivered at greater than 28 weeks gestation and thus had opportunity for GDM screening, and who did not have a prior history of diabetes. Our final sample included 159,373 women (Fig. 1). The study was approved by the Massachusetts Department of Public Health and the Institutional Review Boards of the participating institutions.

Exposure assessment
We used each woman's residential address at the time of delivery to estimate daily PM 2.5 exposure throughout the pregnancy. To create PM 2.5 exposure estimates, we used a hybrid satellite-based spatiotemporal model developed by our group, which has been previously described in detail [13]. This model incorporates aerosol optical depth data from the MODIS Satellite and classic land use regression techniques to generate daily PM 2.5 exposure estimates at the resolution of a 1x1km spatial grid. To create residential exposure estimates, each residential address was linked to the grid in which it was located. The mean daily "out-of-sample" ten-fold cross-validation of the model was excellent (R 2 = 0.88).
To obtain first trimester exposure estimates, we averaged daily exposures from the date of LMP through 12th week of gestation. To obtain second trimester exposure estimates, we averaged daily exposures from the 13th through 24th week of gestation. To create these time windows, we used the birth certificate clinical estimate of gestational age rather than the estimate calculated based on reported LMP because the former has been shown to be a more accurate predictor of gestational age in birth registries [14], and its use is recommended by the American College of Obstetrics and Gynecology for public health research studies [15].
Our dataset included an estimate of residential traffic density [annual average daily traffic (vehicles/day) times length of road (km) within 100 m of the participants' residential address at the time of birth] from the Massachusetts Executive Office of Transportation 2002 road inventory, which we used for consistency with prior analyses of this cohort [16,17]. While our daily estimates of PM 2.5 were temporally and spatially resolved, estimates of neighborhood traffic density were only spatially resolved.

Outcome assessment
Routine clinical screening for GDM is recommended for pregnant women in Massachusetts at the end of the second trimester of gestation (24-28 weeks). If serum glucose one hour after a non-fasting 50 g oral glucose challenge test (GCT) is ≥ 140 mg/dL, women are referred for a 3-hour fasting 100-g oral glucose tolerance test (OGTT). As per American Diabetes Association (ADA) criteria, pregnant women are classified as having GDM if they have ≥ 2 of the following abnormal values on the OGTT: blood glucose > 95 mg/dL at baseline, > 180 mg/dL at 1 h, > 155 mg/dL at 2 h, or > 140 mg/dL at 3 h [18]. For this analysis, we obtained maternal GDM designation from the birth records.

Assessment of covariates
We also retrieved data on other maternal characteristics (age, race/ethnicity, education, smoking status, prenatal insurance, Kotelchuck prenatal care index) as well as infant sex and date of birth directly from birth records. For consistency with prior analyses of this cohort [16,17], we used data from the 2000 United States Census [19] to calculate median annual household income, percent open space, and median value of owneroccupied housing for the census tract associated with each residential address at the time of birth. We abstracted median annual household income and median value of owner-occupied housing directly from the Census, and we calculated percent open space by intersecting 2000 Census tract boundaries with information on land use for recreation and conservation from MassGIS [20].

Statistical analysis
We used logistic regression analyses to evaluate the associations of PM 2.5 exposure and traffic density with GDM. We considered each exposure (first trimester PM 2.5 , second trimester PM 2.5 , and traffic density) in separate models. We first modeled exposures in quartiles to assess for potential non-linearity of exposure-outcome relationships. As we did not see any clear evidence of non-linearity, we also modeled each exposure as a continuous measure and expressed associations per 10-90 percentile increment in exposure. We first fit unadjusted models, followed by full multivariable models for each of the exposures. We included the following covariates potentially associated with air pollution exposure [21,22] and/or GDM [1,23]: maternal age (<20 years, 20-35 years, ≥ 35 years), race/ethnicity (white, black, Asian/ Pacific Islander, Hispanic, other), education (less than high school, high school, some college, bachelor degree, postgraduate degree), smoking habits [never, former, current low (<10 cigarettes/day), current high (>10 cigarettes per day)], and prenatal insurance (public versus private); census tract median household income (continuous), percent open space (continuous), and median value of owner occupied housing (continuous); season of birth [continuous sine and cosine of date as in [24]]; and date of birth (continuous). Additional adjustment for infant sex and Kotelchuck prenatal care index [25] did not change associations; thus, we did not include these variables in our final models.
We also examined associations between air pollution and GDM within strata of maternal age (<20 years, 20-35 years, ≥ 35 years) in light of data suggesting that children and young adults have higher susceptibility to air pollution [e.g. in relation to respiratory outcomes [26]], presumably due to higher respiratory minute volume, activity level, and time spent outdoors [27]. Within strata of maternal age, we adjusted for maternal age as a continuous variable. In addition, due to evidence for a synergistic interaction between air pollution and smoking in relation to outcomes such as respiratory health [28,29] and obesity [30,31], we also examined associations between air pollution and GDM within strata of maternal smoking habits [never, former, current low (<10 cigarettes/day), current high (>10 cigarettes per day)]. We also performed several sensitivity analyses. First, we restricted the sample to women who delivered singleton infants (n = 155,501) for consistency with prior studies that examined the relationship between air pollution and GDM among singleton pregnancies only [8][9][10][11]. Next, as age appeared to be a strong and nonlinear negative confounder in this population, with older women living in less polluted areas but having a substantially increased risk of GDM (Table 1), we represented age as a cubic (age, age 2 , age 3 ) rather than a categorical term in the multivariable model to maximize adjustment for this covariate. In addition, because restricting to nulliparous women decreased our sample size so substantially, we included both nulliparous and multiparous women (n = 358,053) in multivariable analyses adjusted for parity in addition to the covariates listed above, with the knowledge that this larger dataset included non-independent observations. Finally, to further investigate the association we identified between second trimester PM 2.5 exposure and GDM in the youngest mothers, we compared the exposures and sociodemographics of young mothers with versus without GDM. All analyses were conducted using SAS Version 9.4 (SAS Institute Inc, Cary, NC).  (Table 2), consistent with the fact that PM 2.5 emissions in Massachusetts are primarily from regional rather than local traffic sources [32].
Mothers with higher residential PM 2.5 exposure during the second trimester were more likely to be younger, less educated, and nonsmokers. They were also more likely to have public prenatal insurance and live in a census tract with a lower median household income, less open space, and higher median value of owner-occupied housing (Table 1). Women with GDM (vs. those without) were more likely to be Asian/Pacific Islander. They were also more likely to be older and have private prenatal insurance but live in a census tract with lower median household income, less open space, and lower median value of owner-occupied housing ( Table 1). The distribution of characteristics with first trimester air pollution exposure was similar.

Primary analyses
Residential PM 2.5 exposure during the first and second trimesters and neighborhood traffic density were not associated with higher odds of GDM in unadjusted or covariate-adjusted analyses. Unadjusted odds of GDM for women in the highest (Q4) (vs. lowest (Q1)) quartile  (Table 3).

Stratified analyses
When we examined the association between air pollution exposure and GDM within strata of maternal age, we observed an association between second trimester residential PM 2.5 exposure and GDM in mothers who were less than 20 years of age at the time of delivery (Table 4). Within this group of women, those who lived at a residence in the highest (Q4) versus lowest (Q1) quartile of PM 2.5 exposure during the second trimester had 1.97 (95 % CI: 1.17, 3.32) times the odds of developing GDM. Within this stratum, odds of GDM was consistently higher in Q2, Q3, and Q4 versus Q1 of second trimester PM 2.5 exposure, and for each 10-90 percentile increment in exposure, odds of GDM was 1.76 times higher (95 % CI: 1.16, 2.7). In mothers greater than 35 years of age at the time of delivery, those in the highest versus lowest quartile of PM 2.5 exposure during the first trimester had 1.18 (95 % CI: 1.00, 1.39) times the odds of developing GDM, although Q2, Q3, and Q4 versus Q1 odds ratios did not increase monotonically, and, in fact, the Q3 versus Q1 comparison was close to 1. There were no other associations between prenatal pollution and GDM within the age strata (Table 4).
Residential PM 2.5 exposure and traffic density were not consistently associated with GDM within strata of maternal smoking habits. For example, for each 10-90 percentile increment in second trimester PM 2.5 exposure, odds of GDM was 0.97 (95 % CI: 0.89, 1.05) for women who never smoked, 0.90 (95 % CI: 0.70, 1.16) for women who smoked prior to pregnancy, 1.37 (95 % CI:

Sensitivity analyses
When we restricted the sample to women with singleton pregnancies and when we represented age as a cubic rather than a categorical term, there was no change to the pattern of results (data not shown). When we included multiparous women, there was no change to the pattern of results as compared to our primary analyses with the exception of the following: (1) we found no association between first trimester PM 2.5 and GDM in mothers greater than 35 years of age, and (2)

Discussion
In our analysis of Massachusetts birth registry data from 2003-2008, pregnant women with high residential PM 2.5 exposure during the first or second trimester or high neighborhood traffic density had the same odds of developing GDM as women with lower exposures. When we examined this association within strata of maternal age, the youngest mothers (<20 years of age) had increased odds of GDM when exposed to higher residential PM 2.5 during the second trimester. Our findings are consistent with prior populationbased studies that have shown an association between GDM and prenatal exposure to NO x or ozone [9][10][11] but found weaker or no associations with exposure to PM 2.5 [8,10,11] or traffic density [8,9,12]. The overall lack of a consistent association between PM 2.5 exposure and GDM in human observational studies is in contrast to a growing body of epidemiologic literature showing an association between PM 2.5 and type 2 diabetes mellitus, [reviewed in [7]] and rodent studies confirming an association between PM 2.5 and insulin resistance in nonpregnant adults [4,33]. PM 2.5 is primarily thought to lead to insulin resistance through oxidative damage, endothelial dysfunction, and inflammation [4,34], whereas the specific mechanisms of NO x and ozoneinduced insulin resistance are not as well-understood. In one rodent study, ozone-induced insulin resistance was associated with neuronal activation and sympathetic stimulation [35] but not with the increase in circulating inflammatory cytokines that has been observed following PM 2.5 exposure [4,33]. Pregnant women may be more vulnerable to NO x or ozone-specific mechanisms, and the role of each of these pollutants and mechanisms of action in relation to GDM should be studied in pregnant rodent models.
Although prenatal residential PM 2.5 exposure and traffic density were not associated with GDM in analyses of our complete study population, we found higher second trimester PM 2.5 exposure (i.e.-quartiles 2-4 vs. the lowest quartile) to be associated with almost twice the odds of GDM in the subset of women less than 20 years of age. Second trimester exposures are biologically Vehicles/day x km of road within 100 m of residential address relevant, as GDM is a pathologic exacerbation of a physiologic increase in insulin resistance that occurs specifically during the second trimester of pregnancy [36]. Also, the magnitude of our finding is similar to that of other well-known risk factors for GDM. For example, in a meta-analysis, overweight (vs. normal weight) mothers had 1.83 (95 % CI: 1.58, 2.12) times the odds of developing GDM, and those who were obese had 3.52 (95 % CI: 3.24, 3.84) times the odds [37]. Our finding of an association between air pollution and GDM in the youngest mothers may be because maternal age is such a strong risk factor for GDM [1] that older women destined to develop GDM will do so regardless of other factors, whereas young women are more readily influenced by environmental exposures such as air pollution. Additionally, as compared to older adults, children and young adults have been shown to be more susceptible to PM 2.5 -induced health outcomes [26], presumably as a result of higher respiratory minute volume, activity level, and time spent outdoors [27].
We also considered alternative explanations for our finding of an association between PM 2.5 exposure and GDM in the youngest mothers. For example, the possibility that stratification could have reduced negative confounding by age is less likely because PM 2.5 was not associated with GDM in the other age strata, and the association in the full cohort was null even when we included age as a cubic term. Alternatively, our inability to account for pre-pregnancy BMI, which was not recorded on Massachusetts birth certificates until 2011, could have confounded the association in the youngest mothers. However, if there were counfounding by prepregnancy BMI, we would have expected this to have affected all mothers equally, rather than observing it in the youngest age strata only, based on a prior analysis of >170,000 women in which the association between BMI and GDM did not vary by maternal age [38]. Another possible explanation is differential composition of residential PM 2.5 exposure by maternal age (e.g. higher traffic-related and/or ultrafine particle component of PM 2.5 most prevalent in the lower SES neighborhoods of the youngest mothers and also most closely related to GDM). However, there was no association between traffic density and GDM in our complete study population, Vehicles/day x km of road within 100 m of residential address and we found lower rather than higher neighborhood traffic density in the subset of young women with versus without GDM. It is possible that our finding of an association between PM 2.5 and GDM in the youngest mothers may also reflect random chance, particularly given the relatively small number of cases of GDM (n = 179) in the youngest age stratum. Our findings require replication in other populations of young, pregnant women.
In our sensitivity analysis which included multiparous women, first and second trimester PM 2.5 exposure was associated with slightly increased odds of GDM in the subgroup who smoked less than 10 cigarettes per day during pregnancy. Tobacco smoke contains particulates and has previously been shown to operate synergistically with air pollution to increase risk of obesity, another cardiometabolic outcome [30,31]. However, the fact that this association was not present in the heavier smokers (greater than 10 cigarettes per day) makes it difficult to draw a definitive conclusion and suggests the need for replication in other cohorts.
Use of data from the Massachusetts Registry of Vital Records is a unique strength of our study. The registry contains pregnancy data on all Massachusetts residents, includes a large number of cases of GDM, and is free from the selection bias typical of cohort studies. The registry also contains relatively rich data on socioeconomic status, as well as smoking and diabetes history. Also, our PM 2.5 model, which leveraged satellite aerosol optical depth data, calculated estimates at a fine resolution and had a high mean out-of-sample R 2 .
Potential exposure and outcome misclassification are limitations of the present study that may have biased results toward the null. We did not have information on time-activity patterns or residential moves during pregnancy which could have improved accuracy of exposure estimates. Also, traffic density estimates were for 2002 which may not have been relevant for the whole study period. Outcome misclassification could have occurred as a result of underreporting of GDM on the birth certificate (specificity >98 % and sensitivity 46-83 % when compared to medical records) [39] or because women with undiagnosed type 2 diabetes mellitus may have inappropriately been included in the GDM group. In addition, based on limited information available in the birth registry, we were unable to account for every factor that might be related to pollution exposure and GDM risk, such as physical activity, family history of GDM, and maternal pre-pregnancy BMI. Also, census tract covariates utilized 2000 Census tract data which may not have been relevant for the whole study period, and thus residual confounding may exist. Also, we did not have exposure estimates available in this dataset for other criteria pollutants such as NO x or SO 2 . Overall PM 2.5 exposure was low, and it may be that higher levels of exposure are associated with greater risks.
Additional rodent studies are needed to elucidate the extent to which individual criteria pollutants such as PM 2.5 , NO x , SO 2 , black carbon, and ozone, and mixtures of pollutants, are causally linked to development of GDM and to further investigate mechanisms of action. Large population-based studies with information on multiple criteria pollutants and covariates are also needed, particularly in young cohorts with otherwise low risk of GDM.

Conclusions
Although our findings require replication, in our cohort, young women were at increased risk for GDM when exposed to higher residential PM 2.5 during the second trimester of pregnancy. We otherwise found no consistent evidence of an association between first or second trimester residential PM 2.5 exposure or neighborhood traffic density and GDM in pregnant women in Massachusetts with modest levels of air pollution exposure.