findings, we were interested in examining whether methylation of specific CpG sites is 67 associated with expression of non-nearby cis-genes. We thus conducted an expression 68 quantitative trait methylation (eQTM) analysis in nasal airway epithelium from 455 Puerto 69Ricans ages 9 to 20 years, including 219 subjects with asthma (cases) and 236 control subjects. 70 71
Results
72Location of eQTM-methylation probes relative to paired genes 73 By testing associations between methylation probes within 1 Mb of transcription start sites (TSS) 74 of genes and gene expression, we identified 16,867 significant methylation-expression pairs 75 (FDR-P < 0.01, see Methods), comprising 9,103 methylation probes associated with expression 76 of 3,512 genes. We then investigated the position of significant methylation probes in this eQTM 77 analysis in relation to their paired genes. If a methylation probe was associated with expression 78 of multiple genes, we counted such probe for each gene. We found that 11% and 89% of 79 significant eQTM probes were located within and outside genes, respectively -including 4% of 80 eQTM probes in promoter regions (Figure 1a). Most eQTM methylation probes were distant 81 from their target genes (371,840 bp on average) (Figure 1b). 82 83 Because we found distant relationships between methylation and their target genes, we assessed 84 whether eQTM methylation probes in nasal epithelium are enriched in the enhancer regions of 85 their paired genes. For this, we checked the enhancer database for lung tissue 86 (http://enhanceratlas.org/) 3 , since nasal epithelial tissue was not available in that database, and 87 nasal and bronchial epithelial methylation and expression are well correlated 1 . We found that 88 eQTM methylation probes are more likely to be located in enhancer regions of their paired genes 89 than randomly selected control probes (p-value: 1.3x10 -213 ) (Figure 1c, Table 1). 90 91 While most methylation probes near TSS were negatively correlated with gene expression, more 92 distant pairs tended to be positively correlated (Figure 1d). Of the eQTM methylation probes 93 associated with expression of the gene they were located in, 81.9% were negatively correlated 94 with expression (85.3% if in promoter regions) (Figure 1e). In contrast, only 40.2 % of eQTM 95 methylation probes associated with expression of a distant gene (i.e. methylation probes outside 96 of the associated gene) were negatively correlated with expression levels (Figure 1f). 97 98 Most of the top eQTM genes (by eQTM P-value) have been implicated in lung disease (Figure 99 2). PAX8 is associated with bronchodilator response in children with asthma 4 , ECHDC3 is 100 associated with obesity and asthma in children 5 , LSP1 is associated with acute lung inflammation 101 6 , HLA-DQB1 is associated with asthma 7 and total IgE 8 , FRG1B is highly mutated in lung 102 adenocarcinoma 9 , and KANSL1 is associated with pulmonary function 10 . 103 104
Gene Ontology enrichment analysis 105We performed a Gene Ontology enrichment ...