Comparative analysis of transferrin and IgG N-glycosylation in two human populations.

Human plasma transferrin (Tf) N-glycosylation has been mostly studied as a marker for congenital disorders of glycosylation, alcohol abuse, and hepatocellular carcinoma. However, inter-individual variability of Tf N-glycosylation is not known, mainly due to technical limitations of Tf isolation in large-scale studies. Here, we present a highly specific robust high-throughput approach for Tf purification from human blood plasma and detailed characterization of Tf N-glycosylation on the level of released glycans by ultra-high-performance liquid chromatography based on hydrophilic interactions and fluorescence detection (HILIC-UHPLC-FLD), exoglycosidase sequencing, and matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF-MS). We perform a large-scale comparative study of Tf and immunoglobulin G (IgG) N-glycosylation analysis in two human populations and demonstrate that Tf N-glycosylation is associated with age and sex, along with multiple biochemical and physiological traits. Observed association patterns differ compared to the IgG N-glycome corroborating tissue-specific N-glycosylation and specific N-glycans' role in their distinct physiological functions.

Analysis of 29 Targeted Genes for Non-Obstructive Azoospermia: The Relationship between Genetic Testing and Testicular Histology.

PURPOSE: To analyze the presence of potentially pathogenic variants of 29 candidate genes known to cause spermatogenic failure (SPGF) in patients with non-obstructive azoospermia (NOA) who underwent testicular histology. MATERIALS AND METHODS: Forty-eight patients with unexplained NOA referred to the Department of Transfusion Medicine and Transplantation Biology, University Hospital Center Zagreb, Zagreb, Croatia for testicular biopsy. They were divided into three groups: those who had cryptorchidism (n=9), those with varicocele (n=14), and those with idiopathic NOA (n=25). All included patients underwent blood withdrawal for next-generation sequencing (NGS) analysis and gene sequencing. RESULTS: We found a possible genetic cause in 4 patients with idiopathic NOA (16%) and in 2 with cryptorchidism (22%). No pathogenic or possibly pathogenic mutations were identified in patients with varicocele. Variants of undetermined significance (VUS) were found in 11 patients with idiopathic NOA (44%), 3 with cryptorchidism (33%), and 8 patients with varicocele (57%). VUSs of the USP9Y gene were the most frequently as they were found in 14 out of 48 patients (29%). In particular, the VUS USP9Y c.7434+14del was found in 11 patients. They showed varied histological pictures, including Sertoli cell-only syndrome, mixed atrophy, and hypospermatogenesis, regardless of cryptorchidism or varicocele. No direct correlation was found between the gene mutation/variant and the testicular histological picture. CONCLUSIONS: Different mutations of the same gene cause various testicular histological pictures. These results suggest that it is not the gene itself but the type of mutation/variation that determines the testicular histology picture. Based on the data presented above, it remains challenging to design a genetic panel with prognostic value for the outcome of testicular sperm extraction in patients with NOA.

Genetic regulation of post-translational modification of two distinct proteins.

Post-translational modifications diversify protein functions and dynamically coordinate their signalling networks, influencing most aspects of cell physiology. Nevertheless, their genetic regulation or influence on complex traits is not fully understood. Here, we compare the genetic regulation of the same PTM of two proteins - glycosylation of transferrin and immunoglobulin G (IgG). By performing genome-wide association analysis of transferrin glycosylation, we identify 10 significantly associated loci, 9 of which were not reported previously. Comparing these with IgG glycosylation-associated genes, we note protein-specific associations with genes encoding glycosylation enzymes (transferrin - MGAT5, ST3GAL4, B3GAT1; IgG - MGAT3, ST6GAL1), as well as shared associations (FUT6, FUT8). Colocalisation analyses of the latter suggest that different causal variants in the FUT genes regulate fucosylation of the two proteins. Glycosylation of these proteins is thus genetically regulated by both shared and protein-specific mechanisms.

Glycosylation of IgG Associates with Hypertension and Type 2 Diabetes Mellitus Comorbidity in the Chinese Muslim Ethnic Minorities and the Han Chinese.

OBJECTIVES: Hypertension and type 2 diabetes mellitus comorbidity (HDC) is common, which confers a higher risk of cardiovascular disease than the presence of either condition alone. Describing the underlying glycomic changes of immunoglobulin G (IgG) that predispose individuals to HDC may help develop novel protective immune-targeted and anti-inflammatory therapies. Therefore, we investigated glycosylation changes of IgG associated with HDC. METHODS: The IgG N-glycan profiles of 883 plasma samples from the three northwestern Chinese Muslim ethnic minorities and the Han Chinese were analyzed by ultra-performance liquid chromatography instrument. RESULTS: We found that 12 and six IgG N-glycan traits showed significant associations with HDC in the Chinese Muslim ethnic minorities and the Han Chinese, respectively, after adjustment for potential confounders and false discovery rate. Adding the IgG N-glycan traits to the baseline models, the area under the receiver operating characteristic curves (AUCs) of the combined models differentiating HDC from hypertension (HTN), type 2 diabetes mellitus (T2DM), and healthy individuals were 0.717, 0.747, and 0.786 in the pooled samples of Chinese Muslim ethnic minorities, and 0.828, 0.689, and 0.901 in the Han Chinese, respectively, showing improved discriminating performance than both the baseline models and the glycan-based models. CONCLUSION: Altered IgG N-glycan profiles were shown to associate with HDC, suggesting the involvement of inflammatory processes of IgG glycosylation. The alterations of IgG N-glycome, illustrated here for the first time in HDC, demonstrate a biomarker potential, which may shed light on future studies investigating their potential for monitoring or preventing the progression from HTN or T2DM towards HDC.

Replication of 15 loci involved in human plasma protein N-glycosylation in 4802 samples from four cohorts.

Human protein glycosylation is a complex process, and its in vivo regulation is poorly understood. Changes in glycosylation patterns are associated with many human diseases and conditions. Understanding the biological determinants of protein glycome provides a basis for future diagnostic and therapeutic applications. Genome-wide association studies (GWAS) allow to study biology via a hypothesis-free search of loci and genetic variants associated with a trait of interest. Sixteen loci were identified by three previous GWAS of human plasma proteome N-glycosylation. However, the possibility that some of these loci are false positives needs to be eliminated by replication studies, which have been limited so far. Here, we use the largest set of samples so far (4802 individuals) to replicate the previously identified loci. For all but one locus, the expected replication power exceeded 95%. Of the 16 loci reported previously, 15 were replicated in our study. For the remaining locus (near the KREMEN1 gene), the replication power was low, and hence, replication results were inconclusive. The very high replication rate highlights the general robustness of the GWAS findings as well as the high standards adopted by the community that studies genetic regulation of protein glycosylation. The 15 replicated loci present a good target for further functional studies. Among these, eight loci contain genes encoding glycosyltransferases: MGAT5, B3GAT1, FUT8, FUT6, ST6GAL1, B4GALT1, ST3GAL4 and MGAT3. The remaining seven loci offer starting points for further functional follow-up investigation into molecules and mechanisms that regulate human protein N-glycosylation in vivo.

Evaluation of different PNGase F enzymes in immunoglobulin G and total plasma N-glycans analysis.

Glycoproteins, proteins that are co- and posttranslationally modified by sugars (glycans), have significant roles in pathophysiology of many different diseases. One of the main steps in sample preparation for free N-glycan analysis is deglycosylation or glycan removal. The aim of this study was to compare different peptide N-glycosidase F (PNGase F) enzymes (Rapid PNGase F and two recombinant versions) for deglycosylation of total human plasma glycoproteins and different amounts of human immunoglobulin G (IgG). Deglycosylation with different PNGase F enzymes resulted in different IgG and plasma N-glycosylation hydrophilic interaction liquid chromatography ultra-performance liquid chromatography profiles. Additionally, one recombinant version of PNGase F is more efficient in deglycosylation of complex N-glycans compared with Rapid PNGase F and recombinant version of PNGase F from a different manufacturer. In terms of chromatographic peak intensities and coefficient of variation %Area values, all tested versions of PNGase F enzymes were very reproducible and on the similar level when used in optimal conditions. However, care should be taken in terms of which enzyme is used with which protocol, particularly when scaling up.

Global variability of the human IgG glycome.

Immunoglobulin G (IgG) is the most abundant serum antibody which structural characteristics and effector functions are modulated through the attachment of various sugar moieties called glycans. Composition of the IgG N-glycome changes with age of an individual and in different diseases. Variability of IgG glycosylation within a population is well studied and is known to be affected by both genetic and environmental factors. However, global inter-population differences in IgG glycosylation have never been properly addressed. Here we present population-specific N-glycosylation patterns of IgG, analyzed in 5 different populations totaling 10,482 IgG glycomes, and of IgG's fragment crystallizable region (Fc), analyzed in 2,579 samples from 27 populations sampled across the world. Country of residence associated with many N-glycan features and the strongest association was with monogalactosylation where it explained 38% of variability. IgG monogalactosylation strongly correlated with the development level of a country, defined by United Nations health and socioeconomic development indicators, and with the expected lifespan. Subjects from developing countries had low levels of IgG galactosylation, characteristic for inflammation and ageing. Our results suggest that citizens of developing countries may be exposed to environmental factors that can cause low-grade chronic inflammation and the apparent increase in biological age.

Glycosylation of immunoglobulin G is regulated by a large network of genes pleiotropic with inflammatory diseases.

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

Decreased Immunoglobulin G Core Fucosylation, A Player in Antibody-dependent Cell-mediated Cytotoxicity, is Associated with Autoimmune Thyroid Diseases.

Autoimmune thyroid diseases (AITD) are the most common group of autoimmune diseases, associated with lymphocyte infiltration and the production of thyroid autoantibodies, like thyroid peroxidase antibodies (TPOAb), in the thyroid gland. Immunoglobulins and cell-surface receptors are glycoproteins with distinctive glycosylation patterns that play a structural role in maintaining and modulating their functions. We investigated associations of total circulating IgG and peripheral blood mononuclear cells glycosylation with AITD and the influence of genetic background in a case-control study with several independent cohorts and over 3,000 individuals in total. The study revealed an inverse association of IgG core fucosylation with TPOAb and AITD, as well as decreased peripheral blood mononuclear cells antennary α1,2 fucosylation in AITD, but no shared genetic variance between AITD and glycosylation. These data suggest that the decreased level of IgG core fucosylation is a risk factor for AITD that promotes antibody-dependent cell-mediated cytotoxicity previously associated with TPOAb levels.

Heritability of Human Plasma N-Glycome.

The N-glycosylation profile of total human plasma proteins could be a useful biomarker for various pathological states. Reliable high-throughput methods for such profiling have been developed. However, studies of relative importance of genetic and environmental factors in regulating plasma N-glycome are scarce. The aim of our study was to determine the role of genetic factors in phenotypic variation of plasma N-glycan profile through the estimates of its heritability. Thirty-nine total plasma N-glycome traits were analyzed in 2816 individuals from the TwinsUK data set. For the majority of the traits, high heritability estimates (>50%) were obtained pointing at a significant contribution of genetic factors in plasma N-glycome variation, especially for glycans mostly attached to immunoglobulins. We have also found several structures with higher environmental contribution to their variation.

Comprehensive N-glycosylation analysis of immunoglobulin G from dried blood spots.

Immunoglobulin G (IgG) glycans are emerging as a new putative biomarker for biological age and different diseases, requiring a robust workflow for IgG glycome analysis, ideally beginning with a simple and undemanding sampling procedure. Here, we report the first comprehensive study on total N-glycans of IgG isolated from dried blood spots (DBSs), which was performed in a high-throughput mode. We compared the IgG N-glycan profiles originating from DBS with those originating from plasma, compared different media for DBS collection, evaluated analytical variation and assessed IgG N-glycan profile stability for different storage conditions. In conclusion, we show that DBSs are a good and stable source material for a robust IgG N-glycan analysis by ultra-performance liquid chromatography, suitable for blood sampling in conditions where no trained personnel and necessary laboratory equipment are available.

Defining the genetic control of human blood plasma N-glycome using genome-wide association study.

Glycosylation is a common post-translational modification of proteins. Glycosylation is associated with a number of human diseases. Defining genetic factors altering glycosylation may provide a basis for novel approaches to diagnostic and pharmaceutical applications. Here we report a genome-wide association study of the human blood plasma N-glycome composition in up to 3811 people measured by Ultra Performance Liquid Chromatography (UPLC) technology. Starting with the 36 original traits measured by UPLC, we computed an additional 77 derived traits leading to a total of 113 glycan traits. We studied associations between these traits and genetic polymorphisms located on human autosomes. We discovered and replicated 12 loci. This allowed us to demonstrate an overlap in genetic control between total plasma protein and IgG glycosylation. The majority of revealed loci contained genes that encode enzymes directly involved in glycosylation (FUT3/FUT6, FUT8, B3GAT1, ST6GAL1, B4GALT1, ST3GAL4, MGAT3 and MGAT5) and a known regulator of plasma protein fucosylation (HNF1A). However, we also found loci that could possibly reflect other more complex aspects of glycosylation process. Functional genomic annotation suggested the role of several genes including DERL3, CHCHD10, TMEM121, IGH and IKZF1. The hypotheses we generated may serve as a starting point for further functional studies in this research area.

Plasma N-glycome composition associates with chronic low back pain.

BACKGROUND: Low back pain (LBP) is the symptom of a group of syndromes with heterogeneous underlying mechanisms and molecular pathologies, making treatment selection and patient prognosis very challenging. Moreover, symptoms and prognosis of LBP are influenced by age, gender, occupation, habits, and psychological factors. LBP may be characterized by an underlying inflammatory process. Previous studies indicated a connection between inflammatory response and total plasma N-glycosylation. We wanted to identify potential changes in total plasma N-glycosylation pattern connected with chronic low back pain (CLBP), which could give an insight into the pathogenic mechanisms of the disease. METHODS: Plasma samples of 1128 CLBP patients and 760 healthy controls were collected in clinical centers in Italy, Belgium and Croatia and used for N-glycosylation profiling by hydrophilic interaction ultra-performance liquid chromatography (HILIC-UPLC) after N-glycans release, fluorescent labeling and clean-up. Observed N-glycosylation profiles have been compared with a cohort of 126 patients with acute inflammation that underwent abdominal surgery. RESULTS: We have found a statistically significant increase in the relative amount of high-branched (tri-antennary and tetra-antennary) N-glycan structures on CLBP patients' plasma glycoproteins compared to healthy controls. Furthermore, relative amounts of disialylated and trisialylated glycan structures were increased, while high-mannose and glycans containing bisecting N-acetylglucosamine decreased in CLBP. CONCLUSIONS: Observed changes in CLBP on the plasma N-glycome level are consistent with N-glycosylation changes usually seen in chronic inflammation. GENERAL SIGNIFICANCE: To our knowledge, this is a first large clinical study on CLBP patients and plasma N-glycome providing a new glycomics perspective on potential disease pathology.

Profiling and genetic control of the murine immunoglobulin G glycome.

Immunoglobulin G (IgG) glycosylation is essential for function of the immune system, but the genetic and environmental factors that underlie its inter-individual variability are not well defined. The Collaborative Cross (CC) genetic resource harnesses over 90% of the common genetic variation of the mouse. By analyzing the IgG glycome composition of 95 CC strains, we made several important observations: (i) glycome variation between mouse strains was higher than between individual humans, despite all mice having the same environmental influences; (ii) five genetic loci were found to be associated with murine IgG glycosylation; (iii) variants outside traditional glycosylation site motifs affected glycome variation; (iv) bisecting N-acetylglucosamine (GlcNAc) was produced by several strains although most previous studies have reported the absence of glycans containing the bisecting GlcNAc on murine IgGs; and (v) common laboratory mouse strains are not optimal animal models for studying effects of glycosylation on IgG function.

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.

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 glycan patterns are associated with type 2 diabetes in independent European populations.

BACKGROUND: Type 2 diabetes results from interplay between genetic and acquired factors. Glycans on proteins reflect genetic, metabolic and environmental factors. However, associations of IgG glycans with type 2 diabetes have not been described. We compared IgG N-glycan patterns in type 2 diabetes with healthy subjects. METHODS: In the DiaGene study, a population-based case-control study, (1886 cases and 854 controls) 58 IgG glycan traits were analyzed. Findings were replicated in the population-based CROATIA-Korcula-CROATIA-Vis-ORCADES studies (162 cases and 3162 controls), and meta-analyzed. AUCs of ROC-curves were calculated using 10-fold cross-validation for clinical characteristics, IgG glycans and their combination. RESULTS: After correction for extensive clinical covariates, 5 IgG glycans and 13 derived traits significantly associated with type 2 diabetes in meta-analysis (after Bonferroni correction). Adding IgG glycans to age and sex increased the AUC from 0.542 to 0.734. Adding them to the extensive model did not substantially improve the AUC. The AUC for IgG glycans alone was 0.729. CONCLUSIONS: Several IgG glycans and traits firmly associate with type 2 diabetes, reflecting a pro-inflammatory and biologically-aged state. IgG glycans showed limited improvement of AUCs. However, IgG glycans showed good prediction alone, indicating they may capture information of combined covariates. The associations found may yield insights in type 2 diabetes pathophysiology. GENERAL SIGNIFICANCE: This work shows that IgG glycomic changes have biomarker potential and may yield important insights into pathophysiology of complex public health diseases, illustrated here for the first time in type 2 diabetes.

High-throughput analysis of immunoglobulin G glycosylation.

INTRODUCTION: Glycosylation of immunoglobulin G (IgG) is important for its effector functions and was shown to be related to age, sex and disease status of an individual. Adding glycomic information to genome-wide association studies (GWAS) and large clinical trials is enabling insight into the functional relevance of changes in glycosylation, as well as molecular mechanisms behind these changes. Large-scale studies require sensitive, robust and affordable high-throughput methodologies for glycosylation analysis, which are currently available in only a limited number of laboratories. AREAS COVERED: This review focuses on currently used high-throughput approaches for N-glycosylation analysis of IgG, as well as some recent advances in the areas of deglycosylation, trypsin digestion, labeling, purification, derivatization and automation of current workflows. Relevant literature was searched using the PubMed database. Expert commentary: Development, optimization and validation of robust, affordable and simple high-throughput glycosylation analysis methods is essential for discovery and validation of diagnostic and prognostic glycan biomarkers. Although significant advances in glycosylation analysis have been made in recent years, currently used protocols will have to be further optimized to enable subsequent analysis of glycosylation on all levels with the limited initial sample and in the minimal amount of time, which is still a challenging task.