Children's Health Volume 122 | Issue 3 | March 2014
Traffic-Related Air Pollution and Congenital Anomalies in Barcelona
Anna Schembari,1,2,3 Mark J. Nieuwenhuijsen,1,2 Joaquin Salvador,2,4 Audrey de Nazelle,1,2,5 Marta Cirach,1,2 Payam Dadvand,1,2 Rob Beelen,6 Gerard Hoek,6 Xavier Basagaña,1,2 and Martine Vrijheid1,2
Background: A recent meta-analysis suggested evidence for an effect of exposure to ambient air pollutants on risk of certain congenital heart defects. However, few studies have investigated the effects of traffic-related air pollutants with sufficient spatial accuracy.
Objectives: We estimated associations between congenital anomalies and exposure to traffic-related air pollution in Barcelona, Spain.
Method: Cases with nonchromosomal anomalies (n = 2,247) and controls (n = 2,991) were selected from the Barcelona congenital anomaly register during 1994–2006. Land use regression models from the European Study of Cohorts for Air Pollution Effects (ESCAPE), were applied to residential addresses at birth to estimate spatial exposure to nitrogen oxides and dioxide (NOx, NO2), particulate matter with diameter ≤ 10 μm (PM10), 10–2.5 μm (PMcoarse), ≤ 2.5 μm (PM2.5), and PM2.5 absorbance. Spatial estimates were adjusted for temporal trends using data from routine monitoring stations for weeks 3–8 of each pregnancy. Logistic regression models were used to calculate odds ratios (ORs) for 18 congenital anomaly groups associated with an interquartile-range (IQR) increase in exposure estimates.
Results: In spatial and spatiotemporal exposure models, we estimated statistically significant associations between an IQR increase in NO2 (12.2 μg/m3) and coarctation of the aorta (ORspatiotemporal = 1.15; 95% CI: 1.01, 1.31) and digestive system defects (ORspatiotemporal = 1.11; 95% CI: 1.00, 1.23), and between an IQR increase in PMcoarse (3.6 μg/m3) and abdominal wall defects (ORspatiotemporal = 1.93; 95% CI: 1.37, 2.73). Other statistically significant increased and decreased ORs were estimated based on the spatial model only or the spatiotemporal model only, but not both.
Conclusions: Our results overall do not indicate an association between traffic-related air pollution and most groups of congenital anomalies. Findings for coarctation of the aorta are consistent with those of the previous meta-analysis.
Citation: Schembari A, Nieuwenhuijsen MJ, Salvador J, de Nazelle A, Cirach M, Dadvand P, Beelen R, Hoek G, Basagaña X, Vrijheid M. 2014. Traffic-related air pollution and congenital anomalies in Barcelona. Environ Health Perspect 122:317–323; http://dx.doi.org/10.1289/ehp.1306802
Address correspondence to A. Schembari, Center for Research in Environmental Epidemiology (CREAL), Parc de Recerca Biomèdica de Barcelona–PRBB, C. Doctor Aiguader, 88, 08003 Barcelona, Spain. Telephone: 34 93 214 7339. E-mail: email@example.com
This work was jointly supported by the EUROCAT (European Surveillance of Congenital Anomalies) Joint Action, funded by the Public Health Programme 2008–2013 of the European Commission (grant agreement 20102204); the ARIBA (Air Pollution and Pregnancy in Barcelona) project, funded by a Fondo de Investigación Sanitaria (FIS) (grant agreement PI081109); and by the ESCAPE (European Study of Cohorts for Air Pollution Effects) project, funded under the European Union’s Seventh Framework Programme (grant agreement 211250). A.S. is funded by an Instituto de Salud Carlos III-FEDER fellowship (PFIS grant FI 10/00476). P.D. is funded by a Ramón y Cajal fellowship (RYC-2012-10995) awarded by the Spanish Ministry of Economy and Finance.
The authors declare they have no actual or potential competing financial interests.
Received: 14 March 2013
Accepted: 23 December 2013
Advance Publication: 3 January 2014
Final Publication: 1 March 2014
- Supplemental Material (229 KB) PDF
There is a growing body of epidemiologic evidence suggesting that exposure to ambient air pollution may adversely affect the fetus and newborn. Recent studies have found associations between particulate matter (PM) and increased risk of low birth weight, preterm birth, and decrease in birth weight (Dadvand et al. 2013; Pedersen et al. 2013; Sapkota et al. 2010). Major congenital anomalies include structural defects such as heart defects and neural tube defects and chromosomal abnormalities such as Down syndrome, and are diagnosed in 2–4% of births [World Health Organization (WHO) 2012]. They are a main cause of infant mortality and important contributor to childhood and adult morbidity, but their etiology remains largely unknown (Dolk and Vrijheid 2003).
The evidence for an impact of ambient air pollution on congenital anomaly risk is still limited (Agay-Shay et al. 2013; Dadvand et al. 2011a, 2011b; Dolk et al. 2010; Gilboa et al. 2005; Hansen et al. 2009; Hwang and Jaakola 2008; Kim et al. 2007; Marshall et al. 2010; Padula et al. 2013a, 2013b; Rankin et al. 2009; Ritz et al. 2002; Strickland et al. 2009; Vrijheid et al. 2011). Previous studies have focused primarily on the routinely assessed pollutants; only four studies (Agay-Shay et al. 2013; Marshall et al. 2010; Padula et al. 2013a, 2013b) included other specific traffic-related air pollutants such as PM2.5 (PM with diameter ≤ 2.5 μm). Cardiac anomalies or oral clefts were most frequently studied, but available evidence on other anomaly groups, such as defects of the nervous, digestive, or respiratory systems, is scarce (Dolk et al. 2010; Padula et al. 2013a; Rankin et al. 2009). Summary estimates from a recent meta-analysis (Vrijheid et al. 2011) indicated that nitrogen dioxide (NO2) and sulfur dioxide (SO2) were associated with two congenital heart anomalies—coarctation of the aorta and tetralogy of Fallot—and that PM10 (PM with diameter ≤ 10 μm) was associated with atrial septal defects. Except for Dadvand et al. (2011b), exposure assessments in previous studies were based on monitoring data from a limited number of fixed-site stations, and thus did not account for the strong spatial intra-city variation that characterizes traffic-related air pollution (Cyrys et al. 2012; Eeftens et al. 2012b). Dadvand et al. (2011b) estimated spatiotemporal exposure to black smoke and SO2 with higher spatial resolution. The use of more precise spatial air pollution models in urban areas is increasingly recommended to reduce exposure misclassification (Hoek et al. 2008; Jerrett et al. 2004) in studies of traffic-related air pollution and adverse birth outcomes (Aguilera et al. 2009; Ballester et al. 2010; Brauer et al. 2008). Moreover, in studies of congenital anomalies, spatial models need to be combined with temporal adjustments to account for the very specific etiologically relevant time windows of exposure (Ritz and Wilhelm 2008).
Barcelona, Spain, is among the most polluted cities in Europe (Cyrys et al. 2012; Eeftens et al. 2012b). This is partly attributable to its geography; high traffic density (Ajuntament de Barecelona 2007), which is four times higher than London; and large proportion of diesel-powered vehicles, currently 50% (Reche et al. 2011). In Barcelona we developed a refined spatial air pollution exposure metric using a land use regression model (Beelen et al. 2013; Eeftens et al. 2012a). At the same time, Barcelona has a high-quality population-based congenital malformation register, which has collected detailed data on cases of congenital anomalies and control births for > 15 years (Greenlees et al. 2011). Using these two data sources, we investigated associations between specific groups of congenital anomalies and intra-urban traffic-related air pollution exposures during sensitive time windows in pregnancy.
Study population. This study used a population-based case–control design. The source population consisted of congenital anomaly cases and control births from the Barcelona congenital anomalies register, Registro Defectios Congenitos Barcelona (REDCB), which has accrued approximately 13,000 births/year from nearly all maternity units in the city since its start in 1992. The REDCB is part of the Service of Health Information Systems (Servei de Sistemes d’Informació Sanitaria) in the Public Health Agency of Barcelona (Agència de Salut Pública de Barcelona) and a member of the European surveillance of congenital anomalies (EUROCAT) (Greenlees et al. 2011). For the REDCB, full-time nurses actively search for new cases at delivery units, pediatric departments, cytogenetic laboratories, pathology departments, prenatal diagnosis units, and pediatric cardiology services; cases of congenital anomalies suspected but not detected at birth are further followed up until age 1 year. Controls are continuously selected as a 2% representative sample of all live births in Barcelona hospitals, by random date of birth sampling, independent of the cases (nonmatched). Clinical information, including maternal age, is collected from hospital records. As part of the REDCB, information on smoking and education is collected by interview with the mother, using a standardized questionnaire (Salvador et al. 2005). Parental informed consent is completed before the interview, and all consenting cases detected during the maternal period at the hospital and the control mothers are interviewed (35% and 95%, respectively); interviews cannot be held with mothers of cases diagnosed or communicated to the registry authority after the maternity period. In these cases, smoking and maternal education information is abstracted from clinical history records. Further, for the purposes of the current study, the residential addresses of the cases and controls were linked to the MEDEA (Mortalidad en áreas pequeñas Españolas y Desigualdades Socioeconómicas y Ambientales) socioeconomic deprivation index (Dominguez-Berjon et al. 2008) at census-tract level.
The present study was approved by the Clinical Research Ethics Committee of PS–Mar (CEIC project 2008/3115/I and 2006/3394/I) and determined to be exempt from separate informed consent requirements.
Definition of cases and controls. Controls were all live births enrolled in the REDCB during 1994–2006 with a residential address at birth in the city of Barcelona. Of the 3,149 controls included in the REDCB during that time period, 158 were excluded because of a residential address that could not be geocoded or that was outside Barcelona, giving a total of 2,991 (95.0%) controls for analysis.
Cases were live births, stillbirths after 20 weeks of gestation, or terminations of pregnancy after prenatal diagnosis of congenital anomalies, enrolled in the REDCB and born or terminated between 1994 and 2006, with at least one major anomaly in the EUROCAT list of subgroups of congenital anomalies (EUROCAT 2009). Cases with only minor congenital anomalies, such as facial asymmetry (code Q67.0) or talipes or pes calcaneovalgus (Q66.4) were excluded from the study, following the EUROCAT exclusion list (EUROCAT 2009). The range of codes for inclusion was, according to the International Classification of Diseases, 10th Revision (ICD-10) (WHO 1993), the Q chapter. Cases with chromosomal anomalies (Q90–Q92, Q93, Q96–Q99; n = 719), and teratogenic and genetic syndromes (Q86, Q87; n = 72) were excluded. Twins were treated as one outcome (12 pairs) and classified as a case if one or both had a congenital anomaly; siblings were classified as separate outcomes. Cases with residential addresses that could not be geocoded or that were outside Barcelona were excluded (n = 166), as were spontaneous abortion before 20 weeks of gestation (n = 5). In addition, we excluded cases that could not be classified according to the timing of the vulnerable time window of exposure (i.e., 3–8 weeks gestation, when organogenesis occurs) because any information on birth date, gestational age (GA), and date of the last menstrual period (LMP) (n = 74) was available. Of 3,295 cases, 2,247 (68%) were included in analysis. Cases were classified into anomaly subgroups following EUROCAT guidelines (EUROCAT 2009). An a priori decision was made to analyze the following common subgroups of anomalies: neural tube defects (Q00, Q01, Q05), congenital heart disease (Q20–Q26), respiratory system defects (Q30–Q34), orofacial clefts (Q35–Q37), digestive system defects (Q38–Q38, Q402–Q409, Q41–Q45, Q790), abdominal wall defects (Q792, Q793, Q795), urinary system defects (Q60–Q64, Q794), hypospadias (Q54), and limb reduction defects (Q650–652, Q658–Q659, Q660, Q681–Q682, Q688, Q69–Q74). To replicate results from previous studies (Vrijheid et al. 2011), we also considered eight subgroups of cardiac anomalies: transposition of great vessels (Q203), ventricular septal defect (Q210), atrial septal defect (Q211), atrioventricular defect (Q212), tetralogy of Fallot (Q213), tricuspid atresia and stenosis (Q224), pulmonary valve stenosis (Q221), and coarctation of aorta (Q251).
Exposure assessment. This study combined the spatial land use regression (LUR) modeling developed in the European Study of Cohorts for Air Pollution Effects–ESCAPE (ESCAPE 2009) with a temporal adjustment, providing both spatial and spatiotemporal exposure estimates. LUR models are commonly used to explain the small-scale within-city variation of traffic-related air pollutants; temporal adjustment is needed to account for the timing of exposure in relation to organogenesis—during weeks 3 to 8 of pregnancy (Moore et al.1998).
Spatial estimates. LUR models were constructed in Barcelona for the year 2009 as part of the ESCAPE project. Pollutants included NO2, nitrogen oxides (NOx), PM10, PMcoarse (particulate matter with aerodynamic diameter 10–2.5 μm), PM2.5, and PM2.5 absorbance (as a marker of black carbon) (Cyrys et al. 2003). The methodology has been described in detail elsewhere (Beelen et al. 2013; Cyrys et al. 2012; Eeftens et al. 2012a). All addresses were geocoded to assign the spatial exposure estimates. We assumed that the residential address at delivery was constant during the whole pregnancy period and that the city spatial distribution of pollutants and their determinants remained constant over the study period.
Spatiotemporal estimates. The start and end dates of weeks 3–8 of pregnancy were determined from the date of birth and the GA recorded on the register (n = 4,866). When GA was not available, the declared date of LMP was used (n = 261) to determine it. When only the birth date was available, we imputed the GA with the mean GA for the specific type of delivery (live birth, stillbirth, or termination of pregnancy; n = 111).
To estimate exposure during the “critical window” for organogenesis (weeks 3–8), background routine daily ambient air pollution monitoring data were used. NO2 and NOx series were available for the complete study period, from 1994 through 2006 (29% of days with no measurement), and PM10 data from 2000 through 2006 (68% of days with no measurement). On days with no measurement, NO2 and NOx data were estimated using the multiple imputation method developed in the ICE function of the statistical package STATA version 8 (StataCorp, College Station, TX, USA). Because of the large number of days for which PM10 data were missing, we did not impute missing PM10 data. Instead, we restricted our population to those pregnancies with data available on ≥ 2 days per week for at least 5 of the 6 weeks of interest (weeks 3–8 of gestation). This resulted in 63% of the original population. Then, according to the ratio method procedure in the ESCAPE manual for extrapolation back in time (ESCAPE 2009), we calculated the spatiotemporal exposure estimates as follow:
Procedure for extrapolation back in time. The daily spatiotemporal exposure estimates (Cextrapolated) for each pregnancy were the spatial exposure estimation (CESCAPE) multiplied by a daily temporal adjustment (Ratioroutine) that was calculated as the ratio between the daily NO2, NOx, or PM10 measurements at the background routine monitoring station (Cdaily) and its annual 2009 average (Croutine_yearESCAPE):
Cextrapolated = CESCAPE × Ratioroutine, 
Ratioroutine = Cdaily/Croutine_yearESCAPE. 
Daily spatiotemporal estimates were averaged over weeks 3–8 of each pregnancy to obtain the final exposure. We further applied the NOx daily temporal adjustment to PM2.5 absorbance spatial estimates because of their high correlation (Cyrys et al. 2003), and the PM10 daily temporal adjustment to PMcoarse and PM2.5 (Pedersen et al. 2013), to obtain the respective spatiotemporal estimates of exposure to air pollution, following the methodology in ESCAPE (2009).
Statistical analysis. We examined the association between the spatial and spatiotemporal exposure estimates and 18 groups of congenital anomalies using unconditional logistic regression models. Air pollution exposure estimates were entered into the model as continuous covariates. The odds ratios (ORs) were calculated for an interquartile range (IQR) increase in air pollutant concentration measured at the spatial level to ensure consistent exposure contrasts between spatial and spatiotemporal estimates, and comparable contrasts among different pollutants. Models were adjusted for year of birth/termination, season of conception, and maternal age, all as categorical variables, and for the MEDEA socioeconomic deprivation index as a continuous variable (Dadvand et al. 2011a; Ritz et al. 2002).
We tested potential confounding by maternal smoking and maternal education in sensitivity analyses restricted to the populations with data available on these variables (70% and 80% of the full study population, respectively). Information on maternal smoking during each month of pregnancy was used to calculate the proportion of smokers during the first trimester as a two-category variable (yes or quit/no). Maternal education was coded as a four-category variable (primary/secondary/college/university). The assumption of a stable spatial exposure surface over the entire study period was investigated by stratifying the population according to year of birth/termination (before and after 2000). All analyses were performed using the statistical package STATA versions 8 and 10.1 (StataCorp).
The study included 2,247 cases and 2,991 controls. There were no statistically significant differences between cases and controls for season of conception, sex, maternal age, and socioeconomic index (Table 1), but differences were seen for maternal smoking during first trimester and for maternal education: Mothers of cases had less education and smoked more frequently than did mothers of controls. Both these characteristics were derived from the interview with the mother completed by the 35% of cases and 95% of controls; in that subpopulation, the exposure distribution was not differential among cases and controls (see Supplemental Material, Table S1).
Exposure distributions of the air pollutants are presented in Table 2. The IQRs of the spatiotemporal exposure estimates, which were averaged over the 6-week period of organogenesis for each observation, were wider than the IQRs of the spatial estimates, which were averaged over an entire year, because of the short averaging time and seasonal variation. Exposure estimates were somewhat higher among controls than among cases for all exposures. The range of correlations between spatial exposures (r = 0.16–0.93) was much wider than the corresponding range of correlations among the spatiotemporal exposures (r = 0.85–0.91) (Table 3) because of the high seasonal correlation between the pollutants (data not shown) and the use of NOx and PM10 temporal trends to adjust the other air pollutants. Correlations between spatial and spatiotemporal estimates were moderate to low (r = 0.55 and 0.26 for NOx and PM10 respectively).
We observed several statistically significant associations between IQR increases in spatial air pollutant exposure estimates and congenital anomaly groups (see Supplemental Material, Table S2): coarctation of the aorta and NO2 [adjusted (adj)ORspatial = 1.08; 95% CI: 1.00, 1.16], NOx (adjORspatial = 1.06; 95%CI: 1.00, 1.13), and PM2.5 (adjORspatial = 1.25; 95% CI: 1.00, 1.57); transposition of great vessels and PM10 (adjORspatial = 1.27; 95% CI: 1.03, 1.57); digestive system defects and NO2 (adjORspatial = 1.07; 95% CI: 1.02, 1.13) and NOx (adjORspatial = 1.06; 95% CI: 1.01, 1.10); and abdominal wall defects and PMcoarse (adjORspatial = 1.60; 95% CI: 1.04, 2.48). Furthermore, statistically significant negative associations with IQR increases in exposure were observed for all nonchromosomal anomalies and PM2.5 absorbance (adjORspatial = 0.92; 95% CI: 0.86, 0.98), all congenital heart defects and PM2.5 absorbance (adjORspatial = 0.91; 95% CI: 0.83, 0.99), and atrial septal defects and PMcoarse (adjORspatial = 0.73; 95% CI: 0.60, 0.89).
The association between an IQR increase in NO2 (12.2 μg/m3) and coarctation of the aorta was increased in the spatiotemporal exposure model (adjORspatiotemporal = 1.15; 95% CI: 1.01, 1.31) compared with the spatial model, though it was somewhat reduced with an IQR increase in NOx (27.1 μg/m3; adjORspatiotemporal = 1.02; 95% CI: 1.00, 1.04) (Table 4). Similarly, digestive system defects showed a somewhat increased association with NO2 (adjORspatiotemporal = 1.11; 95% CI: 1.00, 1.23) and decreased association with NOx (adjORspatiotemporal = 1.02; 95% CI: 1.00, 1.04). The association between abdominal wall defects and an IQR increase in PMcoarse (3.6 μg/m3) was also higher based on the spatiotemporal model (adjORspatiotemporal = 1.93; 95% CI: 1.37, 2.73) (Table 5) compared with the spatial model. The spatiotemporal model also generated several statistically significant associations that were not observed in the spatial model, including significant positive associations between abdominal wall defects and PM10 and PM2.5, and significant negative associations for ventricular septal defects and PM10, PMcoarse, and PM2.5.
Associations based on spatial models adjusted for maternal smoking (n = 3,752) or maternal education (n = 4,211) were consistent with estimates that were not adjusted for these variables when restricted to the population with available smoking or education data (see Supplemental Material, Tables S3 and S4). However, some associations differed when based on the restricted populations compared with the full study population. In particular, there were more positive associations (ORs > 1) in the restricted population. The direction and magnitude of associations were comparable for anomalies identified from 1994 through 1999 (n = 2,324) and 2000 through 2006 (n = 2,913), consistent with our assumption of a stable exposure surface over the entire study period, though estimates based on the earlier subset were less precise (see Supplemental Material, Table S5).
In this study we used both spatial and spatiotemporal exposure assessment frameworks to evaluate the association between traffic-related air pollution and groups of congenital anomalies. We included pollutants not previously studied in this field: PMcoarse, PM2.5, and PM2.5 absorbance, all of which are characterized by within-city variability (Eeftens et al. 2012b; Jerrett et al. 2005). We found little evidence for an association between most of the traffic-related air pollutants and most of the congenital anomaly groups studied. However, we observed statistically significant positive associations with both spatial and spatiotemporal exposure estimates for NO2 and NOx with coarctation of the aorta and digestive system anomalies, and for PMcoarse and abdominal wall defects.
Given the multiple comparisons involved in testing a large range of congenital anomaly subtypes against six pollution measures, our analyses are likely to have given rise to some chance associations, and results should be interpreted with caution. To avoid errors due to false-positive statistically significant associations (type 1 errors), we emphasized statistically significant associations from the spatial model that persisted after temporal adjustment, in contrast with statistically significant positive or negative associations based on only one of the two exposure models, which we considered more likely to be attributable to chance. Evaluation of consistency with the previous meta-analysis has further been used to interpret the role of chance as an explanation of our main findings.
The positive association between coarctation of the aorta and NO2 is consistent with a recent meta-analysis (Vrijheid et al. 2011) that combined results of four published studies and reported a summary OR for coarctation of the aorta of 1.20 (95% CI: 1.00, 1.44) per 10-ppb increase in NO2 exposure. For comparison, in our study the OR from the spatiotemporal model for a 10-ppb increment in NO2 was 1.23 (95% CI: 1.02, 1.48). Although in our study the number of cases with coarctation of the aorta is small (n = 69), the consistency of findings between spatial and spatiotemporal models, and with the published meta-analysis, appears to strengthen evidence for an association. NO2 and NOx were highly correlated, making it hard to separate their effect, but in the spatiotemporal exposure models NO2 appeared somewhat more strongly associated with coarctation of the aorta than NOx. We also found consistent results between the spatial and the spatiotemporal models when analyzing the association between NO2, NOx, and digestive system anomalies. Only two previous studies evaluated these anomalies: Dolk et al. (2010) reported no association with NO2, but the exposure assessment was based on annual mean levels only; Rankin et al. (2009) evaluated digestive system anomalies in relation to black smoke and also reported no association. In the present study, an association between abdominal wall defects and PMcoarse was observed in both spatial and spatiotemporal models. We are aware of only two previous studies that evaluated associations between air pollutants and abdominal wall defects (Dolk et al. 2010; Padula et al. 2013a), and both reported nonsignificant associations except for the association between PM10 and omphalocele reported by Dolk et al. (2010). Therefore, further research on this outcome is needed.
The exposure assessment method used in this study was a temporally adjusted LUR built to capture and replicate the spatial variability within cities of traffic-related air pollutants (Cyrys et al. 2012; Eeftens et al. 2012b). Additionally, for congenital anomalies, it is of great importance to assess exposure during the critical pregnancy weeks 3–8 (Moore et al. 1998), because any exposure after this period would not contribute to the etiology of a major congenital anomaly. Only recent studies on congenital anomalies and air pollution used a spatiotemporal exposure assessment, but they were not covering the same range of pollutants (Dadvand et al. 2011b) or assessing the exposure at less refined spatial resolution (Agay-Shay et al. 2013; Padula et al. 2013a, 2013b). The spatial LUR model we used for Barcelona resulted in a spatial exposure assessment with similar variability (i.e., similar width of the IQR) to previous studies in this field (Vrijheid et al. 2011). After the temporal adjustment, the exposure variation increased, widening the IQR by > 50% for all pollutants (Table 2), thus giving greater statistical power to detect associations in the spatiotemporal models compared with the spatial models. In particular, NO2 and NOx IQRs (reflecting within-city spatial variation) were larger than for particulate matter, partially explaining the null results for particulate matter in the spatial exposure models. One important assumption in using the LUR is that the city spatial distribution of pollutants and their determinants remained constant over the study period; if this assumption is violated, the ORs may be attenuated (Wu et al. 2011). To test this assumption we split the population by earlier and later years of termination/birth (before 2000 vs. during or after 2000); results were consistent with the overall analyses, although earlier years had wider confidence intervals than the later years. Recently Cesaroni et al. (2012) showed in Rome, a city with characteristics (e.g., climate, traffic) similar to those of Barcelona, a correlation of 96% for the exposure estimates from two LUR models produced 11 years apart. Furthermore, we did not observe long-term trends in NO2 and PM10 (data not shown). A limitation of our exposure assessment method is that we used NOx and PM10 temporal trends to adjust other pollutant exposure estimates. This led to high correlations between the spatiotemporal estimates of these pollutants and thus difficulty in the interpretation of results from these pollutants independently. However, we followed the approach developed in the ESCAPE protocol, which has already been used in the recent publication on air pollution effect on pregnancy outcomes by Pedersen et al. (2013). Another exposure limitation was related to the missing data on the PM10 temporal series, which led to a loss of statistical power in the spatiotemporal models of pollutants temporally adjusted with that series (PM10, PMcoarse, and PM2.5).
Exposure misclassification could be attributable to the time spent in different environments, although we think that it is unlikely to be differential among cases and controls. Exposure misclassification could further have arisen because we estimated outdoor exposure at the residential address at termination/birth as a surrogate of the personal exposure during the first months of pregnancy. In our recent study in Barcelona (Schembari et al. 2013) among 54 pregnant women who carried a personal PM2.5 sampler for 2 days and NOx/NO2 passive badges for 1 week, pregnant women reported time spent at home to be around 60–70% per day, and the correlation between personal exposure and outdoors levels of air pollution ranged from 0.39 for PM2.5 and 0.78 for NOx. This suggests that, particularly for the latter, ambient outdoor levels could act as a good surrogate for personal exposure levels. Residual misclassification may have occurred if women changed residences during pregnancy; but this is unlikely to occur differentially among cases and controls, accordingly with results published by Miller et al. (2010). Residential mobility was estimated to be between 1% and 6% in a recent Spanish study based on four birth cohorts (Estarlich et al. 2011).
We examined a wide range of congenital anomalies that were subtyped following the classification proposed by EUROCAT (2009). The REDCB includes cases identified among live births, stillbirths, and termination of pregnancy following prenatal diagnosis, which is particularly important (Ritz 2010) for some severe congenital anomalies such as neural tube defects or congenital heart disease (27% and 22.6% of anomalies, respectively, in our study). The register follows the EUROCAT guidelines, but ascertainment may be incomplete if the congenital anomaly is not diagnosed prenatally or at birth. Therefore, the prevalence of all congenital anomalies (including live births, fetal deaths, stillbirths, and termination of pregnancy) for Barcelona was lower, during 1994–2006, than the European register—184/10,000 and 233/10,000 respectively (EUROCAT 2012). The register uses active case ascertainment and simultaneous ascertainment of population-based controls; it recodes the address at delivery allowing the geocoding used for the exposure assessment in this study. We adjusted for maternal age, conception season, year of birth, and socioeconomic status, although the latter was classified at census tract level.
Only a subset of participants had data on smoking (70%) and education (80%). The sensitivity analyses on these subpopulations showed little impact of maternal smoking or education on the findings, thus indicating that their confounding effect is likely to be small even in the full study population. Nevertheless, relying on these subpopulation could be misleading: Missing information probably did not occur at random among cases and controls mainly because mothers of cases detected after birth were not interviewed. This may constitute a selective population due to the fact that only specific congenital anomalies are detected after birth; for example, congenital heart defects are frequently detected in the first year of life (Todros et al. 2001). The results of the sensitivity analyses differed for some congenital anomalies/exposure from those on the full population; in particular, PM10 was statistically significantly associated to congenital heart defects in the smoking subpopulation (both adjusted and unadjusted), but not in the main analyses nor in the one adjusted for maternal education. Residual confounding by unmeasured risk factors such as alcohol consumption and dietary factors (e.g., folic acid) could have influenced our findings. However, it is unlikely that such factors would be strongly related to air pollution exposure, and we can partially acount for them by adjusting for social class on area level in the main analyses, and on individual level as sensitivity analysis.
The possible mechanisms and causes underlying the development of congenital anomalies are still unclear, because of their probable multifactorial aetiology (Ritz 2010) and because different types of anomalies are likely to have very distinct etiology. Air pollution–induced oxidative stress during pregnancy has been suggested as one possible mechanism behind pregnancy outcomes (Kannan et al. 2006; Slama et al. 2008), because it regulates the pulmonary and placental inflammation, the hemodynamic responses, and thus the transplacental oxygenation and transportation of nutrients. Oxidative stress may also affect organogenesis and neural crest cell migration and differentiation (Hassler and Moran 1986; Ritz et al. 2002), which play an important role in heart development (Keyte and Redmond Hutson 2012). More recently, van Beynum et al. (2008) suggested a possible gene–environment interaction effect leading to an increased risk of congenital heart defects in mothers exposed to nitric oxide and smoking.
Our results overall do not indicate an association between exposure to traffic-related air pollution and many groups of congenital anomalies in Barcelona, even though the air pollution levels are some of the highest in Europe. The positive association of NO2 and NOx with coarctation of the aorta is consistent with a findings of a meta-analysis of previous studies, and requires further study. Associations of digestive system anomalies with NO2 and NOx, and of abdominal wall defects with PMcoarse, also call for confirmation.
Agay-Shay K, Friger M, Linn S, Peled A, Amitai Y, Peretz C. 2013. Air pollution and congenital heart defects. Environ Res 124:28–34; doi: 10.1016/j.envres.2013.03.005.
Aguilera I, Guxens M, Garcia-Esteban R, Corbella T, Nieuwenhuijsen MJ, Foradada CM, et al. 2009. Association between GIS-based exposure to urban air pollution during pregnancy and birth weight in the INMA Sabadell Cohort. Environ Health Perspect 117:1322–1327; doi: 10.1289/ehp.0800256.
Ajuntament de Barecelona. 2007. 2007 Dades Bàsiques de Mobilitat [in Catalan]. Available: http://w110.bcn.cat/Mobilitat/Continguts/Documents/Fitxers/dadesbasiques2007web.pdf [accessed 13 July 2013].
Ballester F, Estarlich M, Iñiguez C, Llop S, Ramón R, Esplugues A, et al. 2010. Air pollution exposure during pregnancy and reduced birth size: a prospective birth cohort study in Valencia, Spain. Environ Health 9:6; doi: 10.1186/1476-069X-9-6.
Beelen R, Hoek G, Vienneau D, Eeftens M, Dimakopoulou K, Pedeli X, et al. 2013. Development of NO2 and NOx land use regression models for estimating air pollution exposure in 36 study areas in Europe—the ESCAPE project. Atmos Environ 72:10–23.
Brauer M, Lencar C, Tamburic L, Koehoorn M, Demers P, Karr C. 2008. A cohort study of traffic-related air pollution impacts on birth outcomes. Environ Health Perspect 116:680–686; doi: 10.1289/ehp.10952.
Cesaroni G, Porta D, Badaloni C, Stafoggia M, Eeftens M, Meliefste K, et al. 2012. Nitrogen dioxide levels estimated from land use regression models several years apart and association with mortality in a large cohort study. Environ Health 11:48; doi: 10.1186/1476-069X-11-48.
Cyrys J, Eeftens M, Heinrich J, Ampe C, Aarmengaud A, Beelen R, et al. 2012. Variation of NO2 and NOx concentrations between and within 36 European study areas: results from the ESCAPE study. Atmos Environ 62:374–390.
Cyrys J, Heinrich J, Hoek G, Meliefste K, Lewné M, Gehring U, et al. 2003. Comparison between different traffic-related particle indicators: elemental carbon (EC), PM2.5 mass, and absorbance. J Expo Anal Environ Epidemiol 13:134–143.
Dadvand P, Parker J, Bell ML, Bonzini M, Brauer M, Darrow LA, et al. 2013. Maternal exposure to particulate air pollution and term birth weight: a multi-country evaluation of effect and heterogeneity. Environ Health Perspect 121:367–373; doi: 10.1289/ehp.1205575.
Dadvand P, Rankin J, Rushton S, Pless-Mulloli T. 2011b. Association between maternal exposure to ambient air pollution and congenital heart disease: a register-based spatio-temporal analysis. Am J Epidemiol 173:171–182.
Dominguez-Berjon MF, Borrell C, Cano-Serral G, Esnaola S, Nolasco A, Pasarin MI, et al. 2008. Constructing a deprivation index based on census data in large Spanish cities (the MEDEA project). Gaceta Sanitaria 22:179–187.
Eeftens M, Beelen R, Bellander T, Cesaroni G, Cirach M, Declerq C, et al. 2012a. Development of land use regression models for PM2.5, PM2.5 absorbance, PM10 and PMcoarse in 20 European study areas; results of the ESCAPE project. Environ Sci Technol 46:11195–11205.
Eeftens M, Tsai MY, Ampe C, Anwander B, Beelen R, Bellander T, et al. 2012b. Spatial variation of PM2.5, PM10, PM2.5 absorbance and PMcoarse concentrations between and within 20 European study areas and the relationship with NO2—results of the ESCAPE project. Atmos Environ 62:303–317.
ESCAPE. 2009. ESCAPE—European Study of Cohorts for Air Pollution Effects. Available: http://www.escapeproject.eu/index.php [accessed 13 July 2013].
Estarlich M, Ballester F, Aguilera I, Fernandez-Somoano A, Lertxundi A, Llop S, et al. 2011. Residential exposure to outdoor air pollution during pregnancy and anthropometric measures at birth in a multicenter cohort in Spain. Environ Health Perspect 119:1333–1338; doi: 10.1289/ehp.1002918.
EUROCAT. 2009. EUROCAT Guide 1.3. Instructions for the Registration of Congenital Anomalies. Available: http://www.eurocat-network.eu/aboutus/datacollection/guidelinesforregistration/guide1_3instructionmanual [accessed 13 July 2013].
EUROCAT. 2012. Prevalence Tables. Available: http://www.eurocat-network.eu/ACCESSPREVALENCEDATA/PrevalenceTables [accessed 3 April 2012].
Gilboa SM, Mendola P, Olshan AF, Langlois PH, Savitz DA, Loomis D, et al. 2005. Relation between ambient air quality and selected birth defects, seven county study, Texas, 1997–2000. Am J Epidemiol 162:238–252.
Hansen CA, Barnett AG, Jalaludin BB, Morgan GG. 2009. Ambient air pollution and birth defects in Brisbane, Australia. PLoS One 4:e5408; doi: 10.1371/journal.pone.0005408.
Hwang B-F, Jaakkola JJK. 2008. Ozone and other air pollutants and risk of oral clefts. Environ Health Perspect 116:1411–1415; doi: 10.1289/ehp.11311.
Kannan S, Misra DP, Dvonch JT, Krishnakumar A. 2006. Exposures to airborne particulate matter and adverse perinatal outcomes: a biologically plausible mechanistic framework for exploring potential effect modification by nutrition. Environ Health Perspect 114:1636–1642; doi: 10.1289/ehp.9081.
Padula AM, Tager IB, Carmichael SL, Hammond SK, Lurmann F, Shaw GM. 2013a. The association of ambient air pollution and traffic exposures with selected congenital anomalies in the San Joaquin Valley of California. Am J Epidemiol 177(10):1074–1085. doi:10.1093/aje/kws367.
Padula AM, Tager IB, Carmichael SL, Hammond SK, Yang W, Lurmann F, et al. 2013b. Ambient air pollution and traffic exposures and congenital heart defects in the San Joaquin Valley of California. Paediatr Perinat Epidemiol 27(4):329–339. doi:10.1111/ppe.12055.
Reche C, Querol X, Alastuey A, Viana M, Pey J, Moreno T, et al. 2011. New considerations for PM, black carbon and particle number concentration for air quality monitoring across different European cities. Atmos Chem Phys 11:6207–6227.
Salvador J, Borrell A, Lladonosa A. 2005. Increasing detection rates of birth defects by prenatal ultrasound leading to apparent increasing prevalence. Lessons learned from the population-based registry of birth defects of Barcelona. Prenat Diagn 25:991–996.
Slama R, Darrow L, Parker J, Woodruff TJ, Strickland M, Nieuwenhuijsen M, et al. 2008. Meeting report: atmospheric pollution and human reproduction. Environ Health Perspect 116:791–798; doi: 10.1289/ehp.11074.
Strickland MJ, Klein M, Correa A, Reller MD, Mahle WT, Riehle-Colarusso TJ, et al. 2009. Ambient air pollution and cardiovascular malformations in Atlanta, Georgia, 1986–2003. Am J Epidemiol 169:1004–1014.
van Beynum IM, Mooij C, Kapusta L, Heil S, den Heijer M, Blom HJ. 2008. Common 894G>T single nucleotide polymorphism in the gene coding for endothelial nitric oxide synthase (eNOS) and risk of congenital heart defects. Clin Chem Lab Med 46:1369–1375.
Vrijheid M, Martinez D, Manzanares S, Dadvand P, Schembari A, Rankin J, et al. 2011. Ambient air pollution and risk of congenital anomalies: a systematic review and meta-analysis. Environ Health Perspect 119:598–606; doi: 10.1289/ehp.1002946.
WHO (World Health Organization). 2012. Congenital Anomalies. Fact sheet nº370. Available: http://www.who.int/mediacentre/factsheets/fs370/en/ [accessed 13 July 2013].
For a number of years, EHP has posted accepted manuscripts online as what we termed Advance Publications. The concept here was simple—to get the content online as quickly as possible in an effort both to get authors the attention they deserve and to get readers the information they need. This product has served our community well, but it is not without its flaws. … MORE
EHP has launched an updated homepage! This bold, new look will give readers new ways to find our content, and will allow our Editors to highlight important articles that our audience needs to be aware of. … MORE
Featured Children’s Health
Birgit Claus Henn, Adrienne S. Ettinger, Marianne R. Hopkins, Rebecca Jim, Chitra Amarasiriwardena, David C. Christiani, Brent A. Coull, David C. Bellinger, and Robert O. Wright
Diane Gilbert-Diamond, Jennifer A. Emond, Emily R. Baker, Susan A. Korrick, and Margaret R. Karagas
Benjamin B. Green, Margaret R. Karagas, Tracy Punshon, Brian P. Jackson, David J. Robbins, E. Andres Houseman, and Carmen J. Marsit