The Nasal Methylome: A Key to Understanding Allergic Asthma.

Asthma is a complex disease that is heritable (1); however, common genetic variants explain a small portion of asthma risk (2). Epigenetic changes are influenced by the environment (3) and affect the expression of transcription factors that alter the maturation of T lymphocytes (4) and could potentially explain some of the missing heritability in asthma. We have previously reported that allergic asthma–associated DNA methylation changes in nasal epithelia are more substantial in magnitude and have a stronger association with gene expression (5) than we reported for peripheral blood mononuclear cells (PBMCs) (6). Others have reported an IL-13–induced DNA signature in asthmatic airway epithelia (7). However, no study to date has compared asthma-associated DNA methylation marks across different tissues within an individual. The goal of this investigation was to establish proof of concept that nasal epithelium can be used as a proxy for the airway epithelium in studies of allergic asthma. We collected PBMCs, nasal epithelia, and bronchial epithelia from 12 subjects with allergic asthma and 12 control subjects without asthma, all non-Hispanic white nonsmoker adults. There are no significant differences in age (median [low, high]: 37 yr [24, 74 yr] in patients with asthma vs. 29 yr [24, 45 yr] in control subjects; P = 0.16 by Mann-Whitney U test) and sex (33% men in patients with asthma vs. 50% men in control subjects; P = 0.59 by chi-square test) between cases and control subjects; however, serum IgE concentrations are significantly higher in patients with asthma (median [low, high]: cases = 115 [7, 1,085] vs. control subjects = 25.5 [7, 71], P = 0.003 by Mann-Whitney U test). By study design, all patients with asthma were on inhaled corticosteroids for at least 4 weeks before sample collection and not on any nasal steroids. Following standard approaches to quality control (5, 6), we identified single CpG sites (from Illumina 450k array; Illumina, San Diego, CA) or differentially methylated positions (DMPs) associated with allergic asthma after controlling for age, sex, batch effects (array position as identified by principal components analysis), and multiple comparisons (false discovery rate) and after filtering out probe sets with known singlenucleotide polymorphisms at the CpG site in the CEU (Utah residents with ancestry from northern and western Europe) population. To maximize the power for DMP analysis, we first used a mixed-effects linear model for all samples adjusting for tissue and identified 9,598 significant DMPs with false discovery rate–adjusted q value of 0.005. We then ran individual tissue models to identify tissue-specific changes in DNA methylation and determined that, of the 9,598 DMPs, none were significant in PBMCs, 7,484 were significant in nasal epithelia, and 63 were significant in airway epithelia (tissuespecific model q, 0.05). Focusing on the largest effect size (.10% methylation change) resulted in 816 DMPs in nasal epithelia and 41 DMPs in bronchial epithelia (Figure 1). Forty-one DMPs in airway epithelia are associated with 34 unique genes, of which 25 have significant DMPs with greater than 10% change in nasal cells and an additional 7 have significant DMPs with less than 10% change in nasal epithelial cells. To determine the effect of allergic asthma–associated methylation changes on gene expression, we examined array-based (Agilent 83 60k; Agilent, Palo Alto, CA) measures of expression of genes nearest DMPs as a function of methylation change. As in our previous work (5, 6), we identified an enrichment of inverse correlations of DNA methylation and expression in nasal epithelia (P = 2.23 10 for DMPs; Figure 2). This enrichment was not observed in airway epithelia, presumably owing to the small number of data points (Figure 2). Some of the genes with the most prominent changes in DNA methylation and gene expression are consistent in nasal and airway epithelia (CAPN14, CTSC, and FBXL7, among others), whereas we identified other genes with larger changes in nasal epithelia (POSTN, CCL26, GATA2, HLA genes, to name a few). The majority of the genes with the most pronounced methylation–expression relationships were also identified in our previous nasal epithelial study (5), despite differences in age, ethnicity, and potential exposures between the two cohorts. Our study is the first to directly compare DNA methylation and gene expression changes in PBMCs, nasal epithelia, and airway epithelia, and highlights three important points. First, nasal epithelia appear to be a good proxy for airway epithelia in studies of allergic asthma, as almost all allergic asthma–associated methylation changes are captured by nasal epithelia, and many of them have relationships with gene expression in both tissues. This finding is supported by published work demonstrating high concordance of nasal and airway transcriptomes in asthma (8). Second, we identified many more asthmaassociated DNA methylation changes in nasal epithelia than airway epithelia, likely because nasal epithelium is the primary interface with the environment, interacting with allergens and other environmental stimuli. Although it may appear surprising that the asthma airway epithelium biomarker gene POSTN (9) was only identified in the nasal epithelia in our analysis, its methylation change was similar in the two tissues (14.5% hypomethylated in nasal and 13.6% hypomethylated in airway epithelia) but more significant in nasal epithelia, in which it reached genome-wide significant threshold (q = 0.01), whereas this was not the case in airway epithelia (q = 0.12). POSTN is sevenfold up-regulated in nasal and fourfold in airway epithelia. Third, although studies of peripheral blood can shed light on epigenetic and transcriptional changes in immune cells in asthma (6, 10), large sample sizes are needed because of a weaker signal and mixture of cell types that are not specific Supported by NHLBI grant R01-HL101251 and the National Institute for Environmental Health Sciences grant P01-ES18181.