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.

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.

Estimating CDKN2A mutation carrier probability among global familial melanoma cases using GenoMELPREDICT.

BACKGROUND: Although rare in the general population, highly penetrant germline mutations in CDKN2A are responsible for 5%-40% of melanoma cases reported in melanoma-prone families. We sought to determine whether MELPREDICT was generalizable to a global series of families with melanoma and whether performance improvements can be achieved. METHODS: In total, 2116 familial melanoma cases were ascertained by the international GenoMEL Consortium. We recapitulated the MELPREDICT model within our data (GenoMELPREDICT) to assess performance improvements by adding phenotypic risk factors and history of pancreatic cancer. We report areas under the curve (AUC) with 95% confidence intervals (CIs) along with net reclassification indices (NRIs) as performance metrics. RESULTS: MELPREDICT performed well (AUC 0.752, 95% CI 0.730-0.775), and GenoMELPREDICT performance was similar (AUC 0.748, 95% CI 0.726-0.771). Adding a reported history of pancreatic cancer yielded discriminatory improvement (P < .0001) in GenoMELPREDICT (AUC 0.772, 95% CI 0.750-0.793, NRI 0.40). Including phenotypic risk factors did not improve performance. CONCLUSION: The MELPREDICT model functioned well in a global data set of familial melanoma cases. Adding pancreatic cancer history improved model prediction. GenoMELPREDICT is a simple tool for predicting CDKN2A mutational status among melanoma patients from melanoma-prone families and can aid in directing these patients to receive genetic testing or cancer risk counseling.

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.

Gain of power of the general regression model compared to Cochran-Armitage Trend tests: simulation study and application to bipolar disorder.

BACKGROUND: Most genome-wide association studies assumed an additive model of inheritance which may result in significant loss of power when there is a strong departure from additivity. The General Regression Model (GRM), which allows performing an assumption-free test for association by testing for both additive effect and deviation from additive effect, may be more appropriate for association tests. Additionally, GRM allows testing the underlying genetic model. We compared the power of GRM association test to additive and other Cochran-Armitage Trend (CAT) tests through simulations and by applying GRM to a large case/control sample, the bipolar Welcome Trust Case Control Cohort data. Simulations were performed on two sets of case/control samples (1000/1000 and 2000/2000), using a large panel of genetic models. Four association tests (GRM and additive, recessive and dominant CAT tests) were applied to all replicates. RESULTS: We showed that GRM power to detect association was similar or greater than the additive CAT test, in particular in case of recessive inheritance, with up to 67% gain in power. GRM analysis of genome-wide bipolar disorder Welcome Trust Consortium data (1998 cases/3004 controls) showed significant association in the 16p12 region (rs420259; P = 3.4E-7) which has not been identified using the additive CAT test. As expected, rs42025 fitted a non-additive (recessive) model. CONCLUSIONS: GRM provides increased power compared to the additive CAT test for association studies and is easily applicable.

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.

A comprehensive genome-wide analysis of melanoma Breslow thickness identifies interaction between CDC42 and SCIN genetic variants.

Breslow thickness (BT) is a major prognostic factor of cutaneous melanoma (CM), the most fatal skin cancer. The genetic component of BT has only been explored by candidate gene studies with inconsistent results. Our objective was to uncover the genetic factors underlying BT using an hypothesis-free genome-wide approach. Our analysis strategy integrated a genome-wide association study (GWAS) of single nucleotide polymorphisms (SNPs) for BT followed by pathway analysis of GWAS outcomes using the gene-set enrichment analysis (GSEA) method and epistasis analysis within BT-associated pathways. This strategy was applied to two large CM datasets with Hapmap3-imputed SNP data: the French MELARISK study for discovery (966 cases) and the MD Anderson Cancer Center study (1,546 cases) for replication. While no marginal effect of individual SNPs was revealed through GWAS, three pathways, defined by gene ontology (GO) categories were significantly enriched in genes associated with BT (false discovery rate ≤5% in both studies): hormone activity, cytokine activity and myeloid cell differentiation. Epistasis analysis, within each significant GO, identified a statistically significant interaction between CDC42 and SCIN SNPs (pmeta-int =2.2 × 10(-6) , which met the overall multiple-testing corrected threshold of 2.5 × 10(-6) ). These two SNPs (and proxies) are strongly associated with CDC42 and SCIN gene expression levels and map to regulatory elements in skin cells. This interaction has important biological relevance since CDC42 and SCIN proteins have opposite effects in actin cytoskeleton organization and dynamics, a key mechanism underlying melanoma cell migration and invasion.

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.

Unsupervised text mining for assessing and augmenting GWAS results.

Text mining can assist in the analysis and interpretation of large-scale biomedical data, helping biologists to quickly and cheaply gain confirmation of hypothesized relationships between biological entities. We set this question in the context of genome-wide association studies (GWAS), an actively emerging field that contributed to identify many genes associated with multifactorial diseases. These studies allow to identify groups of genes associated with the same phenotype, but provide no information about the relationships between these genes. Therefore, our objective is to leverage unsupervised text mining techniques using text-based cosine similarity comparisons and clustering applied to candidate and random gene vectors, in order to augment the GWAS results. We propose a generic framework which we used to characterize the relationships between 10 genes reported associated with asthma by a previous GWAS. The results of this experiment showed that the similarities between these 10 genes were significantly stronger than would be expected by chance (one-sided p-value<0.01). The clustering of observed and randomly selected gene also allowed to generate hypotheses about potential functional relationships between these genes and thus contributed to the discovery of new candidate genes for asthma.

Integrated pathway and epistasis analysis reveals interactive effect of genetic variants at TERF1 and AFAP1L2 loci on melanoma risk.

Genome-wide association studies (GWASs) have characterized 13 loci associated with melanoma, which only account for a small part of melanoma risk. To identify new genes with too small an effect to be detected individually but which collectively influence melanoma risk and/or show interactive effects, we used a two-step analysis strategy including pathway analysis of genome-wide SNP data, in a first step, and epistasis analysis within significant pathways, in a second step. Pathway analysis, using the gene-set enrichment analysis (GSEA) approach and the gene ontology (GO) database, was applied to the outcomes of MELARISK (3,976 subjects) and MDACC (2,827 subjects) GWASs. Cross-gene SNP-SNP interaction analysis within melanoma-associated GOs was performed using the INTERSNP software. Five GO categories were significantly enriched in genes associated with melanoma (false discovery rate ≤ 5% in both studies): response to light stimulus, regulation of mitotic cell cycle, induction of programmed cell death, cytokine activity and oxidative phosphorylation. Epistasis analysis, within each of the five significant GOs, showed significant evidence for interaction for one SNP pair at TERF1 and AFAP1L2 loci (pmeta-int = 2.0 × 10(-7) , which met both the pathway and overall multiple-testing corrected thresholds that are equal to 9.8 × 10(-7) and 2.0 × 10(-7) , respectively) and suggestive evidence for another pair involving correlated SNPs at the same loci (pmeta-int = 3.6 × 10(-6) ). This interaction has important biological relevance given the key role of TERF1 in telomere biology and the reported physical interaction between TERF1 and AFAP1L2 proteins. This finding brings a novel piece of evidence for the emerging role of telomere dysfunction into melanoma development.

Fine mapping of genetic susceptibility loci for melanoma reveals a mixture of single variant and multiple variant regions.

At least 17 genomic regions are established as harboring melanoma susceptibility variants, in most instances with genome-wide levels of significance and replication in independent samples. Based on genome-wide single nucleotide polymorphism (SNP) data augmented by imputation to the 1,000 Genomes reference panel, we have fine mapped these regions in over 5,000 individuals with melanoma (mainly from the GenoMEL consortium) and over 7,000 ethnically matched controls. A penalized regression approach was used to discover those SNP markers that most parsimoniously explain the observed association in each genomic region. For the majority of the regions, the signal is best explained by a single SNP, which sometimes, as in the tyrosinase region, is a known functional variant. However in five regions the explanation is more complex. At the CDKN2A locus, for example, there is strong evidence that not only multiple SNPs but also multiple genes are involved. Our results illustrate the variability in the biology underlying genome-wide susceptibility loci and make steps toward accounting for some of the "missing heritability."

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.