Molecular and neurodevelopmental benefits to children of closure of a coal burning power plant in China.

Polycyclic aromatic hydrocarbons (PAH) are major toxic air pollutants released during incomplete combustion of coal. PAH emissions are especially problematic in China because of their reliance on coal-powered energy. The prenatal period is a window of susceptibility to neurotoxicants. To determine the health benefits of reducing air pollution related to coal-burning, we compared molecular biomarkers of exposure and preclinical effects in umbilical cord blood to neurodevelopmental outcomes from two successive birth cohorts enrolled before and after a highly polluting, coal-fired power plant in Tongliang County, China had ceased operation. Women and their newborns in the two successive cohorts were enrolled at the time of delivery. We measured PAH-DNA adducts, a biomarker of PAH-exposure and DNA damage, and brain-derived neurotrophic factor (BDNF), a protein involved in neuronal growth, in umbilical cord blood. At age two, children were tested using the Gesell Developmental Schedules (GDS). The two cohorts were compared with respect to levels of both biomarkers in cord blood as well as developmental quotient (DQ) scores across 5 domains. Lower levels of PAH-DNA adducts, higher concentrations of the mature BDNF protein (mBDNF) and higher DQ scores were seen in the 2005 cohort enrolled after closure of the power plant. In the two cohorts combined, PAH-DNA adducts were inversely associated with mBDNF as well as scores for motor (p = 0.05), adaptive (p = 0.022), and average (p = 0.014) DQ. BDNF levels were positively associated with motor (p = 0.018), social (p = 0.001), and average (p = 0.017) DQ scores. The findings indicate that the closure of a coal-burning plant resulted in the reduction of PAH-DNA adducts in newborns and increased mBDNF levels that in turn, were positively associated with neurocognitive development. They provide further evidence of the direct benefits to children's health as a result of the coal plant shut down, supporting clean energy and environmental policies in China and elsewhere.

Introduction

Polycyclic aromatic hydrocarbons (PAH) are a group of compounds formed from the incomplete combustion of coal and other organic material. PAH are ubiquitous in the environment, present in outdoor and indoor air from coal combustion, diesel and other motor vehicle emissions, tobacco smoking and cooking of food [1], [2]. In China, rapid economic expansion coupled with escalating fossil fuel-based energy production has resulted in massive amounts of air pollution. Specifically, coal-fired power plants produce over 70 percent of China's electricity, and new power plants being designed are also run on coal, which contributes to unsafe levels of toxic contaminants, including PAH [3], [4], [5].

Inhalation of airborne PAH leads to formation of DNA adducts, which are considered a valid marker of the biologically effective dose of PAH, reflecting individual variation in metabolism of PAH and DNA repair [6]. PAH are carcinogenic and neurotoxic. Previous studies by the Columbia Center for Children's Environmental Health (CCCEH) in New York City, Krakow and the present Tongliang cohort have shown that the developing fetus is more susceptible than the adult to PAH-DNA adduct formation. They have also shown that prenatal PAH exposure is associated with adverse impacts on child neurodevelopment, including developmental delay, reduced IQ, behavioral problems in childhood [7], [8], [9], [10], [11], [12] and reduced intelligence scores later in childhood [10], [12].

In Tongliang County, the high pollutant concentrations in ambient air prompted the government to shut down the local power plant in May of 2004 to significantly improve community health [4], [13]. This action, announced in advance, provided a unique opportunity to compare air monitoring, biomarker and health outcome data in two successive cohorts of children with and without prenatal exposure to emissions from the coal-fired power plant. In partnership with Chongqing Children's Hospital, the CCCEH carried out two prospective cohort studies between 2002 and 2005, pre-plant shutdown and post-plant shutdown, respectively. We previously reported evidence from the 2002 cohort that prenatal exposure to PAH adversely affected the neurodevelopment of children on the GDS at age 2 [11] and that the later cohort had significantly lower levels of PAH-DNA adducts in cord blood [9]. Additionally, the significant associations observed in the 2002 cohort between PAH-DNA adducts and adverse outcomes on the GDS at age 2 were no longer seen in the 2005 cohort [9].

Here we have analyzed levels of BDNF and their relationship to both adducts and developmental outcomes in the two cohorts. BDNF is critical for neurological survival and cognitive development of the central nervous system (CNS) [14], [15] and is the most widely distributed neurotrophin [16]. BDNF is first released as the proBDNF precursor before being cleaved into mBDNF [14]. Regulation of BDNF depends on its site of release: postsynaptic release is regulated by a Ca2+ influx through ionotropic glutamate receptors and voltage gated Ca2+ channels [17]; release from presynaptic sites is also dependent on mobilizing Ca2+ influx from intracellular stores [17], [18]. Following release, BDNF binds to two different transmembrane proteins: tropomyosin- related TrkB receptor and neurotrophin receptor p75, with higher affinity for the TrkB receptor [17]. Specifically, the binding of TrkB receptor triggers three signaling cascades that ultimately phosphorylate and activate the cAMP responsive element binding protein (CREB) transcription factor to encourage genes transcription essential for neuronal development and synaptic plasticity [17], [19].

Because BDNF and TrkB are expressed in the hippocampus to reinforce and stabilize synaptic connections, BDNF has been widely recognized as a key regulator in long-term potentiation (LTP), one of several phenomena underlying synaptic plasticity, learning and memory [16], [20]. LTP has been divided into early LTP and late LTP, where early LTP (characterized as short term potentiation) is shorter lasting and depends on the modification and translocation of proteins [21], and late LTP lasts several days and depends on de novo protein synthesis [20]. BDNF is an active mediator of neuronal processes in both the developing and mature brain, promoting differentiation, growth, and survival of neurons during development [22]. For late LTP, studies have shown that higher endogenous levels of BDNF are required to activate multiple signaling cascades that may act concertedly to regulate downstream cellular effects for memory formation and maintenance [23]. Therefore, in the mature brain, memory formation and development are associated with increased levels of BDNF and TrkB activation [24].

The goal of this paper is to determine whether mBDNF can provide a potential risk marker in assessing the neurodevelopment effects of exposure to PAH during fetal development, and in documenting the benefits of a regulatory intervention. Because previous studies have shown adverse effects of lower BDNF levels on learning and memory processes [24], we hypothesized that increased levels of PAH-DNA adducts would be associated with decreased levels of mBDNF that, in turn, would mediate the effects on neurodevelopmental outcomes determined by the DQ.

Methods

The overall approach used in this investigation and the comparison of the two cohorts with respect to adducts and DQ scores have been presented elsewhere [9].

Ethics Statement

This study was approved by the Institutional Review Board of Columbia University. All subjects gave informed written consent by completing a form approved by both the Columbia University IRB and by Chongqing Medical University.

Population

The city of Tongliang has a population of approximately 810,000 and is situated in a small basin approximately 3 km in diameter [3]. A coal-fired power plant located south of the town center operated during the dry season from 1 December to 31 May each year prior to 2004 in order to compensate for insufficient hydroelectric power during that time period. This plant was the principal source of local air pollution because in 1995 nearly all domestic heating and cooking units were converted to natural gas, motor vehicles were not a major pollution source, and there were no major coal-burning sources within 20 km of the city [4], [13]. In May 2004, the power plant was closed and replaced by the national grid system of electrical energy.

Study subjects

In the 2002 cohort, as previously reported [9], the subjects were 150 nonsmoking mothers and their newborns enrolled between 4 March 2002 and 19 June 2002 at four hospitals in Tongliang: the Tongliang County Hospital, the Traditional Chinese Medicine Hospital, the Tongliang Maternal Children's Health Hospital, and the Bachuan Hospital. In the 2005 cohort, the subjects were 158 children born at the same hospitals from 2 March 2005 to 23 May 2005 and recruited and studied using the same methods. The women were selected using a screening questionnaire when they checked in for delivery. Eligibility criteria included current nonsmoking status, ≥20 years of age, and residence within 2.5 km of the Tongliang power plant. All but one eligible woman agreed to enter the study. The demographic characteristics of the two cohorts are presented in Table 1.

Table 1

Demographic and Exposure Characteristics of the Cohorts.

		2002 Cohort	2005 Cohort
Mother's Age, years*		25.18±3.152, N = 110	27.909±4.585, N = 107
Mother's Education (%)	<High School	48 (43.64%)	59(55.14%)
	≥High School	62 (56.36%)	48(44.86%)
Gestational Age, days		277.35±11.27, N = 110	276.692±9.188, N = 107
Sex of Newborn, % female		56 (50.91%), N = 110	48(44.86%), N = 107
Birth Head Circumference, cm*		33.77±1.151, N = 110	34.130±1.283, N = 106
Infant's Birth Weight, g		3346.64±398.03, N = 110	3378.505±399.653, N = 107
Infant's Birth Length, cm		50.38±1.690, N = 110	50.212±1.635, N = 104
Heavy Metals	Cord Lead (mg/dL)	3.60±1.59	3.74±1.50
	Cord Mercury (ppb)	6.67±4.43	6.61±2.77
PAH-DNA Adducts, adducts/108 nucleotides*		0.324±0.139, N = 110	.204±.081, N = 107
ETS(hours/day)		0.293±0.586, N = 110	0.297±0.539, N = 107
Income**	Total N	77	107
	<10000	2	22
	10001 to 20000	13	52
	20001 to 30000	30	17
	30001 to 40000	22	6
	40001 to 50000	7	5
	50001 to 60000	3	3
	60001 to 70000	0	0
	>70000	0	2
BDNF, μg/dL*		752.87±463.71, N = 108	1266.57±619.77, N = 107
Gesell Scores	Average	99.42±10.74, N = 110	100.3±7.157, N = 107
	Motor	97.53±11.47, N = 110	97.83±7.821, N = 107
	Adaptive	98.71±14.90, N = 110	101.18±10.96, N = 107
	Language*	102.1±12.83, N = 110	100.47±9.777, N = 107
	Social	99.40±11.79, N = 110	101.83±6.808, N = 107

*significant difference between the two cohorts at alpha  = 0.05 level, using t-test.

** significant difference between the two cohorts at alpha  = 0.05 level, using chi-square test.

Personal interview

A 45-min questionnaire was administered by a trained interviewer after delivery. The questionnaire elicited demographic information, lifetime residential history (location of birth and duration of residence), history of active and passive smoking (including number of household members who smoke), occupational exposure, medication use, alcohol consumption during each trimester of pregnancy, and consumption of PAH-containing meat (frequency of eating fried, broiled, or barbecued meat during the last 2 weeks). Socioeconomic information related to income and education was also collected.

Biological sample collection and biomarker analysis

40–60 mL of umbilical cord blood was collected at delivery and 10 mL of maternal blood (10 mL) was collected within 1 day postpartum. Samples were transported to the field laboratory at the Tongliang County Hospital immediately after collection for processing. Blood samples, the buffy coat, packed red blood cells, and plasma were separated and stored at −70°C. B[a]P is widely used as a representative of PAH because concentrations of individual PAH in the urban setting are highly intercorrelated. Therefore B[a]P –DNA adducts were analyzed as a proxy for PAH-DNA adducts [25]. Details of laboratory methods for analyzing B[a]P-DNA adducts have been described[11]. Immunoassays for plasma levels of BDNF were performed using the BDNF Emax ImmunoAssay System (Promega) according to the manufacturer's instruction.

Measurement of child neurodevelopment using GDS

Two-year-old children in the cohort were administered the Chinese version of the standardized GDS for 0- to 3-year-old children adapted to the Chinese population by the Chinese Pediatric Society from the department of Pediatric Psychiatry in Xinghua Hospital in Shanghai, China. Each child is assigned a DQ in each of four areas: motor, adaptive, language, and social. The standardized mean (±SD) of the DQ is 100±15; a score <85 indicates developmental delay [26]. Testing was conducted by physicians in the same group who were certified in the GDS to maximize reliable assessment and valid interpretation, minimizing both inter-examiner and intra-examiner variability.

Statistical Analysis

We compared characteristics of the two cohorts using the t-test or chi-square test. Multiple linear regression and logistic regression were used to test the relationships among cord adducts, BDNF, and DQ. Covariates in the regression models were selected based on previous findings and included cord blood lead, cord mercury, environmental tobacco smoke (ETS), mother's education, mother's age, gestational age, and gender. Cord adducts were calculated on the original scale, while cord lead was dichotomized at the median of detectable lead in the multiple linear regression models. BDNF values were log transformed to normalize data and remove outliers. DQ scores were handled on the original scale. Cord blood BDNF levels of newborns in the two cohorts were compared using the Mann-Whitney test. We explored whether BDNF was a mediator of the relationship between PAH-DNA adducts and DQ scores by comparing the change in the beta coefficient of cord adducts in two regression models, with and without BDNF as a variable.

Results

The demographics, exposure, and developmental characteristics of each cohort are provided and compared in Table 1. There were no significant differences between the cohorts with respect to these characteristics except in mother's age and income level. The mean level of mBDNF was significantly higher in the 2005 cohort as compared to the 2002 cohort (1266.56 pg/ml vs. 752.87 pg/ml, P<0.05). The mean birth head circumference of the 2005 infants was significantly greater than that of the 2002 cohort (33.766 cm vs. 34.130 cm, P<0.05). The mean PAH-DNA cord adduct level was 0.324 adducts/108 nucleotides in the 2002 cohort and 0.204 adducts/108 nucleotides in the 2005 cohort, also a significant difference (P<0.05). BDNF levels were also significantly higher in the 2005 cohort, with a mean of 752.871 μg/dL in the 2002 cohort compared to 1266.568 μg/dL in the 2005 cohort (P<0.05). Although the differences were not statistically significant in all areas, DQ scores were consistently higher in the 2005 cohort than the 2002 cohort for motor, adaptive, social, and average DQ (Table 1).

In cord blood, PAH-DNA adducts were modestly but significantly inversely correlated with BDNF (r = −0.233, p<0.01, N = 269). The results of multiple regression analysis of the association between biomarkers and DQ scores are shown in Table 2 for levels of PAH-DNA adducts and BDNF protein separately. The regression coefficients (betas) represent the change in DQ score per unit increase in cord adducts or per log unit BDNF. Combining the two cohorts, a unit increase in cord adducts was negatively associated with average (β = −12.113, p = 0.014), motor (β = −10.699, p = 0.050, and adaptive (β = −16.472, p = 0.022) DQ scores after adjusting for cord blood lead, cord mercury, ETS, mother's education, mother's age, gestational age, and gender. A log unit increase in log transformed BDNF level was positively associated with average (β = 2.496, p = 0.017), motor area (β = 2.117, p = 0.018), and social area DQ (β = 3.222, p = 0.001) scores, after adjusting for covariates including cord PAH-DNA adducts, prenatal ETS, cord lead, cord mercury, gender, gestational age, mother's education, mother's age and income. BDNF was deemed to be a potential mediator variable as the betas displayed a change of over 10% after including BDNF as a covariate in the model with adducts as the independent variable and DQs as the dependent variable: the change in average DQ was 21.5%; in motor DQ 33.6%; in adaptive DQ 18.2%; and in social DQ, 37.3%.

Table 2

Regression Analysis for DQ and DNA-Adducts and BDNF.

	Average	Motor	Adaptive	Language	Social
Adducts 1	−12.113, p = .014	−10.699, p = .050	−16.472, p = .022	−11.679, p = .057	−9.544, p = .078
(N = 215)	(−21.786, −2.440)	(−21.409, .010)	(−30.545, −2.398)	(−23.720, .363)	(−20.180, 1.093)
BDNF 2	2.496, p = .017	2.117, p = .018	1.844, p = .086	.368, p = .518	3.222, p = .001
(N = 207)	(.454, 4.539)	(.467, 4.965)	(−.377, 5.604)	(−1.73, 3.416)	(1.694, 6.068)

Adjusting for cord lead (Ln), cord mercury (Ln), ETS (hours/day), mother's education, mother's age, gestational age and gender.

Adjusting for income, cord lead (Ln), cord mercury (Ln), ETS (hours/day), mother's education, mother's age, and gestational age.

Discussion

As hypothesized, comparison of the two cohorts in Tongliang, China, before and after closure of the local power plant, has provided additional evidence of significant measurable benefits on children's development. As previously reported, we found a significant reduction in mean PAH-DNA adduct levels (0.204 adducts/108 nucleotides in the 2005 cohort compared to 0.324 adducts/108 nucleotides in the 2002 cohort). Also, as reported previously, PAH-DNA adducts in cord blood were significantly associated with decrements in motor and language areas, along with overall DQ score in the 2002 cohort [11], consistent with previous evidence that prenatal exposure to B[a]P produces neurodevelopmental effects in the offspring of study animals [27] as well as in children at age 3 years in a New York City cohort [28].

Further, the 2005 cohort had a significantly higher average mBDNF level compared to the 2002 cohort (1266.56 pg/ml vs. 752.87 pg/ml, P<0.05). As PAH-DNA cord adducts were inversely correlated with BDNF levels (r = −0.233, p<0.01), we hypothesized that reduction in BDNF levels as a result of prenatal PAH exposure may have contributed to the adverse neurocognitive effects in the 2002 cohort. Overall, DQ scores were positively associated with BDNF levels. Although the associations were not statistically significant, motor, adaptive, social area, and the average scores were consistently higher in the 2005 cohort compared to the 2002 cohort. Language DQ scores were significantly higher in the 2002 cohort; however language is the least reliable measure for cognitive development because environmental enrichment, such as more parental interaction time, has been shown to reverse potential language deficits [29]. Moreover, our exploratory analyses suggest that BDNF may mediate in part the association between cord adducts and neurocognitive outcomes (average, motor, adaptive and social DQ scores).

The adverse developmental effects observed in the 2002 cohort are consistent with PAH being a maternal contaminant that crosses the placenta to affect neurodevelopment and cognitive function [30]. Furthermore, we have reported that prenatal PAH exposure adversely affects cognitive development [10], [12], [28]. This paradigm of prenatal PAH-induced neurotoxic effects is consistent with results from animal studies, which explained the relationship between downregulation of the gene encoding the receptor tyrosine kinase, MET, a target of B[a]P, which mediates neuronal process growth and synapse formation [31]. In another study, in utero exposure to aerosolized B[a]P altered the expression profile of genes involved in synaptogenesis such as Sp4 transcription factor, as well as glutamate homeostasis and metabolism of reactive oxygen intermediates [32].

These results are consistent with findings of Toledo-Rodriguez et al. that prenatal exposure to ETS (which contains PAH) was associated with increased methylation of the BDNF exon and down-regulation of BDNF[33], and reports of dysfunction in the BDNF system and compromised cognitive function in adult life from fetal exposure to unfavorable intrauterine conditions [34]. This paper is the first to assess BDNF and cognitive development with respect to prenatal exposure to PAH. The mechanisms by which PAH-related reduction of BDNF may cause child developmental problems have not been determined but recent studies have pointed to several potential mechanistic pathways: alteration in the conversion of proBDNF to mBDNF which influences levels of signaling proteins involved in LTP; disruption of activation of TrkB; and induction of methylation changes in BDNF.

First, urban air pollution has been shown to increase levels of plasminogen activator (PA) inhibitor [35], a serine protease necessary to cleave plasminogen into its active form of plasmin [36], [37]. Tissue PA (tPA) is involved in cleaving proBDNF to mBDNF. Defects in tPA-dependent cleavage of proBDNF to mBDNF may downregulate mBDNF, resulting in impaired LTP [36], [37]. Second, the effects of BDNF are believed to be mediated by TrkB, which follows receptor tyrosine kinase activation to initiate three major signaling pathways: phospholipase Cγ (PLCγ), phosphatidylinositol 3-kinase (PI3K), and cascade governed by extracellular signal-regulated kinases (ERK) [38]. Experimentally, prenatal stress decreases hippocampal mBDNF expression and subsequent TrkB- regulated signaling cascades for LTP [39]. The absence of TrkB forces BDNF to selectively bind to p75, resulting in death of the hippocampal neurons [40]. Similarly, higher levels of endogenous proBDNF and low levels of mBDNF were observed in prenatally stressed rats, which influenced synaptic plasticity later in life [37]. Finally, BDNF expression is also regulated by the dynamic methyl CpG binding protein 2 (MeCP2), which may function as a molecular linker between DNA methylation, chromatin remodeling, and ultimately gene repression by recruiting and interacting with histone deacetylases (HDACs) and transcription repressor mSin3A [41], [42]. Experimentally MeCP2 has been shown to bind to methylated CpG sequences in the BDNF promoter III [43] and in BDNF exon IV promoter to induce chromatin compaction and repress gene expression [41]. Prenatal exposure to maternal smoking was associated with enhanced DNA methylation in BDNF exon IV that may lead to long-term deregulation of BDNF [33]. Thus, mBDNF expression is regulated by changes in DNA methylation patterns that have previously been associated with levels of ambient air pollution and prenatal tobacco smoke exposure.

Strengths of the present study are the prospective cohort design, the ability to control for confounders, and the measurement of PAH-DNA adducts as a biological dosimeter for environmental exposure. Additionally, mBDNF is a mechanistically relevant measure of preclinical response/potential risk. A limitation is that we were not able to adjust for postnatal exposure to PAH. However, in this population the major source of PAH was the coal-burning plant emissions prior to its closure.

This study provides insight into the relationship between PAH, PAH-DNA adducts, BDNF, and developmental outcomes. Further studies are needed to determine whether specific enrichment environments can mitigate the damage caused by prenatal environmental exposure to PAH. In conclusion, this provides further evidence that the power plant shutdown in Tongliang, China in May of 2004 directly benefited the health and development of children.