1 Department of Epidemiology, School of Public Health, University of North Carolina, Chapel Hill, NC
2 Human Studies Division, National Health and Environmental Effects Research Laboratory, Environmental Protection Agency, Research Triangle Park, NC
3 Birth Defects Epidemiology and Surveillance Branch, Texas Department of State Health Services, Austin, TX
4 Department of Biostatistics, School of Public Health, University of North Carolina, Chapel Hill, NC
5 Department of Cardiology, Children's Medical Center, Dallas, TX
Reprint requests to Pauline Mendola, Human Studies Division, National Health and Environmental Effects Research Laboratory, Environmental Protection Agency, MD 58A, Research Triangle Park, NC 27711 (e-mail: mendola.pauline{at}epa.gov).
Received for publication October 18, 2004. Accepted for publication March 14, 2005.
![]() |
ABSTRACT |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
abnormalities; air pollution; cleft lip; cleft palate; environment and public health; heart defects, congenital
![]() |
INTRODUCTION |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
A small body of animal toxicology literature suggests that air pollutant exposure can yield adverse reproductive effects (2027
) and potentially heritable gene mutations (28
30
). Other investigators have observed mutations in fetal DNA as a result of exposure to air toxics during pregnancy (31
34
). This suggestive toxicology and the growing epidemiologic evidence for the reproductive toxicity (restricted fetal growth, shortened gestation) of air pollution raise the question of whether air pollution is also an environmental teratogen. One ecologic study found that living in areas with more versus less industrial pollution was associated with higher rates of both spontaneous abortion and congenital anomalies (35
). Similarly, a Ukrainian study that compared mutation rates for 18 different congenital anomalies in three towns representing low, middle, and high pollution levels found higher mutation rates in the more polluted communities (36
). A recent population-based case-control study in southern California that looked at risks of cardiac birth defects and oral clefts associated with exposure to carbon monoxide, nitrogen dioxide, ozone, and particulate matter less than10 µm in aerodynamic diameter (PM10) found associations between 1) carbon monoxide exposure during the second gestational month and an increased risk of isolated ventricular septal defects and 2) ozone exposure during the second gestational month and an elevated risk of isolated aortic artery and valve defects, pulmonary artery and valve defects, and conotruncal defects (37
).
We attempted to corroborate the southern California results by using data from the Birth Defects Epidemiology and Surveillance Branch of the Texas Department of State Health Services. The state maintains a population-based, active surveillance system for major malformations diagnosed before age 1 year. Texas has an extensive network of stationary air pollution monitors placed by the US Environmental Protection Agency. We investigated the association between maternal exposures to carbon monoxide, nitrogen dioxide, ozone, sulfur dioxide, and PM10 during weeks 38 of pregnancy and the risk of selected cardiac birth defects and oral clefts in livebirths and fetal deaths between 1997 and 2000 in seven Texas counties.
![]() |
MATERIALS AND METHODS |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
We analyzed six clinical diagnostic groupings of isolated cardiac and oral cleft birth defects and two clinical diagnostic groupings of multiple cardiac birth defects. The previous study of air quality and birth defects analyzed eight groupings of isolated defects (37); because of small sample sizes, we analyzed multiple rather than isolated conotruncal and endocardial cushion defects. A case with an isolated defect had a cardiac defect or an oral cleft with no other major defects, although the case could have had minor birth defects. Defects were considered minor according to National Birth Defects Prevention Study guidelines (38
). A case with multiple defects had a major cardiac defect or oral cleft and at least one other major defect in either the same organ system or a different organ system. The isolated clinical diagnostic groupings of aortic artery and valve defects, atrial and atrial septal defects, pulmonary artery and valve defects, ventricular septal defects, cleft palate alone, and cleft lip with or without cleft palate, and multiple clinical diagnostic groupings of conotruncal defects and endocardial cushion defects, were created to be comparable with the previous study by Ritz et al. (37
). Table 1 shows the birth defects included in each of the clinical diagnostic groupings, replicating those used in the California analysis. We also devised an alternative classification scheme that was reviewed by a member of our research team (D. E. F.), a pediatric cardiologist. Lastly, we analyzed nine individual birth defects, which included all isolated, multiple, or chromosomal cases of a specific birth defect, in an effort to explore the association in more heterogeneous groupings with greater statistical power. Because we did not use hierarchical birth defect assignments, data on infants and fetuses could be analyzed more than once depending on number of diagnoses.
|
Geocoding of maternal residence
We attempted to geocode maternal residence at delivery of 5,338 cases and 4,574 controls by using the ArcGIS 8.2 and 8.3 mapping program with ESRI StreetMap (Environmental Systems Research Institute Inc., Redlands, California) as the geocoding service. We geocoded 86 percent of cases (n = 4,570) and 80 percent of controls (n = 3,667). We evaluated the characteristics of the geocoded and nongeocoded populations and found no meaningful differences with respect to factors such as maternal age, race/ethnicity, education, cigarette smoking and alcohol consumption during pregnancy, and marital status (data not shown). Cases and controls who could not be geocoded were excluded from the analysis.
Exposure assignment
The Environmental Protection Agency provided raw data for hourly (for gases) or daily (for particulate matter) air pollution concentrations for the seven study counties between January 1, 1996, and December 31, 2000. For gases, we calculated daily monitoring site means and number of hourly measurements contributing to each daily mean. No monitoring data were eliminated nor were any values changed. PM10 and ozone were monitored in all study counties, carbon monoxide and nitrogen dioxide were monitored in six counties (excluding Hidalgo), and sulfur dioxide was monitored in four counties (excluding Bexar, Hidalgo, and Travis).
The distance from each air pollution monitor to each maternal residence at delivery in that county was calculated. Separately for each pollutant, we calculated the average exposure measured at each monitor in the county during weeks 38 of pregnancy. Beginning with the monitor closest to the residence at delivery, we required at least 70 percent of the total number of possible daily means for gases (30/42 daily means) or at least 70 percent of the total number of possible hourly values (
706/1,008 hourly values) to be available for the mean to be assigned to a woman. If these criteria were not met, the next closest monitor in the county was considered until the last monitor measuring the given pollutant was considered. If none of the monitors in the county met the criteria, then no exposure assignment was made. For PM10, the procedure differed slightly because the majority of data was collected only every 6 days; a monitor needed to have at least 50 percent of the total number of possible daily means available (
4/7 daily means) to be assigned to a woman. There were no case-control differences in the proportion of observations assigned exposure estimates: 91 percent for carbon monoxide, 90 percent for nitrogen dioxide, 99 percent for ozone, 55 percent for sulfur dioxide, and 94 percent for PM10.
For each pollutant, at least 85 percent of the calculated means were based on the monitoring data from either the first or the second closest monitor to the maternal residence. The median distances were 8.614.2 km, with minimum distances of 100600 m and maximum distances of 35.554.4 km.
Data analysis
We conducted all analyses with SAS software, version 8.2 (SAS Institute, Inc., Cary, North Carolina). We explored the distribution of covariates and exposure among cases and controls. Considered as potential confounders were alcohol consumption during pregnancy, attendant of delivery (i.e., the person who delivered the baby (physician/nurse-midwife vs. other)), gravidity, marital status, maternal age, maternal education, maternal illness, maternal race/ethnicity, parity, place of delivery, plurality, prenatal care, season of conception, and tobacco use during pregnancy, modeled as shown in tables 2 and 3. These data were obtained from the birth or fetal death certificates for both cases and controls to ensure that data quality and source would be uniform. No covariate was missing for more than 4 percent of the study population. Because we know of only one previous study in this area (37) and no established confounders of these associations, we made no a priori assumptions about which covariates should remain in the multivariable models. The covariates in the final model were determined separately for each combination of an air pollutant and individual birth defect or birth defect grouping. If removal of the covariate resulted in more than a 10 percent change in the estimate of effect, then it was included in the model.
|
|
|
|
![]() |
RESULTS |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Table 4 shows the results for the eight clinical diagnostic birth defect groupings that were reported in the California study (37), and table 5 shows the results of the analyses of the nine individual defects. Odds ratio estimates and 95 percent confidence intervals from unadjusted and covariate-adjusted models were very similar; thus, only adjusted estimates are reported here. Each exposure-outcome model was adjusted for the covariates indicated in the table notes. The majority of analyses did not demonstrate strong evidence of increasing or decreasing risk with increasing exposure.
Clinical diagnostic birth defect groupings
For the clinical diagnostic birth defect groupings (table 4), when we compared the fourth quartile of exposure with the first, carbon monoxide was associated with multiple conotruncal defects (OR = 1.46, 95 percent CI: 1.03, 2.08; p trend = 0.0870), sulfur dioxide was associated with increased risk of isolated ventricular septal defects (OR = 2.16, 95 percent CI: 1.51, 3.09; p trend <0.0001), and PM10 was associated with increased risk of isolated atrial septal defects (OR = 2.27, 95 percent CI: 1.43, 3.60; p trend = 0.0001). PM10 also appeared to be associated with isolated cleft lip with or without cleft palate, although the association was statistically significant in the third quartile of exposure only. Inverse associations were noted for carbon monoxide and risk of isolated atrial septal defects, ozone and risk of isolated ventricular septal defects, sulfur dioxide and risk of isolated atrial septal defects and multiple conotruncal defects, and PM10 and risk of multiple endocardial cushion defects.
Individual birth defects
For individual defects (isolated, multiple, and chromosomal cases combined) (table 5), carbon monoxide exposure was associated with increased risk of tetralogy of Fallot (OR = 2.04, 95 percent CI: 1.26, 3.29; p trend = 0.0017), sulfur dioxide exposure was associated with increased risk of ventricular septal defects (OR = 1.31, 95 percent CI: 1.06, 1.61; p trend = 0.0850), and PM10 exposure was associated with increased risk of atrial septal defects (OR = 1.26, 95 percent CI: 1.03, 1.55; p trend = 0.0096). We also noted several inverse associations: carbon monoxide exposure and atrial septal defects, nitrogen dioxide exposure and ventricular septal defects, and sulfur dioxide exposure and atrial septal defects and ostium secundum. Results were similar, but less precise when we excluded cases with chromosomal anomalies (data not shown).
Simultaneous adjustment for copollutants decreased precision but did not result in meaningful changes in the estimates of the effect of individual pollutants (data not shown). Our analysis did not reveal any evidence of meaningful effect measure modification by infant sex, plurality, maternal education, maternal race, or season of conception (data not shown).
We also modeled these associations by using continuous air pollution variables to calculate a slope per unit change in exposure (data not shown). The results of the continuous models corroborated the quartile findings. For example, a 1-ppm increase in carbon monoxide exposure was associated with increased risk of tetralogy of Fallot (OR = 3.01, 95 percent CI: 1.67, 5.42), and a 10-µg/m3 increase in PM10 exposure was associated with increased risk of isolated atrial septal defects (OR = 1.33, 95 percent CI: 1.11, 1.60).
We conducted two subgroup analyses to assess the robustness of our results to changes in outcome and exposure definitions. First, restricting the analysis to 3,544 heart defect cases (89 percent) who had a cardiac diagnosis confirmed by one of four diagnostic methods (cardiac catheterization, echocardiography, surgery, or autopsy) (37, 38
, 42
) resulted in no remarkable changes in any of the associations with respect to the magnitude, direction, or exposure-response trends (data not shown). The consistency of our results suggests that the group with confirmed cardiac defects did not differ from those with unconfirmed diagnoses. Our second subgroup analysis restricted the study population to those located within 10 km of the assigned monitor for a particular air pollutant. Depending on the defect and pollutant under analysis, this restriction of the study population could result in a sample size reduction of 50 percent or more. Because sample sizes were so much smaller, the associations found in this subgroup analysis were less precise than the results using the full data set, but they were generally consistent with respect to the direction and magnitude of the associations (data not shown).
![]() |
DISCUSSION |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
In a similar study in California, Ritz et al. (37) reported an association between carbon monoxide and isolated ventricular septal defects, which we were unable to replicate. They also reported an association between ozone exposure and increased risk of isolated pulmonary artery and valve defects (OR = 1.36 for quartile 2, OR = 1.42 for quartile 3, and OR = 1.99 for quartile 4). Our odds ratio estimates for this association were of similar magnitude and precision to the California results. Ritz et al. also reported an association between ozone and risk of isolated aortic artery and valve defects (OR = 1.19 for quartile 2, OR = 1.69 for quartile 3, and OR = 2.68 for quartile 4), which we did not see in our data, although the risk of aortic valve stenosis was elevated (OR = 0.97 for quartile 2, OR = 1.13 for quartile 3, and OR = 1.43 for quartile 4). The Ritz et al. California study also noted an association between ozone exposure and conotruncal defects (OR = 1.63 for quartile 2, OR = 1.98 for quartile 3, and OR = 2.50 for quartile 4). We had sufficient data to analyze only the grouping of multiple conotruncal defects (n = 300) and found no association with ozone exposure, although we did find a suggested association with carbon monoxide exposure (table 4). We expected consistency between our ozone results and those from the California study because the ozone levels were comparable. The California study looked at PM10 and reported no associations with any isolated defects. Our study found a noteworthy association between PM10 and isolated atrial septal defects and a suggested association between PM10 and isolated cleft lip with or without cleft palate. The California study did not investigate the effects of sulfur dioxide, providing no source of comparison for our sulfur dioxide results.
A number of factors may explain the differences between our results and those of the California study (37). The carbon monoxide levels in the seven-county region of Texas were much lower and much less variable than those reported in southern California (California's 25th percentile for carbon monoxide exposure was approximately twice as high as the 75th percentile in Texas, and the difference between quartiles 4 and 1 was 1.3 ppm in California and 0.3 ppm in our study). We used a different critical window of exposure in our analyses. Ritz et al. (37
) examined gestational-month-specific exposures, while we calculated a mean exposure for the entire period of susceptibility, from gestational weeks 3 through 8.
Our study had several notable strengths. We conducted a large, population-based analysis using a high-quality birth defects registry with rich air pollution monitoring data. To our knowledge, our study was the first to analyze individual birth defects as well as clinical diagnostic birth defect groupings. We looked at a presumed etiologically relevant exposure window and adjusted for many potential confounders. Despite these efforts, however, residual confounding is possible, either by unmeasured or poorly characterized factors or by other environmental toxicants. A number of unmeasured factors during early pregnancy, such as use during pregnancy of a special prenatal vitamin or a general multivitamin (4350
), residential mobility (51
53
), and time activity patterns such as commuting habits and time spent outside, could all be acting as confounders, effect-measure modifiers, or sources of exposure measurement error.
There were other potential sources of measurement error in our study. Our exposure assignment strategy considered the distance between the home and the monitor and the density of data available during the critical window to develop ambient air pollution exposure estimates. We did not use dispersion models or incorporate meteorologic factors, and we used relatively simple geographic techniques. It is unlikely that the resulting measurement error differed by case-control status; in our data, geographic distributions of cases and controls were similar. National data suggest that people spend about 87 percent of their time indoors, with about 69 percent in a residence (54, 55
). Our exposure assignment method assumed that women spent the majority of their time at or near home and that exposure to ambient air pollutants would be correlated with indoor levels, an assumption supported by the research on exposures to particulate matter (56
59
). We had no data about indoor air pollution or any information on household characteristics such as air conditioner use versus open windows, heating or cooking source, housing quality, environmental tobacco smoke exposure, or external sources of air pollution such as occupational or vehicle-related exposures, any of which could potentially confound our associations.
The study sought to support the findings of a previous study of air quality and birth defects (37) that looked at clinical birth defect groupings and to extend the evaluation by looking at individual birth defects. For comparison with the previous study and to assist future efforts, we reported the results of all quartile analyses. Because of the large number of analyses, we were faced with a significant multiple comparisons problem, and it is likely that some of our statistically significant findings could have resulted from chance. Given the 255 associations presented here (17 outcomes x five pollutants x three quartiles), we had nearly a 100 percent chance of at least one association being significant (1 (1
)n) if all 255 of the individual null hypotheses were true (60
). Instead of applying an overly conservative adjustment (especially because we have correlated hypotheses in which, if quartile 2 of a pollutant exposure is associated with a birth defect, we would expect that quartiles 3 and 4 would be associated as well), we conducted a sequential "step-down" testing procedure that enabled us to control the type I error rate (40
) but present our estimates of association and calculated confidence intervals unadjusted for multiple comparisons.
Our study does not provide strong evidence that air pollution increases the risk of cardiac defects or oral clefts. However, because pathogenesis during cardiac development varies between the different types of congenital heart malformations, we did not expect to find that air pollution exposure increased the risk of all heart defects. Our results supported one of the California study's findings (37), yet we saw more suggestive associations with PM10 exposure than with any other single pollutant, which the California study did not detect. There is regional heterogeneity in the source of the particles and the potentially fetotoxic chemicals that adhere to their surfaces, which could explain why PM10 did not have similar effects in our studies. Our findings are provocative with respect to a proposed biologic mechanism that neural crest cells, essential in both embryonic heart and pharyngeal arch development (61
71
), could be selectively affected by air pollution exposures, leading to conotruncal defects such as tetralogy of Fallot.
Our study contributes to a growing body of epidemiologic literature on the adverse reproductive effects of air pollution exposure. This literature supports the notion that the developing embryo and growing fetus constitute a subpopulation susceptible to air pollution exposure, with adverse effects including reduced fetal growth and shortened gestation. The literature has also suggested that air pollution exposure during pregnancy can cause somatic and heritable gene mutations (2830
, 32
, 33
). Our results do not provide definitive evidence that air pollution exposure is a risk factor for cardiac and oral cleft birth defects, nor do our suggestive findings implicate a single pollutant or a particular developmental process. We suggest that future studies in this area direct more attention to ozone, sulfur dioxide, and particulate matter and focus on individual birth defects in addition to isolated and multiple birth defect groupings. In addition, we suggest that animal models continue to explore how air pollutant exposure may affect neural crest cells.
![]() |
ACKNOWLEDGMENTS |
---|
The authors thank Mary Ethen and Lisa Marengo of the Birth Defects Epidemiology and Surveillance Branch, and Anna Vincent and Rachelle Moore of the Bureau of Vital Statistics, of the Texas Department of State Health Services. The authors acknowledge the guidance of Catherine Harness and Sonya Krogh of Computer Sciences Corporation. They thank Dr. Sid Hunter for his helpful comments on an earlier version of the manuscript.
The information in this document has been funded wholly by the Environmental Protection Agency. It has been subjected to review by the National Health and Environmental Effects Research Laboratory and approved for publication. Approval does not signify that the contents reflect the views of the Agency, nor does mention of trade names or commercial products constitute endorsement or recommendation for use.
Conflict of interest: none declared.
![]() |
References |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
|