Skip to main content

Advertisement

The changes of immunoglobulin G N-glycosylation in blood lipids and dyslipidaemia

Abstract

Background

Alternative N-glycosylation has significant structural and functional consequences on immunoglobulin G (IgG) and can affect immune responses, acting as a switch between pro- and anti-inflammatory IgG functionality. Studies have demonstrated that IgG N-glycosylation is associated with ageing, body mass index, type 2 diabetes and hypertension.

Methods

Herein, we have demonstrated patterns of IgG glycosylation that are associated with blood lipids in a cross-sectional study including 598 Han Chinese aged 20–68 years. The IgG glycome composition was analysed by ultra-performance liquid chromatography.

Results

Blood lipids were positively correlated with glycan peak GP6, whereas they were negatively correlated with GP18 (P < 0.05/57). The canonical correlation analysis indicated that initial N-glycan structures, including GP4, GP6, GP9-12, GP14, GP17, GP18 and GP23, were significantly correlated with blood lipids, including total cholesterol, total triglycerides (TG) and low-density lipoprotein (r = 0.390, P < 0.001). IgG glycans patterns were able to distinguish patients with dyslipidaemia from the controls, with an area under the curve of 0.692 (95% confidence interval 0.644–0.740).

Conclusions

Our findings indicated that a possible association between blood lipids and the observed loss of galactose and sialic acid, as well as the addition of bisecting GlcNAcs, which might be related to the chronic inflammation accompanying with the development and procession of dyslipidaemia.

Background

Glycosylation is an essential posttranslational modification of proteins. It is an integral part of proteins and significantly contributes to their structure and function [1,2,3]. Immunoglobulin G (IgG) plays an important role in the human immune system, and N-glycans attach to the conserved asparagine 297 in the fragment crystallizable (Fc) part of this molecule and act as a switch between pro- and anti-inflammatory IgG functionality [1, 4,5,6]. IgG N-glycosylation impacts the physiology and malfunction of the immune system, and aberrant IgG N-glycosylation is involved in several inflammatory and chronic diseases [7,8,9,10]. Previous studies have revealed that genetic loci associated with variation in IgG glycosylation are also known risk factors for several inflammatory diseases and chronic diseases [11, 12], indicating that IgG glycosylation is not merely a result of complicated enzymatic activities, but is a subtly regulated outcome designed to meet dominant physiological needs.

It is widely recognized that dyslipidaemia is associated with an increased risk of coronary artery disease, stroke, and heart failure [13,14,15,16]. However, the mechanisms by which dyslipidaemia increase the risk of these diseases has not been completely elucidated. Importantly, dyslipidaemia is responsive to changes in systemic inflammation, a recognized causal pathway of cardiovascular and cerebrovascular diseases [17,18,19,20,21,22]. Epidemiological evidence has indicated that blood lipids are associated with several inflammatory markers, such as cytokines and chemokines [17, 23]. Due to the inflammatory role of IgG N-glycosylation, together with its association with ageing [24], obesity [25], type 2 diabetes [26], hypertension [27, 28], ischemic stroke [29], Parkinson’s disease [30] and cancer [31], it can be speculated that blood lipids may be associated with the N-glycosylation of IgG.

In this study, we aimed to determine the association between the IgG N-glycome and blood lipids, including total cholesterol (TC), total triglycerides (TG), high-density lipoprotein (HDL), and low-density lipoprotein (LDL), and explored the use of IgG glycans as biomarkers to diagnose dyslipidaemia to further investigate the role of IgG glycosylation in dyslipidaemia.

Methods

Ethical approval

Written informed consent was obtained from each subject at the beginning of the study, and the study has been approved by the Ethics Committee of the Capital Medical University, Beijing, China. The ethics approval was given in compliance with the Declaration of Helsinki [32].

Participant recruitment

An observational cross-sectional study was conducted during 2012, with a total of 913 participants of Chinese Han ancestry who were recruited from a community-based survey in Beijing [24, 27]. All participants were required to meet the following inclusion criteria: (1) age equal to or greater than 18 years; (2) no medication history during the previous 2 weeks; and (3) signed informed consent prior to participation. Individuals were excluded based on the following criteria: (1) pregnant or lactating women; (2) history of mental illness or infectious disease; and (3) history of malignant tumor, stroke or other cerebrovascular disease, congenital heart disease, acute myocardial infarction, liver disease, renal failure, chronic obstructive pulmonary disease, rheumatoid arthritis or other disease.

Diagnosis of dyslipidaemia

According to the Chinese guidelines for the prevention and control of dyslipidaemia in adults [33], the participants were grouped into cases with TC ≥ 6.2 mmol/L (defined as hypercholesteremia), TG ≥ 2.3 mmol/L (defined as hypertriglyceridemia), HDL < 1.0 mmol/L (defined as the decreased HDL), or LDL ≥ 4.1 mmol/L (defined as the elevated LDL). The controls were those who did not have dyslipidaemia at the time of the study.

Covariates

Each participant was required to complete a thorough health examination, which included anthropometric measurements and an analysis of the physical and chemical properties of a blood sample. Detailed information about subject recruitment and enrolment has been previously described [24, 27]. As in our previous study [24], 12 clinical traits (age, sex, body mass index (BMI), waist–hip ratio (WHR), systolic blood pressure (SBP), diastolic blood pressure (DBP), fasting blood triglycerides (FBG), resting heart rate (RHR), TC, TG, HDL and LDL) were assessed in the current study [24].

To obtain the anthropometric trait data, physical examinations and interviews were carried out by trained nurses and physicians. The BMI was calculated by the formula weight (in kilograms)/height2 (in metres squared). The SBP and DBP were measured twice on the right arm using a standard mercury sphygmomanometer after the subjects had rested for at least 10 min in a sitting position [34]. For the blood parameters (haematology and biochemical traits), fasting blood samples were collected in the morning after an overnight fast by venepuncture. All collected blood samples were processed within 8 h and blood plasma which was separated from whole blood was stored in freezers (− 80 °C) before the measurements. Haematology and biochemical parameters were measured by standard blood chemistry and haematology assays (Hitachi Automatic Analyzer, Model-7600, Tokyo, Japan), which were conducted at the Medical Laboratory of Beijing Xuanwu Hospital, Beijing.

IgG N-glycan analysis

As previously reported, IgG was first isolated from human plasma [24, 35]. The IgG N-glycans were cleaved and analysed by hydrophilic interaction chromatography (HILIC)- ultra-performance liquid chromatography (UPLC) as previously described [35]. In brief, after washing and equilibrating the protein G monolithic plates, 50 µL of plasma was diluted 10× with binding buffer (1× phosphate-buffered saline, pH 7.4), applied to the protein G plates, and immediately washed. The IgGs were eluted with 1 mL of 0.1 M formic acid and immediately neutralized with 1 M ammonium bicarbonate. Next, the IgG N-glycans cleavage and labelling were performed. The released N-glycans were labelled with 2-aminobenzamide, a fluorescent dye used to visualize glycans by UPLC, by multistage mixing with 2-aminobenzamide, dimethyl sulfoxide, glacial acetic acid, and 2-picoline borane. Finally, the IgG N-glycans were analysed by HILIC–UPLC in 24 IgG glycan peaks (GPs). The glycan structures of the glycans per peak were previously reported, and a detailed description is shown in Additional file 1: Table S1 [35]. The normalization and batch correction of the UPLC data was detailed in a previous study [24, 35]. The total area normalization was applied where the amount of glycans in each peak was expressed as a percentage of the total integrated area. The calculation formula is shown in Additional file 1: Table S1.

Among the 24 initial GPs, GP3 was excluded as it did not pass quality control tests. An additional 54 derived traits were calculated for the remaining 23 initial glycans [35]. Of the derived traits, 34 were of interest and were included in the analysis describing the relative abundances of galactosylation, sialylation, bisecting GlcNAc and core fucosylation (Additional file 2: Table S2).

Statistical analysis

The normal distribution of glycans was tested by the Kolmogorov–Smirnov test where P < 0.10 was considered to be statistically significant. Continuous variables underlying the normal distribution were represented as the mean ± standard deviation (SD), otherwise the medians (interquartile ranges) were used. The difference of continuous variables between two groups was tested using Student’s t test or the Wilcoxon rank-sum test. Categorical variables were represented as n (proportion), and the between-group differences were tested by the Chi square test. For the normally distributed variables, including age, BMI, WHR, SBP, DBP, FBG and RHR, t-tests were used to compare the differences between two groups. A generalized linear model (GLM) [36] was fitted to estimate the effect of blood lipids on glycosylation when adjusting for confounding factors, including age, sex, BMI, WHR, SBP, DBP, FBG and RHR, because most of the glycans (dependent variables) were not normally distributed. There were internal associations among the independent variables, which could have induced multi-collinearity in the statistical models. Therefore, principal component analysis (PCA) [37] was used to combine variables to principal components before the regression analysis. Moreover, these principal components were used in the multiple variables association analysis instead of the original variables. The internal associations among the independent variables were analysed using by the R package “corrplot” [38].

Before correlation and regression analysis, the z-score of normalized transformations for IgG N-glycans was applied to add the consistent comparability among IgG N-glycans. Canonical correlation analysis (CCA) [40] was used to determine two sets of variables of the initial glycan structures (x) and the blood lipids (y) and to find the overall correlation between the two sets of variables. The identified variables with a statistically significant impact on the canonical variables were judged by the canonical loadings. Generally, absolute values greater than 0.30 are considered to be significant loadings [39].

To explore IgG N-glycans as biomarkers to diagnose dyslipidaemia, 23 initial glycans were used as predictors. Logistic regression was performed to identify the 23 glycans that were respectively related to cases’ status after adjusting for the effect of age, sex, BMI, WHR, SBP, DBP, FBG and RHR. Considering the internal associations among these variables, we chose to correct sex, Prin.1 and Prin.2 instead of the original variables. The significant glycans were presented as the odds ratio (OR) with a 95% confidence interval (CI) and were graphed by a forest plot, which was generated using the R package “forestplot” [40]. Next, a classification model was established based on the significant glycans. Considering that there were internal associations in the glycans that might have induced multicollinearity in the model, the least absolute shrinkage and selection operator (LASSO) method was used to select glycans to reduce the dimension of data. The LASSO method was carried out by the R package “lars” [41]. The logistic regression model was then used to assess the discrimination of dyslipidaemia by combining the glycans left by the LASSO method. The classification model was only used to assess the ability of IgG N-glycans to be used as diagnostic biomarkers for disease; therefore, in the classification model there was no adjustment for the effect of age, sex, BMI, WHR, SBP, DBP, FBG and RHR. A receiver operating characteristic (ROC) curve was developed for the calculation of the area under the cure (AUC) with a 95% CI.

Data analysis was performed using SPSS Statistics version 21.0 for Windows (IBM Corp., Armonk, NY, USA), SAS software version 9.2 (SAS Institute, Chicago, IL, USA) and R version 3.3.2 (R Core Team 2016). All reported P values were two-sided, and P < 0.05 was considered statistically significant. Particularly, P < 0.05/57 (0.0009) was considered to be statistically significant in the analysis of glycans to correct for multiple comparisons.

Results

Description of clinical traits

Among the 913 recruited participants, 598 (64.50%) had a complete dataset of blood measurements and glycan traits and were included in the further analysis. In total, 8 clinical traits were described and compared between the dyslipidaemia and control groups (Table 1). Most of the traits (6 out of 8) were significantly different between the two groups (P < 0.05). The number of males and the values of age, BMI, WHR, FBG, and SBP in the dyslipidaemia group were significantly higher than those in the control group.

Table 1 Demographic and biochemical characteristics of the participants

Principal component analysis (PCA)

As shown in Additional file 3: Figure S1, the significant correlation coefficients in independent variables ranged from 0.08 to 0.89. PCA was used to combine a few variables to several principal components in order to reduce the multicollinearity of independent variables. The resulting Prin. 1 and Prin. 2 explained up to 54.78% of variance in SBP, DBP, FBG, RHR, age, BMI and WHR (Prin.1 = 0.773 × SBP + 0.783 × DBP + 0.575 × FBG + 0.106 × RHR + 0.352 × Age + 0.688 × BMI + 0.733 × WHR, Prin.2 = 0.361 × SBP + 0.364 × DBP − 0.176 × FBG + 0.739 × RHR − 0.370 × Age − 0.246 × BMI − 0.329 × WHR). The resulting Prin. 3 and Prin. 4 explained up to 85.65% of variance in TC, TG, HDL and LDL (Prin.3 = 0.985 × TC + 0.134 × TG + 0.444 × HDL + 0.923 × LDL, Prin.4 = 0.133 × TC + 0.889 × TG − 0.755 × HDL + 0.092 × LDL). Therefore, in the analysis of the association between IgG N-glycans and blood lipids, the principal components of blood lipids (Prin. 3 and Prin. 4) were used as independent variables instead of blood lipids, and they were adjusted for confounding factors (sex, Prin. 1 and Prin. 2).

The association of IgG N-glycans with blood lipids

The results of the normal distribution tests of glycans are shown in Additional file 4: Table S3. Of the 23 directly measured glycans, 16 were not abnormally distributed (P < 0.10). Therefore, a GLM was used to estimate the effect of the principal components of blood lipids (Prin. 3 and Prin. 4) on glycosylation, adjusting for confounding factors (sex, Prin. 1 and Prin. 2). In total, 2 GPs in the IgG N-glycans were significantly associated with blood lipids after adjusting for the effects of sex, age, BMI, WHR, FBG, RHR, SBP and DBP (P < 0.05/57). As shown in Table 2, TC, TG and LDL were positively correlated with GP6, while they were negatively correlated with GP18. In addition, HDL negatively correlated with GP6 but positively correlated with GP18. Furthermore, 9 initial traits and 7 derived traits in the IgG glycome were significantly associated with blood lipids after adjusting for the effects of sex, age, BMI, WHR, FBG, RHR, SBP and DBP (P < 0.05). Detailed information is shown in the supplemental content (Additional file 5: Table S4).

Table 2 Associations between IgG glycan and blood lipids

Multivariate analyses by canonical correlation analysis (CCA)

The results of CCA showed that there are 4 pairs of canonical variables, with canonical correlations of 0.390 (F = 2.44, P < 0.001), 0.324 (F = 1.88, P < 0.001), 0.240 (F = 1.39, P = 0.053) and 0.197 (F = 1.16, P = 0.282) for each successive pair by the CCA. The first and the second canonical sets were statistically significant, indicating that N-glycan structures were significantly correlated with the blood lipids. As shown in Fig. 1, 10 initial traits (GP4, GP6, GP9-12, GP14, GP17-18 and GP23) tended to be significantly associated with TC, TG and LDL levels in the first canonical set. In addition, the level of GP11 was strongly associated with canonical variables with a loading of 0.645, while the response variable with the highest canonical loading was 0.781 (TC).

Fig. 1
figure1

Canonical structures of the normalized IgG N-glycan and blood lipids in the first canonical set. The absolute value of canonical loadings greater than 0.30 was significant loadings. All of the variables are sorted by the absolute value of their canonical loadings. The positive relationships are represented in black boxes, while negative relationships are showed in red boxes. TC total cholesterol, TG total triglycerides, HDL high-density lipoprotein, LDL low-density lipoprotein

Classification of dyslipidaemia using IgG N-glycans

Among the 23 GPs, 9 GPs (GP1, GP4, GP5, GP6, GP11, GP14, GP18, GP20 and GP21) were found to be significantly different in participants with dyslipidaemia than in controls after adjusting for confounding factors (sex, Prin. 1 and Prin. 2) (Fig. 2 and Additional file 6: Table S5). As shown in Additional file 7: Figure S2, the significant correlation coefficients in independent variables ranged from 0.11 to 0.79. The LASSO method was used to select IgG N-glycans in order to reduce the dimension of data. The remaining glycans used by the LASSO method are listed in Additional file 8: Table S6. The classification model, incorporating significant IgG N-glycans, was able to distinguish dyslipidaemia from controls, and the AUC was 0.692 (95% CI 0.644–0.740) (Fig. 3).

Fig. 2
figure2

Odds ratios (OR) and 95% confidence intervals (95% CI) for the associations of the normalized glycan variables in dyslipidaemia vs controls (adjusted for sex, Prin.1 and Prin.2)

Fig. 3
figure3

Receiver operating characteristic (ROC) curve analysis in regard to binary logistic regression in the prediction of dyslipidaemia. AUC area under the cure; GP4, GP6, GP14, GP18, GP20 and GP21 included in the final model

Discussion

Dyslipidaemia is one of the most important risk factors for atherosclerosis, which triggers the development of various cardiovascular and cerebrovascular diseases [13,14,15,16]. The increased TC, TG and LDL and the decreased HDL could be major public health problems; therefore, the identification of biomarkers that might provide new avenues for the prevention and treatment of dyslipidaemia is urgently needed. To our knowledge, this is the first study to investigate the association of IgG N-glycans with blood lipids and dyslipidaemia.

In the present study, the levels of TC, TG and LDL were positively associated with the levels of GP4 and GP6, while they were negatively correlated with the level of GP18. Furthermore, the results were consistent with those in CCA, in which N-glycan structures were related with blood lipids. GP4 and GP6 are agalactosylated glycans and share the derived trait G0n, which contains nearly the all glycans without galactose. In parallel, GP18 contains 2 galactoses, which was negatively correlated with blood lipids. In addition, GP6 that contains a bisecting GlcNAc was positively associated with blood lipids, while a significant negative correlation existed between blood lipids and GP18, in which 1 sialic acid was found. This is consistent with the results of the association between the derived traits in the IgG glycome and blood lipids. Our findings indicated that the associations between the blood lipids and IgG glycome were independently significant, with a negative association with diagalactosylation and sialylation and a positive association with bisecting agalactosylation and N-acetylgucosylation. Furthermore, the loss of galactose and sialic acid and the addition of a bisecting GlcNAc were consistently observed in patients with dyslipidaemia. The changes of IgG N-glycosylation in dyslipidaemia are consistent with the results of type 2 diabetes, hypertension and ischemic stroke [26, 27, 29, 42].

N-Glycans that are attached to the Fc portion of IgG are important modulators of IgG effector functions [5], and even slight changes in IgG glycosylation can direct pro- and anti-inflammatory actions of immunoglobulins [1, 4]. The decreased IgG galactosylation in rheumatoid arthritis was first reported nearly 32 years ago [43], and recently, the same type of change in IgG glycosylation has been reported in a number of autoimmune and inflammatory diseases [7, 9, 10, 12]. In addition, IgG glycosylation correlates with ageing, obesity, hypertension and cancer [24, 25, 27, 31]. It can be speculated that the loss of galactose is not disease specific but is a general phenomenon that is associated with reducing the anti-inflammatory function of circulating IgG. The addition of sialic acid dramatically changes the physiological role of IgGs, converting them from pro-inflammatory to anti-inflammatory agents [6, 44]. In the interactions between different functional elements of IgG glycosylation, it has been shown that the presence of bisecting GlcNAcs reduce galactosylated IgGs [6, 45, 46]. Accumulating evidence indicates that changes in IgG glycosylation mechanism could be part of the molecular mechanism leading to the promotion of inflammation [6, 8, 44, 47]; therefore, our results indicate inflammation as one of characteristics in dyslipidaemia that can increase the risk of developing other related disorders. This hypothesis is supported by our findings that IgG N-glycosylation is significantly associated with blood lipids and dyslipidaemia.

Previous studies have shown that the inflammatory role of IgG N-glycosylation is associated with the risk factors of dyslipidaemia [48] including aging and obesity [24, 25]. Inflammation, which is characteristic of aging and obesity, is a process associated with different disorders, such as dyslipidaemia [49]. Inflammatory disorders can lead to the activation of several signalling transduction pathways, inflammatory cytokine chemokine production and cell migration, all of which can influence lipid metabolism [49,50,51,52]. Interleukin 6 (IL-6) has been implicated in the pathogenesis of several immune-mediated diseases, and monoclonal antibodies directed against the IL-6 receptor have been developed to treat different inflammatory diseases [53]. IL-6 and interleukin 6 signal transducer (IL6ST) have been identified as dyslipidaemia susceptibility loci [54, 55]. IL6ST has been identified to have the capability to regulate galactosylation of IgG in a genome-wide association study of IgG N-glycosylation [12]. According to the law of Mendelian randomization [56, 57], the genetic loci IL6ST is associated with IgG N-glycosylation and has an effect on dyslipidaemia, indicating that the change of galactosylation and inflammation is the intermediate phenotype between IL6ST and dyslipidaemia and is thus the causal factor of dyslipidaemia. However, this will need to be further explored and validated.

The aberrant glycosylation which induces inflammation may provide exciting insights into the pathogenesis of dyslipidaemia. However, causation is difficult to verify, and the observed changes may be the consequence rather than cause of the disease. Dyslipidaemia accompanied by adipocytes that can produce and secrete cytokines and adipokines [58] may thus affect the structural integrity of IgG glycans, where IgG plays a crucial role in the activation of complement, interacts with Fc receptors and affects antibody-dependent cell mediated cytotoxicity (ADCC) [5]. Dyslipidaemia, as a basic metabolic disease, may trigger changes in IgG glycosylation accompanied by inflammation that can lead to related diseases. Therefore, the casual effect between IgG glycosylation and dyslipidaemia remains unclear.

There are several potential limitations in this study that should be recognized. First, the study was performed with a relatively small sample population. Multiple correction was not used when we selected the initial glycans as diagnostic biomarkers, which may have led to false positive errors. However, to overcome this, the method of reducing the dimensions of data, including PCA and CCA, were applied to examine the association between level of N-glycans and blood lipids. The merit of PCA is that it explores all of the information from principal components, while the merit of CCA is to explore all of the information from two sets of variables; therefore, PCA and CCA help to solve the problem of multicollinearity that is induced by the similarities of the variables. Second, our study is a cross-sectional study, bringing the bias of diagnosis of dyslipidaemia cases. It lacks information regarding the time sequence of events; therefore, we cannot conclude the causal relationship of IgG glycosylation and dyslipidaemia. In addition, the present study is a pilot study to explore the association between N-glycans and blood lipids and dyslipidaemia. Further cohort studies or Mendelian randomization studies [57, 58] with larger sample sizes are needed to provide a more definite explanation about the relationships between N-glycan structures and dyslipidaemia.

Conclusion

In conclusion, the present study showed a possible association between blood lipids and the observed loss of galactose and sialic acid, as well as the addition of bisecting GlcNAcs which might be related to the chronic inflammation accompanying the development of dyslipidaemia. Moreover, IgG N-glycosylation profiles may serve as potential biomarkers for dyslipidaemia, contributing to a move towards personalized medicine. Future studies focused on the mechanisms underlying the causal association of IgG glycosylation and blood lipids or dyslipidaemia are needed.

Abbreviations

ADCC:

antibody dependent cell-mediated cytotoxicity

AUC:

area under the curve

BMI:

body mass index

CCA:

canonical correlation analysis

CI:

confidence intervals

DBP:

diastolic blood pressure

Fc:

fragment crystallizable

FBG:

fasting blood glucose

GLM:

generalized linear model

GP:

glycan peak

GPs:

glycan peaks

HDL:

high-density lipoprotein

HILIC:

hydrophilic interaction chromatography

IgG:

immunoglobulin G

IL-6:

interleukin 6

IL6ST:

interleukin 6 signal transducer

LASSO:

least absolute shrinkage and selection operator

LDL:

low-density lipoprotein

OR:

odds ratio

PCA:

principal component analysis

ROC:

receiver operating characteristic

SBP:

systolic blood pressure

SD:

standard deviation

TC:

total cholesterol

TG:

triglycerides

UPLC:

ultra performance liquid chromatography

WHO:

World Health Organization

WHR:

waist–hip ratio

References

  1. 1.

    Rabinovich GA, van Kooyk Y, Cobb BA. Glycobiology of immune responses. Ann N Y Acad Sci. 2012;1253:1–15.

  2. 2.

    Lauc G, Pezer M, Rudan I, Campbell H. Mechanisms of disease: the human N-glycome. Biochim Biophys Acta. 2016;1860:1574–82.

  3. 3.

    Kristic J, Lauc G. Ubiquitous importance of protein glycosylation. Methods Mol Biol. 2017;1503:1–12.

  4. 4.

    Keser T, Vuckovic F, Barrios C, Zierer J, Wahl A, Akinkuolie AO, et al. Effects of statins on the immunoglobulin G glycome. Biochim Biophys Acta. 2017;1861:1152–8.

  5. 5.

    Hayes JM, Cosgrave EF, Struwe WB, Wormald M, Davey GP, Jefferis R, et al. Glycosylation and Fc receptors. Curr Top Microbiol Immunol. 2014;382:165–99.

  6. 6.

    Biermann MH, Griffante G, Podolska MJ, Boeltz S, Sturmer J, Munoz LE, et al. Sweet but dangerous—the role of immunoglobulin G glycosylation in autoimmunity and inflammation. LUPUS. 2016;25:934–42.

  7. 7.

    Vuckovic F, Kristic J, Gudelj I, Teruel M, Keser T, Pezer M, et al. Association of systemic lupus erythematosus with decreased immunosuppressive potential of the IgG glycome. Arthritis Rheumatol. 2015;67:2978–89.

  8. 8.

    Novokmet M, Lukic E, Vuckovic F, Ethuric Z, Keser T, Rajsl K, et al. Changes in IgG and total plasma protein glycomes in acute systemic inflammation. Sci Rep. 2014;4:4347.

  9. 9.

    Sebastian A, Alzain MA, Asweto CO, Song H, Cui L, Yu X, et al. Glycan biomarkers for rheumatoid arthritis and its remission status in Han Chinese patients. OMICS. 2016;20:343–51.

  10. 10.

    Trbojevic AI, Ventham NT, Theodoratou E, Vuckovic F, Kennedy NA, Kristic J, et al. Inflammatory bowel disease associates with proinflammatory potential of the immunoglobulin G glycome. Inflamm Bowel Dis. 2015;21:1237–47.

  11. 11.

    Zoldos V, Horvat T, Lauc G. Glycomics meets genomics, epigenomics and other high throughput omics for system biology studies. Curr Opin Chem Biol. 2013;17:34–40.

  12. 12.

    Lauc G, Huffman JE, Pucic M, Zgaga L, Adamczyk B, Muzinic A, et al. Loci associated with N-glycosylation of human immunoglobulin G show pleiotropy with autoimmune diseases and haematological cancers. PLoS Genet. 2013;9:e1003225.

  13. 13.

    Abbasi F, Kohli P, Reaven GM, Knowles JW. Hypertriglyceridemia: a simple approach to identify insulin resistance and enhanced cardio-metabolic risk in patients with prediabetes. Diabetes Res Clin Pract. 2016;120:156–61.

  14. 14.

    Kushner PA, Cobble ME. Hypertriglyceridemia: the importance of identifying patients at risk. Postgrad Med. 2016;128:848–58.

  15. 15.

    Maki KC, Guyton JR, Orringer CE, Hamilton-Craig I, Alexander DD, Davidson MH. Triglyceride-lowering therapies reduce cardiovascular disease event risk in subjects with hypertriglyceridemia. J Clin Lipidol. 2016;10:905–14.

  16. 16.

    Gofman JW, Lindgren F. The role of lipids and lipoproteins in atherosclerosis. Science. 1950;111:166–86.

  17. 17.

    Zheng J, Wu J, Chen J, Liu J, Lu Y, Huang C, et al. Therapeutic effects of quercetin on early inflammation in hypertriglyceridemia-related acute pancreatitis and its mechanism. Pancreatology. 2016;16:200–10.

  18. 18.

    Aravindhan V, Madhumitha H. Metainflammation in diabetic coronary artery disease: emerging role of innate and adaptive immune responses. J Diabetes Res. 2016;2016:6264149.

  19. 19.

    Li PY, Wang X, Stetler RA, Chen J, Yu WF. Anti-inflammatory signaling: the point of convergence for medical gases in neuroprotection against ischemic stroke. Med Gas Res. 2016;6:227–31.

  20. 20.

    Miles EA, Rees D, Banerjee T, Cazzola R, Lewis S, Wood R, et al. Age-related increases in circulating inflammatory markers in men are independent of BMI, blood pressure and blood lipid concentrations. Atherosclerosis. 2008;196:298–305.

  21. 21.

    Lubrano V, Del TS, Nicolini G, Di Cecco P, Basta G. Circulating levels of lectin-like oxidized low-density lipoprotein receptor-1 are associated with inflammatory markers. Lipids. 2008;43:945–50.

  22. 22.

    Lockyer S, Rowland I, Spencer JP, Yaqoob P, Stonehouse W. Impact of phenolic-rich olive leaf extract on blood pressure, plasma lipids and inflammatory markers: a randomised controlled trial. Eur J Nutr. 2016;56:1421–32.

  23. 23.

    Li B, Li W, Li X, Zhou H. Inflammation: a novel therapeutic target/direction in atherosclerosis. Curr Pharm Des. 2016;23:1216–27.

  24. 24.

    Yu X, Wang Y, Kristic J, Dong J, Chu X, Ge S, et al. Profiling IgG N-glycans as potential biomarker of chronological and biological ages: a community-based study in a Han Chinese population. Medicine. 2016;95:e4112.

  25. 25.

    Nikolac PM, Pucic BM, Kristic J, Novokmet M, Huffman JE, Vitart V, et al. The association between galactosylation of immunoglobulin G and body mass index. Prog Neuropsychopharmacol Biol Psychiatry. 2014;48:20–5.

  26. 26.

    Lemmers R, Vilaj M, Urda D, Agakov F, Simurina M, Klaric L, et al. IgG glycan patterns are associated with type 2 diabetes in independent European populations. Biochim Biophys Acta. 2017;1861:2240–9.

  27. 27.

    Wang Y, Klaric L, Yu X, Thaqi K, Dong J, Novokmet M, et al. The association between glycosylation of immunoglobulin G and hypertension: a multiple ethnic cross-sectional study. Medicine. 2016;95:e3379.

  28. 28.

    Gao Q, Dolikun M, Stambuk J, Wang H, Zhao F, Yiliham N, et al. Immunoglobulin G N-glycans as potential postgenomic biomarkers for hypertension in the Kazakh population. Omics. 2017;21:380–9.

  29. 29.

    Liu D, Zhao Z, Wang A, Ge S, Wang H, Zhang X, et al. Ischemic stroke is associated with the pro-inflammatory potential of N-glycosylated immunoglobulin G. J Neuroinflamm. 2018;15:123.

  30. 30.

    Russell AC, Simurina M, Garcia MT, Novokmet M, Wang Y, Rudan I, et al. The N-glycosylation of immunoglobulin G as a novel biomarker of Parkinson’s disease. Glycobiology. 2017;27:501–10.

  31. 31.

    Ren S, Zhang Z, Xu C, Guo L, Lu R, Sun Y, et al. Distribution of IgG galactosylation as a promising biomarker for cancer screening in multiple cancer types. Cell Res. 2016;26:963–6.

  32. 32.

    Association WM. World Medical Association Declaration of Helsinki: ethical principles for medical research involving human subjects. J Am Med Assoc. 2004;15:124.

  33. 33.

    Zhao SP. Key points and comments on the 2016 Chinese guideline for the management of dyslipidaemia in adults. Zhonghua Xin Xue Guan Bing Za Zhi. 2016;44:827.

  34. 34.

    WHO. A global brief on hypertension. Geneva: World Health Organization; 2013.

  35. 35.

    Pucic M, Knezevic A, Vidic J, Adamczyk B, Novokmet M, Polasek O, et al. High throughput isolation and glycosylation analysis of IgG-variability and heritability of the IgG glycome in three isolated human populations. Mol Cell Proteomics. 2011;10:M111–10090.

  36. 36.

    Zheng B, Agresti A. Summarizing the predictive power of a generalized linear model. Stat Med. 2000;19:1771–81.

  37. 37.

    Shen D, Zhu H. Spatially weighted principal component regression for high-dimensional prediction. Inf Process Med Imaging. 2015;24:758–69.

  38. 38.

    Hoefle AS, Bangert AM, Stamfort A, Gedrich K, Rist MJ, Lee YM, et al. Metabolic responses of healthy or prediabetic adults to bovine whey protein and sodium caseinate do not differ. J Nutr. 2015;145:467–75.

  39. 39.

    Wilms I, Croux C. Robust sparse canonical correlation analysis. BMC Syst Biol. 2016;10:72.

  40. 40.

    Gordon M, Lumley T. Forestplot: advanced forest plot using ‘grid‘ graphics. http://CRAN.R-project.org/package=forestplot. 2017.

  41. 41.

    Tibshirani R. Regression shrinkage selection via the lasso. J R Stat Soc. 2011;73:273–82.

  42. 42.

    Ercan A, Kohrt WM, Cui J, Deane KD, Pezer M, Yu EW, et al. Estrogens regulate glycosylation of IgG in women and men. JCI Insight. 2017;2:e89703.

  43. 43.

    Parekh RB, Dwek RA, Sutton BJ, Fernandes DL, Leung A, Stanworth D, et al. Association of rheumatoid arthritis and primary osteoarthritis with changes in the glycosylation pattern of total serum IgG. Nature. 1985;316:452–7.

  44. 44.

    Bohm S, Schwab I, Lux A, Nimmerjahn F. The role of sialic acid as a modulator of the anti-inflammatory activity of IgG. Semin Immunopathol. 2012;34:443–53.

  45. 45.

    van de Geijn FE, Wuhrer M, Selman MH, Willemsen SP, de Man YA, Deelder AM, et al. Immunoglobulin G galactosylation and sialylation are associated with pregnancy-induced improvement of rheumatoid arthritis and the postpartum flare: results from a large prospective cohort study. Arthritis Res Ther. 2009;11:R193.

  46. 46.

    Hanzawa K, Suzuki N, Natsuka S. Structures and developmental alterations of N-glycans of zebrafish embryos. Glycobiology. 2016;27(3):228–45.

  47. 47.

    de Jong SE, Selman MH, Adegnika AA, Amoah AS, van Riet E, Kruize YC, et al. IgG1 Fc N-glycan galactosylation as a biomarker for immune activation. Sci Rep. 2016;6:28207.

  48. 48.

    Jakicic JM, Donnelly JE, Jawad AF, Jacobsen DJ, Gunderson SC, Pascale R. Association between blood lipids and different measures of body fat distribution: effects of BMI and age. Int J Obes Relat Metab Disord. 1993;17:131.

  49. 49.

    Esteve E, Ricart W, Fernandez-Real JM. Dyslipidaemia and inflammation: an evolutionary conserved mechanism. Clin Nutr. 2005;24:16–31.

  50. 50.

    Bag-Ozbek A, Giles JT. Inflammation, adiposity, and atherogenic dyslipidaemia in rheumatoid arthritis: is there a paradoxical relationship? Curr Allergy Asthma Rep. 2015;15:497.

  51. 51.

    Rizzo M, Rizvi AA, Rini GB, Berneis K. The therapeutic modulation of atherogenic dyslipidaemia and inflammatory markers in the metabolic syndrome: what is the clinical relevance? Acta Diabetol. 2009;46:1–11.

  52. 52.

    Tang L, Peng H, Xu T, Wang A, Wang G, Tong W, et al. Association of biomarkers of inflammation with dyslipidaemia and its components among Mongolians in China. PLoS ONE. 2014;9:e89023.

  53. 53.

    Yao X, Huang J, Zhong H, Shen N, Faggioni R, Fung M, et al. Targeting interleukin-6 in inflammatory autoimmune diseases and cancers. Pharmacol Therapeut. 2014;141:125–39.

  54. 54.

    Gottardo L, Cosmo SD, Zhang YY, Powers C, Prudente S, Marescotti MC, et al. A polymorphism at the IL6ST (gp130) locus is associated with traits of the metabolic syndrome. Obesity. 2008;16:205–10.

  55. 55.

    Huusko JM, Karjalainen MK, Mahlman M, Haataja R, Kari MA, Andersson S, et al. A study of genes encoding cytokines (IL6, IL10, TNF), cytokine receptors (IL6R, IL6ST), and glucocorticoid receptor (NR3C1) and susceptibility to bronchopulmonary dysplasia. BMC Med Genet. 2014;15:120.

  56. 56.

    Burgess S, Timpson NJ, Ebrahim S, Davey SG. Mendelian randomization: where are we now and where are we going? Int J Epidemiol. 2015;44:379–88.

  57. 57.

    VanderWeele TJ, Tchetgen TE, Cornelis M, Kraft P. Methodological challenges in mendelian randomization. Epidemiology. 2014;25:427–35.

  58. 58.

    Hunter GR, Snyder SW, Kekesszabo T, Nicholson C, Berland L. Intra-abdominal adipose tissue values associated with risk of possessing elevated blood lipids and blood pressure. Obes Res. 2012;2:563–8.

Download references

Authors’ contributions

WYX, LG and WW contributed to the study concept and design. LD, CX, WH and DJ contributed to the acquisition of subjects and/or data. LD, GSQ, ZZY and PHL contributed to the analysis and interpretation of data. LD, SM, WLJ, SMS, GXH, MQ and WYX contributed to the preparation of the manuscript. All authors read and approved the final manuscript.

Acknowledgements

None.

Competing interests

The authors declared that they have no competing interests.

Availability of data and materials

The data are available from the corresponding author.

Consent for publication

Not applicable.

Ethics approval and consent to participate

This study was approved by the ethics committees of Capital Medical University, Beijing, China.

Funding

This work was supported by Grants from National Natural Science Foundation of China (81673247, 81370083, 81773527 and 81530087) and Australian-China Collaborative Grant (NH&MRC-APP1112767–NSFC 81561128020).

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Correspondence to You-Xin Wang.

Additional files

12967_2018_1616_MOESM1_ESM.docx

Additional file 1: Table S1. Structures of the initial IgG glycome.

12967_2018_1616_MOESM2_ESM.docx

Additional file 2: Table S2. The calculation formula of derived glycans.

12967_2018_1616_MOESM3_ESM.docx

Additional file 3: Figure S1. The correlation coefficients in independent variables Statistically significant associations between two variables are shown, while the insignificant correlation coefficients are blank in the boxes. The positive correlations are represented by blue color, while negative correlations are represented by red color. BMI: body mass index; WHR: waist–hip ratio; FBG: fasting blood triglycerides; SBP: systolic blood pressure; DBP: diastolic blood pressure; TC: total cholesterol; TG: triglyceride; HDL: high-density lipoprotein; LDL: low-density lipoprotein; RHR: resting heart rate.

12967_2018_1616_MOESM4_ESM.docx

Additional file 4: Table S3. Description of the IgG glycome.

12967_2018_1616_MOESM5_ESM.docx

Additional file 5: Table S4. Associations between IgG glycan and blood lipid.

12967_2018_1616_MOESM6_ESM.docx

Additional file 6: Table S5. The associations of the normalized glycan variables in dyslipidaemia vs controls.

12967_2018_1616_MOESM7_ESM.docx

Additional file 7: Figure S2. The correlation coefficients in glycans Statistically significant associations between two glycans are shown, while the insignificant correlation coefficients are blank in the boxes. The positive correlations are represented by blue color, while negative correlations are represented by red color.

12967_2018_1616_MOESM8_ESM.docx

Additional file 8: Table S6. The dimension reduction of significant glycans by LASSO method.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Keywords

  • Immunoglobulin G
  • N-Glycosylation
  • Blood lipids
  • Dyslipidaemia