18 Aug 2009:
Correction: Relation of DNA Methylation of 5′-CpG Island of
In a longitudinal cohort of ∼700 children in New York City, the prevalence of asthma (>25%) is among the highest in the US. This high risk may in part be caused by transplacental exposure to traffic-related polycyclic aromatic hydrocarbons (PAHs) but biomarkers informative of PAH-asthma relationships is lacking. We here hypothesized that epigenetic marks associated with transplacental PAH exposure and/or childhood asthma risk could be identified in fetal tissues. Mothers completed personal prenatal air monitoring for PAH exposure determination. Methylation sensitive restriction fingerprinting was used to analyze umbilical cord white blood cell (UCWBC) DNA of 20 cohort children. Over 30 DNA sequences were identified whose methylation status was dependent on the level of maternal PAH exposure. Six sequences were found to be homologous to known genes having one or more 5′-CpG island(s) (5′-CGI). Of these, acyl-CoA synthetase long-chain family member 3 (ACSL3) exhibited the highest concordance between the extent of methylation of its 5′-CGI in UCWBCs and the level of gene expression in matched fetal placental tissues in the initial 20 cohort children. ACSL3 was therefore chosen for further investigation in a larger sample of 56 cohort children. Methylation of the ACSL3 5′-CGI was found to be significantly associated with maternal airborne PAH exposure exceeding 2.41 ng/m3 (OR = 13.8; p<0.001; sensitivity = 75%; specificity = 82%) and with a parental report of asthma symptoms in children prior to age 5 (OR = 3.9; p<0.05). Thus, if validated, methylated ACSL3 5′CGI in UCWBC DNA may be a surrogate endpoint for transplacental PAH exposure and/or a potential biomarker for environmentally-related asthma. This exploratory report provides a new blueprint for the discovery of epigenetic biomarkers relevant to other exposure assessments and/or investigations of exposure-disease relationships in birth cohorts. The results support the emerging theory of early origins of later life disease development.
Citation: Perera F, Tang W-y, Herbstman J, Tang D, Levin L, Miller R, et al. (2009) Relation of DNA Methylation of 5′-CpG Island of ACSL3 to Transplacental Exposure to Airborne Polycyclic Aromatic Hydrocarbons and Childhood Asthma. PLoS ONE 4(2): e4488. doi:10.1371/journal.pone.0004488
Editor: Catherine M. Suter, Victor Chang Cardiac Research Institute, Australia
Received: September 8, 2008; Accepted: December 19, 2008; Published: February 16, 2009
Copyright: © 2009 Perera et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work is in part supported by National Institutes of Health grants: ES009600 (FP), ES008977 (FP), ES015905 (FP, SMH, RM, DT), ES013071 (SMH), ES015584 (SMH), and ES006096 (SMH), ES016817 (WYT), the Department of Defense award DAMD W81XWH-06-1-0373 (WYT), the U.S. Environmental Protection Agency award RD-832141(FP), and the New York Community Trust (FP,RM, DT, JH). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Asthma is the most common chronic childhood disease  and its risk may be strongly influenced by prenatal events . Low income communities experience some of the highest childhood asthma rates in the U.S  and new preventive strategies are lacking due in part to the absence of predictive biomarkers. In our longitudinal cohort study of children residing in urban low-income, minority communities of New York City – the asthma rate exceeds 25% and is among the highest in the nation. Preliminary evidence suggests that transplacental exposure to polycyclic aromatic hydrocarbons (PAHs) derived largely from traffic-related air pollutants may be a risk factor for the early development of asthma-related symptoms in this cohort . However, we could not rule out the possibility that postnatal exposure of children in this cohort to these pollutants may also contribute to the symptoms without additional studies.
According to the hypothesis of the early origins of chronic disease, exposure to agents that are transmissible from mother to fetus has profound influences on the functional capacity and performance of many organs/tissues later in life –. Developmental plasticity allows the fetus to make anticipatory responses to the external environment by altering the course of cellular and organ differentiation in utero in order to gain adaptive advantages for later life challenges. However, a pronounced mismatch between “anticipatory adaptations” made during early life and demands in later life could be a cause of disease . More importantly, environmental insults could “mislead” early organogenesis resulting in serious ailments in later life. The best studied mechanisms underpinning the latter scenario is the induction of aberrant DNA methylation of regulatory sequences in susceptible genes by environmental toxicants during critical developmental periods, leading to inappropriate gene expression and disease pathogenesis in later life –. This body of research suggests that transplacental exposure to high levels of airborne PAHs could cause aberrant DNA methylation changes, leading to dysregulation of gene expression and perhaps childhood asthma. These alterations could occur, not only in target organs/tissues such as fetal lung or the immune system , but also in peripheral tissues such as fetal circulating white blood cells and placental tissues, thus providing opportunities to discover and develop highly sensitive and specific, minimally invasive biomarkers for exposure assessments or clinical indications .
The primary objective of this study was to explore whether transplacental exposure to PAHs in humans induces epigenetic reprogramming involving aberrant DNA methylation of specific genes that might be mechanistically related to childhood asthma or airway inflammation. Study subjects were children born to nonsmoking Dominican and African American women living in Northern Manhattan and the South Bronx. They were participants in the Columbia Center for Children's Environmental Health (CCCEH) cohort study previously described –. Mothers completed prenatal air monitoring and donated a maternal and/or umbilical cord blood sample and, in some cases, fetal placental tissue (FPT) at delivery. Umbilical cord white blood cells (UCWBC) and matched FPT from 20 cohort subjects were used initially for candidate identification and selection, as it is not feasible to obtain lung/immune cells in newborns. UCWBC provide a reasonable surrogate for our target organ/tissue (lung/immune cells) because they contain stem cells that can populate the lungs in later life and also provide a rich source of T cells  which are important producers of cytokines and other asthma mediators. The matched FPT provides ample tissue for obtaining fetal RNA to verify correlative changes in gene expression. A larger sample of 56 cohort subjects was used to validate the markers identified. Methylation of a 5′-CpG island (CGI) in acyl-CoA synthetase long-chain family member 3 (ACSL3) was found to be positively and significantly associated with the level of maternal PAH exposure and with a parental report of asthma symptoms prior to age 5 in this exploratory study.
Identification of candidate genes and selecting the most relevant candidate
Figure 1 shows the work flow used in the discovery of candidates and strategies employed for selecting candidate genes associated with prenatal airborne PAH exposure and with asthma.
Initially, ten umbilical cord white blood cell (UCWBC) samples from children with maternal (prenatal) PAH exposure above and ten with prenatal exposure below the cohort median were selected irrespective of gender and ethnicity. DNA was isolated from the sample was dichotomized according to the median of prenatal PAH exposure obtained from 48-hour personal air monitoring of the mothers during pregnancy (the cohort median is 2.3 ng/m3). Candidates with differentially methylation status between the samples from the two groups were discovered by Methylation Sensitive Restriction Fingerprinting (MSRF) following by subcloning and sequencing. Candidates were further identified by in silico analysis from databases from NCBI and UCSF genomic centre. Candidates which contained CpG island (CGI; CG rich region with >60% GC content and obs/exp ratio >0.6) on its 5′ flanking region were chosen for bisulfite sequencing to confirm the dependency of the methylation status of the CGI in a candidate gene on maternal PAH exposure. Six genes confirmed to have the latter relation were subjected to gene expression analysis in the 14 matched fetal placental tissues (FPT). The ACSL3 was found to have the highest concordance between degree of CGI methylation and level of gene expression in reverse manner. It was therefore selected for the final analyses for correlation between methylation status of its CGI and transplacental PAH exposure and/or parent reported asthma symptoms of childhood asthma up to age 5 in a case-control study comprised of 56 participants. For this part of the study, Methylation specific-PCR (MSPCR) protocol was optimized to analyze the methylation status of R1 in Figure 4A was used as the high throughput method to analyze this larger sample set.
Twenty participants were randomly selected without knowledge of asthma classification: ten with maternal PAH exposure levels below the CCCEH cohort median of 2.3 ng/m3 (“low PAH group”) and ten with prenatal PAH levels above the cohort median (“high PAH group”). Methylation-sensitive restriction fingerprinting ,  was performed to identify candidate sequences displaying differential methylation status between high and low PAH groups. Both hypo- and hyper-methylated candidate sequences in the high PAH group when compared to the low PAH group were selected. In silico analyses revealed that 19 out of 31 sequenced candidates were homologous to known genes (Table 1) and 6 of these aligned to a CGI(s) located in the 5′ flanking region (proximal promoter and/or exon/intron 1) of a known gene (5′-CGI) (Figure 2). These 6 genes therefore were chosen as candidates for further investigation. Unexpectedly, all six candidates were found be hypermethylated sequences in the high PAH group when compared to the low PAH group. A literature search showed that they are expressed in the lung and/or lymphoid tissues and have functions related to inflammatory/or other immune responses (dual specificity phosphatase 22; DUSP22, accession # BC022847, ), DNA damage and repair (RAD21 homolog (S. pombe); RAD21, accession # BC050381, ), fatty acid metabolism (ACSL3, accession # BC041692, ; stearoyl-CoA desaturase 5; SCD5, accession # NM_001037582,  ), gene transcription (Scm-like with four mbt domains 2; SFMBT2, accession # NM_001029880, ; WW domain containing oxidoreductase; WWOX, accession # AF211943,  ). Bisulfite sequencing analyses of the initial 20 UCWBC samples confirmed a good agreement between hypermethylation of all of the candidates' 5′-CGI(s) and either high or low PAH exposure status, consistent with MSRF findings (Table 1).
The CGI(s) (shaded in blue in the genomic DNA sequence) were identified in silico based on a CG content >60%, with an observed/expected ratio of 0.6 according to instructions from MethPrimer. TSS: Transcription start site. The nested PCR-amplified region was indicated by a line with double arrow heads and methylation status of this region was determined by bisulfite genomic sequencing.
We next determined the strength of the relationship between the percent of methylation of a candidate gene's 5′-CGI in UCWBC DNA samples and its transcript expression level in matched FPT RNA samples (Table 2 and 3). The degree of methylation of a region encompassing the entire predicted 5′-CGI for each candidate gene was determined by bisulfite sequencing (Figure 2, the location of the amplicon was marked as BSPCR) and mRNA levels determined by real time RT-PCR. The percentage of promoter methylation (Met%) determined by bisulfite sequencing and relative gene expression levels (RER) determined by real time RT-PCR of all six candidates tested on 14 paired UCWBC DNA and FPT RNA samples are presented in Table 2. The degree of concordance between these two parameters for each candidate is shown in Table 3. We found that the degree of ACSL3 5′-CGI methylation and the level ACSL3 transcript expression exhibited the expected inverse relationship with the highest concordance and significance (τ = −0.45; p<0.01) when compared to the other candidates. These findings indicate that the ACSL3 5′-CGI is likely a regulatory CGI since increased cytosine methylation in this sequence in UCWBC DNA samples was found to be significantly associated with transcriptional silencing of the gene in matched FPTs.
To further demonstrate that the methylation status of the ACSL3 promoter affects gene expression, we used a lung non-small carcinoma cell line H1299 as a model system. Two series of experiments were performed. In the first series, we treated H1299 cells with a DNA demethylation agent, 5-aza-2′-deoxycytidine (5-aza-dC), and found that ACSL3 expression was significantly increased (Figure 3A) along with promoter demethylation (data not shown). These data provide evidence that demethylation of the ACSL3 promoter is associated with enhancement of gene expression. In the second series of experiments, H1299 cells exposed to low doses of benzo[a]pyrene (BaP; 0.01 to 1 nM), a representative PAH, showed a dose-dependent decrease in ACSL3 transcript expression as indicated by real time RT-PCR analysis (Figure 3B) and increased promoter hypermethylation as revealed by bisulfite sequencing (Figure 3C). These data strongly support the notion that the CpG island in the ACSL3 promoter plays a regulatory role in gene expression.
A: Cells were treated with 0.5 and 1.0 µM 5-AZA-dC with DMSO as control every 2 days for a total of 8 days. B: Cells were treated with 0.01, 0.1 and 1.0 nM BaP with DMSO as control every 2 days for a total of 4 days. RNA was isolated, reverse transcribed and underwent real-time PCR. The 2-ΔΔCt method was used to calculate the relative expression level of transcripts normalized to β-actin. *Statistically significant differences between exposed and control was accepted at p<0.05. C: Diagram represents methylation status of ACSL3 promoter of H1299 cells exposed to BaP assayed by bisulfite genomic sequencing. Cells were treated with 0.01, 0.1 and 1.0 nM BaP with DMSO as control every 2 days for a total of 4 days. Open circle: unmethylated CG; closed circle: methylated CG. Putative transcription factor binding sites such as Sp1, AP2, GCF, c-fos and junB are shown in scale with the particular CG sites on the promoter.
Association of methylation status of ACSL3 5′CGI with PAH exposure in the larger sample
ACSL3 was chosen for subsequent assessment of the relationships between prenatal PAH exposure and gene 5′-CGI methylation status and for exploratory analyses of PAH exposure status and asthma classification in a larger sample of 56 children from the CCCEH cohort.
The organization of the 5′-flanking sequence in the promoter region of ACSL3 is shown in Figure 4A. Before we conducted the full investigation on the larger sample of 56 children, the following quality control experiments were performed. First, bisulfite sequencing analyses of an amplicon encompassing the entire predicted ACSL3 5′-CGI (Figure 4A, BSPCR marks the amplicon) were performed on the 20 UWBC DNA samples; five from each of the high or low PAH groups with or without asthma. Bisulfite sequencing data (Figure 4B) showed that hypermethylation of this CGI was found to be strongly associated with maternal exposure to a PAH level higher than the cohort median of 2.3 ng/m3. Asthma status did not seem to affect methylation status of ACSL3 promoter. Second, based on bisulfite sequencing data, we designed and optimized three methylation-sensitive-PCR (MSPCR) protocols ,  to amplify three regions (R1, R2 and R3) within the CGI that exhibited the greatest differential methylation changes between samples from the “high PAH” and the “low PAH” groups regardless of asthma status (Figure 4A, R1–3 are marked). These three regions were closely flanking the transcriptional start site (TSS). We further validated that the methylation status of R1 exhibited the strongest correlation with bisulfite sequencing data from analyses of the entire CGI. Representative data from MSPCR analyses of samples from children with prenatal PAH exposure above (high PAH) or below (low PAH) the cohort median are illustrated in Figure 4C. Thus, MSPCR of R1 was chosen as the high throughput method for investigating the methylation status of the ACSL3 5′-CGI in the larger sample comprised of 56 children.
A) Schematic diagram of the CG content (%) in the 5′ flanking region of the ACSL3 gene. A CGI of 1058 bp was located at the 5′ end of ACSL3 including transcription start site (TSS), untranslated region (UTR) exon1 and intron 1. Individual CG sites are marked as red vertical lines in the genomic DNA sequence. The PCR-amplified regions were indicated by lines and methylation status of this region were determined by BSPCR and MSPCR. B) The diagram represents the methylation status of MSPCR-amplified regions (R1–R3) of 20 samples with high or low PAH exposure with or without asthma assayed by bisulfite genomic sequencing. 6 clones were sequenced from each sample. The PAH level of each sample is shown. Open circle: unmethylated CG; closed circle: methylated CG. MSPCR-primers were then specially designed on particular regions (R1–R3). Overall Met% of ACSL3 promoter is shown. Statistical difference was accepted at p<0.05* when compared with the group of low PAH without asthma. C) Representative results from MSPCR analyses on Region 1 (R1) of samples from the low and high PAH group. M: methylated; and U: unmethylated. The methylation status of ACSL3 5′-CGI of all samples was further analyzed by MS-PCR using the same sets of primers. Results are shown in Table 4 and Table 5. High PAH (≥ the cohort median of 2.3 ng/m3), Low PAH (<the cohort median).
Receiver Operating Curve (ROC) analysis of data from 56 subjects identified the maternal PAH exposure cutpoint value that best distinguished methylated versus unmethylated ACSL3 5′CGI status determined by MSPCR to be 2.41 ng/m3 (Figure 5). The Area Under ROC curve (AUC) value was 0.82 suggesting excellent separability of the subjects with methylated or unmethylated ACSL3 5′CGI. Using the “optimal” cutpoint value, methylation status was predicted with 75% sensitivity, 82% specificity.
Point labels are values of PAH of each children. Statistical software, sas, was used for ROC plot and R was used for determination of cutoff value. The cutoff value corresponds to the desired sensitivity and specificity (or 1-specificity).
Table 4 shows that 30 of the 56 children had a methylated ACSL3 5′-CGI; 81% of these 30 children were born to mothers with PAH exposure equal or exceeding 2.41 ng/m3, the cutpoint value determined by the ROC analysis. Thus, methylation was positively and significantly associated with maternal PAH exposure level. The odds ratio (OR) of having a methylated ACSL3 5′CGI in UCWBC DNA was 13.8, 95% CI = [3.8, 50.2] (p<0.001), given the maternal PAH exposure was ≥2.41 ng/m3, compared to those exposed to <2.41 ng/m3. ACSL3 5′CGI methylation status was not affected by gender or ethnicity.
Relationship between ACSL3 5′ CGI methylation and asthma classification
Table 5 shows the percentage of subjects with methylated or unmethylated ACSL3 5′ CGI according to asthma classification. Methylation of the ACSL3 5′CGI was positively and significantly associated with asthma classification. The odds ratio (OR) relating methylation status and asthma classification among all 56 subjects was 3.9, 95% CI = [1.1, 14.3] (p = 0.03). Eleven of 15 children classified as asthmatic (73%) had UCWBC DNA with a methylated ACSL3 5′-CGI. Table 5 further shows that the median (minimum, maximum) values of maternal PAH exposure differed significantly by methylation status of the ACSL3 5′-CGI: methylation = 3.39 ng/m3 (1.11, 34.5) and unmethylation = 1.70 ng/m3 (0.49, 3.33), respectively (p<0.001). All maternal PAH exposure values except the maximum (34.5 ng/m3) were below 7.7 ng/m3. When the latter analysis was repeated excluding the extreme high PAH value, the p-value was unchanged.
In this exploratory study, we have identified ACSL3 as a candidate biomarker/gene whose 5′-CGI methylation status appears to be related to transplacental PAH exposure and further associated with PAH-associated asthma. Thus, ACSL3 may be the first potential surrogate endpoint for environmentally related childhood asthma. If confirmed, such a biomarker holds promise in assessing PAH exposure and as a clinically relevant predictor for asthma risk in children born to mothers exposed to air pollutants such as traffic-related combustion emissions.
We used a set of stringent criteria for the identification, selection and validation of candidates, which resulted in choosing ACSL3 the final candidate for the 56-sample study. We first used an unbiased method, MSRF, to identify both hypo- and hyper-methylated sequences in the high PAH group compared to the low PAH group. In silico analyses then revealed that 19 out of 31 sequences identified were homologous to known genes and 6 of these aligned to a CGI(s) located in the 5′ flanking region of a known gene (5′-CGI). These 6 genes therefore were chosen as putative candidates for further investigation. By comparing the degree of methylation of their 5′CGI in UCWBC DNA samples and levels of gene expression in FPT RNA samples we identified ACSL3 to have the highest inverse concordance in vivo. We then employed a lung cancer cell line (H1299) to demonstrate that the demethylation of ACSL3 5′-CGI by 5-aza-dC treatment induced enhancement of gene expression in H1299. Conversely, the treatment of H1299 with low doses of BaP, a PAH, was associated with dose-dependent hypermethylaiton of ACSL3 5-CGI and gene silencing. Collectively, these data strongly support methylation of the 5′-CGI of ACSL3 in gene regulation. They also provide evidence that PAH could directly induce ACSL3 promoter methylation. This well though out, step-wise approach allowed us to identify a candidate to be tested in the 56-sample study to determine whether methylation of the ACSL3 5′CGI could be a stand alone biomarker for PAH exposure and/or PAH-associated childhood asthma.
Measuring human exposure to chemicals such as PAHs and determining their relationship to disease outcomes is a challenge for environmental health, especially when children are the study subjects. Although some studies have found a positive correlation between PAH-induced DNA adducts and mutagenesis or cancer development  others have failed to establish clear cause-effect relations between PAH adducts and the induction of these or other disease outcomes , suggesting that other mechanisms may be involved. Moreover, some of the methods for measuring genotoxic mechanisms such as adducts are complex and low-throughput. Therefore, it is of significance that, in this investigative study, we have developed a relatively simple, high throughput assay based on methylation of a sequence in the promoter of the ACSL3 gene that appears to distinguish children born to mothers with high airborne PAH exposure from those who were not (OR = 13.8; p<0.001) with high sensitivity (75%) and specificity (85%). The Area Under ROC curve (AUC) value for this assay was 0.82 suggesting excellent separability of the subjects with the methylated or unmethylated ACSL3sequence. This optimized MSPCR protocol allows for a rapid determination of ACSL3 5′-CGI methylation status in large sample sets, which is a normal requirement in population studies. It also has the added advantage to be developed into a minimally invasive test since only a small amount of DNA from UCWBC is required for the assay. If validated, such an assay is expected to have a high utility in the assessment of PAH exposure in future population-based studies.
Multiple epidemiological studies have suggested that asthma risk is determined by prenatal and/or early life exposure to pollutants and allergens that modify later life lung function and T helper cell allergic phenotype . Epigenetic modifications of regulatory genes crucial to the development of asthma-related pathophysiology have been put forth as a possible mechanistic basis of this disease . Thus, asthma may be caused by epigenetic changes such as DNA methylation, post-translational modifications of histones, dysregulation of non-coding RNAs, and chromatin alterations. At present, the best studied mechanism is DNA methylation-mediated disruption of gene regulation. In this regard, experimental data indicating that DNA methylation of genes critical to the differentiation of T helper cells could skew their response towards or away from an allergic phenotype lend support to this theory –. Additionally, we recently demonstrated that inhaled diesel exposure (a source of PAHs) and intranasal Aspergillus fumigatus induced specific DNA methylation changes in the promoters of interferon-γ and interleukin-4 in T helper cells, respectively, and altered circulating levels of IgE in mice . Despite progress made in experimental studies, fundamental questions related to asthma risk among human populations remain; these can only be addressed by cohort studies using specimens linked to well-annotated clinical and environmental exposure data. Such cohort-driven epigenetic research has the potential to address key questions, such as those concerning the influence of early life factors, environmental exposure, other disease states, and lifestyle factors on susceptibility to asthma. In this proof-of-principle study, we developed an innovative approach to population-based studies of environmental epigenetics, which could easily be implemented in future large-scale investigations.
A range of approaches is available to obtain quantitative and qualitative information on genomic DNA methylation changes . We chose methylation sensitive restriction fingerprinting ( MSRF; ,  ) as a tool for the discovery of PAH-reprogrammable genes in UCWBC due to its high sensitivity, requirement for very small amounts of starting materials (∼100 ng), and the ease in visual identification of differentially methylated sequences among a moderately large number of samples. Using only four arbitrary primers and six permutations of MSRF analyses on UCWBC samples from the high and low PAH exposure groups, we discovered 31 differentially methylated sequences. In future studies more arbitrary primers could be used to cover a larger portion of the genome in order to uncover more candidates . In this study we had overcome another technical hurdle related to the quantification of gene expression in UCWBC. This task is challenging because the amount of cellular components in umbilical cord blood is usually very limited and prohibitive for high quality RNA isolation. To circumvent this problem, in a subset of samples, we assessed the concordance between the extents of methylation in the promoter-CGI(s) of each gene (using UCWBC DNA samples) with its respective transcript expression levels in RNA samples isolated from matched FPT samples. This choice was based on the logic that both UCWBC and FPT are fetal tissues are subjected to similar environmental exposure through the mother. Our results demonstrated good concordance for at least the ACSL3 gene. In future design of such an epigenetic epidemiology experiment, preference should be given to improved techniques for the isolation of RNA and DNA from the UCWBC.
ACSL3 belongs to the acyl-CoA synthetase long chain (ACSL) family of genes which encodes key enzymes in fatty acid metabolism . ACSL3 is expressed in lung and thymic tissue , . Intracellular conversion of long chain fatty acids to acyl-CoA is required for energy production. ACSL catalyzes the initial step in this conversion yielding acyl-CoAs, which are used both in the synthesis of cellular lipids and in degradation via beta-oxidation for energy production , . Moreover, acyl-CoA synthetase can catalyze the ATP-dependent production of arachidonoyl-CoA that is a substrate for lysophosphatidyl acyltransferase involved in arachidonic acid-containing phospholipid remodeling in proliferating T cells . Lipid metabolism influences membrane proteins, including ion channels. In the lung, ACSL3 is found most abundantly in microsomal and cell membrane fractions, with mitochondria as the next richest source . Inhibition of ACS enzyme activity has been shown to cause cell death in p53-defective lung cancer cells . Thus hypermethylation of this gene in T helper cells or lung tissues is expected to diminish fatty acid utilization and beta-oxidation-energy production, and possibly influence membrane phospholipid composition. Whether these functional changes directly affect asthma is unknown. However, several epidemiological studies show that fatty acid composition in milk , , diet , , and umbilical cord blood  affects the development of allergy and other inflammatory diseases including asthma. For example, dietary omega-3 (n-3) fatty acids having a variety of anti-inflammatory and immune-modulating effects may affect allergy and asthma. Further mechanistic studies in in vitro and in vivo experimental models are needed to decipher the role of ACSL3 in childhood asthma. Interestingly, ACSL3 is located in 2q36.1 which has recently been shown to be associated with regions of the asthma susceptibility loci, in specific populations –
Finally, the current finding of a putative epigenetic marker that is associated with PAH exposure and asthma adds to other evidence from the CCCEH cohort that PAHs increase risk of respiratory symptoms and probable asthma. At age 2 years, significantly more difficulty breathing and more probable asthma were reported among children jointly exposed prenatally to PAHs and postnatally to environmental tobacco smoke (ETS) (OR for probable asthma = 7.52, 95% CI 1.71–33.11, p<0.01) . These relationships were reexamined through age 5 years and similar trends were detected, with an observed interaction between high prenatal PAH exposure and postnatal ETS on wheeze and probable asthma at age 5 (OR = 5.48, 95% CI 1.17–25.58, p<0.01) (unpublished data). A parallel cohort study in Krakow, Poland has found a significant relationship between prenatal PAH exposure and respiratory symptoms . In a subset of 333 subjects, PAH exposure measured during pregnancy by personal air monitors was associated with an increased risk for wheezing without cold (RR 3.8; 95% CI: 1.2–12.4), during the course of the infant's first year of life .
In conclusion, no previous studies have examined the effects of prenatal exposure to ambient air pollutants on DNA methylation patterns in genes potentially associated with the asthma phenotype of the offspring. This is an important gap in our understanding of asthma pathogenesis. Given the potential importance of epigenetic mechanisms in the etiology of childhood asthma, research is needed to identify an epigenetic profile related to PAH-associated childhood asthma and determine whether reprogramming events associated with transplacental exposure to PAHs increase risk of childhood asthma. The current report of ACSL3 as a candidate biomarker of PAH exposure and a putative predictor of PAH-associated childhood asthma represents an important first step in this area of research. Finally, the novel design of this study could serve as a prototype approach for future discoveries of additional exposure biomarkers or surrogate disease endpoints in other birth cohort studies.
Materials and Methods
Study subjects were nonsmoking Dominican and African American women and their children residing in Washington Heights, Central Harlem and the South Bronx participating in the CCCEH cohort study previously described –. Written informed consent was obtained from all subjects following procedures approved by the Institutional Review Board of the New-York Presbyterian Medical Center. Full enrollment into the study required completion of prenatal air monitoring and collection of a maternal and/or umbilical cord blood sample at delivery. As of November 1, 2007, 74% of the more than 600 fully enrolled subjects had been retained in the cohort through at least five years. During pregnancy and following delivery, a trained bilingual research team administered repeat questionnaires eliciting environmental and respiratory health histories . There were no differences in demographics of subjects fully enrolled and those lost to follow-up –. Based on parental report of a doctor's diagnosis of asthma or probable asthma before the age of five, children were classified as “asthmatic” or “nonasthmatic.” (“with or without a parental report of asthma”).
Study Sample Selection
All subjects with stored UCWBC and paired FPT were eligible for inclusion in this proof-of-principle study. The population was dichotomized on the median PAH value of the entire cohort (2.3 ng/m3). Half of the 56 subjects in the present study were selected from among the ones above this median (the high PAH exposure group) and half from among those below the cohort median (the lower PAH exposure group). Twenty subjects (10 high and 10 low exposure) were used to identify UCWBC DNA sequences whose methylation status was dependent on maternal PAH exposure (high versus low) and to validate concordant transcript changes in matched FPT samples. All 56 samples were then used to validate the association between PAH exposure and ACSL3 5′CGI methylation and to explore the association between ACSL3 5′CGI methylation and childhood asthma status. Table 6 shows the characteristics of the cohort population and study sample, respectively. The sample was representative of the cohort population in terms of maternal age, child's sex, median PAH exposure, and percentage of children with probable asthma; however, the sample differed with respect to ethnicity. Table 7 showed the PAH distribution and variation in PAH levels in selected and unselected cohorts.
Monitoring and Sample Collection
Personal prenatal air monitoring was conducted for two consecutive days on all participants during the 3rd trimester of pregnancy as previously described –. Levels of pyrene and eight carcinogenic PAHs were analyzed (4). Total PAH levels for each participating mother were computed as the sum of the eight carcinogenic PAHs, where the values below the limit of detection were imputed as 0.125 ng/m3 (the median value between 0 and the limit of detection). Umbilical cord blood (30–60 mL) was collected at delivery. The buffy coat, packed red blood cells, and plasma samples were separated and stored at −80°C. DNA (100 ng–500 ng) was extracted from UCWBCs. For some UCWBC samples, matched FPTs were also available and total RNA was extracted using a TRIZOL® (Invitrogen, CA) according to manufacturer's protocol. All samples were coded.
Methylation Sensitive Restriction Fingerprinting (MSRF)
MSRF was performed in a single-blinded manner on asthma classification, as previously described , . Five hundred ng of UCWBC DNA was restriction-digested using MseI [New England Biolabs, MA] alone or in combination with BstUI [New England Biolabs, MA] and amplified by PCR with various combinations of paired arbitrary primers (Table 8). Following electrophoresis, candidate sequences displaying differential methylation status between the “high PAH group” and “low PAH group” were isolated, extracted from the gel, cloned, identified and aligned with various databases. Sequences matched to known genes and with 5′-CGIs were selected for bisulfite genomic sequencing and confirmation studies of their regulatory function in gene expression. Six candidates [ACSL3, RAD21, DUSP22, SFMBT2, SCD5 and WWOX] were selected.
Treatment of H1299 with 5-aza-dC and BaP
H1299 cell line is a gift from Dr. Carolyn M. Klinge (University of Louisville) and were maintained in RPMI-1640 medium (Invitrogen, CA) supplemented with 10% fetal bovine serum (Hyclone, UT), 10 mM HEPES (Invitrogen, CA), 1 mM sodium pyruvate (Invitrogen, CA) and 4500 mg/L glucose (Invitrogen, CA). For demethylation assays, cells were treated with 0.5 and 1.0 µM 5-aza-dC (Tocris, MO) with DMSO (Sigma, MO) as control every 2 days for a total of 8 days. For BaP treatment, cells were exposed to 0.01, 0.1 and 1.0 nM BaP (Sigma, MO) with DMSO as control every 2 days for a total of 4 days. Medium was replenished with new drug of 5-AZA-dC or BaP every 48 hr. DMSO was added in the control samples. RNA was isolated, reverse transcribed and undergone real-time PCR. DNA was isolated, bisulfite-modified and undergone bisulfite genomic sequencing.
Bisulfite genomic sequencing
Genomic DNA (100 ng) from UCWBC samples or H1299 cells was modified with sodium bisulfite by using the EZ Methylation Modification Kit (Zymo Research, CA) before conducting PCRs. In silico analyses and detailed database searches were used to predict the 5′-CGI(s) in each gene. Primers were designed to amplify ∼0.9 to 1.4 kb fragments encompassing the 5′-CGI(s) of a candidate gene from bisulfite-modified DNA (Table 8). Amplicons were generated from individual UCWBC sample and subcloned into the pCR2.1 vector (Invitrogen, CA); eight or more clones were picked and sequenced (Macrogen, Korea).
Bisulfite-modified genomic DNA was subjected to MSPCR analyses , . Three sets of primers were designed to amplify regions (Figure 3A; R1, R2, and R3) within the 5′-CGI that exhibited the highest degree of differential methylation between the high and low PAH group. For each region, the M primer set amplifies the target sequence when it is methylated; the U primer set only amplifies the unmethylated allele (Table 8). Thirty-five-cycle PCRs were performed with 1 unit Platinum Taq DNA polymerase (Invitrogen, CA) using the following conditions: denature at 94°C for 30 s, anneal at 58°C for 1-min, and extend at 72°C for 1-min, followed by 12-min final extension. PCR products were separated on 2% agarose gel and visualized with ethidium bromide.
Real-time reverse-transcriptase-PCR (real time RT-PCR)
Total RNA was isolated from FPT samples matched to UCWBC samples or from H1299 cells as previously reported , . Steady state mRNA levels of the genes under investigation were quantified by SYBR Green-based real time RT-PCR with normalization to β-actin as previously described , . Primers specific for quantification of transcripts of the selected genes are listed in Table 8. Three separate real time RT-PCRs were performed on each FPT sample to obtain a mean value. The 2-ΔΔCt method was used to calculate the relative expression level of transcripts normalized to β-actin. The mean value from a universal reference RNA (Clontech, CA) was arbitrarily assigned an abundance value of 1.00 for each gene. Mean values from all other samples were normalized against this sample.
The Kendall coefficient of concordance (tau) was calculated to measure the strength of the relationship between the percent of methylation of a gene's promoter CGI in UCWBC DNA samples and its transcript expression level in matched FPT RNA samples. Ninety-five percent confidence intervals and 2-tailed p-values were obtained to determine the expected range and significance of observed tau values.
Chi-square tests were performed to investigate possible differences in the percents of African Americans, males, and methylated subjects between categories of PAH exposure, (<2.41, ≥2.41). Odds ratios and 95% confidence intervals of high PAH exposure for each of the three characteristics were calculated and tested for statistical significance. The cutpoint of 2.41 was determined from an ROC curve analysis following logistic regression of methylation status on continuously measured PAH exposure. This value separated methylated and unmethylated subjects according to ‘high/low’ PAH exposure. We considered MSPCR to be methylated above 2.41 and unmethylated below 2.41.
Logistic regression was performed to evaluate the association between parental report of asthma (dependent variable) and methylation status (independent variable). First, methylation status was modeled alone; then gender and ethnicity were added to the model. The parameter estimate and standard error of methylation status changed less then 15% when adjusted for gender and ethnicity, therefore the unadjusted odds ratio was presented. The equality of the median PAH exposure levels between methylated and unmethylated subjects was tested by the non-parametric Kruskal-Wallis test. For all analyses a 5% alpha level indicated significance unless stated otherwise. Analyses were performed using SAS, Version 9.1 (SAS Institute, Cary, NC).
Conceived and designed the experiments: FP WyT SmH. Performed the experiments: WyT JH DT LL. Analyzed the data: FP WyT JH LL RM SmH. Contributed reagents/materials/analysis tools: FP WyT JH DT LL SmH. Wrote the paper: FP WyT JH DT LL RM SmH.
- 1. Mannino M (1998) in Prevention, CDC. Surveillance summaries in surveillance for asthma-United States, 1960–1995. MMWR 47: 1–28.
- 2. Selgrade MK, Lemanske RF Jr, Gilmour MI, Neas LM, Ward MD, et al. (2006) Induction of asthma and the environment: what we know and need to know. Environ Health Perspect 114: 615–619.
- 3. Nicholas SW, Jean-Louis B, Ortiz B, Northridge M, Shoemaker K, et al. (2005) Addressing the childhood asthma crisis in Harlem: the Harlem Children's Zone Asthma Initiative. Am J Public Health 95: 245–249.
- 4. Perera FP, Rauh V, Tsai WY, Kinney P, Camann D, et al. (2003) Effects of transplacental exposure to environmental pollutants on birth outcomes in a multiethnic population. Environ Health Perspect 111: 201–205.
- 5. Miller RL, Garfinkel R, Horton M, Camann D, Perera FP, et al. (2004) Polycyclic aromatic hydrocarbons, environmental tobacco smoke, and respiratory symptoms in an inner-city birth cohort. Chest 126: 1071–1078.
- 6. Perzanowski MS, Miller RL, Thorne PS, Barr RG, Divjan A, et al. (2006) Endotoxin inner-city homes and the association with atopy and wheeze in the first two years of life. J Allergy Clin Immunology 117: 1082–1089.
- 7. Padmanabhan V (2007) Environment and origin of disease. Rev Endocr Metab Disord 8: 67–69.
- 8. Prins GS, Birch L, Tang WY, Ho SM (2007) Developmental estrogen exposures predispose to prostate carcinogenesis with aging. Reprod Toxicol 23: 374–382.
- 9. Hilakivi-Clarke L, de Assis S (2006) Fetal origins of breast cancer. Trends Endocrinol Metab 17: 340–348.
- 10. Tang WY, Ho SM (2007) Epigenetic reprogramming and imprinting in origins of disease. Rev Endocr Metab Disord 8: 173–182.
- 11. Miller RL, Ho SM (2008) Environmental Epigenetics and Asthma: Current Concepts and Call for Studies. Am J Respir Crit Care Med 177: 567–573.
- 12. Li S, Hursting SD, Davis BJ, McLachlan JA, Barrett JC (2003) Environmental exposure, DNA methylation, and gene regulation: lessons from diethylstilbesterol-induced cancers. Ann NY Acad Sci 983: 161–169.
- 13. Anway MD, Cupp AS, Uzumcu M, Skinner MK (2005) Epigenetic transgenerational actions of endocrine disruptors and male fertility. Science 308: 1466–1469.
- 14. Ho SM, Tang WY, Belmonte dF, Prins GS (2006) Developmental exposure to estradiol and bisphenol A increases susceptibility to prostate carcinogenesis and epigenetically regulates phosphodiesterase type 4 variant 4. Cancer Res 66: 5624–5632.
- 15. Dolinoy DC, Huang D, Jirtle RL (2007) Maternal nutrient supplementation counteracts bisphenol A-induced DNA hypomethylation in early development. Proc Natl Acad Sci USA 104: 13056–13061.
- 16. Tang WY, Newbold R, Mardilovich K, Jefferson W, Cheng RY, et al. (2008) Persistent hypomethylation in the promoter of nucleosomal binding protein 1(Nsbp1) correlates with overexpression of Nsbp1 in mouse uteri neonatally exposed to diethylstilbestrol or genistein. Endocrinology. 2008 Jul 31. [Epub ahead of print].
- 17. Weber P (2001) Role of biomarkers in nutritional science and industry- a commentary. Brit J of Nutrition 86: Suppl. 1S93–S95.
- 18. Stock W, Hoffman R (2000) White blood cells 1: non-malignant disorders. Lancet 355: 1351–1357.
- 19. Huang TH, Laux DE, Hamlin BC, Tran P, Tran H, et al. (1997) Identification of DNA methylation markers for human breast carcinomas using the methylation-sensitive restriction fingerprinting technique. Cancer Res 57: 1030–1034.
- 20. Wu M, Ho SM (2004) PMP24, a gene identified by MSRF, undergoes DNA hypermethylation-associated gene silencing during cancer progression in an LNCaP model. Oncogene 23: 250–259.
- 21. Khaled AR, Bulavin DV, Kittipatarin C, Li WQ, Alvarez M, et al. (2005) Cytokine-driven cell cycling is mediated through Cdc25A. J Cell Biol 169: 755–763.
- 22. Birkenbihl RP, Subramani S (1992) Cloning and characterization of rad21 an essential gene of Schizosaccharomyces pombe involved in DNA double-strand-break repair. Nucleic Acids Res 20: 6605–6011.
- 23. Van Horn CG, Caviglia JM, Li LO, Wang S, Granger DA, et al. (2005) Characterization of recombinant long-chain rat acyl-CoA synthetase isoforms 3 and 6: identification of a novel variant of isoform 6. Biochemistry 44: 1635–42.
- 24. Miyazaki M, Ntambi JM (2003) Role of stearoyl-coenzyme A desaturase in lipid metabolism. Prostaglandins Leukot Essent Fatty Acids 68: 113–121.
- 25. Kuzmin A, Han Z, Golding MC, Mann MR, Latham KE, et al. (2008) The PcG gene Sfmbt2 is paternally expressed in extraembryonic tissues. Gene Expr Patterns 8: 107–116.
- 26. Aqeilan RI, Donati V, Palamarchuk A, et al. (2005) WW domain-containing proteins, WWOX and YAP, compete for interaction with ErbB-4 and modulate its transcriptional function. Cancer Res 65: 6764–6772.
- 27. Ho SM, Tang WY (2007) Techniques used in studies of epigenome dysregulation due to aberrant DNA methylation: an emphasis on fetal-based adult diseases. Reprod Toxicol 23: 267–282.
- 28. Herman JG, Graff JR, Myohanen S, Nelkin BD, Baylin SB (1996) Methylation-specific PCR: a novel PCR assay for methylation status of CpG islands. Proc Natl Acad Sci USA 93: 9821–9826.
- 29. Gyorffy E, Anna L, Kovacs P, Schoket B (2008) Correlation between biomarkers of human exposure to genotoxins with focus on carcinogen-DNA adducts. Mutagenesis 23: 1–18.
- 30. Godschalk RW, Van Schooten FJ, Bartsch H (2003) A critical evaluation of DNA adducts as biological markers for human exposure to polycyclic aromatic compounds. J Biochem Mol Biol 36: 1–11.
- 31. Lee D-U, Agarwal S, Rao A (2002) Th2 lineage commitment and efficient IL-4 production involves extended demethylation of the IL-4 gene. Immunity 16: 649–660.
- 32. Agarwal S, Rao A (1998) Long-range transcriptional regulation of cytokine gene expression. Curr Opin Immunol 10: 345–352.
- 33. Tykocinski LO, Hajkova P, Chang HD, Stamm T, Sozeri O, et al. (2005) A critical control element for interleukin-4 memory expression in T helper lymphocytes. J Biol Chem 280: 28177–28185.
- 34. Jones B, Chen J (2006) Inhibition of IFN-gamma transcription by site-specific methylation during T helper cell development. EMBO J 25: 2443–2452.
- 35. Liu J, Ballaney M, Al Alem U, Quan C, Jin X, et al. (2007) Combined inhaled diesel exhaust particles and allergen exposure alter methylation of T helper genes and IgE production in vivo. Toxicol Sci 102: 76–81.
- 36. Mashek DG, Bornfeldt KE, Coleman RA, Berger J, Bernlohr DA, et al. (2004) Revised nomenclature for the mammalian long-chain acyl-CoA synthetase gene family. J Lipid Res 45: 1958–1961.
- 37. Minekura H, Kang MJ, Inagaki Y, Suzuki H, Sato H, et al. (2001) Genomic organization and transcription units of the human acyl-CoA synthetase 3 gene. Gene 278: 185–192.
- 38. Fujino T, Kang MJ, Suzuki H, Iijima H, Yamamoto T (1996) Molecular characterization and expression of rat acyl-CoA synthetase 3. J Biol Chem 271: 16748–16752.
- 39. Faergeman NJ, Knudsen J (1997) Role of long-chain fatty acyl-CoA esters in the regulation of metabolism and in cell signaling. Biochem J 323(Pt 1): 1–12.
- 40. Coleman RA, Lee DP (2004) Enzymes of triacylglycerol synthesis and their regulation. Prog Lipid Res 43: 134–176.
- 41. Tomita M, Baker RC, Ando S, Santoro TJ (2004) Arachidonoyl-phospholipid remodeling in proliferating murine T cells. Lipids Health Dis 3: 1.
- 42. Das DK (1982) Acyl CoA synthetase from rat lung: purification and properties. Lung 160: 207–218.
- 43. Mashima T, Oh-hara T, Sato S, Mochizuki M, Sugimoto Y, et al. (2005) p53-defective tumors with a functional apoptosome-mediated pathway: a new therapeutic target. J Natl Cancer Inst 97: 765–777.
- 44. Wijga A, Houwelingen AC, Smit HA, Kerkhof M, Vos AP, et al. (2003) Fatty acids in breast milk of allergic and non-allergic mothers: The PIAMA birth cohort study. Pediatr Allergy Immunol 14: 156–162.
- 45. Oddy WH, Pal S, Kusel MM, Vine D, de Klerk NH, et al. (2006) Atopy, eczema and breast milk fatty acids in a high-risk cohort of children followed from birth to 5 yr. Pediatr Allergy Immunol 17: 4–10.
- 46. Mori TA, Beilin LJ (2004) Omega-3 fatty acids and inflammation. Curr Atheroscler Rep 6: 461–467.
- 47. Woods RK, Raven JM, Walters EH, Abramson MJ, Thien FC (2004) Fatty acid levels and risk of asthma in young adults. Thorax 59: 105–110.
- 48. Beck M, Zelczak G, Lentze M-J (2000) Abnormal fatty acid composition in umbilical cord blood of infants at high risk of atopic disease. Acta Paediatr 89: 279–284.
- 49. Bouzigon E, Siroux V, Dizier MH, Lemainque A, Pison C, et al. (2007) Scores of asthma and asthma severity reveal new regions of linkage in EGEA study families. Eur Respir J 30(2): 253–259.
- 50. Choudhry S, Taub M, Mei R, Rodriguez-Santana J, Rodriguez-Cintron W, et al. (2008) Genome-wide screen for asthma in Puerto Ricans: evidence for association with 5q23 region. Hum Genet 123(5): 455–468.
- 51. Jedrychowski W, Galas A, Pac A, Flak E, Camman D, et al. (2005) Prenatal ambient air exposure to polycyclic aromatic hydrocarbons and the occurrence of respiratory symptoms over the first year of life. Eur J Epidemiol 20: 775–782.