Association between arsenic, cadmium, manganese, and lead levels in private wells and birth defects prevalence in North Carolina: a semi-ecologic study
BMC Public Health volume 14, Article number: 955 (2014)
Toxic metals including arsenic, cadmium, manganese, and lead are known human developmental toxicants that are able to cross the placental barrier from mother to fetus. In this population-based study, we assess the association between metal concentrations in private well water and birth defect prevalence in North Carolina.
A semi-ecologic study was conducted including 20,151 infants born between 2003 and 2008 with selected birth defects (cases) identified by the North Carolina Birth Defects Monitoring Program, and 668,381 non-malformed infants (controls). Maternal residences at delivery and over 10,000 well locations measured for metals by the North Carolina Division of Public Health were geocoded. The average level of each metal was calculated among wells sampled within North Carolina census tracts. Individual exposure was assigned as the average metal level of the census tract that contained the geocoded maternal residence. Prevalence ratios (PR) with 95% confidence intervals (CI) were calculated to estimate the association between the prevalence of birth defects in the highest category (≥90th percentile) of average census tract metal levels and compared to the lowest category (≤50th percentile).
Statewide, private well metal levels exceeded the EPA Maximum Contaminant Level (MCL) or secondary MCL for arsenic, cadmium, manganese, and lead in 2.4, 0.1, 20.5, and 3.1 percent of wells tested. Elevated manganese levels were statistically significantly associated with a higher prevalence of conotruncal heart defects (PR: 1.6 95% CI: 1.1-2.5).
These findings suggest an ecologic association between higher manganese concentrations in drinking water and the prevalence of conotruncal heart defects.
The toxic metals arsenic, cadmium, manganese, and lead are known developmental toxicants in humans [1–4]. These metals are able to cross the placental barrier resulting in in utero exposures to the developing fetus [5–7]. Some of the detrimental health outcomes associated with exposure to these metals include spontaneous abortion, stillbirth, low birth weight, preterm birth, reduced fetal growth, impaired neurodevelopment, and congenital malformation [8–15].
Birth defects are a leading cause of infant mortality in the U.S. [16, 17], yet 60-70% have no known cause . In epidemiologic studies, toxic metals including arsenic, cadmium, and lead have been associated with congenital structural malformations (reviewed in ). Exposure to arsenic- and lead-contaminated drinking water has been associated with an increased occurrence of congenital heart defects (CHDs) [19–21] as well as neural tube defects (NTDs) in human populations . While these studies demonstrated elevated associations, some of the reported 95% confidence intervals were consistent with unity. In contrast, the teratogenicity of arsenic, cadmium, manganese, and lead has been demonstrated in animal models, specifically affecting cardiac, limb, musculoskeletal, craniofacial, and central nervous system development [10, 23–26]. It is clear from animal studies that metal exposure can result in structural malformation; however, there have been few population-based human studies that have examined such relationships.
Diet and drinking water are common sources of metal exposure in non-occupationally exposed individuals. In North Carolina, 2.3 million residents use private wells for drinking water , and the water quality of these wells is not federally regulated . Notably, several statewide assessments of North Carolina private wells have indicated geographical regions with naturally-occurring manganese  and arsenic [30–32]. Arsenic levels in some wells can reach up to ~800 ppb a level that is 80 times the current EPA Maximum Contaminant Level (MCL) standard (10 ppb) in public distributions systems and clearly associated with disease risk . These data support that ingestion of metal-contaminated drinking water may be a public health concern in North Carolina. No previous studies in North Carolina have been conducted to investigate the potential association between birth defects and metal concentrations in private well water.
In the present study we examined private well water levels of arsenic, cadmium, manganese, and lead across North Carolina, and used a semi-ecologic study design to estimate the association between metal levels and specific birth defect phenotypes. As individual level data on maternal water consumption were not available, maternal drinking water metal levels were assigned based on geocoded residence within a census tract. The prevalence of twelve specific defect types or groups of defects was compared between census tracts with high versus low average metal levels. This work represents the first effort to assess the association between statewide levels of metals in North Carolina wells with birth defects.
Study design and study population
This study was reviewed by the UNC IRB (#11-0414) and was approved according to regulatory category 5 under the National Institutes of Health description on Research on Human Specimens (44 CFR 46.101(b)). A statewide semi-ecologic study was conducted from a baseline study population that included livebirths between January 2003 and December 2008 that comprised 24,704 infants with birth defects (cases) and 725,690 non-malformed controls in North Carolina. Infants from non-singleton births (n = 25,069), without a geocoded residence at delivery (n = 38,206), or case infants with known chromosomal abnormalities (n = 1,472) were considered ineligible and excluded from this study.
Individual-level data regarding maternal source of drinking water and frequency of consumption were not available. Thus, the census tract corresponding to maternal residence at delivery served as the unit of analysis in this semi-ecologic study. The analysis included 1,563 census tracts comprising 5,271 block groups across North Carolina. Private well water metal data were provided by the North Carolina Department of Health and Human Services (DHHS) within the Division of Public Health. The well water database contained metal concentrations collected from private wells across the state between 1998 and 2010 measured by inductively coupled plasma-mass spectrometry. This database contains historical records of well water samples collected by and analyzed at the State Laboratory of Public Health in Raleigh, North Carolina as we have previously summarized . State or federal laws do not require the testing of existing private wells in North Carolina, though state law requires testing of new wells since 2008. As such, the data in this study represent metal analyses from newly constructed private wells since July 2008 (according to North Carolina rule 15A NCAC 18A.3802) and wells tested between 1998 and 2010 where owners requested testing of existing wells. The total number of wells tested was between 22,000 and 76,000 wells depending on the metal that was analyzed.
Case infants with selected birth defects were identified by the North Carolina Birth Defects Monitoring Program (BDMP). The BDMP is a statewide active surveillance program covering all 100 North Carolina counties . Cases were livebirths between 2003 and 2008 identified through systematic review and abstraction of medical records by BDMP field staff. Diagnoses were confirmed by supporting documentation in the medical record, such as surgical or autopsy reports, medical imaging, and physical exams. Selected birth defects in this analysis included a spectrum of phenotypes comprising NTDs, CHDs, gastrointestinal, genitourinary, and musculoskeletal defects that are of prioritized interest by the Centers for Disease Control and Prevention. Non-malformed controls infants were identified from birth certificate records during the same years.
The following twelve structural defects or groups of defects were included in this study: (1) spina bifida without anencephaly (n = 218); (2) anotia and microtia (n = 94); (3) conotruncal heart defects including common truncus, Tetralogy of Fallot (TOF), and transposition of the great arteries (TGA) (n = 435); (4) atrioventricular septal defects (AVSD) and endocardial cushion defects (ECD) (n = 150); (5) hypoplastic left heart syndrome (HLHS) (n = 142); (6) cleft palate (CP) (n = 351); (7) cleft lip with or without CP (n = 516); (8) esophageal atresia (EA) and tracheo-esophageal fistula (TEF) (n = 140); (9) pyloric stenosis (n = 1,204); (10) reduction defects of the upper and lower limbs (n = 255); (11) gastroschisis (n = 215); and (12) hypospadias (n = 1,994).
Environmental metals of interest included arsenic, cadmium, manganese, and lead. As described above, more than 22,000 measurements of metal concentrations in private wells were collected by the North Carolina DHHS across the state between 1998 and 2010. Only private wells with a corresponding global positioning system (GPS) coordinate or with a complete geocoded address were considered in this study (50-67% of wells sampled statewide depending on the metal analyzed) to minimize the impact of geocoding error. GPS coordinates were standardized to a decimal degree format and private well locations were geocoded according to previously described methods . The number of geocoded wells varied by metal where 46,286 (67%) wells had arsenic measures, 11,381 (50%) wells had cadmium measures, 46,022 (65%) wells had manganese measures, and 46,059 (60%) wells had lead measurements. Geocoded well locations were then assigned to North Carolina census tracts that corresponded to the 2000 census using ESRI ArcGISTM 9.3 (Redlands, CA).
The average level of each of the four metals (arsenic, cadmium, manganese, and lead) was calculated within each census tract across the state. The distributions of average metal levels were mapped using ArcGIS. Census tracts with fewer than ten well measures for each contaminant over the total time period were excluded from further analysis to reduce exposure misclassification in poorly sampled tracts. Approximately 50% of census tracts were excluded due to limited sampling.
Maternal residence at delivery for cases and controls was geocoded by the North Carolina State Center for Health Statistics. Cases and controls were assigned to census tracts in North Carolina based on geocoded maternal residence at delivery using ArcGIS. A dichotomous exposure contrast was developed to compare “exposed” and “unexposed” census tracts at delivery. First, the distribution of average metal concentrations across all tracts with at least ten well water measures was examined, and the 50th and 90th percentiles of average concentration was calculated for each metal. Then, an exposure status was assigned to each census tract for each metal where residences were considered “exposed” if the average concentration of the metal within the corresponding tract was equal to or greater than the 90th percentile of metal concentration across all tracts combined. Likewise, tracts were considered “unexposed” if the average concentration of the metal within the corresponding tract was less than or equal to the 50th percentile of metal concentration across all tracts combined. Thus, the final dichotomous exposure contrast compared exposed census tracts (a tract where metal concentrations exceeded the 90th percentile of the distribution of metal concentrations across the state) and unexposed census tracts. The 50th and 90th percentile of census tract average levels were as follows: 1.29 and 2.54 ppb for arsenic, 0.54 and 1.82 ppb for cadmium, 41.55 and 139.69 ppb for manganese, and 3.52 and 7.28 ppb for lead.
Maternal and infant characteristics were obtained for cases and controls from the North Carolina birth certificate. Covariates were selected a priori using a directed acyclical graph approach. The final set included in the statistical analysis were maternal age (continuous), race/ethnicity (categorized as non-Hispanic white, non-Hispanic Black, Hispanic, or Other), and education (categorized as mother’s highest level of completed education as less high school, high school, or greater than high school).
Crude and adjusted estimates of the association between metal concentrations in drinking water and the prevalence of each birth defect within census tracts were calculated by log-linear regression using SAS 9.3 (SAS Institute Inc., Cary, North Carolina). Prevalence ratios (PRs) with 95% confidence intervals (CIs) were calculated to estimate the association between dichotomized metal exposure status for each metal and the prevalence of birth defects within census tracts, adjusted for maternal age at delivery, race/ethnicity and education status. Any subjects with missing data for one or more covariates were excluded from the adjusted statistical model. To account for multiple testing, Bonferroni-corrected PRs with 99.9% CI estimates (alpha = 0.001) were also calculated.
To examine potential inter-metal correlation, Spearman’s rank correlation was used to examine the relationship between each pairwise metal level averaged to the ecologic unit. Corresponding Spearman’s rank coefficients (r) and p-values were calculated.
Given the correlations observed between levels of arsenic and manganese, a multiplicative interaction effect was estimated for the association between residence in areas of both higher arsenic and manganese levels and the prevalence of each birth defect category. Specifically, the prevalence of birth defects among infants residing in census tracts where both arsenic and manganese levels were greater than or equal to the 90th percentile exposed category were compared to areas where both arsenic and manganese were less than or equal to the 50th percentile unexposed category. Crude and adjusted PRs and corresponding 95% CIs were calculated using log-linear regression as described above.
A sensitivity analysis was conducted using the census block group as the ecologic unit of analysis to determine whether the pattern of findings observed at the tract level was robust. The analysis included 5,271 block groups across North Carolina. Although the smaller block group level as the unit of exposure resulted in a finer geographical exposure assessment, fewer study subjects (roughly 50% of the cohort included with the tract ecologic unit analysis) were included since a larger proportion of the block groups did not have adequate well sampling data (greater than 75%). Maternal residence at delivery and the private well locations were assigned to North Carolina block groups corresponding to the 2000 census. Exposure levels for each metal were calculated for each block group using the same methods as the tract level analysis. The resulting 50th and 90th percentiles of block group average levels were 1.24 and 2.75 ppb for arsenic, 0.50 and 1.90 ppb for cadmium, 41.85 and 163.00 ppb for manganese, and 3.32 and 7.55 ppb for lead. PRs and corresponding 95% CIs were calculated using log-linear regression as for the main analysis.
To further refine the exposure assessment, a secondary sensitivity analysis was conducted using a nested cohort of individuals residing outside of public water supply distribution areas. Using North Carolina public water supply area service data and the spatial overlay tool in ArcGIS, the nested analysis was limited to only those subjects with a geocoded maternal residence at delivery outside of a public water supply area. This selection restricted the analysis to individuals who are more likely to use private drinking water wells (n = 117,203 or 17% of the cohort included in the main analysis). PRs and corresponding 95% CIs were calculated identically to the main analysis using the tract and block group ecologic units.
Study population characteristics
The final eligible study population for analysis included 20,151 case infants and 668,381 non-malformed control infants. Table 1 shows the distribution of demographic characteristics of the mothers and infants. Case and control subjects were similar with respect to age, race, sex, and maternal education status. There were a few notable differences by specific birth defect phenotype. Mothers of infants with anotia/microtia were more likely than controls to be Hispanic, and mothers of infants with gastroschisis were more likely to be young and have lower education compared to controls (data not shown). Prevalence of specific birth defects by maternal race/ethnicity has been reported previously by the North Carolina BDMP .
Metal levels elevated in some private wells
From the well water database records collected between 1998 and 2010, the analysis included between 11,000 and 47,000 wells with geocoded locations (See Methods; Additional file 1: Table S1). Notably, private well metal levels that exceeded the EPA regulatory guidelines for public drinking water were detected in 1,121 (2.4%) wells for arsenic, 16 (0.1%) wells for cadmium, 9,424 (20.5%) wells for manganese, and 1,434 (3.1%) wells for lead (Additional file 1: Table S1).
The geocoded levels of arsenic, cadmium, manganese, and lead in private wells were summarized to average census tract levels (Figure 1). The range of average tract levels was 0.5 to 20.4 ppb for arsenic, 0.5 to 3.0 ppb for cadmium, 15.0 to 1116.3 ppb for manganese, and 2.5 to 1304.2 ppb for lead. The average level of arsenic, manganese, or lead exceeded EPA drinking water guideline levels in at least one census tract (Additional file 1: Table S1). Notable trends of differential metal levels are evident across the state. Both arsenic and manganese levels were elevated around the central area of the state and these levels were correlated with each other (r = 0.21, p < 0.001 for tract-level correlation; Additional file 1: Table S2). Additionally, cadmium levels in wells were elevated in areas near cities (e.g. Raleigh, Charlotte), whereas lead levels appeared to be randomly spatially distributed, had a high coefficient of variation, and were not strongly correlated with other metals (Additional file 1: Table S2). The average levels of lead showed the greatest variation among tracts statewide compared to the variation of the average levels of manganese, arsenic or cadmium (Additional file 1: Table S2).
To examine the spatial occurrence of each birth defect category within census tracts, a dichotomous exposure status was assigned based on 50th and 90th percentiles of average concentration for each metal (see Methods). The patterns of these levels are represented in Figure 1. This dichotomous exposure assignment was developed to sharpen the exposure contrast by emphasizing the tails of the distribution of metal concentrations in our study area and was developed independently from EPA MCL guidelines for metals in public drinking water supplies. Table 2 reports the frequency of case and control infants among the exposed category for each metal. Overall, a similar proportion of case and control infants resided in exposed census tracts for each metal.
Crude and adjusted PRs and 95% CIs were calculated to estimate the association between average concentration of each metal and the prevalence of birth defects within census tracts (Additional file 1: Table S3; Table 3). Table 3 reports the adjusted PR estimates comparing the prevalence of birth defects between the metal-exposed versus unexposed tracts. To facilitate the interpretation of results, this study considered positive effect measures >1.3 or for which the 95% CI exclude the null (1.0) to be noteworthy potential metal-defect associations. The arsenic-exposed category was associated with a statistically significantly decreased prevalence of HLHS (PR: 0.3 95% CI: 0.1-1.0) compared to the unexposed category. Two additional notable elevated PRs were observed for arsenic-exposed versus unexposed and the prevalence of cleft lip with or without cleft palate (PR: 1.3 95% CI: 0.9-1.9) and EA/TEF (PR: 1.3 95% CI: 0.4-4.3), but the observed 95% CIs were consistent with a null result.
The manganese-exposed category was statistically significantly associated with a higher prevalence of conotruncal heart defects (PR: 1.6 95% CI: 1.1-2.5) compared to the unexposed category. An elevated PR was also observed for manganese-exposed versus unexposed categories and the prevalence of anotia/microtia (PR: 1.3 95% CI: 0.5-3.4) and AVSD/ECD (PR: 1.3 95% CI: 0.6-2.5), however the CIs were consistent with unity.
The lead-exposed category was associated with an increased prevalence of HLHS (PR: 1.7 95% CI: 0.9-3.3), but the finding was not statistically significant. The cadmium-exposed group was associated with a significantly decreased prevalence of pyloric stenosis (PR: 0.4 95% CI: 0.3-0.7) compared to the unexposed group. The findings did not suggest patterns of association with cadmium-exposure and an increased prevalence of any birth defect analyzed in this study.
To account for multiple comparisons, Bonferroni correction was applied and 99.9% CIs are reported (Additional file 1: Table S4). The inverse association between cadmium exposure and pyloric stenosis was the only observed association that remained statistically significant after Bonferroni correction (PR: 0.4 99.9% CI: 0.2-1.0).
The potential interaction between joint exposure to the exposed categories of arsenic and manganese and the prevalence of each birth defect was examined. Adjusted PR estimates and 95% CIs showed elevated (though not statistically significant) point estimates for the interaction effect of arsenic and manganese exposure with defect outcomes including anotia/microtia (PR: 2.7 95% CI: 0.7-10.2), conotruncal defects (PR: 1.6 95% CI: 0.8-3.2), AVSD/ECD (PR: 1.5 95% CI: 0.5-4.3) and cleft lip with or without cleft palate (PR: 1.4 95% CI: 0.8-2.6) (Additional file 1: Table S5).
Sensitivity analyses results
To assess the impact of potential exposure misclassification in the ecologic analysis, two sensitivity analyses were conducted (Additional file 1: Tables S6-S8). For both types of sensitivity analyses, the PR estimates had wider CIs due to the reduction in sample size. The PRs and 95% CIs for selected relationships assessed in the sensitivity analyses are shown in Figure 2. The findings suggest patterns of association consistent with the highest category of manganese exposure and conotruncal heart defects (Figure 2A). For example, the association between manganese exposure and conotruncal heart defects was estimated under six conditions including the tract ecologic unit crude PR (95% CIs) of 1.6 (1.1-2.4) (Additional file 1: Table S3) and adjusted PR (95% CI) of 1.6 (1.1-2.5) (Table 1), which increased to 2.0 (0.9-4.3) (Additional file 1: Table S6) when restricted to individuals outside of public supply areas. The same estimates using the block group ecologic unit estimated a crude PR (95% CI) of 1.9 (1.1-3.4), an adjusted PR (95% CI) of 1.9 (1.1-3.4) (Additional file 1: Table S7), which increased to 2.8 (1.2-6.4) (Additional file 1: Table S8) when restricted to individuals outside of public supply areas. Likewise, the association between arsenic exposure and conotruncal heart defects showed a suggestive increasing trend in estimates (Figure 2B; Additional file 1: Tables S6-S8). Other metal-defect associations, such as the relationships between arsenic and cleft lip, arsenic and HLHS, or cadmium and pyloric stenosis, or lead and HLHS, however, show trends suggestive of a null relationship (Figure 2C-F).
In the present study a statewide semi-ecologic analysis was conducted to investigate the potential association between metal concentrations in private well water and birth defects prevalence in North Carolina. Specifically, we examined the association between the prevalence of twelve birth defect groups and levels of arsenic, cadmium, manganese, and lead summarized at the census tract ecologic unit. Taken together, our analyses identified three notable findings: (i) metal levels are elevated in private wells across the state and some metals such as arsenic and manganese co-occur, (ii) residence in manganese-exposed census tracts was statistically significantly associated with increased prevalence of conotruncal heart defects, and (iii) there is suggestive evidence of an interactive effect between residence in arsenic- and manganese-exposed census tracts and the prevalence of conotruncal heart defects.
Elevated environmental metals are spatially correlated in private wells
North Carolina is the fourth largest state population with approximately 2.3 million people relying on private wells for drinking water; the total number of individuals served is larger only in Pennsylvania, California, and Michigan as estimated by the USGS . Arsenic has been previously recognized as a contaminant of concern in certain regions of North Carolina [30–32]. The present study extends our previous findings reported in  to identify census tract and block group levels of average arsenic, manganese, and lead that exceeded respective EPA standards for public distribution systems, and to map spatial trends of these metals in wells across the state.
The levels of metals in private wells is not currently regulated in North Carolina, and routine testing of private well water quality for contaminants is not required by state or federal regulations (state regulation requires a single test only when new wells are constructed). The MCL for arsenic levels in public water supplies is 10 μg/L, and 5 μg/L for cadmium . Lead is regulated according to a Treatment Technique (TT) designed to control corrosiveness that requires that no more than 10% of a system’s water samples measured at consumer taps can exceed 15 μg/L . Manganese has a recommended Secondary Maximum Contaminant Level (SMCL) of 50 μg/L in public water supplies that is based on taste and odor and is not federally enforceable . Based on these standards, 1,121 (2.4%) of arsenic, 16 (0.1%) of cadmium, 1,434 (3.1%) of lead, and 9,424 (20.5%) of manganese samples from individual wells tested between 1998 and 2010 and included in our study exceeded guideline levels. The extent of samples exceeding the SMCL for manganese is substantial, and comparable to a similar proportion of wells (~21%) that exceeded 50 μg/L in a nationwide study . The arsenic data are supported by our previous findings that demonstrated arsenic contamination in areas along a geographic area known as the Carolina terrane . Levels of cadmium or lead in private wells in North Carolina have not been previously reported.
Spatial trends of private well metal contamination were identified across North Carolina. Here we find that manganese levels correlate spatially with arsenic levels in private wells. Levels of arsenic and manganese have been described in North Carolina groundwater [29–32]; however, a spatial correlation has not been previously reported with respect to private well levels. Elevated cadmium levels were evident near major cities including Raleigh and Charlotte. Urbanization and anthropogenic activities are contributing factors to groundwater contamination from metals including cadmium and arsenic . Lead can occur at elevated levels in drinking water by leaching from certain faucets, fittings, or water systems , which may explain the random spatial appearance of levels across the state. Despite the lack of state and federal regulation, it is recommended that individuals relying on private wells for drinking water have their water quality tested annually.
Manganese levels associated with heart defects prevalence
Our findings support a potential association between residence in manganese-exposed census tracts and the prevalence of conotruncal heart defects. The results of the sensitivity analyses showed similar positive trends and statistically significant associations that are suggestive of a potential association between manganese levels and conotruncal defects. To our knowledge, to date no studies have examined the relationship between drinking water manganese and birth defects in humans. A previous study reported that increased placental manganese levels were associated a higher risk of NTDs showed a statistically significant association with spina bifida . In animals, both increased and decreased levels of manganese have been associated with teratogenicity (reviewed in ). Furthermore, manganese exposure has been associated with several adverse newborn outcomes including infant mortality [29, 41], intrauterine growth restriction , and lower birth weight [11, 43]. As manganese can cross the placental barrier and bioaccumulate in fetal and neonatal tissue [5, 7, 44, 45], the potential association between manganese and birth defects is biologically plausible. Manganese exposure is associated with adverse cardiovascular and central nervous system effects in adult populations [3, 46]; however, its role in abnormal prenatal cardiac development requires further study. Dietary, occupational, and industrial exposures, as well as inhalation of manganese aerosols during showering remain important considerations that contribute to total manganese intake in addition to drinking water [3, 47, 48]. These data suggest that additional studies examining the health effects of manganese as a developmental toxicant with potential exposure from drinking water are warranted.
Arsenic levels were modestly associated with birth defects prevalence
The results suggest that residence within arsenic-exposed census tracts was modestly associated with a higher prevalence of cleft lip with or without cleft palate and EA/TEF defects at the census tract ecologic unit. Residence in arsenic-exposed tracts was also significantly associated with decreased prevalence of HLHS at the census tract ecologic unit, although the sensitivity analyses were consistent with a trend towards a null association. Previous population-based studies have shown associations between drinking water arsenic levels and CHDs [19–21], NTDs , and increased prevalence of all birth defects . It should be noted that for two of these studies [21, 22], although the point estimates were elevated in relation to higher concentrations of arsenic in drinking water, the 95% CIs were consistent with the null. While there is strong evidence of the teratogenic effects of arsenic compounds in animal models, currently inorganic arsenic teratogenicity in humans remains debated [24, 26].
Our results show a potential interactive effect between arsenic and manganese on the association with spina bifida, conotruncal heart defects, AVSD heart defects, and cleft lip with or without cleft palate. Specifically, the estimates reveal that the prevalence of these defects among infants residing in census tracts that were considered both arsenic- and manganese-exposed compared to areas that were arsenic- and manganese-unexposed were increased compared to the PR estimates observed for either metal alone. The findings suggest a potential, but not statistically significant, joint effect in areas of arsenic and manganese co-occurrence. These analyses had limited power as reflected in the width of the reported CIs because the total sample size was considerably reduced. This association has not been previously examined in studies of environmental metal exposures and birth defects. There is previous evidence for potential drinking water co-exposure to both manganese and arsenic. In a study of US public wells, the combination of arsenic and manganese detected in the same water source was one of the most commonly occurring pairs of unique mixtures out of over 90 organic and inorganic contaminants that were tested .
Cadmium and lead levels revealed reduced or null associations
The results suggest some associations between residence in areas of the highest cadmium levels and reduced birth defects prevalence. A decreased prevalence in pyloric stenosis was associated with residence in areas of the highest cadmium exposure, a relationship that remained statistically significant after Bonferroni correction. There were no statistically significant positive relationships between birth defect prevalence and residence in areas of the highest cadmium or lead levels. At the time of this study, we are aware of only three studies that examined the relationship between drinking water levels of cadmium and/or lead and birth defects in humans [21, 22, 51]. Two studies reported null associations between drinking water levels and NTDs [22, 51], and the other reported an increased odds of heart defects associated with detectable lead levels in drinking water . Blood lead levels or occupational lead exposure have been associated with CHDs [52, 53], NTDs [54, 55], and cleft lip [55, 56]. Further, animal studies confirm cadmium- and lead-induced malformations in animal models [57–60]. While counterintuitive, there is evidence for protective effects of maternal cigarette smoking, a major source of cadmium, and orofacial clefts, NTDs and conotruncal heart defects , but an inverse relationship with pyloric stenosis prevalence has not been previously reported. The present study could not account for maternal smoking status due to limited sample size, and this may represent a source of unmeasured confounding for some birth defect phenotypes.
Study limitations and strengths
This study had several limitations that are important when interpreting the results. As individual-level environmental data were not available within the birth defects database, it was not possible to evaluate the relationship between personal water consumption patterns and the risk of defects. Data on individual use of well water were unavailable; therefore, any observed associations are limited to the ecologic unit of analysis and the findings are subject to the ecologic fallacy. However, in the absence of individual-level exposure data, ecologic approaches have previously been applied to examine associations between metals in private wells and other health outcomes [62, 63].
Exposure assessment remains a challenge among epidemiological studies of birth defects [64, 65]. For example, there are no existing biomonitoring programs in North Carolina to actively assess environmental exposures during the prenatal period. As a result, research paradigms usually require retrospective exposure assessment that is clearly limited by issues of temporality. Our recent study suggests that prenatal biomonitoring may be warranted in North Carolina , especially among at-risk populations. Such a biomonitoring program would enable a direct prospective assessment of the relationship between environmental exposures and children born with birth defects.
In this study the assigned exposures represent a census tract average over a period of 12 years and BDMP data represent the location of cases at delivery. It is possible that non-differential exposure misclassification was introduced by the databases including a potentially biased sample of tested wells, a large proportion of tracts with limited sampling, and the likelihood of maternal mobility during pregnancy . It is possible that bias may have been introduced by the private well database since it represents a non-random sample of wells for which the owners requested metal analyses. Similarly, the state testing system does not assign unique well IDs, and some wells likely were tested more than once over the 11-year period. There may also be effects of residual confounding related to geographical location such as socioeconomic status and/or residence in rural versus urban areas that were not accounted for in the analysis. We attempt to address some of these limitations by refining exposure assessment, at the expense of statistical precision, by applying two sensitivity analyses. A sensitivity analysis removed individuals residing in areas served by public water supplies. This resulted in approximately 17% of the initial cohort size. In preliminary data from the National Birth Defects Prevention Study (NBDPS), roughly 22% of mothers residing a subset of 19 North Carolina counties over a similar time period reported having a private well tap water source (data not shown). We interpret the NBDPS data with caution, as the subset of 19 counties may not reflect the appropriate proportion of private well users statewide. Given the lack of individual-level data, this ecologic method served as an imperfect proxy for type of water supply as it potentially excludes a fraction of individuals who use private wells but may be within the public distribution service area. Overall, the results of the sensitivity analyses suggest some metal-defect trends. PR estimates using the block group unit and when restricted to individuals outside of public water utility areas were elevated suggesting that exposure misclassification at the tract level may have biased associations towards the null. We note that in both sensitivity analyses, the wide CIs reflect a loss of precision due to decreased sample sizes.
This study had several strengths including an active statewide case ascertainment program, which provided a large sample size of liveborn cases and controls representing six years of active birth defects surveillance in North Carolina. The semi-ecologic design enabled individual-level assessment of newborns with specific defects and covariate adjustment using data from birth certificates. The outcomes we investigated included specific phenotypes or groups of phenotypes based on similar etiologic mechanisms, which is advantageous over previous studies that have grouped multiple defect outcomes. The statewide surveillance program allowed exclusion of case infants born with chromosomal abnormalities, as this is a known risk factor for some structural defects . This study is the first to assess the association between toxic metal levels in drinking water and birth defects in North Carolina, and builds on a limited body of existing literature on the potential relation between toxic metals and risk of birth defects.
Overall, our findings suggest evidence at the ecologic level of possible relationships between elevated concentrations of arsenic, cadmium, manganese, and lead levels in drinking water and specific birth defects. Given that our results support previous evidence of an association between metal exposure and specific heart defects, these relationships should be prioritized for future studies. This is directly in line with ongoing research directed at understanding the biological mechanisms that underlie metal-birth defect associations using animal and cell culture models . Further studies are necessary to validate the potential associations observed in this study and to better understand the complex environmental, genomic, and epigenetic mechanisms involved in potentially preventable environmentally-mediated birth defects.
Atrioventricular septal defect
Congenital heart defect
Department of Health and Human Services
Endocardial cushion defect
Hypoplastic left heart syndrome
Maximum Contaminant Level
Neural tube defect
Secondary maximum contaminant level
Tetralogy of fallot
Transposition of the great arteries
ATSDR: Toxicological Profile for Arsenic. 2007. 2007, Atlanta, GA: CDC, [cited 2011 July 1]; Available from: http://www.atsdr.cdc.gov/ToxProfiles/TP.asp?id=22&tid=3
ATSDR: Toxicological Profile for Lead. 2007, Atlanta, GA: CDC, [cited 2011 July 1]
ATSDR: Toxicological Profile for Manganese. 2008, Atlanta, GA: CDC, [cited 2011 July 1]
ATSDR: Toxicological Profile for Cadmium. 2008, Atlanta, GA: CDC, [cited 2011 July 1]
Rudge CV, Rollin HB, Nogueira CM, Thomassen Y, Rudge MC, Odland JO: The placenta as a barrier for toxic and essential elements in paired maternal and cord blood samples of South African delivering women. J Environ Monit. 2009, 11 (7): 1322-1330. 10.1039/b903805a.
Salpietro CD, Gangemi S, Minciuollo PL, Briuglia S, Merlino MV, Stelitano A, Cristani M, Trombetta D, Saija A: Cadmium concentration in maternal and cord blood and infant birth weight: a study on healthy non-smoking women. J Perinat Med. 2002, 30 (5): 395-399.
Fechter LD: Distribution of manganese in development. Neurotoxicology. 1999, 20 (2–3): 197-201.
Vahter M: Effects of arsenic on maternal and fetal health. Annu Rev Nutr. 2009, 29: 381-399. 10.1146/annurev-nutr-080508-141102.
Wigle DT, Arbuckle TE, Turner MC, Berube A, Yang Q, Liu S, Krewski D: Epidemiologic evidence of relationships between reproductive and child health outcomes and environmental chemical contaminants. J Toxicol Environ Health B Crit Rev. 2008, 11 (5–6): 373-517.
Bellinger DC: Teratogen update: lead and pregnancy. Birth Defects Res A Clin Mol Teratol. 2005, 73 (6): 409-420. 10.1002/bdra.20127.
Grazuleviciene R, Nadisauskiene R, Buinauskiene J, Grazulevicius T: Effects of elevated levels of manganese and iron in drinking water on birth outcomes. Pol J Environ Stud. 2009, 18 (5): 819-825.
Axelrad DA, Bellinger DC, Ryan LM, Woodruff TJ: Dose–response relationship of prenatal mercury exposure and IQ: an integrative analysis of epidemiologic data. Environ Health Perspect. 2007, 115 (4): 609-615. 10.1289/ehp.9303.
Ericson JE, Crinella FM, Clarke-Stewart KA, Allhusen VD, Chan T, Robertson RT: Prenatal manganese levels linked to childhood behavioral disinhibition. Neurotoxicol Teratol. 2007, 29 (2): 181-187. 10.1016/j.ntt.2006.09.020.
Kim Y, Ha EH, Park H, Ha M, Kim Y, Hong YC, Kim EJ, Kim BN: Prenatal lead and cadmium co-exposure and infant neurodevelopment at 6 months of age: The Mothers and Children's Environmental Health (MOCEH) study. Neurotoxicology. 2012, 35C: 15-22.
Yu XD, Yan CH, Shen XM, Tian Y, Cao LL, Yu XG, Zhao L, Liu JX: Prenatal exposure to multiple toxic heavy metals and neonatal neurobehavioral development in Shanghai China. Neurotoxicol Teratol. 2011, 33 (4): 437-443. 10.1016/j.ntt.2011.05.010.
Petrini J, Damus K, Russell R, Poschman K, Davidoff MJ, Mattison D: Contribution of birth defects to infant mortality in the United States. Teratology. 2002, 66 (Suppl 1): S3-S6.
Mathews TJ, MacDorman MF: Infant mortality statistics from the 2005 period linked birth/infant death data set. Natl Vital Stat Rep. 2008, 57: 1-32.
Brent RL: Environmental causes of human congenital malformations: the pediatrician's role in dealing with these complex clinical problems caused by a multiplicity of environmental and genetic factors. Pediatrics. 2004, 113 (4 Suppl): 957-968.
Zierler S, Theodore M, Cohen A, Rothman KJ: Chemical quality of maternal drinking water and congenital heart disease. Int J Epidemiol. 1988, 17 (3): 589-594. 10.1093/ije/17.3.589.
Engel RR, Smith AH: Arsenic in drinking water and mortality from vascular disease: an ecologic analysis in 30 counties in the United States. Arch Environ Health. 1994, 49 (5): 418-427. 10.1080/00039896.1994.9954996.
Aschengrau A, Zierler S, Cohen A: Quality of community drinking water and the occurrence of late adverse pregnancy outcomes. Arch Environ Health. 1993, 48 (2): 105-113. 10.1080/00039896.1993.9938403.
Brender JD, Suarez L, Felkner M, Gilani Z, Stinchcomb D, Moody K, Henry J, Hendricks K: Maternal exposure to arsenic, cadmium, lead, and mercury and neural tube defects in offspring. Environ Res. 2006, 101 (1): 132-139. 10.1016/j.envres.2005.08.003.
Thompson J, Bannigan J: Cadmium: Toxic effects on the reproductive system and the embryo. Reprod Toxicol. 2008, 25 (3): 304-315. 10.1016/j.reprotox.2008.02.001.
Holson JF, Desesso JM, Jacobson CF, Farr CH: Appropriate use of animal models in the assessment of risk during prenatal development: an illustration using inorganic arsenic. Teratology. 2000, 62 (1): 51-71. 10.1002/1096-9926(200007)62:1<51::AID-TERA10>3.0.CO;2-T.
Gilani SH, Alibhai Y: Teratogenicity of metals to chick embryos. J Toxicol Environ Health. 1990, 30 (1): 23-31. 10.1080/15287399009531407.
Shalat SL, Walker DB, Finnell RH: Role of arsenic as a reproductive toxin with particular attention to neural tube defects. J Toxicol Environ Health. 1996, 48 (3): 253-272. 10.1080/009841096161320.
Kenny JF, Barber NL, Hutson SS, Linsey KS, Lovelace JK, Maupin MA: Estimated use of Water use in the United States in 2005. 2009, Reston, VA: U.S. Geological Survey, 1-60.
EPA: Arsenic in Drinking Water. 2006, Washington, DC: USEPA, [cited 2010 May 24]; Available from: http://epa.gov/safewater/arsenic/index.html
Spangler AH, Spangler JG: Groundwater manganese and infant mortality rate by county in North Carolina: an ecological analysis. Ecohealth. 2009, 6 (4): 596-600. 10.1007/s10393-010-0291-4.
Sanders AP, Messier KP, Shehee M, Rudo K, Serre ML, Fry RC: Arsenic in North Carolina: public health implications. Environ Int. 2012, 38 (1): 10-16. 10.1016/j.envint.2011.08.005.
Pippin CG, Butczynski M, Clayton J: Distribution of Total Arsenic in Groundwater in the North Carolina Piedmont, D.o.E.a.N.R.-G. Section. 2003, Mooresville, NC: NCDENR
Kim D, Miranda ML, Tootoo J, Bradey P, Gelfand AE: Spatial modeling for groundwater arsenic levels in North Carolina. Environ Sci Technol. 2011, 45 (11): 4824-4831. 10.1021/es103336s.
SCHS: North Carolina Birth Defects Monitoring Program. 2013, Available from: http://www.schs.state.nc.us/units/bdmp/
NBDPN: State birth defects surveillance program data. BDRA (Part A). 2013, 97 (10): p. S1-S172.
EPA: Drinking Water Contaminants: National Primary Drinking Water Regulations. 2012, [cited 2013 February 14]; Available from: http://water.epa.gov/drink/contaminants/index.cfm
DeSimone LA, Hamilton PA, Gilliom RJ: The quality of our nation's waters - quality of water from domestic wells in principal aquifers of the United States, 1991-2004- overview and major findings. U S Geol Surv Circular. 2009, 1332: 48-
Bhagure GR, Mirgane SR: Heavy metal concentrations in groundwaters and soils of Thane Region of Maharashtra, India. Environ Monit Assess. 2011, 173 (1–4): 643-652.
EPA: Actions You Can Take to Reduce Lead in Drinking Water. 1993, [cited 2013 June 2]; Available from: http://water.epa.gov/drink/info/lead/lead1.cfm
Liu J, Jin L, Zhang L, Li Z, Wang L, Ye R, Zhang Y, Ren A: Placental concentrations of manganese and the risk of fetal neural tube defects. J Trace Elem Med Biol. 2013, 27 (4): 322-325. 10.1016/j.jtemb.2013.04.001.
Gerber GB, Leonard A, Hantson P: Carcinogenicity, mutagenicity and teratogenicity of manganese compounds. Crit Rev Oncol Hematol. 2002, 42 (1): 25-34. 10.1016/S1040-8428(01)00178-0.
Hafeman D, Factor-Litvak P, Cheng Z, van Geen A, Ahsan H: Association between manganese exposure through drinking water and infant mortality in Bangladesh. Environ Health Perspect. 2007, 115 (7): 1107-1112. 10.1289/ehp.10051.
Vigeh M, Yokoyama K, Ramezanzadeh F, Dahaghin M, Fakhriazad E, Seyedaghamiri Z, Araki S: Blood manganese concentrations and intrauterine growth restriction. Reprod Toxicol. 2008, 25 (2): 219-223. 10.1016/j.reprotox.2007.11.011.
Zota AR, Ettinger AS, Bouchard M, Amarasiriwardena CJ, Schwartz J, Hu H, Wright RO: Maternal blood manganese levels and infant birth weight. Epidemiology. 2009, 20 (3): 367-373. 10.1097/EDE.0b013e31819b93c0.
Spencer A: Whole blood manganese levels in pregnancy and the neonate. Nutrition. 1999, 15 (10): 731-734. 10.1016/S0899-9007(99)00144-6.
Takser L, Lafond J, Bouchard M, St-Amour G, Mergler D: Manganese levels during pregnancy and at birth: relation to environmental factors and smoking in a Southwest Quebec population. Environ Res. 2004, 95 (2): 119-125. 10.1016/j.envres.2003.11.002.
Jiang Y, Zheng W: Cardiovascular toxicities upon manganese exposure. Cardiovasc Toxicol. 2005, 5 (4): 345-354. 10.1385/CT:5:4:345.
Bouchard MF: Manganese in drinking water: Bouchard responds. Environ Health Perspect. 2011, 119 (6): A241-10.1289/ehp.1103485R.
Elsner RJ, Spangler JG: Neurotoxicity of inhaled manganese: public health danger in the shower?. Med Hypotheses. 2005, 65 (3): 607-616. 10.1016/j.mehy.2005.01.043.
Kwok RK, Kaufmann RB, Jakariya M: Arsenic in drinking-water and reproductive health outcomes: A study of participants in the Bangladesh integrated nutrition programme. J Health Popul Nutr. 2006, 24 (2): 190-205.
Toccalino PL, Norman JE, Scott JC: Chemical mixtures in untreated water from public-supply wells in the U.S.–occurrence, composition, and potential toxicity. Sci Total Environ. 2012, 431: 262-270.
Macdonell JE, Campbell H, Stone DH: Lead levels in domestic water supplies and neural tube defects in Glasgow. Arch Dis Child. 2000, 82 (1): 50-53. 10.1136/adc.82.1.50.
Jackson LW, Correa-Villasenor A, Lees PS, Dominici F, Stewart PA, Breysse PN, Matanoski G: Parental lead exposure and total anomalous pulmonary venous return. Birth Defects Res A Clin Mol Teratol. 2004, 70 (4): 185-193. 10.1002/bdra.20014.
Bound JP, Harvey PW, Francis BJ, Awwad F, Gatrell AC: Involvement of deprivation and environmental lead in neural tube defects: a matched case–control study. Arch Dis Child. 1997, 76 (2): 107-112. 10.1136/adc.76.2.107.
Cengiz B, Soylemez F, Ozturk E, Cavdar AO: Serum zinc, selenium, copper, and lead levels in women with second-trimester induced abortion resulting from neural tube defects: a preliminary study. Biol Trace Elem Res. 2004, 97 (3): 225-235. 10.1385/BTER:97:3:225.
Irgens A, Kruger K, Skorve AH, Irgens LM: Reproductive outcome in offspring of parents occupationally exposed to lead in Norway. Am J Ind Med. 1998, 34 (5): 431-437. 10.1002/(SICI)1097-0274(199811)34:5<431::AID-AJIM3>3.0.CO;2-T.
Kristensen P, Irgens LM, Daltveit AK, Andersen A: Perinatal outcome among children of men exposed to lead and organic solvents in the printing industry. Am J Epidemiol. 1993, 137 (2): 134-144.
Webster WS, Messerle K: Changes in the mouse neuroepithelium associated with cadmium-induced neural tube defects. Teratology. 1980, 21 (1): 79-88. 10.1002/tera.1420210110.
McClain RM, Becker BA: Teratogenicity, fetal toxicity, and placental transfer of lead nitrate in rats. Toxicol Appl Pharmacol. 1975, 31 (1): 72-82. 10.1016/0041-008X(75)90053-8.
Thompson J, Bannigan J: Effects of cadmium on formation of the ventral body wall in chick embryos and their prevention by zinc pretreatment. Teratology. 2001, 64 (2): 87-97. 10.1002/tera.1050.
Thompson JM, Bannigan JG: Omphalocele induction in the chick embryo by administration of cadmium. J Pediatr Surg. 2007, 42 (10): 1703-1709. 10.1016/j.jpedsurg.2007.05.026.
Grewal J, Carmichael SL, Ma C, Lammer EJ, Shaw GM: Maternal periconceptional smoking and alcohol consumption and risk for select congenital anomalies. Birth Defects Res A Clin Mol Teratol. 2008, 82 (7): 519-526. 10.1002/bdra.20461.
Liu-Mares W, Mackinnon JA, Sherman R, Fleming LE, Rocha-Lima C, Hu JJ, Lee DJ: Pancreatic cancer clusters and arsenic-contaminated drinking water wells in Florida. BMC Cancer. 2013, 13: 111-10.1186/1471-2407-13-111.
Nieder AM, MacKinnon JA, Fleming LE, Kearney G, Hu JJ, Sherman RL, Huang Y, Lee DJ: Bladder cancer clusters in Florida: identifying populations at risk. J Urol. 2009, 182 (1): 46-50. 10.1016/j.juro.2009.02.149. discussion 51
Katz EA, Shaw GM, Schaffer DM: Exposure assessment in epidemiologic studies of birth defects by industrial hygiene review of maternal interviews. Am J Ind Med. 1994, 26 (1): 1-11. 10.1002/ajim.4700260102.
Borja-Aburto VH, Bermudez-Castro O, Lacasana-Navarro M, Kuri P, Bustamante-Montes P, Torres-Meza V: Difficulties of the methods for studying environmental exposure and neural tube defects. Salud Publica Mex. 1999, 41 (Suppl 2): p. S124-p. S131.
Sanders AP, Flood K, Chiang S, Herring AH, Wolf L, Fry RC: Towards prenatal biomonitoring in North Carolina: assessing arsenic, cadmium, mercury, and lead levels in pregnant women. PLoS One. 2012, 7 (3): e31354-10.1371/journal.pone.0031354.
Canfield MA, Ramadhani TA, Langlois PH, Waller DK: Residential mobility patterns and exposure misclassification in epidemiologic studies of birth defects. J Expo Sci Environ Epidemiol. 2006, 16 (6): 538-543. 10.1038/sj.jes.7500501.
Forrester MB, Merz RD: Impact of excluding cases with known chromosomal abnormalities on the prevalence of structural birth defects, Hawaii, 1986–1999. Am J Med Genet A. 2004, 128A (4): 383-388. 10.1002/ajmg.a.30205.
Ahir BK, Sanders AP, Rager JE, Fry RC: Systems biology and birth defects prevention: blockade of the glucocorticoid receptor prevents arsenic-induced birth defects. Environ Health Perspect. 2013, 121 (3): 332-338. 10.1289/ehp.1205659.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2458/14/955/prepub
This research was funded by grants from the NIEHS (T32-ES007018, P42-ES005948-18, P30-ES010126, and R01-ES019315). The authors thank Drs. Thomas Luben at the US Environmental Protection Agency National Center for Environmental Assessment and Pete Weyer at the University of Iowa Center for Health Effects of Environmental Contamination for their assistance in acquisition and analysis of the public water supply data.
The authors declare that they have no competing interests.
RCF, APS, AFO, and REM conceived of and designed the study. APS, TAD, JLW, AHH, and DE carried out the analyses. APS and RCF wrote the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Table S1: Number of geocoded wells for each metal as well as number of census tract and block groups with average metal levels above EPA drinking water standards for public distribution systems. Averages were calculated where more than ten wells were measured in the respected block group or tract. Table S2. Spearman’s rank correlations and coefficient of variation for average (ppb) census tract metal levels and census block group metal levels (noted in parentheses). Table S3. Crude PRs and 95% CIs using the census tract ecologic unit and comparing individuals in the exposed (≥90th percentile) to unexposed (≤50th percentile) groups for average census tract metal levels. Table S4. Bonferroni corrected adjusted PRs and 99.9% CIs using the census tract ecologic unit and comparing individuals in the exposed (≥90th percentile) to unexposed (≤50th percentile) groups for average census tract metal levels. PRs were adjusted for maternal age, race, and education status. Table S5. Crude and adjusted PRs and 95% CIs estimate arsenic-by-manganese interaction using the tract level ecologic unit. PRs were adjusted for maternal age, race, and education status. Table S6. Adjusted PRs and 95% CIs using the census tract ecologic unit and excluding individuals living in public distribution areas. PRs were adjusted for maternal age, race, and education status. Table S7. Crude and adjusted PRs and 95% CIs using the block group ecologic unit and comparing individuals in the exposed (≥90th percentile) to unexposed (≤50th percentile) groups for average census block group metal levels. PRs were adjusted for maternal age, race, and education status. Table S8. Adjusted PRs and CIs at the block group level excluding individuals living in public distribution areas. PRs were adjusted for maternal age, race, and education status. (DOCX 35 KB)
About this article
Cite this article
Sanders, A.P., Desrosiers, T.A., Warren, J.L. et al. Association between arsenic, cadmium, manganese, and lead levels in private wells and birth defects prevalence in North Carolina: a semi-ecologic study. BMC Public Health 14, 955 (2014). https://doi.org/10.1186/1471-2458-14-955