Study population. This study was nested in a panel study with repeated measurements on TRAP and FeNO. From 2011 to 2016, we repeatedly sampled 36 children living in Newark or Elizabeth NJ with mild-to-moderate physician-diagnosed asthma. Each subject was followed for up to 30 days. Power calculations to examine the within-subject association between either exposure and methylation, or methylation and FeNO, were conducted via simulation. Simulations assumed 18 subjects with an average of 5 repeated measures per subject, using a sum across subjects of the fisher-transformed correlations within subject, which follows a normal distribution. Assuming we conducted a two-sided hypothesis test that the correlation was equal to zero, we would have 75.9% power to detect a correlation of 0.35. For within-subject correlations of 0.40 and 0.45, we would have powers of 85.7% and 92.1%, respectively. Therefore, a subset of 18 children (5 buccal samples per child) was randomly selected for the DNA methylation study. The inclusion criteria are: subjects with at least five buccal cell samples; with FeNO measurements, and with less than 20% missing of 24-hour BC data. The Institutional Review Board of Rutgers University approved the epigenetic study protocol. Parents or legal guardians provided informed consent for all subjects.
FeNO measurements. FeNO was measured using the NIOX MINO (Aerocrine, New Providence, NJ) following the manufacturer’s instructions and the standards of the American Thoracic Society and the European Respiratory Society. Concentrations of NO were measured during 10-second exhalations of breath at an exhalation pressure of 10–20 cmH2O to maintain a flow rate of 50 ± 5 ml/second. To control for the possible circadian effects, FeNO measurements were conducted at the same time (around 4:00 p.m.) for all subjects on the same day when collecting buccal cell samples.
Buccal cell samples. Buccal cell samples were collected at the community field site on various weekdays to better capture the variation in DNA methylation on different day-of-week. Each child was provided with two toothbrushes. They were instructed to remove all food from their mouths and used the first toothbrush to brush their teeth before buccal cell collection. Researchers gently brushed buccal mucosa 10 strokes with the second toothbrush and then asked the subject to rinse the mouth with water for 30 seconds. After that, the subjects were told to spit the content in a tube and put the second brush in the expelled water. Isopropyl alcohol (70%) was added to the tube after sampling. In the lab, buccal cell suspensions were centrifuged at 2,500 revolutions per minute (1000 g) for 15 minutes. The pellets were stored frozen at –80°C until used for DNA extraction.
DNA methylation. We studied CpG loci located in promoter regions of NOS1 (6 CpG sites), NOS2A (16 CpG sites), NOS3 (3 CpG sites), ARG1 (8 CpG sites) and ARG2 (13 CpG sites) genes reported from a previous study (13). The location of the gene promoter, amplified regions, and CpG sites that were tested have been published previously (13). Lab technicians who performed DNA methylation analysis were blinded to subject information. Genomic DNA samples were extracted from cell pellets and subjected to bisulfite modification. Samples were then amplified by polymerase chain reaction (PCR) and analyzed using a PSQHS96 Pyrosequencing System (EpigenDx, Hopkinton, MA). Methylation status at each CpG site was measured by calculating the ratio of C (methylated cytosine) relative to T (unmethylated cytosine) (14). For each CpG site tested, the percentage of 5-methylcytosine (%5mC) was calculated and presented as the degree of methylation. For all assays, methylated and unmethylated DNA standards were used as positive and negative control, respectively to validate the method of bisulfite conversion and PCR amplification.
Air pollutants. Personal real-time BC data was recorded by micro-aethalometer (AE51, Aethlabs, Oakland, CA) for up to 30 consecutive days. Subjects carried the monitor on a belt during waking hours, and recharged the monitor at bedside while asleep. NO2 were sampled by passive personal sampler (Ogawa & Co., Pompano Beach, FL) and analyzed by electron absorption spectroscopy. Then the integrated 24-hr NO2 data were obtained for up to 30 consecutive days. Data cleaning algorithms were applied to address artifactual high and low BC readings from the aethalometer (15, 16).
Covariates. Demographic information, including age, gender, and race/ethnicity, were obtained through questionnaires completed by parent or legal guardians, and height and weight were measured, at the beginning of the study.
Statistical analysis. We used SAS software version 9.4 for all analyses. Descriptive analyses were conducted to understand the distributions of air pollution, DNA methylation at each locus and FeNO measurements in general and by subject characteristics. Outcome and exposure were treated as continuous variables in all models. Spearman correlation assessed the associations of percent methylation between loci in the same gene. Since a significantly high correlation of DNA methylations between loci was observed, we applied multivariate analysis of variance (MANOVA) in generalized linear models to study the association between percent methylation and BC exposure while controlling for the correlation of %5mC among loci in the same gene (Supplement Table 1-5).
For genes with statistically significant results in MANOVA tests, we fitted four linear mixed effect models to assess the effects of BC on FeNO at multiple periods prior to the buccal sample collection (Figure 1): 0-6 hours, 7-12 hours, 13-24 hours, 0-24 hours (lag 0 day), 25-48 hours (lag 1 day) 49-72 hours (lag 2 days) 73-96 hours (lag 3 days). The four models include: unadjusted model (model 1); model adjusted for week number and day-of-week when collecting the outcomes (model 2); model further adjusted for age, gender and race/ethnicity (model 3) and model further adjusted for NO2 (model 4). BC and all covariates were incorporated as fixed-effect terms while a random-effect term for each subject was added to account for the correlations of repeated measurements collected from the same subject. To better understand the role of epigenetic changes in response to TRAP at varying time lag periods, we examined the association between BC and percent DNA methylation in models 1, 2, and 3 at all lag periods.
All tests were 2-sided at a 5% significance level. The effect estimates of DNA methylation were presented as the change and its 95% confidence interval (CIs) in %5mC per log-transformed and lag-specific interquartile range (IQR) increase in BC. The effect estimates of FeNO were presented as the relative change and its 95% CIs per log-transformed and lag-specific IQR increase in BC.
Missing data problems are usually addressed by including only subjects without missing data in any variables required for an analysis. Such method is subject to bias and loss of information (17). We conducted a sensitivity analysis by using multiple imputation for missing BC and NO2 data and fitted the same models to examine the lag patterns of imputed BC, DNA methylation and FeNO levels (Supplement Table 5-7).