Network propagation for GWAS analysis:a practical guide to leveraging molecular networks for disease gene discovery

MOTIVATION: Genome-wide association studies (GWAS) have enabled large-scale analysis of the role of genetic variants in human disease. Despite impressive methodological advances, subsequent clinical interpretation and application remains challenging when GWAS suffer from a lack of statistical power. In recent years, however, the use of information diffusion algorithms with molecular networks has led to fruitful insights on disease genes.RESULTS: We present an overview of the design choices and pitfalls that prove crucial in the application of network propagation methods to GWAS summary statistics. We highlight general trends from the literature, and present benchmark experiments to expand on these insights selecting as case study three diseases and five molecular networks. We verify that the use of gene-level scores based on GWAS P-values offers advantages over the selection of a set of 'seed' disease genes not weighted by the associated P-values if the GWAS summary statistics are of sufficient quality. Beyond that, the size and the density of the networks prove to be important factors for consideration. Finally, we explore several ensemble methods and show that combining multiple networks may improve the network propagation approach.</p

Network propagation for GWAS analysis:a practical guide to leveraging molecular networks for disease gene discovery

MOTIVATION: Genome-wide association studies (GWAS) have enabled large-scale analysis of the role of genetic variants in human disease. Despite impressive methodological advances, subsequent clinical interpretation and application remains challenging when GWAS suffer from a lack of statistical power. In recent years, however, the use of information diffusion algorithms with molecular networks has led to fruitful insights on disease genes.RESULTS: We present an overview of the design choices and pitfalls that prove crucial in the application of network propagation methods to GWAS summary statistics. We highlight general trends from the literature, and present benchmark experiments to expand on these insights selecting as case study three diseases and five molecular networks. We verify that the use of gene-level scores based on GWAS P-values offers advantages over the selection of a set of 'seed' disease genes not weighted by the associated P-values if the GWAS summary statistics are of sufficient quality. Beyond that, the size and the density of the networks prove to be important factors for consideration. Finally, we explore several ensemble methods and show that combining multiple networks may improve the network propagation approach.</p

Using eQTL weights to improve power for genome-wide association studies: a genetic study of childhood asthma

Increasing evidence suggests that single nucleotide polymorphisms (SNPs) associated with complex traits are more likely to be expression quantitative trait loci (eQTLs). Incorporating eQTL information hence has potential to increase power of genome-wide association studies (GWAS). In this paper, we propose using eQTL weights as prior information in SNP based association tests to improve test power while maintaining control of the family-wise error rate (FWER) or the false discovery rate (FDR). We apply the proposed methods to the analysis of a GWAS for childhood asthma consisting of 1296 unrelated individuals with German ancestry. The results confirm that eQTLs are enriched for previously reported asthma SNPs. We also find that some SNPs are insignificant using procedures without eQTL weighting, but become significant using eQTL-weighted Bonferroni or Benjamini–Hochberg procedures, while controlling the same FWER or FDR level. Some of these SNPs have been reported by independent studies in recent literature. The results suggest that the eQTL-weighted procedures provide a promising approach for improving power of GWAS. We also report the results of our methods applied to the large-scale European GABRIEL consortium data

EGEA Collection: A Biobank Devoted to Asthma and Asthma-related Phenotypes

The Epidemiological study on the Genetics and Environment of Asthma (EGEA) is a 20-year cohort including cases with asthma (n = 388), their first degree relatives (n = 1244) and population-based controls (n = 415) recruited in the early 90’s from five French cities. Participants were extensively characterized regarding environment and respiratory health, and a specific effort was made to trace and bank biological samples (ISO 9001 certification since 2006). Access to data and samples is opened to researchers wishing to develop new scientific collaborative programs. The survey has already led to more than a hundred papers with almost one third in collaboration with national and international teams.</p

Identification of OCA2 as a novel locus for the co-morbidity of asthma-plus-eczema

Background: Numerous genes have been associated with the three most common allergic diseases (asthma, allergic rhinitis or eczema) but these genes explain only a part of the heritability. In the vast majority of genetic studies, complex phenotypes such as co- morbidity of two of these diseases, have not been considered. This may partly explain missing heritability. Objective: To identify genetic variants specifically associated with the co-morbidity of asthma-plus-eczema. Methods: We first conducted a meta-analysis of four GWAS (Genome-Wide Association Study) of the combined asthma-plus-eczema phenotype (total of 8807 European-ancestry subjects of whom 1208 subjects had both asthma and eczema). To assess whether the association with SNP(s) was specific to the co- morbidity, we also conducted a meta-analysis of homogeneity test of association according to disease status (“asthma-plus-eczema” vs. the presence of only one disease “asthma only or eczema only”). We then used a joint test by combining the two test statistics from the co-morbidity-SNP association and the phenotypic heterogeneity of SNP effect meta-analyses. Results: Seven SNPs were detected for specific association to the asthma-plus-eczema co-morbidity, two with significant and five with suggestive evidence using the joint test after correction for multiple testing. The two significant SNPs are located in the OCA2 gene (Oculocutaneous Albinism II), a new locus never detected for significant evidence of association with any allergic disease. This gene is a promising candidate gene, because of its link to skin and lung diseases, and to epithelial barrier and immune mechanisms. Conclusion: Our study underlines the importance of studying sub-phenotypes as co-morbidities to detect new susceptibility genes

Analysis of GWAS-nominated loci for lung cancer and COPD revealed a new asthma locus

Background: Asthma, lung cancer (LC) and chronic obstructive pulmonary disease (COPD) are three respiratory diseases characterized by complex mechanisms underlying and genetic predispositions, with asthma having the highest calculated heritability. Despite efforts deployed in the last decades, only a small part of its heritability has been elucidated. It was hypothesized that shared genetic factors by these three diseases could help identify new asthma loci. Methods: GWAS-nominated LC and COPD loci were selected among studies performed in Caucasian cohorts using the GWAS Catalog. Genetic analyses were carried out for these loci in the Saguenay–Lac-Saint-Jean (SLSJ) asthma familial cohort and then replicated in two independent cohorts (the Canadian Cohort Obstructive Lung Disease [CanCOLD] and the Epidemiological Study of the Genetics and Environment of Asthma [EGEA]). Results: Analyses in the SLSJ cohort identified 2851 and 4702 genetic variants to be replicated in the CanCOLD and EGEA cohorts for LC and COPD loci respectively. Replication and meta-analyses allowed the association of one new locus with asthma, 2p24.3, from COPD studies. None was associated from LC studies reported. Conclusions: The approach used in this study contributed to better understand the heritability of asthma with shared genetic backgrounds of respiratory diseases

Genome-wide interaction study of early-life smoking exposure on time-to-asthma onset in childhood

Background: Asthma, a heterogeneous disease with variable age of onset, results from the interplay between genetic and environmental factors. Early-life tobacco smoke (ELTS) exposure is a major asthma risk factor. Only a few genetic loci have been reported to interact with ELTS exposure in asthma. Objective: Our aim was to identify new loci interacting with ELTS exposure on time-to-asthma onset (TAO) in childhood.Methods: We conducted genome-wide interaction analyses of ELTS exposure on time-to-asthma onset in childhood in five European-ancestry studies (totaling 8,273 subjects) using Cox proportional-hazard model. The results of all five genome-wide analyses were meta-analyzed.Results: The 13q21 locus showed genome-wide significant interaction with ELTS exposure (P=4.3x10-8 for rs7334050 within KLHL1 with consistent results across the five studies). Suggestive interactions (P&lt;5x10-6) were found at three other loci: 20p12 (rs13037508 within MACROD2; P=4.9x10-7), 14q22 (rs7493885 near NIN; P=2.9x10-6) and 2p22 (rs232542 near CYP1B1; P=4.1x10-6). Functional annotations and the literature showed that the lead SNPs at these four loci influence DNA methylation in the blood and are located nearby CpG sites reported to be associated with exposure to tobacco smoke components, which strongly support our findings.Conclusion and Clinical Relevance: We identified novel candidate genes interacting with ELTS exposure on time-to-asthma onset in childhood. These genes have plausible biological relevance related to tobacco smoke exposure. Further epigenetic and functional studies are needed to confirm these findings and to shed light on the underlying mechanisms

Identification of a new locus at 16q12 associated with time-to-asthma onset

International audienceBackground: Asthma is a heterogeneous disease in which age-of-onset plays an important role.Objective: We sought to identify the genetic variants associated with time-to-asthma onset.Methods: We conducted a large-scale meta-analysis of nine genome-wide association studies of time-to-asthma onset (total of 5,462 asthmatics with a broad range of age-of-asthma onset and 8,424 controls of European ancestry) performed using survival analysis techniques.Results: We detected five regions associated with time-to-asthma onset at genome-wide significant level (P<5x10-8). We evidenced a new locus in 16q12 region (near cylindromatosis turban tumor syndrome gene (CYLD)) and confirmed four asthma risk regions: 2q12 (IL1RL1), 6p21 (HLA-DQA1), 9p24 (IL33) and 17q12-q21 (ZPBP2-GSDMA). Conditional analyses identified two distinct signals at 9p24 (both upstream of IL33) and at 17q12-q21 (near ZPBP2 and within GSDMA). These seven distinct loci explained together 6.0% of the variance in time-to-asthma onset. In addition, we showed that genetic variants at 9p24 and 17q12-q21 were strongly associated with an earlier onset of childhood asthma (P≤0.002) whereas 16q12 SNP was associated with a later asthma onset (P=0.04). A high burden of disease risk alleles at these loci was associated with earlier age-of-asthma onset (4 years versus 9-12 years, P=10-4).Conclusion: The new susceptibility region for time-to-asthma onset at 16q12 harbors variants that correlate with the expression of CYLD and NOD2 (nucleotide-binding oligomerization domain 2), two strong candidates for asthma. This study demonstrates that incorporating the variability of age-of-asthma onset in asthma modeling is a helpful approach in the search for disease susceptibility genes

Replication of Association between ADAM33 Polymorphisms and Psoriasis

Polymorphisms in ADAM33, the first gene identified in asthma by positional cloning, have been recently associated with psoriasis. No replication study of this association has been published so far. Data available in the French EGEA study (Epidemiological study on Genetics and Environment of Asthma, bronchial hyperresponsivensess and Atopy) give the opportunity to attempt to replicate the association between ADAM33 and psoriasis in 2002 individuals. Psoriasis (n = 150) has been assessed by questionnaire administered by an interviewer and a sub-sample of subjects with early-onset psoriasis (n = 74) has been identified based on the age of the subjects at time of interview (<40 years). Nine SNPs in ADAM33 and 11 SNPs in PSORS1 were genotyped. Association analysis was conducted by using two methods, GEE regression-based method and a likelihood-based method (LAMP program). The rs512625 SNP in ADAM33 was found associated with psoriasis at p = 0.01, the usual threshold required for replication (OR [95% CI] for heterozygotes compared to the reference group of homozygotes for the most frequent allele = 0.61 [0.42;0.89]). The rs628977 SNP, which was not in linkage disequilibrium with rs512625, was significantly associated with early-onset psoriasis (p = 0.01, OR [95% CI] for homozygotes for the minor allele compared to the reference group = 2.52 [1.31;4.86]). Adjustment for age, sex, asthma and a PSORS1 SNP associated with psoriasis in the EGEA data did not change the significance of these associations. This suggests independent effects of ADAM33 and PSORS1 on psoriasis. This is the first study that replicates an association between genetic variants in ADAM33 and psoriasis. Interestingly, the 2 ADAM33 SNPs associated with psoriasis in the present analysis were part of the 3-SNPs haplotypes showing the strongest associations in the initial study. The identification of a pleiotropic effect of ADAM33 on asthma and psoriasis may contribute to the understanding of these common immune-mediated diseases