Abstract
Free full text
Immunoglobulin G glycome composition in transition from premenopause to postmenopause
Summary
Gonadal hormones affect immunoglobulin G (IgG) glycosylation, and the more proinflammatory IgG glycome composition might be one of the molecular mechanisms behind the increased proinflammatory phenotype in perimenopause. Using ultra-high-performance liquid chromatography, we analyzed IgG glycome composition in 5,080 samples from 1940 pre-, peri-, and postmenopausal women. Statistically significant decrease in galactosylation and sialylation was observed in postmenopausal women. Furthermore, during the transition from pre- to postmenopausal period, the rate of increase in agalactosylated structures (0.051/yr; 95%CI = 0.043–0.059, p < 0.001) and decrease in digalactosylated (−0.043/yr; 95%CI = −0.050 to −0.037, p < 0.001) and monosialylated glycans (−0.029/yr; 95%CI = −0.034 to −0.024, p < 0.001) were significantly higher than in either pre- or postmenopausal periods. The conversion to the more proinflammatory IgG glycome and the resulting decrease in the ability of IgG to suppress low-grade chronic inflammation may be an important molecular mechanism mediating the increased health risk in perimenopause and postmenopause.
Introduction
Glycosylation of immunoglobulin G (IgG) is an important regulator of the immune system (Martinić Kavur et al., 2021). In addition to directly affecting the effector functions of IgG by promoting its binding to different Fc receptors (Nimmerjahn and Ravetch, 2008a), glycans attached to IgG have numerous other roles in the regulation of the immune system (Ipsen-Escobedo and Nimmerjahn, 2018; Seeling et al., 2017). Glycans are also important mediators of the anti-inflammatory activity of IgG (Karsten et al., 2012; Schwab and Nimmerjahn, 2013) as well as the driving force for one of the molecular mechanisms underlying immunosuppressive activity of intravenous IgG (IVIG) therapy used to treat a variety of immunological disorders (Nimmerjahn and Ravetch, 2008b). Changes in IgG glycosylation are also believed to be an important contributor to aging at the molecular level through a process called inflammaging (Franceschi et al., 2018). Furthermore, hyposialylated IgG, which associates with ageing, was shown to contribute to the development of hypertension in animal models through FcγRII on endothelial cells (Peng et al., 2019; Sundgren et al., 2015).
IgG glycome composition is tightly regulated and, despite the absence of a direct genetic template, is heritable to a significant extent (Menni et al., 2013; Pučić et al., 2011). Alternative glycosylation (the attachment of different glycan structures to the same glycosylation site) is structurally and functionally analogous to coding mutations, but instead of being defined by sequence variation in individual genes, glycome composition is inherited as a complex trait encoded by multiple genes (Krištić et al., 2018). Genome-wide association studies identified a network of over 30 genes that associate with the IgG glycome composition (Klarić et al., 2020; Lauc et al., 2013; Shen et al., 2017). One of the regulators of IgG glycosylation are sex hormones, and the association between estradiol and the IgG glycome composition was confirmed to be causal (Ercan et al., 2017). A recent randomized placebo-controlled clinical study demonstrated that the deprivation of gonadal hormones resulted in an increase of biological age measured by IgG glycans, which was completely prevented by estradiol supplementation (Juric et al., 2020).
Large population studies revealed that IgG glycome composition changes with age and that in females this change is particularly pronounced in the time preceding the average age of menopause (Krištić et al., 2014; Stambuk et al., 2020). However, the association between changes in the IgG glycome composition and menopause has not yet been validated. Perimenopause is accompanied by different symptoms and is very often misdiagnosed due to changing hormone levels and symptom overlap with other conditions. Menopause associates with an increased risk of developing a range of diseases affecting metabolism, bones, cardiovascular system, brain, and even cancer (Lobo et al., 2014). Early perimenopause diagnosis would enable timely interventions to ease the symptoms, decrease risks of developing accompanying diseases, and timely monitoring of any changes in women's health status due to hormonal changes. To fill this knowledge gap and explore the association of IgG glycome changes with perimenopause, we analyzed IgG N-glycome composition in 5,080 samples from 1,940 females in multiple time points during their transition from premenopause to postmenopause, along with 274 samples from 113 males. As an additional control, samples from 301 females in an independent validation cohort were analyzed.
Results
IgG glycome composition in premenopause and menopause
IgG glycome composition was analyzed in multiple samples from 1,940 females and 113 males from the UK registry of adult twins. The 1,087 samples from 732 individual twins (500 families) were collected prior to the onset of menopause, whereas 3,993 samples from 1,678 individual twins (969 families) were collected after the onset of menopause (Table S1). In addition, IgG glycome was analyzed in 113 men (274 samples) as an additional control. A representative chromatogram with detailed structures of individual IgG glycans is shown in Figure S1. Mixed modelling was used to determine average levels of individual IgG glycans in premenopausal women, menopausal women, and men. Large and statistically significant differences were observed in multiple glycans, mainly reflecting decreased galactosylation and sialylation of glycans in menopausal women (Figure 1). Levels of individual IgG glycans are available in Figure S2.
Average yearly change of IgG glycans in perimenopause
It is known that IgG galactosylation and sialylation decrease with age (Krištić et al., 2014), thus we attempted to evaluate whether the changes in glycans associated with the transition to menopause were more extensive than the age-related changes. Because we had multiple samples from the same individual, we were able to calculate the rate of change of individual IgG glycans over time. A subset of women in our cohort (n = 379) entered menopause between two sampling time points, thus for them, we were also able to calculate the rate of changes in individual IgG glycans during the perimenopausal period. The comparison of the rate of age-related IgG glycome changes (when there was no change of the menopause status) with the rate of IgG glycome changes during perimenopause (i.e., between time points when the transition to menopause occurred) revealed statistically significant differences for a number of glycans (Figure S3). The most prominent differences in the perimenopause period were the significantly higher rates of increase in the agalactosylated structures (G0) and the decrease of digalactosylated (G2) and monosialylated (S1) glycans (Figure 2). During perimenopause women had a significant higher rate of increase in agalactosylated structures (0.051/yr; 95%CI = 0.043–0.059, p < 0.001) and decrease in digalactosylated (−0.043/yr; 95%CI = −0.050 to −0.037, p < 0.001) and monosialylated glycans (−0.029/yr; 95%CI = −0.034 to −0.024, p < 0.001), compared with premenopausal women.
Prediction of perimenopause using IgG glycans
Using single-point measurements of IgG N-glycans we attempted to predict menopause status. The lasso model was fitted with menopausal status (TRUE for those who entered menopause and FALSE for all others) as dependent variable and age and all directly measured IgG glycans as predictors. The model was fitted on a randomly chosen one twin per family (training set) in the subset of women with a two time point median age between 45 and 55 years. The remaining twins were used for testing the fitted model. The same train and test subset model with only mean age as a predictor were used to fit and test the model for comparison (Figure 3A). IgG N-glycome in combination with age measured in a single time point has shown a good classification performance, with the AUC of 0.853 (95% CI 0.823–0.880) outperforming only age as a predictor of perimenopause (AUC of 0.818, 95% CI 0.783–0.849). Because changes in the IgG glycome composition were the most prominent in the period of perimenopause, we used the rate of changes in glycans to predict whether given women already entered the period of perimenopause or not (Figure 3B). The addition of the information about changes in the IgG glycan levels to a prediction model notably increases AUC value compared with using only age (0.694, 95% CI 0.637–0.751 compared with 0.507, 95% CI 0.443–0.570, respectively).
IgG glycome composition in replication cohorts
The sampling in replication cohorts was cross-sectional, and there was only one time point for each sample, which did not allow to infer any causal relationships such as monitoring changes in IgG N-glycome over time. Because of that limitation, this was not a perfect replication cohort for our model validation; therefore, it was not possible to assess whether someone is perimenopausal or not. What we could replicate are the consistent differences between premenopausal and menopausal females of any age between 45 and 55 years (Figure 4).
Discussion
By analyzing IgG glycome in multiple samples from the same individuals, we have shown the association between the period of perimenopause and the extensive changes in the IgG glycome composition. Aging associates with the transition of the IgG glycome from an inflammation-suppressive to a proinflammatory composition, but this change is much faster in the period when women are progressing from a regular cycle to menopause, as was shown in several population studies (Krištić et al., 2014; Stambuk et al., 2020). The most extensive change in the IgG glycome during perimenopause was a decrease in galactosylation and sialylation, a change that is generally associated with aging (Martinić Kavur et al., 2021). The depletion of estrogen during perimenopause and menopause is an important component of aging, as low hormone levels increase vulnerability to diseases in hormone-responsive tissues, including the brain, bone, and the cardiovascular system (Lobo et al., 2014). Menopause is, therefore, more complex than just being a natural condition resulting in symptoms; it needs to be considered as a long-term hormone deficiency potentially leading to higher health risks if not appropriately managed with lifestyle adjustments or hormone replacement therapy (HRT) in consonance with a gynecologist and/or general practitioner (Gambacciani et al., 2019).
The age-related gradual decrease in the level of galactosylated IgG and the increase in agalactosylated IgG is known to exacerbate low-grade chronic inflammation (De Martinis et al., 2005; Franceschi et al., 2000; Monti et al., 2017) and acts as an effector of proinflammatory pathological changes (Dall’Olio et al., 2013; Franceschi et al., 2007). Decreased levels of IgG galactosylation and/or sialylation have been associated with a wide range of autoimmune conditions, such as systemic lupus erythematosus, inflammatory bowel disease, and many others. Changes in these IgG glycosylation traits are also associated with disease progression, disease activity, symptom severity, and response to treatment. In many cases, these changes occur before the onset of symptoms (Gudelj et al., 2018a).
The rate of change from “young” to “old” IgG glycome increased significantly in the perimenopausal period, suggesting that the loss of immunosuppressive functions of IgG may be an underlying molecular mediator of at least some health risks in this period. The results of this study also indicate the potential of the IgG glycome composition as a biomarker for perimenopause. Perimenopause can last up to 15 years and is difficult to diagnose due to highly irregular hormonal cycles (National Institute for Health and Care Excellence, 2015). As a result of poor awareness and inappropriate use of hormonal tests, women are often misdiagnosed with conditions such as fibromyalgia, migraines, depression, or chronic fatigue syndrome and are frequently prescribed antidepressants despite there being no evidence to support their use to improve the low mood associated with perimenopause or menopause (Leonhardt, 2019). The importance of early diagnosis of perimenopause is supported by the research that has shown that early use of HRT leads to a greater reduction of future disease incidence as well as minimizing the duration of symptoms (Ward and Deneris, 2018).
Recent placebo-controlled randomized trials demonstrated the importance of estrogens in the regulation of IgG glycosylation (Ercan et al., 2017; Juric et al., 2020), thus it is reasonable to assume that the decrease in estrogen levels is an important driver of changes in IgG glycosylation in perimenopause. Contrary to estrogens that have a very short half-life in circulation, IgG has a half-life of 3 weeks (Vidarsson et al., 2014), which implies that the current composition of the IgG glycome is affected by average concentrations of estrogen in the past 3–4 weeks. Because changes in the IgG glycome during perimenopause are so extensive, they could potentially be used as a diagnostic tool. In fact, IgG glycome composition could potentially be used to determine long-term average concentrations of sex hormones in the same way glycated hemoglobin (HbA1c) is used to determine long-term concentrations of blood glucose.
Conclusions
Changes in IgG glycosylation from “young” to “old” glycome associate with many health risks that accompany menopause (Martinić Kavur et al., 2021). In some diseases such as rheumatoid arthritis (Ercan et al., 2010; Gudelj et al., 2018b) and cardiovascular diseases (Kifer et al., 2021; Menni et al., 2018; Wittenbecher et al., 2020), this change was shown to occur years before disease onset, suggesting the causal role of IgG glycans in disease development. Proinflammatory IgG glycome and the resulting decrease in the ability of IgG to suppress low-grade chronic inflammation may be an important molecular mechanism that mediates the increased health risk in perimenopause, thus this topic should be studied in more detail.
Limitations of the study
The main limitation of this study is the relatively long time between individual sampling time points (6.9 ± 2.6 years) that was not tailored for the occurrence of menopause. Therefore, in some cases, perimenopause was just a short part of the period in which the change in the IgG N-glycome occurred, whereas in others sampling might have occurred within the perimenopause period, with part of the perimenopause being before and another part after a specific sampling. Hopefully, a future study that would have more frequent sampling during the perimenopause period would enable an even more accurate assessment of IgG glycome association with perimenopause.
STARMethods
Key resources table
REAGENT or RESOURCE | SOURCE | IDENTIFIER |
---|---|---|
Biological samples | ||
Human serum samples | TwinsUK Registry | http://www.twinsuk.ac.uk |
Human plasma samples, CROATIA-Korčula | “10001 Dalmatians” – Croatian National Biobank (Krištić et al., 2014) | http://www.mefst.unist.hr/znanost/istrazivacke-skupine-i-laboratoriji/10-001-dalmatinac-hrvatska-biobanka/5035 |
Human plasma samples, CROATIA-Vis | “10001 Dalmatians” – Croatian National Biobank (Krištić et al., 2014) | http://www.mefst.unist.hr/znanost/istrazivacke-skupine-i-laboratoriji/10-001-dalmatinac-hrvatska-biobanka/5035 |
Chemicals, peptides, and recombinant proteins | ||
Protein G monolithic 96-well plate | BIA Separations | Pucic et al. 2011 |
Acetonitrile, LC/MS grade | Honeywell | CAS#75-05-8, EC#200-835-2 |
Formic acid | Merck | Cat#100264, CAS#64-18-6, EC#200-579-1 |
Ammonium bicarbonate | Sigma Aldrich | Cat#09830, CAS#1066-33-7, EC#213-911-5 |
2-aminobenzamide (2-AB) | Sigma Aldrich | CAS# 88-68-6, EC# 201-851-2 |
Dimethyl sulfoxide (DMSO) | Sigma Aldrich | CAS# 67-68-5, EC# 200-664-3 |
Glacial acetic acid | Merck | CAS# 64-19-7, EC#200-580-7 |
1xPBS (phosphate buffer saline); prepared from: | prepared in-house | |
NaCI (Sodium chloride) | Alkaloid Skopje | CAS#7647-14-5, EC#231-598-3 |
Na2HPO4 (Disodium phosphate) | BIOCHEM Chemopharma | CAS#7558-79-4, EC#231-448-7 |
KH₂PO₄ (Monopotassium phosphate) | BIOCHEM Chemopharma | CAS#7778-77-0, EC#231-913-4 |
KCl (Potassium chloride) | Sigma Aldrich | CAS#7447-40-7, EC#231-211-8 |
Critical commercial assays | ||
GlycoWorks RapiFluor-MS N-Glycan Kit | Waters Corporation | Cat#176003910 |
Software and algorithms | ||
Empower 3 software, build 3471 | Waters Corporation | https://www.waters.com/waters/promotionDetail.htm?id=10184996&locale=en_US |
R programming language | R Core Team | https://www.r-project.org/ |
R package ‘sva’ | (Leek et al., 2012) | https://bioconductor.org/packages/release/bioc/html/sva.html |
R package ‘lme4’ | (Bates et al., 2015) | https://cran.r-project.org/web/packages/lme4/index.html |
R package ‘emmeans’ | (Lenth, 2020) | https://cran.r-project.org/web/packages/emmeans/index.html |
R package ‘caret’ | (Kuhn, 2020) | https://cran.r-project.org/web/packages/caret/vignettes/caret.html |
R package ‘glmnet’ | (Friedman et al., 2010) | https://cran.r-project.org/web/packages/glmnet/index.html |
R package ‘pROC’ | (Robin et al., 2011) | https://cran.r-project.org/web/packages/pROC/index.html |
Resource availability
Lead contact
The lead contact for further information and requests is Gordan Lauc ([email protected]).
Materials availability
The study did not generate new unique reagents.
Experimental model and subject details
Subject details
Study subjects were individuals enrolled in the TwinsUK registry, a national register of adult twins recruited as volunteers without selecting for any particular disease or traits (Verdi et al., 2019). Data on menopause was gathered through the multiple-choice questionnaire asking women at each visit: “What is your menopausal status?” with possible answers “Premenopausal”, “Going through the menopause”, “Postmenopausal” and “Don’t know” and “How old were you when you became postmenopausal (when you stopped having periods for one year or more)?”. Consequently, a sample was defined as menopausal if the questionnaire resulted in a “Postmenopausal” answer or age greater than the reported age of becoming postmenopausal. Similarly, a sample was defined as premenopausal if the questionnaire resulted with a “Premenopausal” answer or age at sampling was less than the reported age of becoming postmenopausal reduced by 1 year. Any inconsistent (i.e. “Premenopausal” answer at age after becoming menopausal) or unreliable data (i.e. answered as “Going through the menopause”) were dropped from the study (Table S1). The TwinsUK study was approved by NRES Committee London–Westminster, and all twins provided informed written consent.
Validation cohort
The validation cohort consisted of two different population studies, CROATIA-Vis and CROATIA-Korcula, described elsewhere in detail (Polašek et al., 2009; Rudan et al., 2009; Vitart et al., 2008). These studies were part of a larger “10001 Dalmatians” study of Croatian island isolates, in this case, islands of Vis and Korcula, within the scope of the Croatian National Biobank. All human participants included in this study signed informed consent and the study was approved by the appropriate Ethics Committee of the University of Split Medical School.
Method details
Isolation of IgG
IgG isolation from serum samples was performed on protein G monolithic 96-well plates (Trbojević-Akmačić et al., 2017), 100 μL of serum was diluted 1:7 (v/v) with 1xPBS (phosphate buffer saline) followed by filtration on a 0.45 μm GHP filter plate (Pall Corporation, USA). Diluted and filtered serum samples were then transferred to the protein G monolithic plate (BIA Separations, Slovenia), which was then washed three times with 2 ml 1xPBS. Elution of IgG was performed with 1 mL of 0.1 mol/L formic acid (Merck, Germany) followed by immediate neutralization of the mixture with ammonium bicarbonate (Acros Organics, USA) to pH 7.0. IgG-containing eluate was aliquoted into a PCR plate (average mass of 15 μg of IgG per sample) and dried in a vacuum centrifuge.
Deglycosylation, labelling and purification of IgG N-glycans
Deglycosylation, RapiFluor-MS labelling and purification of IgG N-glycans were performed using the GlycoWorks RapiFluor-MS N-Glycan Kit obtained from Waters Corporation (USA). The entire procedure was done following Waters’ protocol (Waters Corporation, 2017, https://www.waters.com/webassets/cms/support/docs/715004793en.pdf), which was further adapted to suit the high-throughput analysis in the 96-well PCR plate format (Deriš et al., 2021). Isolated and dried IgG was first resuspended in 10.8 μL ultrapure water, swiftly denatured with 3 μL 5 % RapiGest SF solution in a 3 min reaction at 99°C and enzymatically deglycosylated with 1.2 μL GlycoWorks Rapid PNGase F in a 5 min reaction at 50°C. Released N-glycans were then labeled with the RapiFluor-MS label. RapiFluor-MS dye contains a rapid tagging function group, N-hydroxysuccinimide carbamate, which rapidly modifies glycosylamine-bearing N-glycans after their enzymatic release, yielding a stable urea linkage (Lauber et al., 2015). Labeled N-glycans were purified by hydrophilic interaction chromatography-solid phase extraction (HILIC-SPE) clean-up. In the end, the samples of released and labelled IgG glycans were stored at -20 °C until further use.
HILIC-UHPLC-FLR analysis of RapiFlour-MS labelled IgG N-glycans
RapiFluor-MS labeled IgG N-glycans were analysed using ultra-high-performance liquid chromatography based on hydrophilic interactions with fluorescence detection (HILIC-UHPLC-FLR) on Waters Acquity UPLC H-class instruments. The instruments were controlled and monitored with the Empower 3 software, build 3471 (Waters, USA). Chromatographic separation of glycan structures was performed on Waters UPLC Glycan bridged ethylene hybrid (BEH) Amide chromatographic columns (130 Å, 1.7 μm BEH particles, 2.1 × 100 mm). Solvent A was 50 mmol/L ammonium formate, pH 4.4, while solvent B was 100% LC-MS grade acetonitrile. A linear gradient of 75–61.5 % acetonitrile (v/v) was used at a flow rate of 0.4 mL/min over 30 min in a 42-min analytical run. Each chromatogram was separated into 22 glycan peaks containing IgG glycan structures (Figure S1.), by automated integration (Agakova et al., 2017) resulting in relative quantification of IgG N-glycans (total area normalization). Specific N-glycan structures found in each glycan peak according to Keser et al. (2018) are presented in Table S2.
Validation cohort
IgG was immobilized in a block of sodium dodecyl sulfate–polyacrylamide gel and N-glycans were released by digestion with PNGase F (ProZyme, Hayward, CA), in a 96-well microtiter plate to achieve the best throughput. Released N-glycans were then labeled with 25 μL fluorescent 2-aminobenzamide (2-AB) dye dissolved in dimethyl sulfoxide (DMSO) and glacial acetic acid mixture (85:15, v/v). The excess label was removed by solid-phase extraction using Whatman 3MM chromatography paper (Royle et al., 2008). Labeled glycans were eluted with water and stored at −20°C until further use. UHPLC analysis was performed as described previously (Krištić et al., 2014); 2-AB labeled IgG N-glycans were separated on a Waters UPLC Glycan BEH Amide chromatographic columns with 100 mmol/L ammonium formate, pH 4.4, as solvent A and 100% LC-MS grade acetonitrile as solvent B. Separation method used a linear gradient of 75–62% acetonitrile (v/v) at a flow rate of 0.4mL/min in a 25-min analytical run. Chromatograms were separated into 24 glycan peaks.
Quantification and statistical analysis
Statistical analysis
After UHPLC analysis, each measured glycan peak was normalized by total chromatogram area. Observed relative glycan abundances were log transformed and batch corrected using ComBat method (R package ‘sva’ (Leek et al., 2012)). After batch correction, values were back-transformed to relative abundance and derived traits were calculated according to formulas presented in Table S3.
Overlapping measured N-glycomes and information available on menopause found in TwinsUK register resulted in 5354 samples from 1940 females and 113 males.
For each glycan trait, the effect of menopause on the abundance was estimated by mixed modelling (R package ‘lme4’ (Bates et al., 2015)). Logit transformed relative abundance was used as dependent variable and sex, menopausal status (nested within sex), age and age-menopausal status interaction as fixed factors, as well as family ID and individual ID (nested within family ID) as random intercepts and age as random slope. Using the model, mean values were estimated for each group and pairwise compared using t-test (post-hoc) (R package ‘emmeans’ (Lenth, 2020)). Given multiple testing for a large number of glycans, the control of false discovery rate (FDR) was done by adjusting p values according to Benjamini-Hochberg’s method.
For the analysis of the rate of change of glycan abundances, average yearly changes were calculated from two time points according to the formula: rate = (logit(glycan2) – logit(glycan1))/(age2 – age1), where glycan is relative abundance resulted from normalized chromatogram, and age is represented in years, assuming age2 > age1. Depending on menopausal status in the first and second time point measurements were classified into three categories: premenopausal (both timepoints before menopause), perimenopausal (first timepoint before, and second in the menopause), and menopausal (both timepoints in menopause). The diagram clarifying which time points were used to calculate the rate of change in IgG-N-glycome composition is shown in Figure S4. Mean rates were estimated by the mixed model, in which rate was dependent variable, and sex, menopausal group (nested within menopausal status) and mean age (calculated as (age2+age1)/2) were set as fixed factors, as well as family ID and individual ID (nested within family ID) as random intercepts. Estimated mean rates for each menopausal group and males were pairwise compared using a t-test (post-hoc). FDR control was done by adjusting p values according to Benjamini-Hochberg’s method.
N-glycome (observed in one timepoint) potential for diagnosis of perimenopause was analysed in a subset of women aged between 45 and 55 years by logistic regression with the L1 regularization technique (R package ‘caret’ and ‘glmnet’ (Friedman et al., 2010; Kuhn, 2020)). The subset was split into training (randomly chosen one measurement per family) and test (everything else) subset. The logistic model was defined with menopausal status (TRUE/FALSE) as the dependent variable, and age and all directly measured glycans as independent variables. To avoid overfitting, coefficients of the logistic model were estimated by 10-fold cross validation on a training subset data, with hyperparameter lambda set to lambdamin. The model was tested on a test subset; results are presented by receiver operating characteristic (ROC) curve. In parallel, the null model was fitted in the same way but using only age as a predictor. Observed ROC curves of the two models were compared by a bootstrap test (B = 2000) (R package ‘pROC’ (Robin et al., 2011)).
The rate of N-glycome change (based on two time points) potential for diagnosis of perimenopause was analysed in a subset of women with median two-timepoint age between 45 and 55 years in the same way as diagnosis of menopause with IgG glycome analysed in one timepoint. The perimenopausal category was defined as TRUE, while others (postmenopausal and premenopausal) were defined as FALSE. All statistical analyses were done using R (R Core Team, 2020).
Validation cohort
For the following statistical analysis, only female samples of the matching ages, 45 to 55 years of age, were selected, which included 62 premenopausal females in the Vis study and 86 in the Korcula study and 70 and 83 menopausal females in the Vis and Korcula study respectively, as presented in the Table S4. Menopausal status was projected from a survey in which women indicated the age of becoming menopausal or answered the question (YES / NO) whether they had a regular menstrual cycle. Also, all women who use hormone replacement therapy (HRT) were excluded from the analysis. For each glycan trait, the effect of menopause on the abundance was estimated by mixed modelling (R package ‘lme4’ (Bates et al., 2015)). Logit transformed relative abundance was used as the dependent variable and menopausal status and age as fixed factors, as well as cohort ID as a random intercept. Using the model, mean values were estimated for each menopausal status group and compared using t-test (post-hoc) (R package ‘emmeans’ (Lenth, 2020)). Given multiple testing for the number of glycans, the control of false discovery rate (FDR) was done by adjusting p values according to Benjamini-Hochberg’s method.
Acknowledgements
This study has been supported by the European Structural and Investment Funds IRI “CardioMetabolic” grant (#KK.01.2.1.02.0321), Center of Competence in Molecular Diagnostics grant (#KK.01.2.2.03.0006), and Croatian National Center of Research Excellence in Personalized Healthcare grant (#KK.01.1.1.01.0010) and the Human Glycome Project. Equipment and products from Waters and New England Biolabs were used for this research. The Department of Twin Research receives support from grants from the Wellcome Trust (212904/Z/18/Z) and the Medical Research Council (MRC)/British Heart Foundation (BHF) Ancestry and Biological Informative Markers for Stratification of Hypertension (AIM-HY; MR/M016560/1), European Union, Chronic Disease Research Foundation (CDRF), Zoe Global Ltd., the NIHR Clinical Research Facility and Biomedical Research Centre (based at Guy’|'s and St Thomas' NHS Foundation Trust in partnership with King’|'s College London). C.M. is funded by the Chronic Disease Research Foundation and by the MRC AIM-HY project grant.
Author contributions
D.K.: Data analysis and interpretation, Visualization, Writing—Original draft preparation, Writing—Review and editing; H.D.: Quantification of IgG N-glycans, Data interpretation, Writing—Review and editing; A.Ci.: Quantification of IgG N-glycans, Data interpretation, Writing—Review and editing; T.P.: Quantification of IgG N-glycans, Data interpretation, Writing—Review and editing; A.Cv.: Quantification of IgG N-glycans, Data interpretation, Writing—Original draft preparation, Writing—Review and editing; I.T.-A.: Quantification of IgG N-glycans, Data interpretation, Writing—Review and editing; I.K.: Data interpretation, Writing—Review and editing; O.P.: Data interpretation, Writing—Review and editing; L.N.: Data interpretation, Writing—Review and editing; T.S.; Conceptualization, Data interpretation, Writing—Review and editing; C.M.: Conceptualization, Data interpretation, Writing—Review and editing; G.L.: Conceptualization, Supervision, Data interpretation, Writing—Original draft preparation, Writing—Review and editing.
Declaration of interests
A patent covering all the main aspects of the use of IgG glycome as a predictor of menopause and perimenopause has been filed by Genos Ltd (application number: P20210509A). The application is currently pending. D.K., C.M., and G.L. are named as co-inventors on the patent application. G.L. is the founder and owner of Genos Ltd, a biotech company specializing in high-throughput glycomics that also has several patents in the field. H.D., A.Ci. T.P. and I.T.-A. are employees of Genos Ltd. The remaining authors declare no competing interests.
Footnotes
Supplemental information can be found online at https://doi.org/10.1016/j.isci.2022.103897.
Data and code availability
- • The raw data reported in this study cannot be deposited in a public repository because the data are confidential records. To request access to data from the Twins UK Registry contact data access manager Victoria Vazquez (https://twinsuk.ac.uk/resources-for-researchers/access-our-data/) and to access the data from Vis and Korčula studies, contact Ozren Polašek ([email protected]).
- • This paper does not report original code.
- • Any additional information required to reanalyze the data reported in this paper may be available from the lead contact upon request.
References
- Agakova A., Vučković F., Klarić L., Lauc G., Agakov F. Automated integration of a UPLC glycomic profile. Methods Mol. Biol. 2017;1503:217–233. [Abstract] [Google Scholar]
- Bates D., Mächler M., Bolker B., Walker S. Fitting linear mixed-effects models using lme4. J. Stat. Softw. 2015;67:1–48. [Google Scholar]
- Dall’Olio F., Vanhooren V., Chen C.C., Slagboom P.E., Wuhrer M., Franceschi C. N-glycomic biomarkers of biological aging and longevity: a link with inflammaging. Ageing Res. Rev. 2013;12:685–698. [Abstract] [Google Scholar]
- De Martinis M., Franceschi C., Monti D., Ginaldi L. Inflamm-ageing and lifelong antigenic load as major determinants of ageing rate and longevity. FEBS Lett. 2005;579:2035–2039. [Abstract] [Google Scholar]
- Deriš H., Cindrić A., Lauber M., Petrović T., Bielik A., Taron C.H., Van Wingerden M., Lauc G., Trbojević-Akmačić I. Robustness and repeatability of GlycoWorks RapiFluor-MS IgG N-glycan profiling in a long-term high-throughput glycomic study. Glycobiology. 2021;31:1062–1067. [Abstract] [Google Scholar]
- Ercan A., Cui J., Chatterton D.E., Deane K.D., Hazen M.M., Brintnell W., O’Donnell C.I., Derber L.A., Weinblatt M.E., Shadick N.A., et al. Aberrant IgG galactosylation precedes disease onset, correlates with disease activity, and is prevalent in autoantibodies in rheumatoid arthritis. Arthritis Rheum. 2010;62:2239–2248. [Europe PMC free article] [Abstract] [Google Scholar]
- Ercan A., Kohrt W.M., Cui J., Deane K.D., Pezer M., Yu E.W., Hausmann J.S., Campbell H., Kaiser U.B., Rudd P.M., et al. Estrogens regulate glycosylation of IgG in women and men. JCI Insight. 2017;2:e89703. [Europe PMC free article] [Abstract] [Google Scholar]
- Franceschi C., Bonafe M., Valensin S., Olivieri F., De Luca M., Ottaviani E., De Benedictis G. Inflamm-aging. An evolutionary perspective on immunosenescence. Ann. N. Y. Acad. Sci. 2000;908:244–254. [Abstract] [Google Scholar]
- Franceschi C., Capri M., Monti D., Giunta S., Olivieri F., Sevini F., Panourgia M.P., Invidia L., Celani L., Scurti M., et al. Inflammaging and anti-inflammaging: a systemic perspective on aging and longevity emerged from studies in humans. Mech. Ageing Dev. 2007;128:92–105. [Abstract] [Google Scholar]
- Franceschi C., Garagnani P., Parini P., Giuliani C., Santoro A. Inflammaging: a new immune–metabolic viewpoint for age-related diseases. Nat. Rev. Endocrinol. 2018;14:576–590. [Abstract] [Google Scholar]
- Friedman J., Hastie T., Tibshirani R. Regularization paths for generalized linear models via coordinate descent. J. Stat. Softw. 2010;33:1–22. [Europe PMC free article] [Abstract] [Google Scholar]
- Gambacciani M., Cagnacci A., Lello S. Hormone replacement therapy and prevention of chronic conditions. Climacteric. 2019;22:303–306. [Abstract] [Google Scholar]
- Gudelj I., Lauc G., Pezer M. Immunoglobulin G glycosylation in aging and diseases. Cell. Immunol. 2018;333:65–79. 10.1016/j.cellimm.2018.07.009. [Abstract] [CrossRef] [Google Scholar]
- Gudelj I., Salo P.P., Trbojević-Akmačić I., Albers M., Primorac D., Perola M., Lauc G. Low galactosylation of IgG associates with higher risk for future diagnosis of rheumatoid arthritis during 10 years of follow-up. Biochim. Biophys. Acta. 2018;1864:2034–2039. [Abstract] [Google Scholar]
- Ipsen-Escobedo A., Nimmerjahn F. More rules, still exceptions: understanding immunomodulatory antibody activity in vivo. Cancer Cell. 2018;33:545–546. [Abstract] [Google Scholar]
- Juric J., Kohrt W.M., Kifer D., Gavin K.M., Pezer M., Nigrovic P.A., Lauc G. Effects of estradiol on biological age measured using the glycan age index. AGING-US. 2020;12:19756–19765. [Europe PMC free article] [Abstract] [Google Scholar]
- Karsten C.M., Pandey M.K., Figge J., Kilchenstein R., Taylor P.R., Rosas M., McDonald J.U., Orr S.J., Berger M., Petzold D., et al. Anti-inflammatory activity of IgG1 mediated by Fc galactosylation and association of FcgammaRIIB and dectin-1. Nat. Med. 2012;18:1401–1406. [Europe PMC free article] [Abstract] [Google Scholar]
- Keser T., Pavic T., Lauc G., Gornik O. Comparison of 2-aminobenzamide, procainamide and RapiFluor-MS as derivatizing agents for high-throughput HILIC-UPLC-FLR-MS N-glycan analysis. Front. Chem. 2018;6:1–12. [Europe PMC free article] [Abstract] [Google Scholar]
- Kifer D., Louca P., Cvetko A., Deriš H., Cindrić A., Grallert H., Peters A., Polašek O., Gornik O., Mangino M., et al. N-glycosylation of immunoglobulin G predicts incident hypertension. J. Hypertens. 2021;39:2527–2533. [Europe PMC free article] [Abstract] [Google Scholar]
- Klarić L., Tsepilov Y.A., Stanton C.M., Mangino M., Sikka T.T., Esko T., Pakhomov E., Salo P., Deelen J., McGurnaghan S.J., et al. Glycosylation of immunoglobulin G is regulated by a large network of genes pleiotropic with inflammatory diseases. Sci. Adv. 2020;6:eaax0301. [Europe PMC free article] [Abstract] [Google Scholar]
- Krištić J., Vučković F., Menni C., Klarić L., Keser T., Beceheli I., Pučić-Baković M., Novokmet M., Mangino M., Thaqi K., et al. Glycans are a novel biomarker of chronological and biological ages. J. Gerontol. A Biol. Sci. Med. Sci. 2014;69:779–789. [Europe PMC free article] [Abstract] [Google Scholar]
- Krištić J., Zaytseva O.O., Ram R., Nguyen Q., Novokmet M., Vučković F., Vilaj M., Trbojević-Akmačić I., Pezer M., Davern K.M., et al. Profiling and genetic control of the murine immunoglobulin G glycome. Nat. Chem. Biol. 2018;14:516–524. [Abstract] [Google Scholar]
- Kuhn M. CRAN Repository; 2020. [Google Scholar]
- Lauber M.A., Yu Y.Q., Brousmiche D.W., Hua Z., Koza S.M., Magnelli P., Guthrie E., Taron C.H., Fountain K.J. Rapid preparation of released N -glycans for HILIC analysis using a labeling reagent that facilitates sensitive fluorescence and ESI-MS detection. Anal. Chem. 2015;87:5401–5409. [Abstract] [Google Scholar]
- Lauc G., Huffman J.E., Pučić M., Zgaga L., Adamczyk B., Mužinić A., Novokmet M., Polašek O., Gornik O., Krištić J., et al. Loci associated with N-glycosylation of human immunoglobulin G show pleiotropy with autoimmune diseases and haematological cancers. PLoS Genet. 2013;9:e1003225. [Europe PMC free article] [Abstract] [Google Scholar]
- Leek J.T., Johnson W.E., Parker H.S., Jaffe A.E., Storey J.D. The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics. 2012;28:882–883. [Europe PMC free article] [Abstract] [Google Scholar]
- Lenth R. CRAN Repository; 2020. [Google Scholar]
- Leonhardt M. Low mood and depressive symptoms during perimenopause – should General Practitioners prescribe hormone replacement therapy or antidepressants as the first-line treatment? Post Reprod. Heal. 2019;25:124–130. [Abstract] [Google Scholar]
- Lobo R.A., Davis S.R., De Villiers T.J., Gompel A., Henderson V.W., Hodis H.N., Lumsden M.A., Mack W.J., Shapiro S., Baber R.J. Prevention of diseases after menopause. Climacteric. 2014;17:540–556. [Abstract] [Google Scholar]
- Martinić Kavur M., Lauc G., Pezer M. Reference Module in Chemistry, Molecular Sciences and Chemical Engineering. Elsevier; 2021. [Google Scholar]
- Menni C., Gudelj I., MacDonald-Dunlop E., Mangino M., Zierer J., Bešić E., Joshi P.K., Trbojević-Akmačić I., Chowienczyk P.J., Spector T.D., et al. Glycosylation profile of immunoglobulin G is cross-sectionally associated with cardiovascular disease risk score and subclinical atherosclerosis in two independent cohorts. Circ. Res. 2018;122:1555–1564. [Europe PMC free article] [Abstract] [Google Scholar]
- Menni C., Keser T., Mangino M., Bell J.T., Erte I., Akmačić I., Vučković F., Baković M.P., Gornik O., McCarthy M.I., et al. Glycosylation of immunoglobulin G: role of genetic and epigenetic influences. PLoS One. 2013;8:e82558. [Europe PMC free article] [Abstract] [Google Scholar]
- Monti D., Ostan R., Borelli V., Castellani G., Franceschi C. Inflammaging and human longevity in the omics era. Mech. Ageing Dev. 2017;165:129–138. [Abstract] [Google Scholar]
- National Institute for Health and Care Excellence . Nice; 2015. [Abstract] [Google Scholar]
- Nimmerjahn F., Ravetch J.V. Fcgamma receptors as regulators of immune responses. Nat. Rev. Immunol. 2008;8:34–47. [Abstract] [Google Scholar]
- Nimmerjahn F., Ravetch J.V. Anti-inflammatory actions of intravenous immunoglobulin. Annu. Rev. Immunol. 2008;26:513–533. [Abstract] [Google Scholar]
- Peng J., Vongpatanasin W., Sacharidou A., Kifer D., Yuhanna I.S., Banerjee S., Tanigaki K., Polasek O., Chu H., Sundgren N.C., et al. Supplementation with the sialic acid precursor N-Acetyl-D-Mannosamine breaks the link between obesity and hypertension. Circulation. 2019;140:2005–2018. [Europe PMC free article] [Abstract] [Google Scholar]
- Polašek O., Marušić A., Rotim K., Hayward C., Vitart V., Huffman J., Campbell S., Janković S., Boban M., Biloglav Z., et al. Genome-wide association study of anthropometric traits in Korcula Island, Croatia. Croat. Med. J. 2009;50:7–16. [Europe PMC free article] [Abstract] [Google Scholar]
- Pučić M., Knežević A., Vidič J., Adamczyk B., Novokmet M., Polašek O., Gornik O., Šupraha-Goreta S., Wormald M.R., Redžic I., 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.010090. [Europe PMC free article] [Abstract] [Google Scholar]
- R Core Team . R Foundation for Statistical Computing; 2020. [Google Scholar]
- Robin X., Turck N., Hainard A., Tiberti N., Lisacek F., Sanchez J.-C., Müller M. pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinformatics. 2011;12:77. [Europe PMC free article] [Abstract] [Google Scholar]
- Royle L., Campbell M.P., Radcliffe C.M., White D.M., Harvey D.J., Abrahams J.L., Kim Y.G., Henry G.W., Shadick N.A., Weinblatt M.E., et al. HPLC-based analysis of serum N-glycans on a 96-well plate platform with dedicated database software. Anal. Biochem. 2008;376:1–12. [Abstract] [Google Scholar]
- Rudan I., Marušić A., Janković S., Rotim K., Boban M., Lauc G., Grković I., Dogaš Z., Zemunik T., Vatavuk Z., et al. 10 001 Dalmatians:” Croatia launches its national biobank. Croat. Med. J. 2009;50:4. [Europe PMC free article] [Abstract] [Google Scholar]
- Schwab I., Nimmerjahn F. Intravenous immunoglobulin therapy: how does IgG modulate the immune system? Nat. Rev. Immunol. 2013;13:176–189. [Abstract] [Google Scholar]
- Seeling M., Brückner C., Nimmerjahn F. Differential antibody glycosylation in autoimmunity: sweet biomarker or modulator of disease activity? Nat. Rev. Rheumatol. 2017;13:621–630. [Abstract] [Google Scholar]
- Shen X., Klarić L., Sharapov S., Mangino M., Ning Z., Wu D., Trbojević-Akmačić I., Pučić-Baković M., Rudan I., Polašek O., et al. Multivariate discovery and replication of five novel loci associated with Immunoglobulin G N-glycosylation. Nat. Commun. 2017;8:447. [Europe PMC free article] [Abstract] [Google Scholar]
- Stambuk J., Nakic N., Vuckovic F., Pucic-Bakovic M., Razdorov G., Trbojevic-Akmacic I., Novokmet M., Keser T., Vilaj M., Stambuk T., et al. Global variability of the human IgG glycome. AGING-US. 2020;12:15222–15259. [Europe PMC free article] [Abstract] [Google Scholar]
- Sundgren N.C., Vongpatanasin W., Boggan B.M.D., Tanigaki K., Yuhanna I.S., Chambliss K.L., Mineo C., Shaul P.W. IgG receptor FcγRIIB plays a key role in obesity-induced hypertension. Hypertension. 2015;65:456–462. [Europe PMC free article] [Abstract] [Google Scholar]
- Trbojević-Akmačić I., Ugrina I., Lauc G. Comparative analysis and validation of different steps in glycomics studies. Methods Enzymol. 2017;586:37–55. [Abstract] [Google Scholar]
- Verdi S., Abbasian G., Bowyer R.C.E., Lachance G., Yarand D., Christofidou P., Mangino M., Menni C., Bell J.T., Falchi M., et al. TwinsUK: the UK adult twin registry update. Twin Res. Hum. Genet. 2019;22:523–529. [Abstract] [Google Scholar]
- Vidarsson G., Dekkers G., Rispens T. IgG subclasses and allotypes: from structure to effector functions. Front. Immunol. 2014;5:520. [Europe PMC free article] [Abstract] [Google Scholar]
- Vitart V., Rudan I., Hayward C., Gray N.K., Floyd J., Palmer C.N.A., Knott S.A., Kolcic I., Polasek O., Graessler J., et al. SLC2A9 is a newly identified urate transporter influencing serum urate concentration, urate excretion and gout. Nat. Genet. 2008;40:437–442. [Abstract] [Google Scholar]
- Ward K., Deneris A. An update on menopause management. J. Midwifery Women’s Heal. 2018;63:168–177. [Abstract] [Google Scholar]
- Wittenbecher C., Štambuk T., Kuxhaus O., Rudman N., Vučković F., Štambuk J., Schiborn C., Rahelić D., Dietrich S., Gornik O., et al. Plasma N-glycans as emerging biomarkers of cardiometabolic risk: a prospective investigation in the EPIC-Potsdam cohort study. Diabetes Care. 2020;43:661–668. [Abstract] [Google Scholar]
Articles from iScience are provided here courtesy of Elsevier
Full text links
Read article at publisher's site: https://doi.org/10.1016/j.isci.2022.103897
Read article for free, from open access legal sources, via Unpaywall: http://www.cell.com/article/S2589004222001675/pdf
Citations & impact
Impact metrics
Citations of article over time
Alternative metrics
Discover the attention surrounding your research
https://www.altmetric.com/details/122856922
Article citations
Glycan clock of ageing-analytical precision and time-dependent inter- and i-individual variability.
Geroscience, 46(6):5781-5796, 14 Jun 2024
Cited by: 1 article | PMID: 38877341 | PMCID: PMC11494675
Immunoglobulin G N-glycan markers of accelerated biological aging during chronic HIV infection.
Nat Commun, 15(1):3035, 10 Apr 2024
Cited by: 0 articles | PMID: 38600088 | PMCID: PMC11006954
Anecdotal Evidence Elucidates the Aging Process.
Curr Aging Sci, 17(3):175-179, 01 Jan 2024
Cited by: 0 articles | PMID: 37424351
Impact of estrogen on IgG glycosylation and serum protein glycosylation in a murine model of healthy postmenopause.
Front Endocrinol (Lausanne), 14:1243942, 11 Sep 2023
Cited by: 0 articles | PMID: 37766692 | PMCID: PMC10519799
Polyfunctional antibodies: a path towards precision vaccines for vulnerable populations.
Front Immunol, 14:1183727, 27 Jun 2023
Cited by: 4 articles | PMID: 37600816 | PMCID: PMC10433199
Review Free full text in Europe PMC
Go to all (12) article citations
Similar Articles
To arrive at the top five similar articles we use a word-weighted algorithm to compare words from the Title and Abstract of each citation.
IgG N-glycome changes during the course of severe COVID-19: An observational study.
EBioMedicine, 81:104101, 27 Jun 2022
Cited by: 17 articles | PMID: 35773089 | PMCID: PMC9234382
Association of systemic lupus erythematosus with decreased immunosuppressive potential of the IgG glycome.
Arthritis Rheumatol, 67(11):2978-2989, 01 Nov 2015
Cited by: 133 articles | PMID: 26200652 | PMCID: PMC4626261
Inflammatory bowel disease associates with proinflammatory potential of the immunoglobulin G glycome.
Inflamm Bowel Dis, 21(6):1237-1247, 01 Jun 2015
Cited by: 103 articles | PMID: 25895110 | PMCID: PMC4450892
Can we suppress chronic systemic inflammation and postpone age-related diseases by targeting the IgG glycome?
Expert Opin Ther Targets, 28(6):491-499, 02 Nov 2023
Cited by: 0 articles | PMID: 37897176
Review
Funding
Funders who supported this work.
British Heart Foundation
Medical Research Council (1)
MICA: Ancestry and biological Informative Markers for stratification of HYpertension: The AIM HY study
Professor Phil Chowienczyk, King's College London
Grant ID: MR/M016560/1
Wellcome Trust (1)
Grant ID: 212904/Z/18/Z