Toma Keser

Glycans Are a Novel Biomarker of Chronological and Biological Ages

Fine structural details of glycans attached to the conserved N-glycosylation site significantly not only affect function of individual immunoglobulin G (IgG) molecules but also mediate inflammation at the systemic level. By analyzing IgG glycosylation in 5,117 individuals from four European populations, we have revealed very complex patterns of changes in IgG glycosylation with age. Several IgG glycans (including FA2B, FA2G2, and FA2BG2) changed considerably with age and the combination of these three glycans can explain up to 58% of variance in chronological age, significantly more than other markers of biological age like telomere lengths. The remaining variance in these glycans strongly correlated with physiological parameters associated with biological age. Thus, IgG glycosylation appears to be closely linked with both chronological and biological ages. Considering the important role of IgG glycans in inflammation, and because the observed changes with age promote inflammation, changes in IgG glycosylation also seem to represent a factor contributing to aging. Significance Statement Glycosylation is the key posttranslational mechanism that regulates function of immunoglobulins, with multiple systemic repercussions to the immune system. Our study of IgG glycosylation in 5,117 individuals from four European populations has revealed very extensive and complex changes in IgG glycosylation with age. The combined index composed of only three glycans explained up to 58% of variance in age, considerably more than other biomarkers of age like telomere lengths. The remaining variance in these glycans strongly correlated with physiological parameters associated with biological age; thus, IgG glycosylation appears to be closely linked with both chronological and biological ages. The ability to measure human biological aging using molecular profiling has practical applications for diverse fields such as disease prevention and treatment, or forensics.

Association of systemic lupus erythematosus with decreased immunosuppressive potential of the IgG glycome.

Glycans 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.

Changes in IgG and total plasma protein glycomes in acute systemic inflammation

Recovery after cardiac surgery is a complex process that has to compensate for both individual variability and extensive tissue damage in the context of systemic inflammation. Protein glycosylation is essential in many steps of the inflammatory cascade, but due to technological limitations the role of individual variation in glycosylation in systemic inflammation has not been addressed until now. We analysed composition of the total plasma and IgG N-glycomes in 107 patients undergoing cardiac surgery. In nearly all individuals plasma N-glycome underwent the same pattern of changes in the first 72 h, revealing a general mechanism of glycosylation changes. To the contrary, changes in the IgG glycome were very individualized. Bi-clustering analysis revealed the existence of four distinct patterns of changes. One of them, characterized by a rapid increase in galactosylated glycoforms, was associated with nearly double mortality risk measured by EuroSCORE II. Our results indicate that individual variation in IgG glycosylation changes during acute systemic inflammation associates with increased mortality risk and indicates new avenues for the development of personalized diagnostic and therapeutic approach.

Glycosylation of Immunoglobulin G: Role of Genetic and Epigenetic Influences

Objective To determine the extent to which genetic and epigenetic factors contribute to variations in glycosylation of immunoglobulin G (IgG) in humans. Methods 76 N-glycan traits in circulating IgG were analyzed by UPLC in 220 monozygotic and 310 dizygotic twin pairs from TwinsUK. A classical twin study design was used to derive the additive genetic, common and unique environmental components defining the variance in these traits. Epigenome-wide association analysis was performed using the Illumina 27k chip. Results 51 of the 76 glycan traits studied have an additive genetic component (heritability, h 2)≥ 0.5. In contrast, 12 glycan traits had a low genetic contribution (h2<0.35). We then tested for association between methylation levels and glycan levels (P<2 x10-6). Among glycan traits with low heritability probe cg08392591 maps to a CpG island 5’ from the ANKRD11 gene, a p53 activator on chromosome 16. Probe cg26991199 maps to the SRSF10 gene involved in regulation of RNA splicing and particularly in regulation of splicing of mRNA precursors upon heat shock. Among those with high heritability we found cg13782134 (mapping to the NRN1L gene) and cg16029957 mapping near the QPCT gene to be array-wide significant. The proportion of array-wide epigenetic associations was significantly larger (P<0.005) among glycans with low heritability (42%) than in those with high heritability (6.2%). Conclusions Glycome analyses might provide a useful integration of genetic and non-genetic factors to further our understanding of the role of glycosylation in both normal physiology and disease.

Glycosylation Profile of IgG in Moderate Kidney Dysfunction

Glycans constitute the most abundant and diverse form of the post-translational modifications, and animal studies have suggested the involvement of IgG glycosylation in mechanisms of renal damage. Here, we explored the associations between IgG glycans and renal function in 3274 individuals from the TwinsUK registry. We analyzed the correlation between renal function measured as eGFR and 76 N-glycan traits using linear regressions adjusted for covariates and multiple testing in the larger population. We replicated our results in 31 monozygotic twin pairs discordant for renal function. Results from both analyses were then meta-analyzed. Fourteen glycan traits were associated with renal function in the discovery sample (P<6.5×10(-4)) and remained significant after validation. Those glycan traits belong to three main glycosylation features: galactosylation, sialylation, and level of bisecting N-acetylglucosamine of the IgG glycans. These results show the role of IgG glycosylation in kidney function and provide novel insight into the pathophysiology of CKD and potential diagnostic and therapeutic targets.

Effects of statins on the immunoglobulin G glycome

BACKGROUND Statins are among the most widely prescribed medications worldwide and usually many individuals involved in clinical and population studies are on statin therapy. Immunoglobulin G (IgG) glycosylation has been associated with numerous cardiometabolic risk factors. METHODS The aim of this study was to investigate the possible association of statin use with N-glycosylation of IgG. The association was analyzed in two large population cohorts (TwinsUK and KORA) using hydrophilic interaction liquid chromatography (HILIC-UPLC) in the TwinsUK cohort and reverse phase liquid chromatography coupled with electrospray mass spectrometry (LC-ESI-MS) in the KORA cohort. Afterwards we investigated the same association for only one statin (rosuvastatin) in a subset of individuals from the randomized double-blind placebo-controlled JUPITER study using LC-ESI-MS for IgG glycome and HILIC-UPLC for total plasma N-glycome. RESULTS In the TwinsUK population, the use of statins was associated with higher levels of core-fucosylated biantennary glycan structure with bisecting N-acetylglucosamine (FA2B) and lower levels of core-fucosylated biantennary digalactosylated monosialylated glycan structure (FA2G2S1). The association between statin use and FA2B was replicated in the KORA cohort. In the JUPITER trial we found no statistically significant differences between the randomly allocated placebo and rosuvastatin groups. CONCLUSIONS In the TwinsUK and KORA cohorts, statin use was associated with a small increase of pro-inflammatory IgG glycan, although this finding was not confirmed in a subset of participants from the JUPITER trial. GENERAL SIGNIFICANCE Even if the association between IgG N-glycome and statins exists, it is not large enough to pose a problem for glycomic studies.

The Association Between Low Back Pain and Composition of IgG Glycome

Low back pain (LBP) is a common debilitating condition which aetiology and pathogenesis are poorly understood. We carried out a first so far analysis of associations between LBP and plasma IgG N-glycome in a sample of 4511 twins from TwinsUK database assessed for LBP, lumbar disc degeneration (LDD) as its possible cause, and IgG-glycan levels. Using weighted correlation network analysis, we established a correlation between LBP and glycan modules featured by glycans that either promote or block antibody-dependent cell-mediated cytotoxicity (ADCC). The levels of four glycan traits representing two of those modules were statistically significantly different in monozygotic twins discordant for LBP. Also, the trend to higher prevalence of systemic inflammatory disorders was shown for twins with low level of fucosylated glycans and high level of non-fucosylated glycans. Core fucosylation of IgG is a “safety switch” reducing ADCC, thus our results suggest the involvement of ADCC and associated inflammation in pathogenesis of LBP. No correlation between LDD scores and glycans was found assuming that the inflammation may not be a part of LDD. These data provide a new insight into understanding the complex pathophysiology of LBP and suggest glycan levels as a possible biomarker for inflammation-related subtypes of LBP.

IgG glycosylation and DNA methylation are interconnected with smoking

BACKGROUND Glycosylation is one of the most common post-translation modifications with large influences on protein structure and function. The effector function of immunoglobulin G (IgG) alters between pro- and anti-inflammatory, based on its glycosylation. IgG glycan synthesis is highly complex and dynamic. METHODS With the use of two different analytical methods for assessing IgG glycosylation, we aim to elucidate the link between DNA methylation and glycosylation of IgG by means of epigenome-wide association studies. In total, 3000 individuals from 4 cohorts were analyzed. RESULTS The overlap of the results from the two glycan measurement panels yielded DNA methylation of 7 CpG-sites on 5 genomic locations to be associated with IgG glycosylation: cg25189904 (chr.1, GNG12); cg05951221, cg21566642 and cg01940273 (chr.2, ALPPL2); cg05575921 (chr.5, AHRR); cg06126421 (6p21.33); and cg03636183 (chr.19, F2RL3). Mediation analyses with respect to smoking revealed that the effect of smoking on IgG glycosylation may be at least partially mediated via DNA methylation levels at these 7 CpG-sites. CONCLUSION Our results suggest the presence of an indirect link between DNA methylation and IgG glycosylation that may in part capture environmental exposures. GENERAL SIGNIFICANCE An epigenome-wide analysis conducted in four population-based cohorts revealed an association between DNA methylation and IgG glycosylation patterns. Presumably, DNA methylation mediates the effect of smoking on IgG glycosylation.

Network inference from glycoproteomics data reveals new reactions in the IgG glycosylation pathway.

Immunoglobulin G (IgG) is a major effector molecule of the human immune response, and aberrations in IgG glycosylation are linked to various diseases. However, the molecular mechanisms underlying protein glycosylation are still poorly understood. We present a data-driven approach to infer reactions in the IgG glycosylation pathway using large-scale mass-spectrometry measurements. Gaussian graphical models are used to construct association networks from four cohorts. We find that glycan pairs with high partial correlations represent enzymatic reactions in the known glycosylation pathway, and then predict new biochemical reactions using a rule-based approach. Validation is performed using data from a GWAS and results from three in vitro experiments. We show that one predicted reaction is enzymatically feasible and that one rejected reaction does not occur in vitro. Moreover, in contrast to previous knowledge, enzymes involved in our predictions colocalize in the Golgi of two cell lines, further confirming the in silico predictions.IgG glycosylation is an important factor in immune function, yet the molecular details of protein glycosylation remain poorly understood. The data-driven approach presented here uses large-scale plasma IgG mass spectrometry measurements to infer new biochemical reactions in the glycosylation pathway.

Publisher Correction: Network inference from glycoproteomics data reveals new reactions in the IgG glycosylation pathway

Correction to:Nature Communications (2017) 8:1231. doi:10.1038/s41467-017-01525-0