|Home | About | Journals | Submit | Contact Us | Français|
Hazardous air pollutants are plausible candidate exposures for autism spectrum disorders. They have been explored in recent studies for their role in the development of these disorders.
We used a prevalent case-control design to screen perinatal exposure to 35 hazardous air pollutants for further investigation in autism etiology. We included 383 children with autism spectrum disorders and, as controls, 2829 children with speech and language impairment. All participants were identified from the records-based surveillance of 8-year-old children conducted by the Autism and Developmental Disabilities Monitoring Network in North Carolina (for children born in 1994 and 1996) and West Virginia (born in 1992 and 1994). Exposures to ambient concentrations of metal, particulate, and volatile organic air pollutants in the census tract of the child’s birth residence were assigned from the 1996 National Air Toxics Assessment annual-average model. We estimated odds ratios (ORs) for autism spectrum disorders and corresponding 95% confidence intervals (CIs), comparing across the 20th and 80th percentiles of log-transformed hazardous air pollutant concentration among the selected controls, using semi-Bayes logistic models and adjusting for sampling variables (surveillance year and state), a priori demographic confounders from the birth certificate and census, and covarying air pollutants.
We estimated many near-null ORs, including those for metals, established human neurodevelopmental toxicants, and several pollutants that were elevated in a similar study in California. Hazardous air pollutants with more precise and elevated OR estimates included methylene chloride, 1.4 (95% CI = 0.7–2.5), quinoline, 1.4 (1.0–2.2), and styrene, 1.8 (1.0–3.1).
Our screening design was limited by exposure misclassification of air pollutants and the use of an alternate developmental disorder as the control group, both of which may have biased results toward the null. Despite these limitations, methylene chloride, quinoline, and styrene emerged (based on this analysis and prior epidemiologic evidence) as candidates that warrant further investigation for a possible role in autism etiology.
Hazardous air pollutants include hundreds of metal, particulate, and volatile organic compounds known to harm human health. Many are plausible candidate exposures in autism etiology because they have neurotoxic or immunotoxic properties,1 and theories of autism etiology include damage to the developing nervous system and immune perturbation.2–4 Airborne chemical exposures may be especially important in contributing to neurodevelopmental disorders because inhaled particles or metals have been found to be delivered directly to the brain through the olfactory bulb.5,6
Although hazardous air pollutants are not routinely monitored in the United States, emissions data have been used to model annual-average census-tract ambient concentrations in the National Air Toxics Assessment (NATA) program. A study in California studied hazardous air pollutants and autism,7 using NATA output to assign exposures. Children with autism spectrum disorders from a state developmental service agency and a large health maintenance organization were compared with children from birth certificate rosters. Hazardous air pollutants were examined individually and combined into composite scores based on chemical structure or mechanism (eg, endocrine disruptors). From a selected list of 25 air pollutants, positive associations were found for the chlorinated solvent group, the metal group, and several individual pollutants.
Despite calls for data to fill gaps in our understanding of associations between environmental agents and autism,8,9 this California study remains one of the few studies contributing actual data. Intense investigation of many environmental agents is time- and cost-prohibitive. The NATA model presents an opportunity to screen a large number of pollutants with biologic plausibility for a role in autism. Results of such screenings can be considered in combination with toxicologic literature to identify the candidates that deserve more indepth research. Semi-Bayes hierarchical methods are useful in multiple-comparisons situations in which investigators must choose which comparisons to investigate further. This is true especially when there is substantial cost associated with future investigations,10,11 as there is with autism and environmental exposures.
We conducted a screening analysis of hazardous air pollutants and autism spectrum disorders using a prevalent case-control design in North Carolina and West Virginia. As in the California study,7 exposure was assigned using output from the NATA model, but we additionally included some design improvements. We included a more complete group of autism spectrum disorders cases by reviewing records from schools in addition to other sources12 and verified residency of our control group at age 8. To account for the fact that hazardous air pollutants are highly correlated and that more than one of them might have effects on autism, we simultaneously adjusted for measured air pollutants and stabilized the estimates using semi-Bayes models.
Cases with autism spectrum disorders and controls with speech and language impairment were identified from the Autism and Developmental Disabilities Monitoring Network (ADDM) in North Carolina and West Virginia.12 The monitoring network screens developmental records of children in health and educational settings through the 8th year of life, an age at which most affected children have been identified.13 Autism spectrum disorders diagnoses are fairly stable.14 The North Carolina site evaluated children living in the counties of Alamance, Chatham, Davidson, Durham, Forsyth, Guilford, Orange, and Randolph in 2002 and 2004 (born in 1994 and 1996). The West Virginia site evaluated children living in the entire state in 2000 and 2002 (born in 1992 and 1994). We limited participants in this study to children who had resided in the surveillance region at the time of birth, to obtain in-state birth certificate data and to preclude a potential selection bias induced by differential case versus control in-migration from regions with different concentrations of air pollutants.
Cases were all children with developmental records documenting characteristics and behaviors that met a standardized definition for autism spectrum disorders based on the Diagnostic and Statistical Manual of the American Psychiatric Association (DSM-IV-TR),15 regardless of whether a previous diagnosis of autism spectrum disorders had been documented.
Controls were defined as all children in the surveillance system in North Carolina and West Virginia with a school designation of speech and language impairment without documentation of other serious developmental problems (eg, autism spectrum disorders, intellectual disabilities etc). We assumed that speech and language impairments would not be appreciably affected by air pollution. We selected this control group because (1) the children had equivalent access to developmental evaluations, and (2) we could confirm their residency at age 8. Family residential mobility is common, and has been associated with lifestyle factors such as maternal age and income.16 Children who out-migrate may have different air pollution exposures, necessitating the selection of a control group meeting the same residency requirements as cases.
To obtain the residential address in early life and demographic characteristics, we linked each child to corresponding birth certificate data. In West Virginia, we randomly sampled one-third of all children identified with speech and language impairment within each surveillance year to reduce the burden on the vital statistics department performing birth certificate matches.
We assigned exposures during the perinatal period because it is a period of susceptibility to exogenous agents in autism etiology.17 Ambient hazardous air pollutant concentrations were not directly measured in North Carolina or West Virginia in the early to mid-1990s, so we assigned exposure using modeled concentrations from the NATA program. NATA was developed as a tool to prioritize geographic areas and individual pollutants requiring study or intervention.18 We selected NATA 1996 estimates over those from 1990 or 1999 because the 1996 model had improved data input compared with earlier models and was closer in time to the birth-years of our cohort.19
NATA-1996 included diesel particulate matter, components of vehicle exhaust of concern to the Office of Transportation and Air Quality (ethylbenzene, hexane, methyl-tert-butyl ether, propionaldehyde, styrene, toluene, and xylenes), and 32 priority pollutants of concern in urban areas (selected based on toxicity-weighted emission levels). Our analytic method required independent exposure groups. Therefore, we selected the combined group of 7 polycyclic aromatic hydrocarbons (PAH) from NATA instead of the encompassing group of 16 compounds termed polycyclic organic matter because previous studies suggesting human neurodevelopmental toxicity have focused on PAH.20
The NATA model assumes Gaussian air dispersion and uses overlapping spatial grids to predict annual-average ambient concentrations of 41 hazardous air pollutants for each census tract.21 The primary source of inputs into the model is the National Emissions Inventory, which includes the location and rate of pollutant release from point sources such as manufacturing, power-generating, and waste incineration facilities, area sources such as dry cleaners and gas stations, and motor vehicles. Inputs also include meteorologic and secondary-pollutant formation data. We used total ambient concentrations obtained by adding modeled concentrations to an invariant clean background level which is known only for a few hazardous air pollutants.
We assigned individual exposures using the modeled concentrations corresponding to the census tract of the residential address from the birth certificate. Census tracts were determined using a variety of methods, including the US Census Internet-based address search feature22 (n = 2321), geocoding in ArcGIS 9.2 (n = 155), Internet delivery databases of the US Postal Service23 (n = 214), identifying rural zip codes that were enclosed by census tracts (n = 156), and assigning to a random location within a zip code when other methods were unsuccessful (n = 366). Some methods provided year 2000 census tracts, which were translated to 1990 boundaries using census tract relationship files.24
We examined distributions of hazardous air pollutant concentrations by state and by level of urbanicity. We examined pairwise Spearman rank correlations for all pollutants.
We estimated prevalence odds ratios for each air pollutant and autism spectrum disorders by calculating odds ratios (ORs) and 95% confidence intervals (CIs) using logistic regression. ORs from categorical coding of air pollutants and from quadratic splines of pollutants and autism spectrum disorders suggested monotonic trends. Air pollutant concentrations were log-transformed because resultant ORs were similar in magnitude to ORs with categorical coding of air pollutants but more precise. We calculated ORs corresponding to the odds of autism spectrum disorders for high air pollutant concentrations (80th percentile) over odds estimated for low concentrations (20th percentile). Percentile values were derived from controls from both states combined. All ORs were adjusted for study design factors of surveillance year and state to account for the stratified sampling and the different control sampling fraction in West Virginia.
We used a directed acyclic graph to evaluate covariates for the potential to confound the ORs.25 Adjusted models included suspected autism spectrum disorders risk factors (maternal age26 and smoking in pregnancy27) social class variables potentially related to air pollutants and unknown environmental causes of autism spectrum disorders (marital status, maternal education, race, and census tract median household income) and urbanicity (because air pollution and prevalence of autism28 are higher in urban areas).
Hazardous air pollutants share sources and are therefore correlated with each other. To accurately estimate ORs attributable to a given pollutant, and also to account for multiple comparisons, we simultaneously included all hazardous air pollutants in one semi-Bayes hierarchical model, which improves the plausibility and stability of estimates.29 The semi-Bayes model adjusted the beta coefficient for each hazardous air pollutant toward the mean of its exchangeability group, with the magnitude of shrinkage dependent on the precision of its conventional likelihood estimate and a pre-specified variance of the assumed normal distribution for that parameter (τ2). We selected a τ2 of 0.209 corresponding to expectations that 95% of ORs would fall within a 6-fold range. Semi-Bayes models were fit using SAS IML code developed by Witte et al.30 In this 2-stage hierarchical model, coefficients and covariance matrices from a conventional logistical model were regressed on second-stage data including exchangeability variables and the value for τ2.
We included exchangeability variables expected to predict the magnitude of ORs, considering both determinants of the true OR and determinants of the observed OR, such as exposure misclassification. Hazardous air pollutants were grouped by the following categorical exchangeability variables: an indicator of demonstrated human developmental neurotoxicity,20,31,32 any evidence of placental transfer,33 similarity between NATA model estimates and measured ambient concentrations categorized as low, medium, and high,34 and an indicator that outdoor pollutant concentrations differ from indoor concentrations.35,36 However, including these exchangeability predictors produced results that were similar to those assuming that all estimates arose from one exchangeability group.
We evaluated modification by state because emission inventories are created at the state level and because polluting sources and sample characteristics exhibited variation between North Carolina and West Virginia. We used likelihood ratio tests of cross-products between pollutants and state, and between pollutants and urbanicity, with an alpha level of 0.10 adjusting for a priori demographic confounders in single-pollutant logistic models to evaluate whether estimated prevalence ratios should be combined across these strata. We used likelihood ratio tests and inspection of stratified ORs to examine whether boys or girls were more susceptible to the effects of exposure. We also repeated analyses limiting cases to those with autism spectrum disorders and comorbid intellectual disabilities (n = 161), and autism spectrum disorders without evidence of this comorbidity (n = 234), in case one subgroup was more susceptible to effects of exposure than the other. Finally, we performed a sensitivity analysis of potentially greater misclassification of exposures for children born prior to NATA-1996, by performing analyses separately for North Carolina children born in 1996 and 1994.
In North Carolina, 220 of the 311 children identified with autism spectrum disorders had a matching North Carolina birth certificate. We determined a census tract for all 220, of whom 206 were born in the surveillance counties. In West Virginia, 189 of the 257 children identified with autism spectrum disorders had a West Virginian birth certificate, and we determined a census tract for 177. In North Carolina, 1947 of the 2584 children with speech and language impairment had a matching North Carolina birth certificate. We determined a census tract for 1901, of whom 1733 were born in the surveillance counties. In West Virginia, 4287 children with speech and language impairment were identified. Of these, 1420 were randomly selected for birth certificate matching, 1146 were found, and we determined a census tract for 1096. For North Carolina children with autism spectrum disorders we had evidence that 86% without a matching birth certificate were born out of state, suggesting that most children excluded from further analyses because of a lack of birth certificate match were not part of our intended sample because they were born out of state.
Children with autism spectrum disorders were more likely to be boys, to be firstborn, to have mothers with higher education, and to reside in urban areas, compared with children with speech and language impairment (Table 1). In North Carolina, most children lived in 100% urban or mixed-urbanicity census tracts. In West Virginia, most children lived in 100% rural or mixed census tracts.
We excluded the following air pollutants from further analysis because they had no detectable variation above background levels: carbon tetrachloride, ethylene dibromide, ethylene dichloride, hexachlorobenzene, and polychlorinated biphenyls (Table 2). Coke-oven emissions were present only in West Virginia. Concentrations of all pollutants were higher in urban than nonurban areas. Concentrations were higher in North Carolina than West Virginia for most pollutants even when comparing regions of similar urbanicity. Variability of most pollutants was greater in nonurban areas. Most air pollutants were correlated with all other air pollutants, with Spearman rank correlations of at least 0.5 and often much higher—except for coke oven emissions, which were not correlated with any other pollutant. Air pollutants known to be emitted from mobile sources37 were especially correlated (>0.8).
In unadjusted analyses, ORs were elevated above the null for many pollutants (Fig.). After adjustment for a priori demographic confounders, all estimates were attenuated, so that most adjusted estimates were near-null. For example, the OR for beryllium compounds was 1.5 (95% CI = 1.1–2.0) before adjustment and 1.0 (0.6 –1.4) after adjustment. For xylenes, the unadjusted was 1.5 (1.2–2.0) and the adjusted was 1.0 (0.6–1.6). After adjusting for other air pollutants in semi-Bayes models, ORs were less precise and many remained near-null, but some were shifted upwards or downwards. Estimates shifted upward included those for chromium compounds, dichloropropene, methylene chloride, propionaldehyde, quinoline, styrene, and tetrachloroethane. Estimates shifted downward below the null included ethylbenzene and toluene.
Associations between many pollutants and autism spectrum disorders differed by both state and level of urbanicity when tested in single-pollutant adjusted models. This modification was less pronounced after adjusting for other pollutants in semi-Bayes models (Table 3). Level of urbanicity and state were associated in our sample, with North Carolina having a greater proportion of children living in urban areas than West Virginia. There was a suggestion that urbanicity partly accounted for differences in the magnitude of estimates stratified by state for some pollutants. For example, the ORs for methylene chloride, quinoline, and styrene were more elevated in North Carolina than in West Virginia, and in urban compared with nonurban areas.
Most stable ORs from the fully-adjusted semi-Bayes models were near-null overall, in both states, and within levels of urbanicity (Table 3). Pollutants that had more-precise elevated associations overall and within states and levels of urbanicity were methylene chloride, quinoline, and styrene, whereas hydrazine and lead had negative associations.
Sex of the child satisfied our screening criteria for interaction with many pollutants. After adjusting for other pollutants in semi-Bayes models, confidence intervals of ORs stratified by sex were imprecise for girls and exhibited substantial overlap with those for boys, so that differences could be due to chance. Point estimates were not consistently greater for either boys or girls. Some estimates that showed differences were butadiene (boys, 1.5 [0.7–3.4] and girls, 0.9 [0.3–2.9]); mercury compounds (boys, 1.0 [0.5–1.9] and girls, 1.8 [0.7–4.4]); styrene (boys, 2.2 [1.1–4.1] and girls, 0.7 [0.2–2.6]); and tetrachloroethane (boys, 1.6 [0.5–4.6] and girls, 0.7 [0.1–3.9]).
When models were fit separately for autism spectrum disorders with and without comorbid intellectual disabilities, resultant ORs from semi-Bayes models were mostly within the confidence intervals of the estimate from the other case group, with no clear patterns of increased point estimates for one or the other case group. Some estimates with greater differences were coke oven emissions (with intellectual disabilities, 1.9 [1.1–3.2] and without, 1.0 [0.6–1.6]), dichloropropene (with comorbidity, 1.0 [0.3–3.1] and without, 2.2 [0.8–6.2]), mercury compounds (with comorbidity, 1.7 [0.8–3.6] and without, 0.9 [0.5–1.8]), and methylene chloride (with comorbidity, 1.7 [0.8–3.7] and without, 1.0 [0.5–2.1]).
When North Carolina models were fit separately for children born in 1994 and 1996, ORs were very unstable (results not shown).
We used semi-Bayes analyses to screen perinatal exposure to 35 hazardous air pollutants for a role in contributing to autism spectrum disorders at age 8. These methods can outperform conventional maximum-likelihood-estimation techniques for prediction,38 and are useful in settings where multiple exposures are evaluated for prioritization in future studies.10 After accounting for family and census-tract social characteristics, urbanicity, and potential confounding by correlated air pollutants, we estimated mostly near-null associations, including pollutants that are established human developmental neurotoxicants: polycyclic aromatic hydrocarbons,20 arsenic, lead, manganese, mercury, and toluene.31 Some pollutants stood out with estimates that were precise but elevated compared with those for all other air pollutants—especially quinoline and styrene. Our results are mostly divergent from those from a previous study of hazardous air pollutants and autism spectrum disorders in California,7 except for methylene chloride.
Methylene chloride is a solvent and propellant widely used in manufacturing, and is emitted from point and area sources. Ambient airborne exposure to methylene chloride was positively associated with autism spectrum disorders prevalence in the previous California study.7 We found an elevated OR for methylene chloride, but with a confidence interval including the null. Chronic exposure to methylene chloride has established central nervous system effects.33 A study examining maternal exposure to solvents including methylene chloride found decrements in child cognitive, language, and behavioral functioning.39 Methylene chloride is a pollutant that should be further studied for a role in autism spectrum disorders etiology, using improved exposure assessment.
Unlike the previous study, we estimated elevated associations for quinoline and styrene. The California study did not examine quinoline because of a lack of evidence of developmental toxicity or endocrine-disrupting characteristics of quinoline. Styrene was examined only as part of a composite score that included 5 aromatic solvents, with a resultant null estimate. It is possible that any effect of styrene in the previous study was diluted by its inclusion in a group with other pollutants, if the other pollutants had null or inverse associations with autism spectrum disorders. Neurodevelopmental effects of quinoline and styrene have not been previously reported. Quinoline is a chemical intermediate used in organic processing such as petroleum refining and coal mining, with emissions primarily from area sources. Traffic is the primary source of styrene emissions.
We estimated near-null associations for several pollutants that had positive associations in previous studies: cadmium, mercury, nickel, trichloroethylene, and vinyl chloride.7,40,41 Estimates for cadmium, mercury, and nickel were elevated prior to adjustment in our study but were attenuated upon adjustment. This suggests that the earlier positive associations for these metals may have been due to confounding by urbanicity, other pollutants, or unknown factors associated with urbanicity or pollutant concentrations. In contrast, trichloroethylene and vinyl chloride estimates were near-null both before and after adjustment for confounding in our study. Low-exposure levels in our study sample could have contributed to null findings, especially for vinyl chloride, for which outdoor concentrations were 10 times higher in the California sample. Overall our study areas were characterized by outdoor pollutant concentrations 2–3 times lower than those in the California study. The range of ambient concentrations of these chemicals that may impair neurodevelopment is largely unknown for these pollutants, limiting our ability to discern the impact of studying potentially low levels. In contrast, null findings could reflect that inhalation exposure to these chemicals from ambient sources is not as important as other routes of exposure for these chemicals. Several of the included pollutants have substantial sources of human exposure from diet, tobacco smoke, and indoor environments.
Differences in results between our study and the California study may also have been influenced by the different control groups. Our control group had similar access to developmental testing and met the same residential requirements as the cases, thereby providing an estimate of exposure for the population source that gave rise to the case group. The age 8 address or access to developmental testing services could not be confirmed for children included from the birth certificate rosters used in the previous study.7 Speech and language impairment was chosen over other disorders included in the ADDM surveillance network because, under this protocol, children designated with speech and language impairment but no other developmental problems never met surveillance criteria for autism spectrum disorders. In addition, speech and language impairments are relatively mild and more common than autism spectrum disorders, and therefore, provide a sufficiently-sized control sample.42
We assumed that isolated speech and language impairments were not caused by exposure to any of the included pollutants. This assumption is consistent with one study examining maternal occupational exposure and speech and language impairment that estimated near-null associations with solvents, pesticides, dusts, and metals.43 However, environmental contributions to speech and language development have been scarcely studied,43 and it remains possible that impairments may be influenced by air pollutants. The below-null estimates for lead and toluene in our study could be due to a positive association between these chemicals and speech and language impairment to a greater extent than for autism spectrum disorders. Furthermore, ORs from our study may underestimate the true association with autism spectrum disorders for any pollutant positively associated with speech and language impairment.
It is possible that only an unknown subgroup of children was susceptible to a given air pollutant. We had little information to define such subgroups, and for the few variables we did have (gender and comorbid intellectual disabilities), we had limited power to detect modification.
Diagnosis of autism spectrum disorders is complex, and detection may vary geographically with access to health or educational services.40 This can complicate studies of autism spectrum disorders when the geographic gradient of environmental exposures may correlate with factors promoting case ascertainment. Our study is less likely to have been affected by an ascertainment bias, in part because our case group was not required to have a prior diagnosis and because case and control selection required similar access to developmental testing (as evidenced by their designation of speech and language impairment). We additionally controlled for urbanicity, which may serve as a proxy of increased access to health services, and is related to autism spectrum disorders prevalence and hazardous air pollutant concentrations.
We lacked information on personal exposure to pollutants during early pregnancy, and so assigned individual exposure to each pollutant using output from a temporally and spatially averaged air pollution model (NATA). NATA addresses hazardous air pollutants of highest health concern and greatest mass of emissions. NATA uses extensive national and local emissions inputs and consequently is more complete than exposure-assignment methods that use only proximity to roads or large polluting sources. Our design took advantage of the ability of NATA to distinguish persons with high exposures from those with low exposures, relying on spatial contrasts more than temporal contrasts.
NATA 1996 output has been found to underpredict monitored ambient concentrations, but to preserve rank order of pollutant level.21,35,44–46 Model accuracy varies by pollutant.46 Discrepancies between model output and measured ambient concentrations are thought to be due primarily to the completeness of emissions inventory inputs. For example, not all businesses report to the Toxic Release Inventory, which is an important source of industrial emissions model inputs. However, inputs for the 1996 model are more complete than those for previous years.19
NATA model performance varies regionally across the United States.21 Both North Carolina and West Virginia were ranked as having the highest level of participation in reporting to the National Toxics Inventory with the lowest proportion (<10%) having uncertain location.46 Still, conclusions about accuracy for our sample are limited by the lack of monitoring of pollutants in North Carolina and West Virginia for the relevant time period, precluding validation studies in these states.
Even where NATA reflects outdoor concentration rankings, there are additional levels of uncertainties in using these outputs to characterize individual exposure. Our design assigned the average concentration for the mother’s residential census tract. This does not account for daily activity patterns that include leaving the tract, geocoding errors, or spatial variability in pollutant concentrations within the tract. All of these can be important determinants of individual exposure.47 Greater errors will arise in assigning exposures for persons residing in high-pollutant areas, especially for pollutants dominated by point sources.46 The ability of model output to serve as a surrogate of personal exposure has been evaluated directly and deemed reasonable for most air pollutants.35 Ability to reflect personal exposure is especially good for mobile source volatile organic compounds (eg, benzene) and less accurate for volatile organic compounds with substantial indoor sources (eg, chloroform).35
Exposure misclassification also arises from temporal averaging. NATA output pertains to an average for 1996, whereas the period of susceptibility may be a small window during gestation. Furthermore, our design assumed that the rank order of air pollutants remained stable from 1992 through 1996, the periods of birth for our subjects. Overall, concentrations of air pollutants during this period were decreasing due to regulations that address point sources and on-road vehicle emissions.48 Concentrations could have changed during these periods in North Carolina or West Virginia, for example, due to the opening or closing of a large polluting source or changes in mountain-top coal removal in West Virginia. We lacked the statistical power to determine if our results were sensitive to potentially greater misclassification for earlier birth years.
We attempted to account for uncertainties in the ability of NATA-model estimates to reflect personal exposure by including 2 exchangeability variables in our semi-Bayes models: similarity between modeled and measured concentrations, and similarity between outdoor and indoor concentrations for a given pollutant. Inclusion of these variables did not influence results. However, these exchangeability variables had limitations. For example, accuracy rankings focused on the ability to predict pollutant concentrations at the mean instead of the ability to differentiate high-exposure from low-exposure areas.35,36
The use of NATA to assign exposure in our study is certainly the largest limitation. However, because autism spectrum disorder cases and speech and language impairment controls met similar residency requirements and were derived from the same data sources, it is unlikely that this misclassification differed by outcome status. Hence, exposure misclassification has likely attenuated the reported estimates toward the null, so that true effects may be larger. In addition, exposure misclassification limited our ability to control for confounding due to correlated pollutants. Simultaneously adjusting for multiple pollutants measured at the ecologic level may have produced unexpected results. For example, outdoor measurements of some criteria pollutants are more correlated with personal exposure to another pollutant than the measured pollutant.49 If this were true in these data, it is possible that we may have inadvertently adjusted a given pollutant for a proxy of itself, thereby biasing estimates toward the null.
In summary, several facets of our study may have produced estimates biased toward the null, including the use of a control group that may have associations with the exposure, exposure misclassification, and our use of semi-Bayes shrinkage estimators. For these reasons our design may have failed to screen in chemicals that could potentially have a role in autism etiology. Despite these limitations, we estimated elevated associations for some hazardous air pollutants. Our results are consistent with previous findings regarding methylene chloride7 and suggest quinoline and styrene as possible targets for follow-up in future investigation of autism etiology.
We thank Barbara Becker-Cottrill, Julie O’Malley, and Tom Leonard for their help with West Virginia surveillance and birth certificate data, Cynthia Cassell and Bob Meyer for assistance with North Carolina birth certificate data, Ted Palma for information regarding the NATA program, and Kristen Rappazzo for geocoding assistance.
Supported in part by grants from the Centers for Disease Control Prevention and the National Institute of Environmental Health Sciences (P30ES10126 and T32 ES007018).
The findings and conclusions in this report are those of the authors and do not necessarily represent the official position of the Centers for Disease Control and Prevention.