USP15 regulates type I interferon response and is required for pathogenesis of neuroinflammation.

Genes and pathways in which inactivation dampens tissue inflammation present new opportunities for understanding the pathogenesis of common human inflammatory diseases, including inflammatory bowel disease, rheumatoid arthritis and multiple sclerosis. We identified a mutation in the gene encoding the deubiquitination enzyme USP15 (Usp15L749R) that protected mice against both experimental cerebral malaria (ECM) induced by Plasmodium berghei and experimental autoimmune encephalomyelitis (EAE). Combining immunophenotyping and RNA sequencing in brain (ECM) and spinal cord (EAE) revealed that Usp15L749R-associated resistance to neuroinflammation was linked to dampened type I interferon responses in situ. In hematopoietic cells and in resident brain cells, USP15 was coexpressed with, and functionally acted together with the E3 ubiquitin ligase TRIM25 to positively regulate type I interferon responses and to promote pathogenesis during neuroinflammation. The USP15-TRIM25 dyad might be a potential target for intervention in acute or chronic states of neuroinflammation.

Diffusion imaging genomics provides novel insight into early mechanisms of cerebral small vessel disease.

Cerebral small vessel disease (cSVD) is a leading cause of stroke and dementia. Genetic risk loci for white matter hyperintensities (WMH), the most common MRI-marker of cSVD in older age, were recently shown to be significantly associated with white matter (WM) microstructure on diffusion tensor imaging (signal-based) in young adults. To provide new insights into these early changes in WM microstructure and their relation with cSVD, we sought to explore the genetic underpinnings of cutting-edge tissue-based diffusion imaging markers across the adult lifespan. We conducted a genome-wide association study of neurite orientation dispersion and density imaging (NODDI) markers in young adults (i-Share study: N = 1 758, (mean[range]) 22.1[18-35] years), with follow-up in young middle-aged (Rhineland Study: N = 714, 35.2[30-40] years) and late middle-aged to older individuals (UK Biobank: N = 33 224, 64.3[45-82] years). We identified 21 loci associated with NODDI markers across brain regions in young adults. The most robust association, replicated in both follow-up cohorts, was with Neurite Density Index (NDI) at chr5q14.3, a known WMH locus in VCAN. Two additional loci were replicated in UK Biobank, at chr17q21.2 with NDI, and chr19q13.12 with Orientation Dispersion Index (ODI). Transcriptome-wide association studies showed associations of STAT3 expression in arterial and adipose tissue (chr17q21.2) with NDI, and of several genes at chr19q13.12 with ODI. Genetic susceptibility to larger WMH volume, but not to vascular risk factors, was significantly associated with decreased NDI in young adults, especially in regions known to harbor WMH in older age. Individually, seven of 25 known WMH risk loci were associated with NDI in young adults. In conclusion, we identified multiple novel genetic risk loci associated with NODDI markers, particularly NDI, in early adulthood. These point to possible early-life mechanisms underlying cSVD and to processes involving remyelination, neurodevelopment and neurodegeneration, with a potential for novel approaches to prevention.

Genomic variants associated with age at diagnosis of childhood-onset type 1 diabetes.

Age at diagnosis (AAD) of Type 1 diabetes (T1D) is determined by the age at onset of the autoimmune attack and by the rate of beta cell destruction that follows. Twin studies found that T1D AAD is strongly influenced by genetics, notably in young children. In young UK, Finnish, Sardinian patients AAD-associated genomic variants were previously identified, which may vary across populations and with time. In 1956 children of European ancestry born in mainland France in 1980-2008 who declared T1D before 15 years, we tested 94 T1D-associated SNPs for their association with AAD using nonparametric Kruskal-Wallis test. While high-risk HLA genotypes were not found to be associated with AAD, fourteen SNPs located in 12 non-HLA loci showed a strong association (2.9 × 10-12 < P < 1.4 × 10-3 after FDR correction). Four of these loci have been associated with AAD in previous cohorts (GSDMB, IL2, TNFAIP3, IL1), supporting a partially shared genetic influence on AAD of T1D in the studied European populations. In contrast, the association of 8 new loci CLEC16A, TYK2, ERBB3, CCR7, FCRL3, DNAH2, FGF3/4, and HPSE2 with AAD is novel. The 12 protein-coding genes located within these loci are involved in major immune pathways or in predisposition to other autoimmune diseases, which suggests a prominent role for these genes in the early immune mechanisms of beta cell destruction.

Contribution of genetic variants in the development of familial premature coronary artery disease in a cohort of cardiac patients.

Coronary artery disease (CAD), the most prevalent cardiovascular disease, is the leading cause of death worldwide. Heritable factors play a significant role in the pathogenesis of CAD. It has been proposed that approximately one-third of patients with CAD have a positive family history, and individuals with such history are at ~1.5-fold increased risk of CAD in their lifespans. Accordingly, the long-recognized familial clustering of CAD is a strong risk factor for this disease. Our study aimed to identify candidate genetic variants contributing to CAD by studying a cohort of 60 large Iranian families with at least two members in different generations afflicted with premature CAD (PCAD), defined as established disease at ≤45 years in men and ≤55 years in women. Exome sequencing was performed for a subset of the affected individuals, followed by prioritization and Sanger sequencing of candidate variants in all available family members. Subsequently, apparently healthy carriers of potential risk variants underwent coronary computed tomography angiography (CCTA), followed by co-segregation analysis of the combined data. Putative causal variants were identified in seven genes, ABCG8, CD36, CYP27A1, PIK3C2G, RASSF9, RYR2, and ZFYVE21, co-segregating with familial PCAD in seven unrelated families. Among these, PIK3C2G, RASSF9, and ZFYVE21 are novel candidate CAD susceptibility genes. Our findings indicate that rare variants in genes identified in this study are involved in CAD development.

Common Genetic Variants Contribute to Risk of Transposition of the Great Arteries.

RATIONALE: Dextro-transposition of the great arteries (D-TGA) is a severe congenital heart defect which affects approximately 1 in 4,000 live births. While there are several reports of D-TGA patients with rare variants in individual genes, the majority of D-TGA cases remain genetically elusive. Familial recurrence patterns and the observation that most cases with D-TGA are sporadic suggest a polygenic inheritance for the disorder, yet this remains unexplored. OBJECTIVE: We sought to study the role of common single nucleotide polymorphisms (SNPs) in risk for D-TGA. METHODS AND RESULTS: We conducted a genome-wide association study in an international set of 1,237 patients with D-TGA and identified a genome-wide significant susceptibility locus on chromosome 3p14.3, which was subsequently replicated in an independent case-control set (rs56219800, meta-analysis P=8.6x10-10, OR=0.69 per C allele). SNP-based heritability analysis showed that 25% of variance in susceptibility to D-TGA may be explained by common variants. A genome-wide polygenic risk score derived from the discovery set was significantly associated to D-TGA in the replication set (P=4x10-5). The genome-wide significant locus (3p14.3) co-localizes with a putative regulatory element that interacts with the promoter of WNT5A, which encodes the Wnt Family Member 5A protein known for its role in cardiac development in mice. We show that this element drives reporter gene activity in the developing heart of mice and zebrafish and is bound by the developmental transcription factor TBX20. We further demonstrate that TBX20 attenuates Wnt5a expression levels in the developing mouse heart. CONCLUSIONS: This work provides support for a polygenic architecture in D-TGA and identifies a susceptibility locus on chromosome 3p14.3 near WNT5A. Genomic and functional data support a causal role of WNT5A at the locus.

Dyslipidemia, inflammation, calcification, and adiposity in aortic stenosis: a genome-wide study.

AIMS: Although highly heritable, the genetic etiology of calcific aortic stenosis (AS) remains incompletely understood. The aim of this study was to discover novel genetic contributors to AS and to integrate functional, expression, and cross-phenotype data to identify mechanisms of AS. METHODS AND RESULTS: A genome-wide meta-analysis of 11.6 million variants in 10 cohorts involving 653 867 European ancestry participants (13 765 cases) was performed. Seventeen loci were associated with AS at P ≤ 5 × 10-8, of which 15 replicated in an independent cohort of 90 828 participants (7111 cases), including CELSR2-SORT1, NLRP6, and SMC2. A genetic risk score comprised of the index variants was associated with AS [odds ratio (OR) per standard deviation, 1.31; 95% confidence interval (CI), 1.26-1.35; P = 2.7 × 10-51] and aortic valve calcium (OR per standard deviation, 1.22; 95% CI, 1.08-1.37; P = 1.4 × 10-3), after adjustment for known risk factors. A phenome-wide association study indicated multiple associations with coronary artery disease, apolipoprotein B, and triglycerides. Mendelian randomization supported a causal role for apolipoprotein B-containing lipoprotein particles in AS (OR per g/L of apolipoprotein B, 3.85; 95% CI, 2.90-5.12; P = 2.1 × 10-20) and replicated previous findings of causality for lipoprotein(a) (OR per natural logarithm, 1.20; 95% CI, 1.17-1.23; P = 4.8 × 10-73) and body mass index (OR per kg/m2, 1.07; 95% CI, 1.05-1.9; P = 1.9 × 10-12). Colocalization analyses using the GTEx database identified a role for differential expression of the genes LPA, SORT1, ACTR2, NOTCH4, IL6R, and FADS. CONCLUSION: Dyslipidemia, inflammation, calcification, and adiposity play important roles in the etiology of AS, implicating novel treatments and prevention strategies.

Gene-mapping study of extremes of cerebral small vessel disease reveals TRIM47 as a strong candidate.

Cerebral small vessel disease is a leading cause of stroke and a major contributor to cognitive decline and dementia, but our understanding of specific genes underlying the cause of sporadic cerebral small vessel disease is limited. We report a genome-wide association study and a whole-exome association study on a composite extreme phenotype of cerebral small vessel disease derived from its most common MRI features: white matter hyperintensities and lacunes. Seventeen population-based cohorts of older persons with MRI measurements and genome-wide genotyping (n = 41 326), whole-exome sequencing (n = 15 965), or exome chip (n = 5249) data contributed 13 776 and 7079 extreme small vessel disease samples for the genome-wide association study and whole-exome association study, respectively. The genome-wide association study identified significant association of common variants in 11 loci with extreme small vessel disease, of which the chr12q24.11 locus was not previously reported to be associated with any MRI marker of cerebral small vessel disease. The whole-exome association study identified significant associations of extreme small vessel disease with common variants in the 5' UTR region of EFEMP1 (chr2p16.1) and one probably damaging common missense variant in TRIM47 (chr17q25.1). Mendelian randomization supports the causal association of extensive small vessel disease severity with increased risk of stroke and Alzheimer's disease. Combined evidence from summary-based Mendelian randomization studies and profiling of human loss-of-function allele carriers showed an inverse relation between TRIM47 expression in the brain and blood vessels and extensive small vessel disease severity. We observed significant enrichment of Trim47 in isolated brain vessel preparations compared to total brain fraction in mice, in line with the literature showing Trim47 enrichment in brain endothelial cells at single cell level. Functional evaluation of TRIM47 by small interfering RNAs-mediated knockdown in human brain endothelial cells showed increased endothelial permeability, an important hallmark of cerebral small vessel disease pathology. Overall, our comprehensive gene-mapping study and preliminary functional evaluation suggests a putative role of TRIM47 in the pathophysiology of cerebral small vessel disease, making it an important candidate for extensive in vivo explorations and future translational work.

Shared genomic segment analysis with equivalence testing.

An important aspect of disease gene mapping is replication, that is, a putative finding in one group of individuals is confirmed in another set of individuals. As it can happen by chance that individuals share an estimated disease position, we developed a statistical approach to determine the p-value for multiple individuals or families to share a possibly small number of candidate susceptibility variants. Here, we focus on candidate variants for dominant traits that have been obtained by our previously developed heterozygosity analysis, and we are testing the sharing of candidate variants obtained for different individuals. Our approach allows for multiple pathogenic variants in a gene to contribute to disease, and for estimated disease variant positions to be imprecise. Statistically, the method developed here falls into the category of equivalence testing, where the classical null and alternative hypotheses of homogeneity and heterogeneity are reversed. The null hypothesis situation is created by permuting genomic locations of variants for one individual after another. We applied our methodology to the ALSPAC data set of 1,927 whole-genome sequenced individuals, where some individuals carry a pathogenic variant for the BRCA1 gene, but no two individuals carry the same variant. Our shared genomic segment analysis found significant evidence for BRCA1 pathogenic variants within ±5 kb of a given DNA variant.

Correction: Whole exome sequencing study identifies novel rare and common Alzheimer's-Associated variants involved in immune response and transcriptional regulation.

A correction to this paper has been published and can be accessed via a link at the top of the paper.

Whole exome sequencing study identifies novel rare and common Alzheimer's-Associated variants involved in immune response and transcriptional regulation.

The Alzheimer's Disease Sequencing Project (ADSP) undertook whole exome sequencing in 5,740 late-onset Alzheimer disease (AD) cases and 5,096 cognitively normal controls primarily of European ancestry (EA), among whom 218 cases and 177 controls were Caribbean Hispanic (CH). An age-, sex- and APOE based risk score and family history were used to select cases most likely to harbor novel AD risk variants and controls least likely to develop AD by age 85 years. We tested ~1.5 million single nucleotide variants (SNVs) and 50,000 insertion-deletion polymorphisms (indels) for association to AD, using multiple models considering individual variants as well as gene-based tests aggregating rare, predicted functional, and loss of function variants. Sixteen single variants and 19 genes that met criteria for significant or suggestive associations after multiple-testing correction were evaluated for replication in four independent samples; three with whole exome sequencing (2,778 cases, 7,262 controls) and one with genome-wide genotyping imputed to the Haplotype Reference Consortium panel (9,343 cases, 11,527 controls). The top findings in the discovery sample were also followed-up in the ADSP whole-genome sequenced family-based dataset (197 members of 42 EA families and 501 members of 157 CH families). We identified novel and predicted functional genetic variants in genes previously associated with AD. We also detected associations in three novel genes: IGHG3 (p = 9.8 × 10-7), an immunoglobulin gene whose antibodies interact with ß-amyloid, a long non-coding RNA AC099552.4 (p = 1.2 × 10-7), and a zinc-finger protein ZNF655 (gene-based p = 5.0 × 10-6). The latter two suggest an important role for transcriptional regulation in AD pathogenesis.

Independent risk factors for simvastatin-related myopathy and relevance to different types of muscle symptom.

AIMS: Statins are widely used to prevent cardiovascular events, but little is known about the impact of different risk factors for statin-related myopathy or their relevance to reports of other types of muscle symptom. METHODS AND RESULTS: An observational analysis was undertaken of 171 clinically adjudicated cases of myopathy (defined as unexplained muscle pain or weakness with creatine kinase >10× upper limit of normal) and, separately, of 15 208 cases of other muscle symptoms among 58 390 individuals with vascular disease treated with simvastatin for a mean of 3.4 years. Cox proportional hazards models were used to identify independent predictors of myopathy. The rate of myopathy was low: 9 per 10 000 person-years of simvastatin therapy. Independent risk factors for myopathy included: simvastatin dose, ethnicity, sex, age, body mass index, medically treated diabetes, concomitant use of niacin-laropiprant, verapamil, beta-blockers, diltiazem and diuretics. In combination, these risk factors predicted more than a 30-fold risk difference between the top and bottom thirds of a myopathy risk score (hazard ratio : 34.35, 95% CI: 12.73-92.69, P across thirds = 9·1 × 10-48). However, despite the strong association with myopathy, this score was not associated with the other reported muscle symptoms (P across thirds = 0.93). Likewise, although SLCO1B1 genotype was associated with myopathy, it was not associated with other muscle symptoms. CONCLUSIONS: The absolute risk of simvastatin-related myopathy is low, but individuals at higher risk can be identified to help guide patient management. The lack of association of the myopathy risk score with other muscle symptoms reinforces randomized placebo-controlled evidence that statins do not cause the vast majority of reported muscle symptoms.

Dominant gut Prevotella copri in gastrectomised non-obese diabetic Goto-Kakizaki rats improves glucose homeostasis through enhanced FXR signalling.

AIMS/HYPOTHESIS: Drug and surgical-based therapies in type 2 diabetes are associated with altered gut microbiota architecture. Here we investigated the role of the gut microbiome in improved glucose homeostasis following bariatric surgery. METHODS: We carried out gut microbiome analyses in gastrectomised (by vertical sleeve gastrectomy [VSG]) rats of the Goto-Kakizaki (GK) non-obese model of spontaneously occurring type 2 diabetes, followed by physiological studies in the GK rat. RESULTS: VSG in the GK rat led to permanent improvement of glucose tolerance associated with minor changes in the gut microbiome, mostly characterised by significant enrichment of caecal Prevotella copri. Gut microbiota enrichment with P. copri in GK rats through permissive antibiotic treatment, inoculation of gut microbiota isolated from gastrectomised GK rats, and direct inoculation of P. copri, resulted in significant improvement of glucose tolerance, independent of changes in body weight. Plasma bile acids were increased in GK rats following inoculation with P. copri and P. copri-enriched microbiota from VSG-treated rats; the inoculated GK rats then showed increased liver glycogen and upregulated expression of Fxr (also known as Nr1h4), Srebf1c, Chrebp (also known as Mlxipl) and Il10 and downregulated expression of Cyp7a1. CONCLUSIONS: Our data underline the impact of intestinal P. copri on improved glucose homeostasis through enhanced bile acid metabolism and farnesoid X receptor (FXR) signalling, which may represent a promising opportunity for novel type 2 diabetes therapeutics.

ZBTB7B (ThPOK) Is Required for Pathogenesis of Cerebral Malaria and Protection against Pulmonary Tuberculosis.

We used a genome-wide screen in N-ethyl-N-nitrosourea (ENU)-mutagenized mice to identify genes in which recessive loss-of-function mutations protect against pathological neuroinflammation. We identified an R367Q mutation in the ZBTB7B (ThPOK) protein in which homozygosity causes protection against experimental cerebral malaria (ECM) caused by infection with Plasmodium berghei ANKA. Zbtb7bR367Q homozygous mice show a defect in the lymphoid compartment expressed as severe reduction in the number of single-positive CD4 T cells in the thymus and in the periphery, reduced brain infiltration of proinflammatory leukocytes in P. berghei ANKA-infected mice, and reduced production of proinflammatory cytokines by primary T cells ex vivo and in vivo Dampening of proinflammatory immune responses in Zbtb7bR367Q mice is concomitant to increased susceptibility to infection with avirulent (Mycobacterium bovis BCG) and virulent (Mycobacterium tuberculosis H37Rv) mycobacteria. The R367Q mutation maps to the first DNA-binding zinc finger domain of ThPOK and causes loss of base contact by R367 in the major groove of the DNA, which is predicted to impair DNA binding. Global immunoprecipitation of ThPOK-containing chromatin complexes coupled to DNA sequencing (ChIP-seq) identified transcriptional networks and candidate genes likely to play key roles in CD4+ CD8+ T cell development and in the expression of lineage-specific functions of these cells. This study highlights ThPOK as a global regulator of immune function in which alterations may affect normal responses to infectious and inflammatory stimuli.

Dysregulation at multiple points of the kynurenine pathway is a ubiquitous feature of renal cancer: implications for tumour immune evasion.

BACKGROUND: Indoleamine 2,3-dioxygenase (IDO), the first step in the kynurenine pathway (KP), is upregulated in some cancers and represents an attractive therapeutic target given its role in tumour immune evasion. However, the recent failure of an IDO inhibitor in a late phase trial raises questions about this strategy. METHODS: Matched renal cell carcinoma (RCC) and normal kidney tissues were subject to proteomic profiling. Tissue immunohistochemistry and gene expression data were used to validate findings. Phenotypic effects of loss/gain of expression were examined in vitro. RESULTS: Quinolate phosphoribosyltransferase (QPRT), the final and rate-limiting enzyme in the KP, was identified as being downregulated in RCC. Loss of QPRT expression led to increased potential for anchorage-independent growth. Gene expression, mass spectrometry (clear cell and chromophobe RCC) and tissue immunohistochemistry (clear cell, papillary and chromophobe), confirmed loss or decreased expression of QPRT and showed downregulation of other KP enzymes, including kynurenine 3-monoxygenase (KMO) and 3-hydroxyanthranilate-3,4-dioxygenase (HAAO), with a concomitant maintenance or upregulation of nicotinamide phosphoribosyltransferase (NAMPT), the key enzyme in the NAD+ salvage pathway. CONCLUSIONS: Widespread dysregulation of the KP is common in RCC and is likely to contribute to tumour immune evasion, carrying implications for effective therapeutic targeting of this critical pathway.

Mutations in SULT2B1 Cause Autosomal-Recessive Congenital Ichthyosis in Humans.

Ichthyoses are a clinically and genetically heterogeneous group of genodermatoses associated with abnormal scaling of the skin over the whole body. Mutations in nine genes are known to cause non-syndromic forms of autosomal-recessive congenital ichthyosis (ARCI). However, not all genetic causes for ARCI have been discovered to date. Using whole-exome sequencing (WES) and multigene panel screening, we identified 6 ARCI-affected individuals from three unrelated families with mutations in Sulfotransferase family 2B member 1 (SULT2B1), showing their causative association with ARCI. Cytosolic sulfotransferases form a large family of enzymes that are involved in the synthesis and metabolism of several steroids in humans. We identified four distinct mutations including missense, nonsense, and splice site mutations. We demonstrated the loss of SULT2B1 expression at RNA and protein levels in keratinocytes from individuals with ARCI by functional analyses. Furthermore, we succeeded in reconstructing the morphologic skin alterations in a 3D organotypic tissue culture model with SULT2B1-deficient keratinocytes and fibroblasts. By thin layer chromatography (TLC) of extracts from these organotypic cultures, we could show the absence of cholesterol sulfate, the metabolite of SULT2B1, and an increased level of cholesterol, indicating a disturbed cholesterol metabolism of the skin upon loss-of-function mutation in SULT2B1. In conclusion, our study reveals an essential role for SULT2B1 in the proper development of healthy human skin. Mutation in SULT2B1 leads to an ARCI phenotype via increased proliferation of human keratinocytes, thickening of epithelial layers, and altered epidermal cholesterol metabolism.

Association of variants in HTRA1 and NOTCH3 with MRI-defined extremes of cerebral small vessel disease in older subjects.

We report a composite extreme phenotype design using distribution of white matter hyperintensities and brain infarcts in a population-based cohort of older persons for gene-mapping of cerebral small vessel disease. We demonstrate its application in the 3C-Dijon whole exome sequencing (WES) study (n = 1924, nWESextremes = 512), with both single variant and gene-based association tests. We used other population-based cohort studies participating in the CHARGE consortium for replication, using whole exome sequencing (nWES = 2,868, nWESextremes = 956) and genome-wide genotypes (nGW = 9924, nGWextremes = 3308). We restricted our study to candidate genes known to harbour mutations for Mendelian small vessel disease: NOTCH3, HTRA1, COL4A1, COL4A2 and TREX1. We identified significant associations of a common intronic variant in HTRA1, rs2293871 using single variant association testing (Pdiscovery = 8.21 × 10-5, Preplication = 5.25 × 10-3, Pcombined = 4.72 × 10-5) and of NOTCH3 using gene-based tests (Pdiscovery = 1.61 × 10-2, Preplication = 3.99 × 10-2, Pcombined = 5.31 × 10-3). Follow-up analysis identified significant association of rs2293871 with small vessel ischaemic stroke, and two blood expression quantitative trait loci of HTRA1 in linkage disequilibrium. Additionally, we identified two participants in the 3C-Dijon cohort (0.4%) carrying heterozygote genotypes at known pathogenic variants for familial small vessel disease within NOTCH3 and HTRA1. In conclusion, our proof-of-concept study provides strong evidence that using a novel composite MRI-derived phenotype for extremes of small vessel disease can facilitate the identification of genetic variants underlying small vessel disease, both common variants and those with rare and low frequency. The findings demonstrate shared mechanisms and a continuum between genes underlying Mendelian small vessel disease and those contributing to the common, multifactorial form of the disease.

Heterozygosity mapping for human dominant trait variants.

Homozygosity mapping is a well-known technique to identify runs of homozygous variants that are likely to harbor genes responsible for autosomal recessive disease, but a comparable method for autosomal dominant traits has been lacking. We developed an approach to map dominant disease genes based on heterozygosity frequencies of sequence variants in the immediate vicinity of a dominant trait. We demonstrate through theoretical analysis that DNA variants surrounding an inherited dominant disease variant tend to have increased heterozygosity compared with variants elsewhere in the genome. We confirm existence of this phenomenon in sequence data with known dominant pathogenic variants obtained on family members and in unrelated population controls. A computer-based approach to estimating empirical significance levels associated with our test statistics shows genome-wide p-values smaller than 0.05 for many but not all of the individuals carrying a pathogenic variant.

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.

Detecting multi-way epistasis in family-based association studies.

The era of genome-wide association studies (GWAS) has led to the discovery of numerous genetic variants associated with disease. Better understanding of whether these or other variants interact leading to differential risk compared with individual marker effects will increase our understanding of the genetic architecture of disease, which may be investigated using the family-based study design. We present M-TDT (the multi-locus transmission disequilibrium test), a tool for detecting family-based multi-locus multi-allelic effects for qualitative or quantitative traits, extended from the original transmission disequilibrium test (TDT). Tests to handle the comparison between additive and epistatic models, lack of independence between markers and multiple offspring are described. Performance of M-TDT is compared with a multifactor dimensionality reduction (MDR) approach designed for investigating families in the hypothesis-free genome-wide setting (the multifactor dimensionality reduction pedigree disequilibrium test, MDR-PDT). Other methods derived from the TDT or MDR to investigate genetic interaction in the family-based design are also discussed. The case of three independent biallelic loci is illustrated using simulations for one- to three-locus alternative hypotheses. M-TDT identified joint-locus effects and distinguished effectively between additive and epistatic models. We showed a practical example of M-TDT based on three genes already known to be implicated in malaria susceptibility. Our findings demonstrate the value of M-TDT in a hypothesis-driven context to test for multi-way epistasis underlying common disease etiology, whereas MDR-PDT-based methods are more appropriate in a hypothesis-free genome-wide setting.