Selective immunoglobulin A deficiency (IgAD) is the most common primary immunodeficiency in Europeans. Our GWAS meta-analysis of 1,635 IgAD patients and 4,852 controls identified four new significant (P < 5x10 −8 ) loci and association with a rare IFIH1 variant (Ile923Val). Peak novel variants (PVT1 P = 4.3x10 −11 , ATG13-AMBRA1 P = 6.7x10 −10 , AHI1 P = 8.4x10 −10 and CLEC16A P = 1.4x10 −9 ) overlapped with autoimmune markers (3/4) and correlated with 21 putative regulatory variants, including eQTLs for AHI1 and DEXI and DNase hypersensitivity in FOXP3+ T regulatory cells. A pathway analysis of the meta-analysis results showed a striking association with the KEGG pathway for IgA production (pathway P < 0.0001): where 22 of 30 annotated pathway genes contained at least one variant with a P-value ≤0.05 in the IgAD metaanalysis. These data suggest that a complex network of genetic effects, including genes known to influence the biology of IgA production, contribute to IgAD.IgAD is the most common primary immunodeficiency and is defined by serum IgA level <0.07 g/L.1 The prevalence in Europeans is 1:600.1 Secretory IgA (sIgA) is important for mucosal immunity and gut commensalism2,3, and clinical features of IgAD include recurrent mucosal infections. In IgAD, B cells fail to terminally differentiate into IgA+ plasma cells, however IL-21 can restore B cell IgA production in vitro4,5. IgAD is strongly associated with HLA6 and aggregates in families with autoimmunity.7 The prevalence of Celiac disease is 35 times higher in IgAD patients whereas the prevalence of systemic lupus erythematosus (SLE) and type 1 diabetes (T1D) is 10 times higher.8Ferreira et al.6 previously imputed HLA-B, HLA-DRB1 and HLA-DQB1 in IgAD cases and controls (>2,700 individuals).9 The primary signal mapped to HLA-DQB1*02 (OR = 2.8, P = 7.7x10 −57 ), due to combined independent effects of the HLA-B*08:01-DRB1*03:01-DQB1*02 and HLA-DRB1*07:01-DQB1*02 haplotypes.9 There was a secondary signal at HLA-DRB1*01:02 (OR = 4.28, P = 5.86x10 −17 ) and a protective effect for HLA-DRB1*15:01 (OR = 0.13, P = 2.24x10 −35 ).9 HLA-DQB1*02:01 is protective for IgA nephropathy (OR = 0.71, P = 2.61x10 −13 ). 10 None of the non-MHC IgAD loci overlapped with IgA nephropathy.10 A previous GWAS identified the common allele IFIH1 Thr946Ala (OR = 0.62, control allele frequency = 0.39) as the first and, to date, only non-HLA genome-wide significant IgAD To expand our understanding of IgAD risk, we studied four new IgAD cohorts, and performed a GWAS meta-analysis of ~9.5M SNPs in 1,635 cases and 4,852 controls (Table 1). Genotypes for untyped markers were imputed for each cohort separately (1000 Genomes Project) and genotypes for variants fully typed in the entire cohort. Up to four controls per new case were iteratively selected based on ancestry eigenvectors14 to minimize population substructure (Online Methods). Genomic inflation factor values were minimal ( Table 1), indicating that population substructure was adequately addressed. Association analyses were c...