Complement genes contribute sex-biased vulnerability in diverse disorders.

Many common illnesses, for reasons that have not been identified, differentially affect men and women. For instance, the autoimmune diseases systemic lupus erythematosus (SLE) and Sjögren's syndrome affect nine times more women than men1, whereas schizophrenia affects men with greater frequency and severity relative to women2. All three illnesses have their strongest common genetic associations in the major histocompatibility complex (MHC) locus, an association that in SLE and Sjögren's syndrome has long been thought to arise from alleles of the human leukocyte antigen (HLA) genes at that locus3-6. Here we show that variation of the complement component 4 (C4) genes C4A and C4B, which are also at the MHC locus and have been linked to increased risk for schizophrenia7, generates 7-fold variation in risk for SLE and 16-fold variation in risk for Sjögren's syndrome among individuals with common C4 genotypes, with C4A protecting more strongly than C4B in both illnesses. The same alleles that increase risk for schizophrenia greatly reduce risk for SLE and Sjögren's syndrome. In all three illnesses, C4 alleles act more strongly in men than in women: common combinations of C4A and C4B generated 14-fold variation in risk for SLE, 31-fold variation in risk for Sjögren's syndrome, and 1.7-fold variation in schizophrenia risk among men (versus 6-fold, 15-fold and 1.26-fold variation in risk among women, respectively). At a protein level, both C4 and its effector C3 were present at higher levels in cerebrospinal fluid and plasma8,9 in men than in women among adults aged between 20 and 50 years, corresponding to the ages of differential disease vulnerability. Sex differences in complement protein levels may help to explain the more potent effects of C4 alleles in men, women's greater risk of SLE and Sjögren's syndrome and men's greater vulnerability to schizophrenia. These results implicate the complement system as a source of sexual dimorphism in vulnerability to diverse illnesses.

Genetic Inhibition of APOL1 Pore-Forming Function Prevents APOL1-Mediated Kidney Disease.

SIGNIFICANCE STATEMENT: African Americans are at increased risk of CKD in part due to high-risk (HR) variants in the apolipoprotein L1 ( APOL1 ) gene, termed G1/G2. A different APOL1 variant, p.N264K , reduced the risk of CKD and ESKD among carriers of APOL1 HR variants to levels comparable with individuals with APOL1 low-risk variants in an analysis of 121,492 participants of African ancestry from the Million Veteran Program (MVP). Functional genetic studies in cell models showed that APOL1 p.N264K blocked APOL1 pore-forming function and ion channel conduction and reduced toxicity of APOL1 HR mutations. Pharmacologic inhibitors that mimic this mutation blocking APOL1 -mediated pore formation may be able to prevent and/or treat APOL1 -associated kidney disease. BACKGROUND: African Americans are at increased risk for nondiabetic CKD in part due to HR variants in the APOL1 gene. METHODS: We tested whether a different APOL1 variant, p.N264K , modified the association between APOL1 HR genotypes (two copies of G1/G2) and CKD in a cross-sectional analysis of 121,492 participants of African ancestry from the MVP. We replicated our findings in the Vanderbilt University Biobank ( n =14,386) and National Institutes of Health All of Us ( n =14,704). Primary outcome was CKD and secondary outcome was ESKD among nondiabetic patients. Primary analysis compared APOL1 HR genotypes with and without p.N264K . Secondary analyses included APOL1 low-risk genotypes and tested for interaction. In MVP, we performed sequential logistic regression models adjusting for demographics, comorbidities, medications, and ten principal components of ancestry. Functional genomic studies expressed APOL1 HR variants with and without APOL1 p.N264K in cell models. RESULTS: In the MVP cohort, 15,604 (12.8%) had two APOL1 HR variants, of which 582 (0.5%) also had APOL1 p.N264K . In MVP, 18,831 (15%) had CKD, 4177 (3%) had ESKD, and 34% had diabetes. MVP APOL1 HR, without p.N264K , was associated with increased odds of CKD (odds ratio [OR], 1.72; 95% confidence interval [CI], 1.60 to 1.85) and ESKD (OR, 3.94; 95% CI, 3.52 to 4.41). In MVP, APOL1 p.N264K mitigated the renal risk of APOL1 HR, in CKD (OR, 0.43; 95% CI, 0.28 to 0.65) and ESKD (OR, 0.19; CI 0.07 to 0.51). In the replication cohorts meta-analysis, APOL1 p.N264K mitigated the renal risk of APOL1 HR in CKD (OR, 0.40; 95% CI, 0.18 to 0.92) and ESKD (OR, 0.19; 95% CI, 0.05 to 0.79). In the mechanistic studies, APOL1 p.N264K blocked APOL1 pore-forming function and ion channel conduction and reduced toxicity of APOL1 HR variants. CONCLUSIONS: APOL1 p.N264K is associated with reduced risk of CKD and ESKD among carriers of APOL1 HR to levels comparable with individuals with APOL1 low-risk genotypes.

Rare protein-altering variants in ANGPTL7 lower intraocular pressure and protect against glaucoma.

Protein-altering variants that are protective against human disease provide in vivo validation of therapeutic targets. Here we use genotyping data from UK Biobank (n = 337,151 unrelated White British individuals) and FinnGen (n = 176,899) to conduct a search for protein-altering variants conferring lower intraocular pressure (IOP) and protection against glaucoma. Through rare protein-altering variant association analysis, we find a missense variant in ANGPTL7 in UK Biobank (rs28991009, p.Gln175His, MAF = 0.8%, genotyped in 82,253 individuals with measured IOP and an independent set of 4,238 glaucoma patients and 250,660 controls) that significantly lowers IOP (ß = -0.53 and -0.67 mmHg for heterozygotes, -3.40 and -2.37 mmHg for homozygotes, P = 5.96 x 10-9 and 1.07 x 10-13 for corneal compensated and Goldman-correlated IOP, respectively) and is associated with 34% reduced risk of glaucoma (P = 0.0062). In FinnGen, we identify an ANGPTL7 missense variant at a greater than 50-fold increased frequency in Finland compared with other populations (rs147660927, p.Arg220Cys, MAF Finland = 4.3%), which was genotyped in 6,537 glaucoma patients and 170,362 controls and is associated with a 29% lower glaucoma risk (P = 1.9 x 10-12 for all glaucoma types and also protection against its subtypes including exfoliation, primary open-angle, and primary angle-closure). We further find three rarer variants in UK Biobank, including a protein-truncating variant, which confer a strong composite lowering of IOP (P = 0.0012 and 0.24 for Goldman-correlated and corneal compensated IOP, respectively), suggesting the protective mechanism likely resides in the loss of interaction or function. Our results support inhibition or down-regulation of ANGPTL7 as a therapeutic strategy for glaucoma.

Genome-wide association study identifies kallikrein 5 in type 2 inflammation-low asthma.

BACKGROUND: Clinical studies of type 2 (T2) cytokine-related neutralizing antibodies in asthma have identified a substantial subset of patients with low levels of T2 inflammation who do not benefit from T2 cytokine neutralizing antibody treatment. Non-T2 mechanisms are poorly understood in asthma but represent a redefined unmet medical need. OBJECTIVE: We sought to gain a better understanding of genetic contributions to T2-low asthma. METHODS: We utilized an unbiased genome-wide association study of patients with moderate to severe asthma stratified by T2 serum biomarker periostin. We also performed additional expression and biological analysis for the top genetic hits. RESULTS: We identified a novel protective single nucleotide polymorphism at chr19q13.41, which is selectively associated with T2-low asthma and establishes Kallikrein-related peptidase 5 (KLK5) as the causal gene mediating this association. Heterozygous carriers of the single nucleotide polymorphisms have reduced KLK5 expression. KLK5 is secreted by human bronchial epithelial cells and elevated in asthma bronchial alveolar lavage. T2 cytokines IL-4 and IL-13 downregulate KLK5 in human bronchial epithelial cells. KLK5, dependent on its catalytic function, induces epithelial chemokine/cytokine expression. Finally, overexpression of KLK5 in airway or lack of an endogenous KLK5 inhibitor, SPINK5, leads to spontaneous airway neutrophilic inflammation. CONCLUSION: Our data identify KLK5 to be the causal gene at a novel locus at chr19q13.41 associated with T2-low asthma.

Genome-wide association study and functional validation implicates JADE1 in tauopathy.

Primary age-related tauopathy (PART) is a neurodegenerative pathology with features distinct from but also overlapping with Alzheimer disease (AD). While both exhibit Alzheimer-type temporal lobe neurofibrillary degeneration alongside amnestic cognitive impairment, PART develops independently of amyloid-ß (Aß) plaques. The pathogenesis of PART is not known, but evidence suggests an association with genes that promote tau pathology and others that protect from Aß toxicity. Here, we performed a genetic association study in an autopsy cohort of individuals with PART (n = 647) using Braak neurofibrillary tangle stage as a quantitative trait. We found some significant associations with candidate loci associated with AD (SLC24A4, MS4A6A, HS3ST1) and progressive supranuclear palsy (MAPT and EIF2AK3). Genome-wide association analysis revealed a novel significant association with a single nucleotide polymorphism on chromosome 4 (rs56405341) in a locus containing three genes, including JADE1 which was significantly upregulated in tangle-bearing neurons by single-soma RNA-seq. Immunohistochemical studies using antisera targeting JADE1 protein revealed localization within tau aggregates in autopsy brains with four microtubule-binding domain repeats (4R) isoforms and mixed 3R/4R, but not with 3R exclusively. Co-immunoprecipitation in post-mortem human PART brain tissue revealed a specific binding of JADE1 protein to four repeat tau lacking N-terminal inserts (0N4R). Finally, knockdown of the Drosophila JADE1 homolog rhinoceros (rno) enhanced tau-induced toxicity and apoptosis in vivo in a humanized 0N4R mutant tau knock-in model, as quantified by rough eye phenotype and terminal deoxynucleotidyl transferase dUTP nick end-labeling (TUNEL) in the fly brain. Together, these findings indicate that PART has a genetic architecture that partially overlaps with AD and other tauopathies and suggests a novel role for JADE1 as a modifier of neurofibrillary degeneration.

Paired Immunoglobulin-like Type 2 Receptor Alpha G78R variant alters ligand binding and confers protection to Alzheimer's disease.

Paired Immunoglobulin-like Type 2 Receptor Alpha (PILRA) is a cell surface inhibitory receptor that recognizes specific O-glycosylated proteins and is expressed on various innate immune cell types including microglia. We show here that a common missense variant (G78R, rs1859788) of PILRA is the likely causal allele for the confirmed Alzheimer's disease risk locus at 7q21 (rs1476679). The G78R variant alters the interaction of residues essential for sialic acid engagement, resulting in >50% reduced binding for several PILRA ligands including a novel ligand, complement component 4A, and herpes simplex virus 1 (HSV-1) glycoprotein B. PILRA is an entry receptor for HSV-1 via glycoprotein B, and macrophages derived from R78 homozygous donors showed significantly decreased levels of HSV-1 infection at several multiplicities of infection compared to homozygous G78 macrophages. We propose that PILRA G78R protects individuals from Alzheimer's disease risk via reduced inhibitory signaling in microglia and reduced microglial infection during HSV-1 recurrence.

A Crohn's disease variant in Atg16l1 enhances its degradation by caspase 3.

Crohn's disease is a debilitating inflammatory bowel disease (IBD) that can involve the entire digestive tract. A single-nucleotide polymorphism (SNP) encoding a missense variant in the autophagy gene ATG16L1 (rs2241880, Thr300Ala) is strongly associated with the incidence of Crohn's disease. Numerous studies have demonstrated the effect of ATG16L1 deletion or deficiency; however, the molecular consequences of the Thr300Ala (T300A) variant remains unknown. Here we show that amino acids 296-299 constitute a caspase cleavage motif in ATG16L1 and that the T300A variant (T316A in mice) significantly increases ATG16L1 sensitization to caspase-3-mediated processing. We observed that death-receptor activation or starvation-induced metabolic stress in human and murine macrophages increased degradation of the T300A or T316A variants of ATG16L1, respectively, resulting in diminished autophagy. Knock-in mice harbouring the T316A variant showed defective clearance of the ileal pathogen Yersinia enterocolitica and an elevated inflammatory cytokine response. In turn, deletion of the caspase-3-encoding gene, Casp3, or elimination of the caspase cleavage site by site-directed mutagenesis rescued starvation-induced autophagy and pathogen clearance, respectively. These findings demonstrate that caspase 3 activation in the presence of a common risk allele leads to accelerated degradation of ATG16L1, placing cellular stress, apoptotic stimuli and impaired autophagy in a unified pathway that predisposes to Crohn's disease.

Genetics of rheumatoid arthritis contributes to biology and drug discovery.

A major challenge in human genetics is to devise a systematic strategy to integrate disease-associated variants with diverse genomic and biological data sets to provide insight into disease pathogenesis and guide drug discovery for complex traits such as rheumatoid arthritis (RA). Here we performed a genome-wide association study meta-analysis in a total of >100,000 subjects of European and Asian ancestries (29,880 RA cases and 73,758 controls), by evaluating ∼10 million single-nucleotide polymorphisms. We discovered 42 novel RA risk loci at a genome-wide level of significance, bringing the total to 101 (refs 2 - 4). We devised an in silico pipeline using established bioinformatics methods based on functional annotation, cis-acting expression quantitative trait loci and pathway analyses--as well as novel methods based on genetic overlap with human primary immunodeficiency, haematological cancer somatic mutations and knockout mouse phenotypes--to identify 98 biological candidate genes at these 101 risk loci. We demonstrate that these genes are the targets of approved therapies for RA, and further suggest that drugs approved for other indications may be repurposed for the treatment of RA. Together, this comprehensive genetic study sheds light on fundamental genes, pathways and cell types that contribute to RA pathogenesis, and provides empirical evidence that the genetics of RA can provide important information for drug discovery.

Previously reported placebo-response-associated variants do not predict patient outcomes in inflammatory disease Phase III trial placebo arms.

In clinical trials, a placebo response refers to improvement in disease symptoms arising from the psychological effect of receiving a treatment rather than the actual treatment under investigation. Previous research has reported genomic variation associated with the likelihood of observing a placebo response, but these studies have been limited in scope and have not been validated. Here, we analyzed whole-genome sequencing data from 784 patients undergoing placebo treatment in Phase III Asthma or Rheumatoid Arthritis trials to assess the impact of previously reported variation on patient outcomes in the placebo arms and to identify novel variants associated with the placebo response. Contrary to expectations based on previous reports, we did not observe any statistically significant associations between genomic variants and placebo treatment outcome. Our findings suggest that the biological origin of the placebo response is complex and likely to be variable between disease areas.

A mutation in APP protects against Alzheimer's disease and age-related cognitive decline.

The prevalence of dementia in the Western world in people over the age of 60 has been estimated to be greater than 5%, about two-thirds of which are due to Alzheimer's disease. The age-specific prevalence of Alzheimer's disease nearly doubles every 5 years after age 65, leading to a prevalence of greater than 25% in those over the age of 90 (ref. 3). Here, to search for low-frequency variants in the amyloid-ß precursor protein (APP) gene with a significant effect on the risk of Alzheimer's disease, we studied coding variants in APP in a set of whole-genome sequence data from 1,795 Icelanders. We found a coding mutation (A673T) in the APP gene that protects against Alzheimer's disease and cognitive decline in the elderly without Alzheimer's disease. This substitution is adjacent to the aspartyl protease ß-site in APP, and results in an approximately 40% reduction in the formation of amyloidogenic peptides in vitro. The strong protective effect of the A673T substitution against Alzheimer's disease provides proof of principle for the hypothesis that reducing the ß-cleavage of APP may protect against the disease. Furthermore, as the A673T allele also protects against cognitive decline in the elderly without Alzheimer's disease, the two may be mediated through the same or similar mechanisms.

Identifying and mitigating batch effects in whole genome sequencing data.

BACKGROUND: Large sample sets of whole genome sequencing with deep coverage are being generated, however assembling datasets from different sources inevitably introduces batch effects. These batch effects are not well understood and can be due to changes in the sequencing protocol or bioinformatics tools used to process the data. No systematic algorithms or heuristics exist to detect and filter batch effects or remove associations impacted by batch effects in whole genome sequencing data. RESULTS: We describe key quality metrics, provide a freely available software package to compute them, and demonstrate that identification of batch effects is aided by principal components analysis of these metrics. To mitigate batch effects, we developed new site-specific filters that identified and removed variants that falsely associated with the phenotype due to batch effect. These include filtering based on: a haplotype based genotype correction, a differential genotype quality test, and removing sites with missing genotype rate greater than 30% after setting genotypes with quality scores less than 20 to missing. This method removed 96.1% of unconfirmed genome-wide significant SNP associations and 97.6% of unconfirmed genome-wide significant indel associations. We performed analyses to demonstrate that: 1) These filters impacted variants known to be disease associated as 2 out of 16 confirmed associations in an AMD candidate SNP analysis were filtered, representing a reduction in power of 12.5%, 2) In the absence of batch effects, these filters removed only a small proportion of variants across the genome (type I error rate of 3%), and 3) in an independent dataset, the method removed 90.2% of unconfirmed genome-wide SNP associations and 89.8% of unconfirmed genome-wide indel associations. CONCLUSIONS: Researchers currently do not have effective tools to identify and mitigate batch effects in whole genome sequencing data. We developed and validated methods and filters to address this deficiency.

Two functional lupus-associated BLK promoter variants control cell-type- and developmental-stage-specific transcription.

Efforts to identify lupus-associated causal variants in the FAM167A/BLK locus on 8p21 are hampered by highly associated noncausal variants. In this report, we used a trans-population mapping and sequencing strategy to identify a common variant (rs922483) in the proximal BLK promoter and a tri-allelic variant (rs1382568) in the upstream alternative BLK promoter as putative causal variants for association with systemic lupus erythematosus. The risk allele (T) at rs922483 reduced proximal promoter activity and modulated alternative promoter usage. Allelic differences at rs1382568 resulted in altered promoter activity in B progenitor cell lines. Thus, our results demonstrated that both lupus-associated functional variants contribute to the autoimmune disease association by modulating transcription of BLK in B cells and thus potentially altering immune responses.

Genetic variants near TNFAIP3 on 6q23 are associated with systemic lupus erythematosus.

Systemic lupus erythematosus (SLE) is an autoimmune disease influenced by genetic and environmental factors. We carried out a genome-wide association scan and replication study and found an association between SLE and a variant in TNFAIP3 (rs5029939, meta-analysis P = 2.89 x 10(-12), OR = 2.29). We also found evidence of two independent signals near TNFAIP3 associated with SLE, including one previously associated with rheumatoid arthritis (RA). These results establish that variants near TNFAIP3 contribute to differential risk of SLE and RA.

Unraveling multiple MHC gene associations with systemic lupus erythematosus: model choice indicates a role for HLA alleles and non-HLA genes in Europeans.

We have performed a meta-analysis of the major-histocompatibility-complex (MHC) region in systemic lupus erythematosus (SLE) to determine the association with both SNPs and classical human-leukocyte-antigen (HLA) alleles. More specifically, we combined results from six studies and well-known out-of-study control data sets, providing us with 3,701 independent SLE cases and 12,110 independent controls of European ancestry. This study used genotypes for 7,199 SNPs within the MHC region and for classical HLA alleles (typed and imputed). Our results from conditional analysis and model choice with the use of the Bayesian information criterion show that the best model for SLE association includes both classical loci (HLA-DRB1(∗)03:01, HLA-DRB1(∗)08:01, and HLA-DQA1(∗)01:02) and two SNPs, rs8192591 (in class III and upstream of NOTCH4) and rs2246618 (MICB in class I). Our approach was to perform a stepwise search from multiple baseline models deduced from a priori evidence on HLA-DRB1 lupus-associated alleles, a stepwise regression on SNPs alone, and a stepwise regression on HLA alleles. With this approach, we were able to identify a model that was an overwhelmingly better fit to the data than one identified by simple stepwise regression either on SNPs alone (Bayes factor [BF] > 50) or on classical HLA alleles alone (BF > 1,000).

High-density SNP mapping of the HLA region identifies multiple independent susceptibility loci associated with selective IgA deficiency.

Selective IgA deficiency (IgAD; serum IgA<0.07 g/l) is the most common form of human primary immune deficiency, affecting approximately 1∶600 individuals in populations of Northern European ancestry. The polygenic nature of IgAD is underscored by the recent identification of several new risk genes in a genome-wide association study. Among the characterized susceptibility loci, the association with specific HLA haplotypes represents the major genetic risk factor for IgAD. Despite the robust association, the nature and location of the causal variants in the HLA region remains unknown. To better characterize the association signal in this region, we performed a high-density SNP mapping of the HLA locus and imputed the genotypes of common HLA-B, -DRB1, and -DQB1 alleles in a combined sample of 772 IgAD patients and 1,976 matched controls from 3 independent European populations. We confirmed the complex nature of the association with the HLA locus, which is the result of multiple effects spanning the entire HLA region. The primary association signal mapped to the HLA-DQB1*02 allele in the HLA Class II region (combined P = 7.69×10(-57); OR = 2.80) resulting from the combined independent effects of the HLA-B*0801-DRB1*0301-DQB1*02 and -DRB1*0701-DQB1*02 haplotypes, while additional secondary signals were associated with the DRB1*0102 (combined P = 5.86×10(-17); OR = 4.28) and the DRB1*1501 (combined P = 2.24×10(-35); OR = 0.13) alleles. Despite the strong population-specific frequencies of HLA alleles, we found a remarkable conservation of these effects regardless of the ethnic background, which supports the use of large multi-ethnic populations to characterize shared genetic association signals in the HLA region. We also provide evidence for the location of association signals within the specific extended haplotypes, which will guide future sequencing studies aimed at characterizing the precise functional variants contributing to disease pathogenesis.

Using gene expression to improve the power of genome-wide association analysis.

BACKGROUND/AIMS: Genome-wide association (GWA) studies have reported susceptible regions in the human genome for many common diseases and traits; however, these loci only explain a minority of trait heritability. To boost the power of a GWA study, substantial research endeavors have been focused on integrating other available genomic information in the analysis. Advances in high through-put technologies have generated a wealth of genomic data and made combining SNP and gene expression data become feasible. RESULTS: In this paper, we propose a novel procedure to incorporate gene expression information into GWA analysis. This procedure utilizes weights constructed by gene expression measurements to adjust p values from a GWA analysis. RESULTS from simulation analyses indicate that the proposed procedures may achieve substantial power gains, while controlling family-wise type I error rates at the nominal level. To demonstrate the implementation of our proposed approach, we apply the weight adjustment procedure to a GWA study on serum interferon-regulated chemokine levels in systemic lupus erythematosus patients. The study results can provide valuable insights for the functional interpretation of GWA signals. AVAILABILITY: The R source code for implementing the proposed weighting procedure is available at http://www.biostat.umn.edu/∼yho/research.html.

Lupus nephritis susceptibility loci in women with systemic lupus erythematosus.

Lupus nephritis is a manifestation of SLE resulting from glomerular immune complex deposition and inflammation. Lupus nephritis demonstrates familial aggregation and accounts for significant morbidity and mortality. We completed a meta-analysis of three genome-wide association studies of SLE to identify lupus nephritis-predisposing loci. Through genotyping and imputation, >1.6 million markers were assessed in 2000 unrelated women of European descent with SLE (588 patients with lupus nephritis and 1412 patients with lupus without nephritis). Tests of association were computed using logistic regression adjusting for population substructure. The strongest evidence for association was observed outside the MHC and included markers localized to 4q11-q13 (PDGFRA, GSX2; P=4.5×10(-7)), 16p12 (SLC5A11; P=5.1×10(-7)), 6p22 (ID4; P=7.4×10(-7)), and 8q24.12 (HAS2, SNTB1; P=1.1×10(-6)). Both HLA-DR2 and HLA-DR3, two well established lupus susceptibility loci, showed evidence of association with lupus nephritis (P=0.06 and P=3.7×10(-5), respectively). Within the class I region, rs9263871 (C6orf15-HCG22) had the strongest evidence of association with lupus nephritis independent of HLA-DR2 and HLA-DR3 (P=8.5×10(-6)). Consistent with a functional role in lupus nephritis, intra-renal mRNA levels of PDGFRA and associated pathway members showed significant enrichment in patients with lupus nephritis (n=32) compared with controls (n=15). Results from this large-scale genome-wide investigation of lupus nephritis provide evidence of multiple biologically relevant lupus nephritis susceptibility loci.

A common haplotype of interferon regulatory factor 5 (IRF5) regulates splicing and expression and is associated with increased risk of systemic lupus erythematosus.

Systemic lupus erythematosus (SLE) is a complex autoimmune disease characterized by activation of the type I interferon (IFN) pathway. Here we convincingly replicate association of the IFN regulatory factor 5 (IRF5) rs2004640 T allele with SLE in four independent case-control cohorts (P = 4.4 x 10(-16)) and by family-based transmission disequilibrium test analysis (P = 0.0006). The rs2004640 T allele creates a 5' donor splice site in an alternate exon 1 of IRF5, allowing expression of several unique IRF5 isoforms. We also identify an independent cis-acting variant associated with elevated expression of IRF5 and linked to the exon 1B splice site. Haplotypes carrying the variant associated with elevated expression and lacking the exon 1B donor site do not confer risk of SLE. Thus, a common IRF5 haplotype driving elevated expression of multiple unique isoforms of IRF5 is an important genetic risk factor for SLE, establishing a causal role for type I IFN pathway genes in human autoimmunity.

Transferability of tag SNPs in genetic association studies in multiple populations.

A general question for linkage disequilibrium-based association studies is how power to detect an association is compromised when tag SNPs are chosen from data in one population sample and then deployed in another sample. Specifically, it is important to know how well tags picked from the HapMap DNA samples capture the variation in other samples. To address this, we collected dense data uniformly across the four HapMap population samples and eleven other population samples. We picked tag SNPs using genotype data we collected in the HapMap samples and then evaluated the effective coverage of these tags in comparison to the entire set of common variants observed in the other samples. We simulated case-control association studies in the non-HapMap samples under a disease model of modest risk, and we observed little loss in power. These results demonstrate that the HapMap DNA samples can be used to select tags for genome-wide association studies in many samples around the world.

Association of NCF2, IKZF1, IRF8, IFIH1, and TYK2 with systemic lupus erythematosus.

Systemic lupus erythematosus (SLE) is a complex trait characterised by the production of a range of auto-antibodies and a diverse set of clinical phenotypes. Currently, ~8% of the genetic contribution to SLE in Europeans is known, following publication of several moderate-sized genome-wide (GW) association studies, which identified loci with a strong effect (OR>1.3). In order to identify additional genes contributing to SLE susceptibility, we conducted a replication study in a UK dataset (870 cases, 5,551 controls) of 23 variants that showed moderate-risk for lupus in previous studies. Association analysis in the UK dataset and subsequent meta-analysis with the published data identified five SLE susceptibility genes reaching genome-wide levels of significance (P(comb)<5×10(-8)): NCF2 (P(comb) = 2.87×10(-11)), IKZF1 (P(comb) = 2.33×10(-9)), IRF8 (P(comb) = 1.24×10(-8)), IFIH1 (P(comb) = 1.63×10(-8)), and TYK2 (P(comb) = 3.88×10(-8)). Each of the five new loci identified here can be mapped into interferon signalling pathways, which are known to play a key role in the pathogenesis of SLE. These results increase the number of established susceptibility genes for lupus to ~30 and validate the importance of using large datasets to confirm associations of loci which moderately increase the risk for disease.