Motivation Glycosylation elaborates the structures and functions of glycoproteins; glycoproteins are common post-translationally modified proteins and are heterogeneous and non-deterministically synthesized as an evolutionarily driven mechanism that elaborates the functions of glycosylated gene products. Glycoproteins, accounting for approximately half of all proteins, require specialized proteomics data analysis methods due to micro- and macro-heterogeneities as a given glycosite can be divided into several glycosylated forms, each of which must be quantified. Sampling of heterogeneous glycopeptides is limited by mass spectrometer speed and sensitivity, resulting in missing values. In conjunction with the low sample size inherent to glycoproteomics, a specialized toolset is needed to determine if observed changes in glycopeptide abundances are biologically significant or due to data quality limitations. Results We developed an R package, Relative Assessment of m/z Identifications by Similarity (RAMZIS), that uses similarity metrics to guide researchers to a more rigorous interpretation of glycoproteomics data. RAMZIS uses a permutation test to generate contextual similarity, which assesses the quality of mass spectral data and outputs a graphical demonstration of the likelihood of finding biologically significant differences in glycosylation abundance datasets. Investigators can assess dataset quality, holistically differentiate glycosites, and identify which glycopeptides are responsible for glycosylation pattern change. RAMZIS is validated by theoretical cases and a proof-of-concept application. RAMZIS enables comparison between datasets too stochastic, small, or sparse for interpolation while acknowledging these issues in its assessment. Using this tool, researchers will be able to rigorously define the role of glycosylation and the changes that occur during biological processes. Availability and implementation https://github.com/WillHackett22/RAMZIS.

Glioblastoma (GBM) is the most common and malignant primary brain tumor. The extracellular matrix, also known as the matrisome, helps determine glioma invasion, adhesion, and growth. Little attention, however, has been paid to glycosylation of the extracellular matrix components that constitute the majority of glycosylated protein mass and presumed biological properties. To acquire a comprehensive understanding of the biological functions of the matrisome and its components, including proteoglycans (PGs) and glycosaminoglycans (GAGs), in GBM tumorigenesis, and to identify potential biomarker candidates, we studied the alterations of GAGs, including heparan sulfate (HS) and chondroitin sulfate (CS), the core proteins of PGs, and other glycosylated matrisomal proteins in GBM subtypes versus control human brain tissue samples. We scrutinized the proteomics data to acquire in-depth site-specific glycoproteomic profiles of the GBM subtypes that will assist in identifying specific glycosylation changes in GBM. We observed an increase in CS 6-O sulfation and a decrease in HS 6-O sulfation, accompanied by an increase in unsulfated CS and HS disaccharides in GBM versus control samples. Several core matrisome proteins, including PGs (decorin, biglycan, agrin, prolargin, glypican-1, and chondroitin sulfate proteoglycan 4), tenascin, fibronectin, hyaluronan link protein 1 and 2, laminins, and collagens, were differentially regulated in GBM versus controls. Interestingly, a higher degree of collagen hydroxyprolination was also observed for GBM versus controls. Further, two PGs, chondroitin sulfate proteoglycan 4 and agrin, were significantly lower, about 6-fold for isocitrate dehydrogenase-mutant, compared to the WT GBM samples. Differential regulation of O-glycopeptides for PGs, including brevican, neurocan, and versican, was observed for GBM subtypes versus controls. Moreover, an increase in levels of glycosyltransferase and glycosidase enzymes was observed for GBM when compared to control samples. We also report distinct protein, peptide, and glycopeptide features for GBM subtypes comparisons. Taken together, our study informs understanding of the alterations to key matrisomal molecules that occur during GBM development. (Data are available via ProteomeXchange with identifier PXD028931, and the peaks project file is available at Zenodo with DOI 10.5281/zenodo.5911810).

Complex protein glycosylation occurs through biosynthetic steps in the secretory pathway that create macro- and microheterogeneity of structure and function. Required for all life forms, glycosylation diversifies and adapts protein interactions with binding partners that underpin interactions at cell surfaces and pericellular and extracellular environments. Because these biological effects arise from heterogeneity of structure and function, it is necessary to measure their changes as part of the quest to understand nature. Quite often, however, the assumption behind proteomics that posttranslational modifications are discrete additions that can be modeled using the genome as a template does not apply to protein glycosylation. Rather, it is necessary to quantify the glycosylation distribution at each glycosite and to aggregate this information into a population of mature glycoproteins that exist in a given biological system. To date, mass spectrometric methods for assigning singly glycosylated peptides are well-established. But it is necessary to quantify glycosylation heterogeneity accurately in order to gauge the alterations that occur during biological processes. The task is to quantify the glycosylated peptide forms as accurately as possible and then apply appropriate bioinformatics algorithms to the calculation of micro- and macro-similarities. In this review, we summarize current approaches for protein quantification as they apply to this glycoprotein similarity problem.

The spike protein of SARS-CoV-2, the virus responsible for the global pandemic of COVID-19, is an abundant, heavily glycosylated surface protein that plays a key role in receptor binding and host cell fusion, and is the focus of all current vaccine development efforts. Variants of concern are now circulating worldwide that exhibit mutations in the spike protein. Protein sequence and glycosylation variations of the spike may affect viral fitness, antigenicity, and immune evasion. Global surveillance of the virus currently involves genome sequencing, but tracking emerging variants should include quantitative measurement of changes in site-specific glycosylation as well. In this work, we used data-dependent acquisition (DDA) and data-independent acquisition (DIA) mass spectrometry to quantitatively characterize the five N-linked glycosylation sites of the glycoprotein standard alpha-1-acid glycoprotein (AGP), as well as the 22 sites of the SARS-CoV-2 spike protein. We found that DIA compared favorably to DDA in sensitivity, resulting in more assignments of low-abundance glycopeptides. However, the reproducibility across replicates of DIA-identified glycopeptides was lower than that of DDA, possibly due to the difficulty of reliably assigning low-abundance glycopeptides confidently. The differences in the data acquired between the two methods suggest that DIA outperforms DDA in terms of glycoprotein coverage but that overall performance is a balance of sensitivity, selectivity, and statistical confidence in glycoproteomics. We assert that these analytical and bioinformatics methods for assigning and quantifying glycoforms would benefit the process of tracking viral variants as well as for vaccine development.

Protein glycosylation that mediates interactions among viral proteins, host receptors, and immune molecules is an important consideration for predicting viral antigenicity. Viral spike proteins, the proteins responsible for host cell invasion, are especially important to be examined. However, there is a lack of consensus within the field of glycoproteomics regarding identification strategy and false discovery rate (FDR) calculation that impedes our examinations. As a case study in the overlap between software, here as a case study, we examine recently published SARS-CoV-2 glycoprotein datasets with four glycoproteomics identification software with their recommended protocols: GlycReSoft, Byonic, pGlyco2, and MSFragger-Glyco. These software use different Target-Decoy Analysis (TDA) forms to estimate FDR and have different database-oriented search methods with varying degrees of quantification capabilities. Instead of an ideal overlap between software, we observed different sets of identifications with the intersection. When clustering by glycopeptide identifications, we see higher degrees of relatedness within software than within glycosites. Taking the consensus between results yields a conservative and non-informative conclusion as we lose identifications in the desire for caution; these non-consensus identifications are often lower abundance and, therefore, more susceptible to nuanced changes. We conclude that present glycoproteomics softwares are not directly comparable, and that methods are needed to assess their overall results and FDR estimation performance. Once such tools are developed, it will be possible to improve FDR methods and quantify complex glycoproteomes with acceptable confidence, rather than potentially misleading broad strokes.

Wild-type transthyretin-associated (ATTRwt) amyloidosis is an age-related disease that causes heart failure in older adults. This disease frequently features cardiac amyloid fibril deposits that originate from dissociation of the tetrameric protein, transthyretin (TTR). Unlike hereditary TTR (ATTRm) amyloidosis, where amino acid replacements destabilize the native protein, in ATTRwt amyloidosis, amyloid-forming TTR lacks protein sequence alterations. The initiating cause of fibril formation in ATTRwt amyloidosis is unclear, and thus, it seems plausible that other factors are involved in TTR misfolding and unregulated accumulation of wild-type TTR fibrils. We believe that clusterin (CLU, UniProtKB P10909), a plasma circulating glycoprotein, plays a role in the pathobiology of ATTRwt amyloidosis. Previously, we have suggested a role for CLU in ATTRwt amyloidosis based on our studies showing that (1) CLU codeposits with non-native TTR in amyloid fibrils from ATTRwt cardiac tissue, (2) CLU interacts only with non-native (monomeric and aggregated) forms of TTR, and (3) CLU serum levels in patients with ATTRwt are significantly lower compared to healthy controls. In the present study, we provide comprehensive detail of compositional findings from mass spectrometry analyses of amino acid and glycan content of CLU purified from ATTRwt and control sera. The characterization of oligosaccharide content in serum CLU derived from patients with ATTRwt amyloidosis is novel data. Moreover, results comparing CLU oligosaccharide variations between patient and healthy controls are original and provide further evidence for the role of CLU in ATTRwt pathobiology, possibly linked to disease-specific structural features that limit the chaperoning capacity of CLU.

Influenza A virus (IAV) mutates rapidly, resulting in antigenic drift and poor year-to-year vaccine effectiveness. One challenge in designing effective vaccines is that genetic mutations frequently cause amino acid variations in IAV envelope protein hemagglutinin (HA) that create new N-glycosylation sequons; resulting N-glycans cause antigenic shielding, allowing viral escape from adaptive immune responses. Vaccine candidate strain selection currently involves correlating antigenicity with HA protein sequence among circulating strains, but quantitative comparison of site-specific glycosylation information may likely improve the ability to design vaccines with broader effectiveness against evolving strains. However, there is poor understanding of the influence of glycosylation on immunodominance, antigenicity, and immunogenicity of HA, and there are no well-tested methods for comparing glycosylation similarity among virus samples. Here, we present a method for statistically rigorous quantification of similarity between two related virus strains that considers the presence and abundance of glycopeptide glycoforms. We demonstrate the strength of our approach by determining that there was a quantifiable difference in glycosylation at the protein level between WT IAV HA from A/Switzerland/9715293/2013 (SWZ13) and a mutant strain of SWZ13, even though no N-glycosylation sequons were changed. We determined site-specifically that WT and mutant HA have varying similarity at the glycosylation sites of the head domain, reflecting competing pressures to evade host immune response while retaining viral fitness. To our knowledge, our results are the first to quantify changes in glycosylation state that occur in related proteins of considerable glycan heterogeneity. Our results provide a method for understanding how changes in glycosylation state are correlated with variations in protein sequence, which is necessary for improving IAV vaccine strain selection. Understanding glycosylation will be especially important as we find new expression vectors for vaccine production, as glycosylation state depends greatly on the host species.