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.

Identification of novel genes influencing eosinophil-specific protein levels in asthma families.

BACKGROUND: Eosinophils play a key role in the asthma allergic response by releasing cytotoxic molecules such as eosinophil cationic protein (ECP) and eosinophil-derived neurotoxin (EDN) that generate epithelium damages. OBJECTIVE: We sought to identify genetic variants influencing ECP and EDN levels in asthma-ascertained families. METHODS: We performed univariate and bivariate genome-wide association analyses of ECP and EDN levels in 1018 subjects from the EGEA study with follow-up in 153 subjects from the Saguenay-Lac-Saint-Jean study and combined the results of these 2 studies through meta-analysis. We then conducted Bayesian statistical fine mapping together with quantitative trait locus and functional annotation analyses to identify the most likely functional genetic variants and candidate genes. RESULTS: We identified 5 genome-wide significant loci (P < 5 × 10^-8) including 7 distinct signals associated with ECP and/or EDN levels. The genes targeted by our fine mapping and functional search include RNASE2 and RNASE3 (14q11), which encode EDN and ECP, respectively, and 4 other genes that regulate ECP and EDN levels. These 4 genes were JAK1 (1p31), a transcription factor that plays a key role in the immune response and acts as a potential therapeutic target for eosinophilic asthma; ARHGAP25 (2p13), which is involved in leukocyte recruitment to inflammatory sites; NDUFA4 (7p21), which encodes a component of the mitochondrial respiratory chain and is involved in cellular response to stress; and CTSL (9q22), which is involved in immune response, extracellular remodeling, and allergic inflammation. CONCLUSION: Analysis of specific phenotypes produced by eosinophils allows the identification of genes that play a major role in allergic response and inflammation, and offers potential therapeutic targets for asthma.

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.

Associations between specific IgE sensitization to 26 respiratory allergen molecules and HLA class II alleles in the EGEA cohort.

BACKGROUND: Allergy, the most frequent immune disorder affecting 30% of the world's population, is the consequence of immunoglobin E (IgE) sensitization to allergens. Among the genetic factors suspected to be involved in allergy, the HLA class-II genomic region is a strong candidate. OBJECTIVE: To assess the association between HLA class-II alleles and specific IgE (sIgE) sensitization to a large number of respiratory allergen molecules. METHODS: The analysis relied on 927 participants of the EGEA cohort, including 497 asthmatics. The study focuses on 26 aeroallergens recognized by sIgE in at least 5% of the study population (determined with the MEDALL chip with sIgE ≥ 0.3 ISU) and 23 imputed HLA class-II alleles. For each sIgE sensitization and HLA class-II allele, we fitted a logistic regression model accounting for familial dependence and adjusted for gender, age, and genetic principal components. p-values were corrected for multiple comparisons (False Discovery Rate). RESULTS: Most of the 19 statistically significant associations observed regard pollen allergens (mugwort Art v 1, olive tree Ole e 1, timothy grass Phl p 2, Phl p 5 and plantain Pla l 1), three were mold allergen (Alternaria Alt a 1), and a single one regards house dust mite allergen (Der p 7). No association was observed with pet allergens. The strongest associations were found with mugwort Art v 1 (OR = 5.42 (95%CI, 3.30; 8.88), 4.14 (2.65; 6.47), 3.16 (1.88; 5.31) with DQB1*05:01, DQA1*01:01 and DRB1*01:01, respectively). CONCLUSION: Our results support the important role of HLA class-II alleles as immune response genes predisposing their carriers for sensitization to various major pollen allergens.

Interactive effect between ATPase-related genes and early-life tobacco smoke exposure on bronchial hyper-responsiveness detected in asthma-ascertained families.

BACKGROUND: A positional cloning study of bronchial hyper-responsiveness (BHR) at the 17p11 locus in the French Epidemiological study on the Genetics and Environment of Asthma (EGEA) families showed significant interaction between early-life environmental tobacco smoke (ETS) exposure and genetic variants located in DNAH9. This gene encodes the heavy chain subunit of axonemal dynein, which is involved with ATP in the motile cilia function.Our goal was to identify genetic variants at other genes interacting with ETS in BHR by investigating all genes belonging to the 'ATP-binding' and 'ATPase activity' pathways which include DNAH9, are targets of cigarette smoke and play a crucial role in the airway inflammation. METHODS: Family-based interaction tests between ETS-exposed and unexposed BHR siblings were conducted in 388 EGEA families. Twenty single-nucleotide polymorphisms (SNP) showing interaction signals (p≤5.10-3) were tested in the 253 Saguenay-Lac-Saint-Jean (SLSJ) families. RESULTS: One of these SNPs was significantly replicated for interaction with ETS in SLSJ families (p=0.003). Another SNP reached the significance threshold after correction for multiple testing in the combined analysis of the two samples (p=10-5). Results were confirmed using both a robust log-linear test and a gene-based interaction test. CONCLUSION: The SNPs showing interaction with ETS belong to the ATP8A1 and ABCA1 genes, which play a role in the maintenance of asymmetry and homeostasis of lung membrane lipids.

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 (totalling 8273 subjects) using Cox proportional-hazard model. The results of all five genome-wide analyses were meta-analysed. RESULTS: The 13q21 locus showed genome-wide significant interaction with ELTS exposure (P = 4.3 × 10-8 for rs7334050 within KLHL1 with consistent results across the five studies). Suggestive interactions (P < 5 × 10-6 ) were found at three other loci: 20p12 (rs13037508 within MACROD2; P = 4.9 × 10-7 ), 14q22 (rs7493885 near NIN; P = 2.9 × 10-6 ) and 2p22 (rs232542 near CYP1B1; P = 4.1 × 10-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. CONCLUSIONS 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.

A novel role for ciliary function in atopy: ADGRV1 and DNAH5 interactions.

BACKGROUND: Atopy, an endotype underlying allergic diseases, has a substantial genetic component. OBJECTIVE: Our goal was to identify novel genes associated with atopy in asthma-ascertained families. METHODS: We implemented a 3-step analysis strategy in 3 data sets: the Epidemiological Study on the Genetics and Environment of Asthma (EGEA) data set (1660 subjects), the Saguenay-Lac-Saint-Jean study data set (1138 subjects), and the Medical Research Council (MRC) data set (446 subjects). This strategy included a single nucleotide polymorphism (SNP) genome-wide association study (GWAS), the selection of related gene pairs based on statistical filtering of GWAS results, and text-mining filtering using Gene Relationships Across Implicated Loci and SNP-SNP interaction analysis of selected gene pairs. RESULTS: We identified the 5q14 locus, harboring the adhesion G protein-coupled receptor V1 (ADGRV1) gene, which showed genome-wide significant association with atopy (rs4916831, meta-analysis P value = 6.8 × 10-9). Statistical filtering of GWAS results followed by text-mining filtering revealed relationships between ADGRV1 and 3 genes showing suggestive association with atopy (P ≤ 10-4). SNP-SNP interaction analysis between ADGRV1 and these 3 genes showed significant interaction between ADGRV1 rs17554723 and 2 correlated SNPs (rs2134256 and rs1354187) within the dynein axonemal heavy chain 5 (DNAH5) gene (Pmeta-int = 3.6 × 10-5 and 6.1 × 10-5, which met the multiple-testing corrected threshold of 7.3 × 10-5). Further conditional analysis indicated that rs2134256 alone accounted for the interaction signal with rs17554723. CONCLUSION: Because both DNAH5 and ADGRV1 contribute to ciliary function, this study suggests that ciliary dysfunction might represent a novel mechanism underlying atopy. Combining GWAS and epistasis analysis driven by statistical and knowledge-based evidence represents a promising approach for identifying new genes involved in complex traits.

SigMod: an exact and efficient method to identify a strongly interconnected disease-associated module in a gene network.

MOTIVATION: Apart from single marker-based tests classically used in genome-wide association studies (GWAS), network-assisted analysis has become a promising approach to identify a set of genes associated with disease. To date, most network-assisted methods aim at finding genes connected in a background network, whatever the density or strength of their connections. This can hamper the findings as sparse connections are non-robust against noise from either the GWAS results or the network resource. RESULTS: We present SigMod, a novel and efficient method integrating GWAS results and gene network to identify a strongly interconnected gene module enriched in high association signals. Our method is formulated as a binary quadratic optimization problem, which can be solved exactly through graph min-cut algorithms. Compared to existing methods, SigMod has several desirable properties: (i) edge weights quantifying confidence of connections between genes are taken into account, (ii) the selection path can be computed rapidly, (iii) the identified gene module is strongly interconnected, hence includes genes of high functional relevance, and (iv) the method is robust against noise from either the GWAS results or the network resource. We applied SigMod to both simulated and real data. It was found to outperform state-of-the-art network-assisted methods in identifying disease-associated genes. When SigMod was applied to childhood-onset asthma GWAS results, it successfully identified a gene module enriched in consistently high association signals and made of functionally related genes that are biologically relevant for asthma. AVAILABILITY AND IMPLEMENTATION: An R package SigMod is available at: https://github.com/YuanlongLiu/SigMod. CONTACT: yuanlong.liu@inserm.fr. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Vitamin D levels and susceptibility to asthma, elevated immunoglobulin E levels, and atopic dermatitis: A Mendelian randomization study.

BACKGROUND: Low circulating vitamin D levels have been associated with risk of asthma, atopic dermatitis, and elevated total immunoglobulin E (IgE). These epidemiological associations, if true, would have public health importance, since vitamin D insufficiency is common and correctable. METHODS AND FINDINGS: We aimed to test whether genetically lowered vitamin D levels were associated with risk of asthma, atopic dermatitis, or elevated serum IgE levels, using Mendelian randomization (MR) methodology to control bias owing to confounding and reverse causation. The study employed data from the UK Biobank resource and from the SUNLIGHT, GABRIEL and EAGLE eczema consortia. Using four single-nucleotide polymorphisms (SNPs) strongly associated with 25-hydroxyvitamin D (25OHD) levels in 33,996 individuals, we conducted MR studies to estimate the effect of lowered 25OHD on the risk of asthma (n = 146,761), childhood onset asthma (n = 15,008), atopic dermatitis (n = 40,835), and elevated IgE level (n = 12,853) and tested MR assumptions in sensitivity analyses. None of the four 25OHD-lowering alleles were associated with asthma, atopic dermatitis, or elevated IgE levels (p ≥ 0.2). The MR odds ratio per standard deviation decrease in log-transformed 25OHD was 1.03 (95% confidence interval [CI] 0.90-1.19, p = 0.63) for asthma, 0.95 (95% CI 0.69-1.31, p = 0.76) for childhood-onset asthma, and 1.12 (95% CI 0.92-1.37, p = 0.27) for atopic dermatitis, and the effect size on log-transformed IgE levels was -0.40 (95% CI -1.65 to 0.85, p = 0.54). These results persisted in sensitivity analyses assessing population stratification and pleiotropy and vitamin D synthesis and metabolism pathways. The main limitations of this study are that the findings do not exclude an association between the studied outcomes and 1,25-dihydoxyvitamin D, the active form of vitamin D, the study was underpowered to detect effects smaller than an OR of 1.33 for childhood asthma, and the analyses were restricted to white populations of European ancestry. This research has been conducted using the UK Biobank Resource and data from the SUNLIGHT, GABRIEL and EAGLE Eczema consortia. CONCLUSIONS: In this study, we found no evidence that genetically determined reduction in 25OHD levels conferred an increased risk of asthma, atopic dermatitis, or elevated total serum IgE, suggesting that efforts to increase vitamin D are unlikely to reduce risks of atopic disease.

A common 16p11.2 inversion underlies the joint susceptibility to asthma and obesity.

The prevalence of asthma and obesity is increasing worldwide, and obesity is a well-documented risk factor for asthma. The mechanisms underlying this association and parallel time trends remain largely unknown but genetic factors may be involved. Here, we report on a common ~0.45 Mb genomic inversion at 16p11.2 that can be accurately genotyped via SNP array data. We show that the inversion allele protects against the joint occurrence of asthma and obesity in five large independent studies (combined sample size of 317 cases and 543 controls drawn from a total of 5,809 samples; combined OR = 0.48, p = 5.5 × 10(-6)). Allele frequencies show remarkable worldwide population stratification, ranging from 10% in East Africa to 49% in Northern Europe, consistent with discordant and extreme genetic drifts or adaptive selections after human migration out of Africa. Inversion alleles strongly correlate with expression levels of neighboring genes, especially TUFM (p = 3.0 × 10(-40)) that encodes a mitochondrial protein regulator of energy balance and inhibitor of type 1 interferon, and other candidates for asthma (IL27) and obesity (APOB48R and SH2B1). Therefore, by affecting gene expression, the ~0.45 Mb 16p11.2 inversion provides a genetic basis for the joint susceptibility to asthma and obesity, with a population attributable risk of 39.7%. Differential mitochondrial function and basal energy balance of inversion alleles might also underlie the potential selection signature that led to their uneven distribution in world populations.

DNA methylation within melatonin receptor 1A (MTNR1A) mediates paternally transmitted genetic variant effect on asthma plus rhinitis.

BACKGROUND: Asthma and allergic rhinitis (AR) are common allergic comorbidities with a strong genetic component in which epigenetic mechanisms might be involved. OBJECTIVE: We aimed to identify novel risk loci for asthma and AR while accounting for parent-of-origin effect. METHODS: We performed a series of genetic analyses, taking into account the parent-of-origin effect in families ascertained through asthma: (1) genome-wide linkage scan of asthma and AR in 615 European families, (2) association analysis with 1233 single nucleotide polymorphisms (SNPs) covering the significant linkage region in 162 French Epidemiological Study on the Genetics and Environment of Asthma families with replication in 154 Canadian Saguenay-Lac-Saint-Jean asthma study families, and (3) association analysis of disease and significant SNPs with DNA methylation (DNAm) at CpG sites in 40 Saguenay-Lac-Saint-Jean asthma study families. RESULTS: We detected a significant paternal linkage of the 4q35 region to asthma and allergic rhinitis comorbidity (AAR; P = 7.2 × 10(-5)). Association analysis in this region showed strong evidence for the effect of the paternally inherited G allele of rs10009104 on AAR (P = 1.1 × 10(-5), reaching the multiple-testing corrected threshold). This paternally inherited allele was also significantly associated with DNAm levels at the cg02303933 site (P = 1.7 × 10(-4)). Differential DNAm at this site was found to mediate the identified SNP-AAR association. CONCLUSION: By integrating genetic and epigenetic data, we identified that a differentially methylated CpG site within the melatonin receptor 1A (MTNR1A) gene mediates the effect of a paternally transmitted genetic variant on the comorbidity of asthma and AR. This study provides a novel insight into the role of epigenetic mechanisms in patients with allergic respiratory diseases.

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

BACKGROUND: 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 (TAO). METHODS: We conducted a large-scale meta-analysis of 9 genome-wide association studies of TAO (total of 5462 asthmatic patients with a broad range of age of asthma onset and 8424 control subjects of European ancestry) performed by using survival analysis techniques. RESULTS: We detected 5 regions associated with TAO at the genome-wide significant level (P < 5 × 10-8). We evidenced a new locus in the 16q12 region (near cylindromatosis turban tumor syndrome gene [CYLD]) and confirmed 4 asthma risk regions: 2q12 (IL-1 receptor-like 1 [IL1RL1]), 6p21 (HLA-DQA1), 9p24 (IL33), and 17q12-q21 (zona pellucida binding protein 2 [ZPBP2]-gasdermin A [GSDMA]). Conditional analyses identified 2 distinct signals at 9p24 (both upstream of IL33) and 17q12-q21 (near ZPBP2 and within GSDMA). Together, these 7 distinct loci explained 6.0% of the variance in TAO. In addition, we showed that genetic variants at 9p24 and 17q12-q21 were strongly associated with an earlier onset of childhood asthma (P ≤ .002), whereas the 16q12 single nucleotide polymorphism was associated with later asthma onset (P = .04). A high burden of disease risk alleles at these loci was associated with earlier age of asthma onset (4 vs 9-12 years, P = 10-4). CONCLUSION: The new susceptibility region for TAO at 16q12 harbors variants that correlate with the expression of CYLD and nucleotide-binding oligomerization domain 2 (NOD2), 2 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.

Interaction between the DNAH9 gene and early smoke exposure in bronchial hyperresponsiveness.

A previous genome-wide linkage scan of bronchial hyperresponsiveness (BHR) in the French Epidemiological study on the Genetics and Environment of Asthma (EGEA) families, performed in the presence of a gene×early-life environmental tobacco smoke (ETS) exposure interaction, showed the strongest interaction in the 17p11 region where linkage was detected only among unexposed siblings. Our goal was to conduct fine-scale mapping of 17p11 to identify single nucleotide polymorphisms (SNPs) interacting with ETS that influence BHR.Analyses were performed in 388 French EGEA asthmatic families, using a two-step strategy: 1) selection of SNPs displaying family-based association test (FBAT) association signals (p≤0.01) with BHR in unexposed siblings, and 2) a FBAT homogeneity test between exposed and unexposed siblings plus a robust log-linear interaction test.A single SNP reached the threshold (p≤3×10(-3)) for significant interaction with ETS using both interaction tests, after accounting for multiple testing. Results were replicated in 253 French-Canadian families, but not in 341 UK families, probably due in part to differences in phenotypic features between datasets.The SNP showing significant interaction with ETS belongs toDNAH9(dynein, axonemal, heavy chain 9), a promising candidate gene involved in respiratory cilia mobility and associated with primary ciliary dyskinesia, a disease associated with abnormalities of pulmonary function.

The nuclear factor I/A (NFIA) gene is associated with the asthma plus rhinitis phenotype.

BACKGROUND: A previous genome-wide linkage scan in 295 families of the French Epidemiological Study on the Genetics and Environment of Asthma (EGEA) showed strong evidence of linkage of the 1p31 region to the combined asthma plus allergic rhinitis (AR) phenotype. OBJECTIVE: Our purpose was to conduct fine-scale mapping of the 1p31 linkage region to identify the genetic variants associated with asthma plus AR. METHODS: Association analyses with the asthma plus rhinitis phenotype were first conducted in the EGEA family sample using the family-based association method (FBAT) and logistic regression. The test of homogeneity of association between asthma plus AR versus asthma alone or AR alone was also applied. Replication of EGEA findings was sought in French-Canadian and United Kingdom family samples. RESULTS: We found a significant association between asthma plus rhinitis and a 1p31 genetic variant (P = 2 × 10(-5) for rs12122228, which reached the multiple testing-corrected threshold) in EGEA using FBAT. There was evidence of heterogeneity of association between asthma plus AR versus asthma alone or AR alone (P = .03). A Meta-analysis of FBAT results from EGEA and French-Canadian families improved evidence for both association and heterogeneity (P = 5 × 10(-6) and P = .008, respectively), whereas a meta-analysis of EGEA, French-Canadian, and United Kingdom samples based on logistic regression slightly increased the evidence for heterogeneity. CONCLUSION: The single nucleotide polymorphism specifically associated to asthma plus rhinitis is located in the flanking 5' untranslated region of the nuclear factor I/A (NFIA) gene, a strong candidate gene for asthma and AR.

Fraction of exhaled nitric oxide values in childhood are associated with 17q11.2-q12 and 17q12-q21 variants.

BACKGROUND: The fraction of exhaled nitric oxide (Feno) value is a biomarker of eosinophilic airway inflammation and is associated with childhood asthma. Identification of common genetic variants associated with childhood Feno values might help to define biological mechanisms related to specific asthma phenotypes. OBJECTIVE: We sought to identify the genetic variants associated with childhood Feno values and their relation with asthma. METHODS: Feno values were measured in children age 5 to 15 years. In 14 genome-wide association studies (N = 8,858), we examined the associations of approximately 2.5 million single nucleotide polymorphisms (SNPs) with Feno values. Subsequently, we assessed whether significant SNPs were expression quantitative trait loci in genome-wide expression data sets of lymphoblastoid cell lines (n = 1,830) and were related to asthma in a previously published genome-wide association data set (cases, n = 10,365; control subjects: n = 16,110). RESULTS: We identified 3 SNPs associated with Feno values: rs3751972 in LYR motif containing 9 (LYRM9; P = 1.97 × 10(-10)) and rs944722 in inducible nitric oxide synthase 2 (NOS2; P = 1.28 × 10(-9)), both of which are located at 17q11.2-q12, and rs8069176 near gasdermin B (GSDMB; P = 1.88 × 10(-8)) at 17q12-q21. We found a cis expression quantitative trait locus for the transcript soluble galactoside-binding lectin 9 (LGALS9) that is in linkage disequilibrium with rs944722. rs8069176 was associated with GSDMB and ORM1-like 3 (ORMDL3) expression. rs8069176 at 17q12-q21, but not rs3751972 and rs944722 at 17q11.2-q12, were associated with physician-diagnosed asthma. CONCLUSION: This study identified 3 variants associated with Feno values, explaining 0.95% of the variance. Identification of functional SNPs and haplotypes in these regions might provide novel insight into the regulation of Feno values. This study highlights that both shared and distinct genetic factors affect Feno values and childhood asthma.

Genetic heterogeneity of asthma phenotypes identified by a clustering approach.

The aim of the study was to identify genetic variants associated with refined asthma phenotypes enabling multiple features of the disease to be taken into account. Latent class analysis (LCA) was applied in 3001 adults ever having asthma recruited in the frame of three epidemiological surveys (the European Community Respiratory Health Survey (ECRHS), the Swiss Study on Air Pollution and Lung Disease in Adults (SAPALDIA) and the Epidemiological Study on the Genetics and Environment of Asthma (EGEA)). 14 personal and phenotypic characteristics, gathered from questionnaires and clinical examination, were used. A genome-wide association study was conducted for each LCA-derived asthma phenotype, compared to subjects without asthma (n=3474). The LCA identified four adult asthma phenotypes, mainly characterised by disease activity, age of asthma onset and atopic status. Associations of genome-wide significance (p<1.25 × 10(-7)) were observed between "active adult-onset nonallergic asthma" and rs9851461 flanking CD200 (3q13.2) and between "inactive/mild nonallergic asthma" and rs2579931 flanking GRIK2 (6q16.3). Borderline significant results (2.5 × 10(-7) < p <8.2 × 10(-7)) were observed between three single nucleotide polymorphisms (SNPs) in the ALCAM region (3q13.11) and "active adult-onset nonallergic asthma". These results were consistent across studies. 15 SNPs identified in previous genome-wide association studies of asthma have been replicated with at least one asthma phenotype, most of them with the "active allergic asthma" phenotype. Our results provide evidence that a better understanding of asthma phenotypic heterogeneity helps to disentangle the genetic heterogeneity of asthma.

Human leukocyte antigen class II variants and adult-onset asthma: does occupational allergen exposure play a role?

Recently, a locus centred on rs9273349 in the HLA-DQ region emerged from genome-wide association studies of adult-onset asthma. We aimed to further investigate the role of human leukocyte antigen (HLA) class II in adult-onset asthma and a possible interaction with occupational exposures. We imputed classical HLA-II alleles from 7579 single-nucleotide polymorphisms in 6025 subjects (1202 with adult-onset asthma) from European cohorts: ECRHS, SAPALDIA, EGEA and B58C, and from surveys of bakers and agricultural workers. Based on an asthma-specific job-exposure matrix, 2629 subjects had ever been exposed to high molecular weight (HMW) allergens. We explored associations between 23 common HLA-II alleles and adult-onset asthma, and tested for gene-environment interaction with occupational exposure to HMW allergens. Interaction was also tested for rs9273349. Marginal associations of classical HLA-II alleles and adult-onset asthma were not statistically significant. Interaction was detected between the DPB1*03:01 allele and exposure to HMW allergens (p = 0.009), in particular to latex (p = 0.01). In the unexposed group, the DPB1*03:01 allele was associated with adult-onset asthma (OR 0.67, 95%CI 0.53-0.86). HMW allergen exposures did not modify the association of rs9273349 with adult-onset asthma. Common classical HLA-II alleles were not marginally associated with adult-onset asthma. The association of latex exposure and adult-onset asthma may be modified by DPB1*03:01.

Determinants of immunoglobulin G responses to respiratory syncytial virus and rhinovirus in children and adults.

Introduction: Exposure to respiratory viruses is a significant cause of morbidity and affects virus-specific antibody levels. Little is known about determinants associated with immune response to these viruses. We aimed to investigate the determinants of respiratory syncytial virus (RSV)- and rhinovirus (RV)- specific IgG responses in both children and adults. Methods: The study is based on the EGEA cohort, composed of 530 samples of children in EGEA1 (1991-95) and 1241 samples of adults in EGEA2 (2003-07). Cumulative RV-specific IgG levels (species A, B and C) and IgG levels to RSV-G protein were measured by using micro-array technoloy. Multiple linear mixed models (random effect to account for familial dependence) were performed to assess associations between age, sex, body mass index (BMI), tobacco smoke exposure and season of blood sampling with RSV-and RV-specific IgG levels. Results: In children (11.1 ± 2.8 years old, 57% boys), higher RV-specific IgG levels were associated with older age (only for RV-B), female sex and lower BMI, while only older age was associated with higher RSV-specific IgG levels. In adults (43.5 ± 16.7 years old, 48% men), younger age, female sex, lower BMI, active smoking and all seasons except summer were associated with higher RV-specific IgG levels. Older age, active smoking and all seasons except summer were associated with higher RSV-specific IgG levels. Conclusion: Personal and seasonal determinants of RSV- and RV-specific IgG levels seem to vary according to the respiratory virus type and between children and adults, suggesting different patterns of responses along the life course.

Trait selection strategy in multi-trait GWAS: Boosting SNP discoverability.

Since the first genome-wide association studies (GWASs), thousands of variant-trait associations have been discovered. However, comprehensively mapping the genetic determinant of complex traits through univariate testing can require prohibitive sample sizes. Multi-trait GWAS can circumvent this issue and improve statistical power by leveraging the joint genetic architecture of human phenotypes. Although many methodological hurdles of multi-trait testing have been solved, the strategy to select traits has been overlooked. In this study, we conducted multi-trait GWAS on approximately 20,000 combinations of 72 traits using an omnibus test as implemented in the Joint Analysis of Summary Statistics. We assessed which genetic features of the sets of traits analyzed were associated with an increased detection of variants compared with univariate screening. Several features of the set of traits, including the heritability, the number of traits, and the genetic correlation, drive the multi-trait test gain. Using these features jointly in predictive models captures a large fraction of the power gain of the multi-trait test (Pearson's r between the observed and predicted gain equals 0.43, p < 1.6 × 10-60). Applying an alternative multi-trait approach (Multi-Trait Analysis of GWAS), we identified similar features of interest, but with an overall 70% lower number of new associations. Finally, selecting sets based on our data-driven models systematically outperformed the common strategy of selecting clinically similar traits. This work provides a unique picture of the determinant of multi-trait GWAS statistical power and outlines practical strategies for multi-trait testing.