ObjectiveGlycans attached to the Fc portion of IgG are important modulators of IgG effector functions. Interindividual differences in IgG glycome composition are large and they associate strongly with different inflammatory and autoimmune diseases. IKZF1, HLA–DQ2A/B, and BACH2 genetic loci that affect IgG glycome composition show pleiotropy with systemic lupus erythematosus (SLE), indicating a potentially causative role of aberrant IgG glycosylation in SLE. We undertook this large multicenter case–control study to determine whether SLE is associated with altered IgG glycosylation.MethodsUsing ultra‐performance liquid chromatography analysis of released glycans, we analyzed the composition of the IgG glycome in 261 SLE patients and 247 matched controls of Latin American Mestizo origin (the discovery cohort) and in 2 independent replication cohorts of different ethnicity (108 SLE patients and 193 controls from Trinidad, and 106 SLE patients and 105 controls from China).ResultsMultiple statistically significant differences in IgG glycome composition were observed between patients and controls. The most significant changes included decreased galactosylation and sialylation of IgG (which regulate proinflammatory and antiinflammatory actions of IgG) as well as decreased core fucose and increased bisecting N‐acetylglucosamine (which affect antibody‐dependent cell‐mediated cytotoxicity).ConclusionThe IgG glycome in SLE patients is significantly altered in a way that decreases immunosuppressive action of circulating immunoglobulins. The magnitude of observed changes is associated with the intensity of the disease, indicating that aberrant IgG glycome composition or changes in IgG glycosylation may be an important molecular mechanism in SLE.

Effector functions of immunoglobulin G (IgG) are regulated by the composition of a glycan moiety, thus affecting activity of the immune system. Aberrant glycosylation of IgG has been observed in many diseases, but little is understood about the underlying mechanisms. We performed a genome-wide association study of IgG N-glycosylation (N = 8090) and, using a data-driven network approach, suggested how associated loci form a functional network. We confirmed in vitro that knockdown of IKZF1 decreases the expression of fucosyltransferase FUT8, resulting in increased levels of fucosylated glycans, and suggest that RUNX1 and RUNX3, together with SMARCB1, regulate expression of glycosyltransferase MGAT3. We also show that variants affecting the expression of genes involved in the regulation of glycoenzymes colocalize with variants affecting risk for inflammatory diseases. This study provides new evidence that variation in key transcription factors coupled with regulatory variation in glycogenes modifies IgG glycosylation and has influence on inflammatory diseases.

Joint modeling of a number of phenotypes using multivariate methods has often been neglected in genome-wide association studies and if used, replication has not been sought. Modern omics technologies allow characterization of functional phenomena using a large number of related phenotype measures, which can benefit from such joint analysis. Here, we report a multivariate genome-wide association studies of 23 immunoglobulin G (IgG) N-glycosylation phenotypes. In the discovery cohort, our multi-phenotype method uncovers ten genome-wide significant loci, of which five are novel (IGH, ELL2, HLA-B-C, AZI1, FUT6-FUT3). We convincingly replicate all novel loci via multivariate tests. We show that IgG N-glycosylation loci are strongly enriched for genes expressed in the immune system, in particular antibody-producing cells and B lymphocytes. We empirically demonstrate the efficacy of multivariate methods to discover novel, reproducible pleiotropic effects.