Genome-wide association studies (GWAS) with intermediate phenotypes, like changes in metabolite and protein levels, provide functional evidence to map disease associations and translate them into clinical applications. However, although hundreds of genetic variants have been associated with complex disorders, the underlying molecular pathways often remain elusive. Associations with intermediate traits are key in establishing functional links between GWAS-identified risk-variants and disease end points. Here we describe a GWAS using a highly multiplexed aptamer-based affinity proteomics platform. We quantify 539 associations between protein levels and gene variants (pQTLs) in a German cohort and replicate over half of them in an Arab and Asian cohort. Fifty-five of the replicated pQTLs are located in trans. Our associations overlap with 57 genetic risk loci for 42 unique disease end points. We integrate this information into a genome-proteome network and provide an interactive web-tool for interrogations. Our results provide a basis for novel approaches to pharmaceutical and diagnostic applications.

The Immunoglobulin G (IgG) glycome is well known for its heterogeneity and shows a significant degree of variation within populations. IgG glycome composition is influenced both by genes and by environment, making it an excellent biomarker of a person's general health state, i.e. biological age. IgG glycosylation appears to be highly regulated, both during homeostasis and in cases of its disturbance. Changes in IgG glycosylation patterns have been observed in aging and in various diseases. Differential IgG glycosylation is known to modulate IgG effector functions and is involved in disease development and progression, representing both a predisposition and a functional mechanism involved in disease pathology. This makes IgG glycosylation analysis a promising add-on to improve existing disease biomarkers.

The immunologic potency of IgG is modulated by glycosylation, but mechanisms regulating this process are undefined. A role for sex hormones is suggested by differences in IgG glycans between women and men, most prominently with respect to galactose. We therefore assessed IgG galactosylation in 713 healthy adults from 2 cohorts as well as in 159 subjects from 4 randomized controlled studies of endocrine manipulation: postmenopausal women receiving conjugated estrogens, raloxifene, or placebo; premenopausal women deprived of gonadal hormones with leuprolide and treated with estradiol or placebo; men deprived of gonadal hormones with goserelin and given testosterone or placebo; and men deprived of gonadal hormones with goserelin and given testosterone or placebo together with anastrozole to block conversion of testosterone to estradiol. Menopause was associated with an increase in agalactosylated IgG glycans, particularly in the most abundant fucosylated nonbisected (G0F) glycoform. Conjugated estrogens and raloxifene reduced G0F glycans in postmenopausal women, while in premenopausal women leuprolide increased G0F glycans in a manner reversed by estradiol. Among men, goserelin increased G0F glycans, an effect blocked by testosterone through conversion to estradiol. These results establish estrogens as an in vivo modulator of IgG galactosylation in both women and men, defining a pathway by which sex modulates immunity.

The pathogenesis of autoimmune thyroid diseases (AITD) is poorly understood and the association between different immune features and the germline variants involved in AITD are yet unclear. We previously observed systemic depletion of IgG core fucosylation and antennary α1,2 fucosylation in peripheral blood mononuclear cells in AITD, correlated with anti-thyroid peroxidase antibody (TPOAb) levels. Fucose depletion is known to potentiate strong antibody-mediated NK cell activation and enhanced target antigen-expressing cell killing. In autoimmunity, this may translate to autoantibody-mediated immune cell recruitment and attack of self-antigen expressing normal tissues. Hence, we investigated the crosstalk between immune cell traits, secreted proteins, genetic variants and the glycosylation patterns of serum IgG, in a multi-omic and cross-sectional study of 622 individuals from the TwinsUK cohort, 172 of whom were diagnosed with AITD. We observed associations between two genetic variants (rs505922 and rs687621), AITD status, the secretion of Desmoglein-2 protein, and the profile of two IgG N-glycan traits in AITD, but further studies need to be performed to better understand their crosstalk in AITD. On the other side, enhanced afucosylated IgG was positively associated with activatory CD335^- CD314⁺ CD158b⁺ NK cell subsets. Increased levels of the apoptosis and inflammation markers Caspase-2 and Interleukin-1α positively associated with AITD. Two genetic variants associated with AITD, rs1521 and rs3094228, were also associated with altered expression of the thyrocyte-expressed ligands known to recognize the NK cell immunoreceptors CD314 and CD158b. Our analyses reveal a combination of heightened Fc-active IgG antibodies, effector cells, cytokines and apoptotic signals in AITD, and AITD genetic variants associated with altered expression of thyrocyte-expressed ligands to NK cell immunoreceptors. Together, TPOAb responses, dysregulated immune features, germline variants associated with immunoactivity profiles, are consistent with a positive autoreactive antibody-dependent NK cell-mediated immune response likely drawn to the thyroid gland in AITD.