Improving Statistical Certainty of Glycosylation Similarity between Influenza A Virus Variants Using Data-Independent Acquisition Mass Spectrometry.

Amino acid sequences of immunodominant domains of hemagglutinin (HA) on the surface of influenza A virus (IAV) evolve rapidly, producing viral variants. HA mediates receptor recognition, binding and cell entry, and serves as the target for IAV vaccines. Glycosylation, a post-translational modification that places large branched polysaccharide molecules on proteins, can modulate the function of HA and shield antigenic regions allowing for viral evasion from immune responses. Our previous work showed that subtle changes in the HA protein sequence can have a measurable change in glycosylation. Thus, being able to quantitatively measure glycosylation changes in variants is critical for understanding how HA function may change throughout viral evolution. Moreover, understanding quantitatively how the choice of viral expression systems affects glycosylation can help in the process of vaccine design and manufacture. Although IAV vaccines are most commonly expressed in chicken eggs, cell-based vaccines have many advantages, and the adoption of more cell-based vaccines would be an important step in mitigating seasonal influenza and protecting against future pandemics. Here, we have investigated the use of data-independent acquisition (DIA) mass spectrometry for quantitative glycoproteomics. We found that DIA improved the sensitivity of glycopeptide detection for four variants of A/Switzerland/9715293/2013 (H3N2): WT and mutant, each expressed in embryonated chicken eggs and Madin-Darby canine kidney cells. We used the Tanimoto similarity metric to quantify changes in glycosylation between WT and mutant and between egg-expressed and cell-expressed virus. Our DIA site-specific glycosylation similarity comparison of WT and mutant expressed in eggs confirmed our previous analysis while achieving greater depth of coverage. We found that sequence variations and changing viral expression systems affected distinct glycosylation sites of HA. Our methods can be applied to track glycosylation changes in circulating IAV variants to bolster genomic surveillance already being done, for a more complete understanding of IAV evolution.

In-Depth Matrisome and Glycoproteomic Analysis of Human Brain Glioblastoma Versus Control Tissue.

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

Calculating Glycoprotein Similarities From Mass Spectrometric Data.

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.

Measuring Site-specific Glycosylation Similarity between Influenza a Virus Variants with Statistical Certainty.

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.

Data-independent acquisition mass spectrometry for site-specific glycoproteomics characterization of SARS-CoV-2 spike protein.

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.

The Need for Community Standards to Enable Accurate Comparison of Glycoproteomics Algorithm Performance.

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.

Glycosylation of Serum Clusterin in Wild-Type Transthyretin-Associated (ATTRwt) Amyloidosis: A Study of Disease-Associated Compositional Features Using Mass Spectrometry Analyses.

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.