Epigenome-wide association study reveals longitudinally stable DNA methylation differences in CD4+ T cells from children with IgE-mediated food allergy

Food allergy is mediated by a combination of genetic and environmental risk factors, potentially mediated by epigenetic mechanisms. CD4+ T-cells are key drivers of the allergic response, and may therefore harbor epigenetic variation in association with the disease phenotype. Here we retrospectively examined genome-wide DNA methylation profiles (~450?000 CpGs) from CD4+ T-cells on a birth cohort of 12 children with IgE-mediated food allergy diagnosed at 12-months, and 12 non-allergic controls. DNA samples were available at two time points, birth and 12-months. Case:control comparisons of CD4+ methylation profiles identified 179 differentially methylated probes (DMP) at 12-months and 136 DMP at birth (FDR-adjusted P value < 0.05, delta β > 0.1). Approximately 30% of DMPs were coincident with previously annotated SNPs. A total of 96 allergy-associated non-SNP DMPs were present at birth when individuals were initially disease-free, potentially implicating these loci in the causal pathway. Pathway analysis of differentially methylated genes identified several MAP kinase signaling molecules. Mass spectrometry was used to validate 15 CpG sites at 3 candidate genes. Combined analysis of differential methylation with gene expression profiles revealed gene expression differences at some but not all allergy associated differentially methylated genes. Thus, dysregulation of DNA methylation at MAPK signaling-associated genes during early CD4+ T-cell development may contribute to suboptimal T-lymphocyte responses in early childhood associated with the development of food allergy.     

Genetic Control of Individual Differences in Gene-Specific Methylation in Human Brain

We have observed extensive interindividual differences in DNA methylation of 8590 CpG sites of 6229 genes in 153 human adult cerebellum samples, enriched in CpG island “shores” and at further distances from CpG islands. To search for genetic factors that regulate this variation, we performed a genome-wide association study (GWAS) mapping of methylation quantitative trait loci (mQTLs) for the 8590 testable CpG sites. cis association refers to correlation of methylation with SNPs within 1 Mb of a CpG site. 736 CpG sites showed phenotype-wide significant cis association with 2878 SNPs (after permutation correction for all tested markers and methylation phenotypes). In trans analysis of methylation, which tests for distant regulation effects, associations of 12 CpG sites and 38 SNPs remained significant after phenotype-wide correction. To examine the functional effects of mQTLs, we analyzed 85 genes that were with genetically regulated methylation we observed and for which we had quality gene expression data. Ten genes showed SNP-methylation-expression three-way associations—the same SNP simultaneously showed significant association with both DNA methylation and gene expression, while DNA methylation was significantly correlated with gene expression. Thus, we demonstrated that DNA methylation is frequently a heritable continuous quantitatively variable trait in human brain. Unlike allele-specific methylation, genetic polymorphisms mark both cis- and trans-regulatory genetic sites at measurable distances from their CpG sites. Some of the genetically regulated DNA methylation is directly connected with genetically regulated gene expression variation.     

DNA hypermethylation and DNA hypomethylation is present at different loci in chronic kidney disease

Genetic risk factors for chronic kidney disease (CKD) are being identified through international collaborations. By comparison, epigenetic risk factors for CKD have only recently been considered using population-based approaches. DNA methylation is a major epigenetic modification that is associated with complex diseases, so we investigated methylome-wide loci for association with CKD. A total of 485,577 unique features were evaluated in 255 individuals with CKD (cases) and 152 individuals without evidence of renal disease (controls). Following stringent quality control, raw data were quantile normalized and β values calculated to reflect the methylation status at each site. The difference in methylation status was evaluated between cases and controls with resultant P values adjusted for multiple testing. Genes with significantly increased and decreased levels of DNA methylation were considered for biological relevance by functional enrichment analysis using KEGG pathways in Partek Genomics Suite. Twenty-three genes, where more than one CpG per loci was identified with P_adjusted < 10^?8, demonstrated significant methylation changes associated with CKD and additional support for these associated loci was sought from published literature. Strong biological candidates for CKD that showed statistically significant differential methylation include CUX1, ELMO1, FKBP5, INHBA-AS1, PTPRN2, and PRKAG2 genes; several genes are differentially methylated in kidney tissue and RNA-seq supports a functional role for differential methylation in ELMO1 and PRKAG2 genes. This study reports the largest, most comprehensive, genome-wide quantitative evaluation of DNA methylation for association with CKD. Evidence confirming methylation sites influence development of CKD would stimulate research to identify epigenetic therapies that might be clinically useful for CKD.     

Comparative Anatomy of Chromosomal Domains with Imprinted and Non-Imprinted Allele-Specific DNA Methylation

Allele-specific DNA methylation (ASM) is well studied in imprinted domains, but this type of epigenetic asymmetry is actually found more commonly at non-imprinted loci, where the ASM is dictated not by parent-of-origin but instead by the local haplotype. We identified loci with strong ASM in human tissues from methylation-sensitive SNP array data. Two index regions (bisulfite PCR amplicons), one between the C3orf27 and RPN1 genes in chromosome band 3q21 and the other near the VTRNA2-1 vault RNA in band 5q31, proved to be new examples of imprinted DMRs (maternal alleles methylated) while a third, between STEAP3 and C2orf76 in chromosome band 2q14, showed non-imprinted haplotype-dependent ASM. Using long-read bisulfite sequencing (bis-seq) in 8 human tissues we found that in all 3 domains the ASM is restricted to single differentially methylated regions (DMRs), each less than 2kb. The ASM in the C3orf27-RPN1 intergenic region was placenta-specific and associated with allele-specific expression of a long non-coding RNA. Strikingly, the discrete DMRs in all 3 regions overlap with binding sites for the insulator protein CTCF, which we found selectively bound to the unmethylated allele of the STEAP3-C2orf76 DMR. Methylation mapping in two additional genes with non-imprinted haplotype-dependent ASM, ELK3 and CYP2A7, showed that the CYP2A7 DMR also overlaps a CTCF site. Thus, two features of imprinted domains, highly localized DMRs and allele-specific insulator occupancy by CTCF, can also be found in chromosomal domains with non-imprinted ASM. Arguing for biological importance, our analysis of published whole genome bis-seq data from hES cells revealed multiple genome-wide association study (GWAS) peaks near CTCF binding sites with ASM.     

Gene expression analysis uncovers novel hedgehog interacting protein (HHIP) effects in human bronchial epithelial cells

Hedgehog interacting protein (HHIP) was implicated in chronic obstructive pulmonary disease (COPD) by genome-wide association studies (GWAS). However, it remains unclear how HHIP contributes to COPD pathogenesis. To identify genes regulated by HHIP, we performed gene expression microarray analysis in a human bronchial epithelial cell line (Beas-2B) stably infected with HHIP shRNAs. HHIP silencing led to differential expression of 296 genes; enrichment for variants nominally associated with COPD was found. Eighteen of the differentially expressed genes were validated by real-time PCR in Beas-2B cells. Seven of 11 validated genes tested in human COPD and control lung tissues demonstrated significant gene expression differences. Functional annotation indicated enrichment for extracellular matrix and cell growth genes. Network modeling demonstrated that the extracellular matrix and cell proliferation genes influenced by HHIP tended to be interconnected. Thus, we identified potential HHIP targets in human bronchial epithelial cells that may contribute to COPD pathogenesis.     

Epigenome-wide association study reveals longitudinally stable DNA methylation differences in CD4+ T cells from children with IgE-mediated food allergy

Food allergy is mediated by a combination of genetic and environmental risk factors, potentially mediated by epigenetic mechanisms. CD4+ T-cells are key drivers of the allergic response, and may therefore harbor epigenetic variation in association with the disease phenotype. Here we retrospectively examined genome-wide DNA methylation profiles (~450 000 CpGs) from CD4+ T-cells on a birth cohort of 12 children with IgE-mediated food allergy diagnosed at 12-months, and 12 non-allergic controls. DNA samples were available at two time points, birth and 12-months. Case:control comparisons of CD4+ methylation profiles identified 179 differentially methylated probes (DMP) at 12-months and 136 DMP at birth (FDR-adjusted P value < 0.05, delta β > 0.1). Approximately 30% of DMPs were coincident with previously annotated SNPs. A total of 96 allergy-associated non-SNP DMPs were present at birth when individuals were initially disease-free, potentially implicating these loci in the causal pathway. Pathway analysis of differentially methylated genes identified several MAP kinase signaling molecules. Mass spectrometry was used to validate 15 CpG sites at 3 candidate genes. Combined analysis of differential methylation with gene expression profiles revealed gene expression differences at some but not all allergy associated differentially methylated genes. Thus, dysregulation of DNA methylation at MAPK signaling-associated genes during early CD4+ T-cell development may contribute to suboptimal T-lymphocyte responses in early childhood associated with the development of food allergy.     

Tuberculosis infection and lung adenocarcinoma: Mendelian randomization and pathway analysis of genome-wide association study data from never-smoking Asian women

We investigated whether genetic susceptibility to tuberculosis (TB) influences lung adenocarcinoma development among never-smokers using TB genome-wide association study (GWAS) results within the Female Lung Cancer Consortium in Asia. Pathway analysis with the adaptive rank truncated product method was used to assess the association between a TB-related gene-set and lung adenocarcinoma using GWAS data from 5512 lung adenocarcinoma cases and 6277 controls. The gene-set consisted of 31 genes containing known/suggestive associations with genetic variants from previous TB-GWAS. Subsequently, we followed-up with Mendelian Randomization to evaluate the association between TB and lung adenocarcinoma using three genome-wide significant variants from previous TB-GWAS in East Asians. The TB-related gene-set was associated with lung adenocarcinoma (p = 0.016). Additionally, the Mendelian Randomization showed an association between TB and lung adenocarcinoma (OR = 1.31, 95% CI: 1.03, 1.66, p = 0.027). Our findings support TB as a causal risk factor for lung cancer development among never-smoking Asian women.     

Genome-scale analysis of DNA methylation in colorectal cancer using Infinium HumanMethylation450 BeadChips

Illumina’s Infinium HumanMethylation450 BeadChip arrays were used to examine genome-wide DNA methylation profiles in 22 sample pairs from colorectal cancer (CRC) and adjacent tissues and 19 colon tissue samples from cancer-free donors. We show that the methylation profiles of tumors and healthy tissue samples can be clearly distinguished from one another and that the main source of methylation variability is associated with disease status. We used different statistical approaches to evaluate the methylation data. In general, at the CpG-site level, we found that common CRC-specific methylation patterns consist of at least 15,667 CpG sites that were significantly different from either adjacent healthy tissue or tissue from cancer-free subjects. Of these sites, 10,342 were hypermethylated in CRC, and 5,325 were hypomethylated. Hypermethylated sites were common in the maximum number of sample pairs and were mostly located in CpG islands, where they were significantly enriched for differentially methylated regions known to be cancer-specific. In contrast, hypomethylated sites were mostly located in CpG shores and were generally sample-specific. Despite the considerable variability in methylation data, we selected a panel of 14 highly robust candidates showing methylation marks in genes SND1, ADHFE1, OPLAH, TLX2, C1orf70, ZFP64, NR5A2, and COL4A. This set was successfully cross-validated using methylation data from 209 CRC samples and 38 healthy tissue samples from The Cancer Genome Atlas consortium (AUC = 0.981 [95% CI: 0.9677–0.9939], sensitivity = 100% and specificity = 82%). In summary, this study reports a large number of loci with novel differential methylation statuses, some of which may serve as candidate markers for diagnostic purposes.     

An epigenome-wide association meta-analysis of prenatal maternal stress in neonates: A model approach for replication

Prenatal maternal stress exposure has been associated with neonatal differential DNA methylation. However, the available evidence in humans is largely based on candidate gene methylation studies, where only a few CpG sites were evaluated. The aim of this study was to examine the association between prenatal exposure to maternal stress and offspring genome-wide cord blood methylation using different methods. First, we conducted a meta-analysis and follow-up pathway analyses. Second, we used novel region discovery methods [i.e., differentially methylated regions (DMRs) analyses]. To this end, we used data from two independent population-based studies, the Generation R Study (n = 912) and the Avon Longitudinal Study of Parents and Children (ALSPAC, n = 828), to (i) measure genome-wide DNA methylation in cord blood and (ii) extract a prenatal maternal stress composite. The meta-analysis (n_total = 1,740) revealed no epigenome-wide (meta P <1.00e-07) associations of prenatal maternal stress exposure with neonatal differential DNA methylation. Follow-up analyses of the top hits derived from our epigenome-wide meta-analysis (meta P <1.00e-04) indicated an over-representation of the methyltransferase activity pathway. We identified no Bonferroni-corrected (P <1.00e-06) DMRs associated with prenatal maternal stress exposure. Combining data from two independent population-based samples in an epigenome-wide meta-analysis, the current study indicates that there are no large effects of prenatal maternal stress exposure on neonatal DNA methylation. Such replication efforts are essential in the search for robust associations, whether derived from candidate gene methylation or epigenome-wide studies.     

Intronic Parent-of-Origin Dependent Differential Methylation at the Actn1 Gene Is Conserved in Rodents but Is Not Associated with Imprinted Expression

Parent-of-origin differential DNA methylation has been associated with regulation of the preferential expression of paternal or maternal alleles of imprinted genes. Based on this association, recent studies have searched for parent-of-origin dependent differentially methylated regions in order to identify new imprinted genes in their vicinity. In a previous genome-wide analysis of mouse brain DNA methylation, we found a novel differentially methylated region in a CpG island located in the last intron of the alpha 1 Actinin (Actn1) gene. In this region, preferential methylation of the maternal allele was observed; however, there were no reports of imprinted expression of Actn1. Therefore, we have tested if differential methylation of this region is common to other tissues and species and affects the expression of Actn1. We have found that Actn1 differential methylation occurs in diverse mouse tissues. Moreover, it is also present in other murine rodents (rat), but not in the orthologous human region. In contrast, we have found no indication of an imprinted effect on gene expression of Actn1 in mice: expression is always biallelic regardless of sex, tissue type, developmental stage or isoform. Therefore, we have identified a novel parent-of-origin dependent differentially methylated region that has no apparent association with imprinted expression of the closest genes. Our findings sound a cautionary note to genome-wide searches on the use of differentially methylated regions for the identification of imprinted genes and suggest that parent-of-origin dependent differential methylation might be conserved for functions other that the control of imprinted expression.     

Inflammation-associated DNA methylation patterns in epithelium of ulcerative colitis

Aberrant DNA methylation patterns have been reported in inflamed tissues and may play a role in disease. We studied DNA methylation and gene expression profiles of purified intestinal epithelial cells from ulcerative colitis patients, comparing inflamed and non-inflamed areas of the colon. We identified 577 differentially methylated sites (false discovery rate <0.2) mapping to 210 genes. From gene expression data from the same epithelial cells, we identified 62 differentially expressed genes with increased expression in the presence of inflammation at prostate cancer susceptibility genes PRAC1 and PRAC2. Four genes showed inverse correlation between methylation and gene expression; ROR1, GXYLT2, FOXA2, and, notably, RARB, a gene previously identified as a tumor suppressor in colorectal adenocarcinoma as well as breast, lung and prostate cancer. We highlight targeted and specific patterns of DNA methylation and gene expression in epithelial cells from inflamed colon, while challenging the importance of epithelial cells in the pathogenesis of chronic inflammation.     

Genome-wide methylation profiling identifies novel methylated genes in neuroblastoma tumors

Neuroblastoma is a very heterogeneous tumor of childhood. The clinical spectra range from very aggressive metastatic disease to spontaneous regression, even without therapy. Aberrant DNA methylation pattern is a common feature of most cancers. For neuroblastoma, it has been demonstrated both for single genes as well as genome-wide, where a so-called methylator phenotype has been described. Here, we present a study using Illumina 450K methylation arrays on 60 neuroblastoma tumors. We show that aggressive tumors, characterized by International Neuroblastoma Risk Group (INRG) as stage M, are hypermethylated compared to low-grade tumors. On the contrary, INRG stage L tumors display more non-CpG methylation. The genes with the highest number of hypermethylated CpG sites in INRG M tumors are TERT, PCDHGA4, DLX5, and DLX6-AS1. Gene ontology analysis showed a representation of neuronal tumor relevant gene functions among the differentially methylated genes. For validation, we used a set of independent tumors previously analyzed with the Illumina 27K methylation arrays, which confirmed the differentially methylated sites. Top candidate genes with aberrant methylation were analyzed for altered gene expression through the R2 platform (http://r2.amc.nl), and for correlations between methylation and gene expression in a public dataset. Altered expression in nonsurvivors was found for the genes B3GALT4 and KIAA1949, CLIC5, DLX6-AS, TERT, and PIRT, and strongest correlations were found for TRIM36, KIAA0513, and PIRT. Our data indicate that methylation profiling can be used for patient stratification and informs on epigenetically deregulated genes with the potential of increasing our knowledge about the underlying mechanisms of tumor development.     

Epigenetic modification of retinoic acid-treated human embryonic stem cells

Epigenetic modification of the genome through DNA methylation is the key to maintaining the differentiated state of human embryonic stem cells (hESCs), and it must be reset during differentiation by retinoic acid (RA) treatment. A genome-wide methylation/gene expression assay was performed in order to identify epigenetic modifications of RA-treated hESCs. Between undifferentiated and RA-treated hESCs, 166 differentially methylated CpG sites and 2,013 differentially expressed genes were discovered. Combined analysis of methylation and expression data revealed that 19 genes (STAP2, VAMP8, C10orf26, WFIKKN1, ELF3, C1QTNF6, C10orf10, MRGPRF, ARSE, LSAMP, CENTD3, LDB2, POU5F1, GSPT2, THY1, ZNF574, MSX1, SCMH1, and RARB) were highly correlated with each other. The results provided in this study will facilitate future investigations into the interplay between DNA methylation and gene expression through further functional and biological studies.     

Two-marker association tests yield new disease associations for coronary artery disease and hypertension

It has been postulated that multiple-marker methods may have added ability, over single-marker methods, to detect genetic variants associated with disease. The Wellcome Trust Case Control Consortium (WTCCC) provided the first successful large genome-wide association studies (GWAS) which included single-marker association analyses for seven common complex diseases. Of those signals detected, only one was associated with coronary artery disease (CAD), and none were identified for hypertension (HTN). Our objective was to find additional genetic associations and pathways for cardiovascular disease by examining the WTCCC data for variants associated with CAD and HTN using two-marker testing methods. We applied two-marker association testing to the WTCCC dataset, which includes ~2,000 affected individuals with each disorder, and a shared pool of ~3,000 controls, all genotyped using Affymetrix GeneChip 500 K arrays. For CAD, we detected single nucleotide polymorphisms (SNP) pairs in three genes showing genome-wide significance: HFE2, STK32B, and DIPC2. The most notable SNP pairs in a non-protein-coding region were at 9p21, a known major CAD-associated region. For HTN, we detected SNP pairs in five genes: GPR39, XRCC4, MYO6, ZFAT, and MACROD2. Four further associated SNP pair regions were at least 70 kb from any known gene. We have shown that novel, multiple-marker, statistical methods can be of use in finding variants in GWAS. We describe many new, associated variants for both CAD and HTN and describe their known genetic mechanisms.     

Refining Susceptibility Loci of Chronic Obstructive Pulmonary Disease with Lung eqtls

Chronic obstructive pulmonary disease (COPD) is the fourth leading cause of mortality worldwide. Recent genome-wide association studies (GWAS) have identified robust susceptibility loci associated with COPD. However, the mechanisms mediating the risk conferred by these loci remain to be found. The goal of this study was to identify causal genes/variants within susceptibility loci associated with COPD. In the discovery cohort, genome-wide gene expression profiles of 500 non-tumor lung specimens were obtained from patients undergoing lung surgery. Blood-DNA from the same patients were genotyped for 1,2 million SNPs. Following genotyping and gene expression quality control filters, 409 samples were analyzed. Lung expression quantitative trait loci (eQTLs) were identified and overlaid onto three COPD susceptibility loci derived from GWAS; 4q31 (HHIP), 4q22 (FAM13A), and 19q13 (RAB4B, EGLN2, MIA, CYP2A6). Significant eQTLs were replicated in two independent datasets (n = 363 and 339). SNPs previously associated with COPD and lung function on 4q31 (rs1828591, rs13118928) were associated with the mRNA expression of HHIP. An association between mRNA expression level of FAM13A and SNP rs2045517 was detected at 4q22, but did not reach statistical significance. At 19q13, significant eQTLs were detected with EGLN2. In summary, this study supports HHIP, FAM13A, and EGLN2 as the most likely causal COPD genes on 4q31, 4q22, and 19q13, respectively. Strong lung eQTL SNPs identified in this study will need to be tested for association with COPD in case-control studies. Further functional studies will also be needed to understand the role of genes regulated by disease-related variants in COPD.     

An approach based on a genome-wide association study reveals candidate loci for narcolepsy

Narcolepsy is a sleep disorder characterized by excessive daytime sleepiness, cataplexy, and a pathological manifestation of rapid eye movement during sleep. Narcoleptic pathogenesis is triggered by both genetic and environmental factors. Recently, development of genome-wide association studies (GWAS) has identified new genetic factors, with many more susceptibility genes yet to be elucidated. Using a new approach that consists of a combination of GWAS and an extensive database search for candidate genes, we picked up 202 candidate genes and performed a replication study in 222 narcoleptic patients and 380 controls. Statistical analysis indicated that six genes, NFATC2, SCP2, CACNA1C, TCRA, POLE, and FAM3D, were associated with narcolepsy (P < 0.001). Some of these associations were further supported by gene expression analyses and an association study in essential hypersomnia (EHS), CNS hypersonia similar to narcolepsy. This novel approach will be applicable to other GWAS in the search of disease-related susceptibility genes.     

Epigenetic Impact of Long-Term Shiftwork: Pilot Evidence From Circadian Genes and Whole-Genome Methylation Analysis

Epigenetic association studies have demonstrated differential promoter methylation in the core circadian genes in breast cancer cases relative to cancer-free controls. The current pilot study aims to investigate whether epigenetic changes affecting breast cancer risk could be caused by circadian disruption through exposure to light at night. Archived DNA samples extracted from whole blood of 117 female subjects from a prospective cohort conducted in Denmark were included in this study. A polymerase chain reaction (PCR)-based method was used for detection of gene-promoter methylation, whereas genome-wide methylation analysis was performed using the Illumina Infinium Methylation Chip. Long-term shiftwork resulted in the same promoter hypomethylation of CLOCK and hypermethylation of CRY2, as was previously observed in breast cancer case-control studies. Genome-wide methylation analysis further discovered widespread methylation alterations in shiftworkers, including changes in many methylation- and cancer-relevant genes. Pathway analysis of the genes with altered methylation patterns revealed several cancer-related pathways. One of the top three networks generated was designated as “DNA replication, recombination, and repair, gene expression, behavior” with ESR1 (estrogen receptor α) featured most prominently in the network, underscoring the potential breast cancer relevance of the genes differentially methylated in long-term shiftworkers. These results, although exploratory, demonstrate the first evidence of the cancer-relevant epigenetic effects of night shiftwork, which warrant further investigation. Considering there are millions of shiftworkers worldwide, understanding the effects of this exposure may lead to novel strategies for cancer prevention and new policies regulating shiftwork. (Author correspondence: yong.zhu@yale.edu)     

GWAS loci associated with Chagas cardiomyopathy influences DNA methylation levels

A recent genome-wide association study (GWAS) identified a locus in chromosome 11 associated with the chronic cardiac form of Chagas disease. Here we aimed to elucidate the potential functional mechanism underlying this genetic association by analyzing the correlation among single nucleotide polymorphisms (SNPs) and DNA methylation (DNAm) levels as cis methylation quantitative trait loci (cis-mQTL) within this region. A total of 2,611 SNPs were tested against 2,647 DNAm sites, in a subset of 37 chronic Chagas cardiomyopathy patients and 20 asymptomatic individuals from the GWAS. We identified 6,958 significant cis-mQTLs (False Discovery Rate [FDR]<0.05) at 1 Mb each side of the GWAS leading variant, where six of them potentially modulate the expression of the SAC3D1 gene, the reported gene in the previous GWAS. In addition, a total of 268 cis-mQTLs showed differential methylation between chronic Chagas cardiomyopathy patients and asymptomatic individuals. The most significant cis-mQTLs mapped in the gene bodies of POLA2 (FDR = 1.04x10^-11), PLAAT3 (FDR = 7.22x10^-03), and CCDC88B (FDR = 1.89x10^-02) that have been associated with cardiovascular and hematological traits in previous studies. One of the most relevant interactions correlated with hypermethylation of CCDC88B. This gene is involved in the inflammatory response, and its methylation and expression levels have been previously reported in Chagas cardiomyopathy. Our findings support the functional relevance of the previously associated genomic region, highlighting the regulation of novel genes that could play a role in the chronic cardiac form of the disease.     

Lung eQTLs to Help Reveal the Molecular Underpinnings of Asthma

Genome-wide association studies (GWAS) have identified loci reproducibly associated with pulmonary diseases; however, the molecular mechanism underlying these associations are largely unknown. The objectives of this study were to discover genetic variants affecting gene expression in human lung tissue, to refine susceptibility loci for asthma identified in GWAS studies, and to use the genetics of gene expression and network analyses to find key molecular drivers of asthma. We performed a genome-wide search for expression quantitative trait loci (eQTL) in 1,111 human lung samples. The lung eQTL dataset was then used to inform asthma genetic studies reported in the literature. The top ranked lung eQTLs were integrated with the GWAS on asthma reported by the GABRIEL consortium to generate a Bayesian gene expression network for discovery of novel molecular pathways underpinning asthma. We detected 17,178 cis- and 593 trans- lung eQTLs, which can be used to explore the functional consequences of loci associated with lung diseases and traits. Some strong eQTLs are also asthma susceptibility loci. For example, rs3859192 on chr17q21 is robustly associated with the mRNA levels of GSDMA (P = 3.55×10⁻¹⁵¹). The genetic-gene expression network identified the SOCS3 pathway as one of the key drivers of asthma. The eQTLs and gene networks identified in this study are powerful tools for elucidating the causal mechanisms underlying pulmonary disease. This data resource offers much-needed support to pinpoint the causal genes and characterize the molecular function of gene variants associated with lung diseases.