Sulfation of the N-linked oligosaccharides of influenza virus hemagglutinin: temporal relationships and localization of sulfotransferases

Application of network smoothing to glycan LC-MS profiling

Motivation

Glycosylation is one of the most heterogeneous and complex protein post-translational modifications. Liquid chromatography coupled mass spectrometry (LC-MS) is a common high throughput method for analyzing complex biological samples. Accurate study of glycans require high resolution mass spectrometry. Mass spectrometry data contains intricate sub-structures that encode mass and abundance, requiring several transformations before it can be used to identify biological molecules, requiring automated tools to analyze samples in a high throughput setting. Existing tools for interpreting the resulting data do not take into account related glycans when evaluating individual observations, limiting their sensitivity.

Results

We developed an algorithm for assigning glycan compositions from LC-MS data by exploring biosynthetic network relationships among glycans. Our algorithm optimizes a set of likelihood scoring functions based on glycan chemical properties but uses network Laplacian regularization and optionally prior information about expected glycan families to smooth the likelihood and thus achieve a consistent and more representative solution. Our method was able to identify as many, or more glycan compositions compared to previous approaches, and demonstrated greater sensitivity with regularization. Our network definition was tailored to N-glycans but the method may be applied to glycomics data from other glycan families like O-glycans or heparan sulfate where the relationships between compositions can be expressed as a graph.

Availability and implementation Built Executable

http://www.bumc.bu.edu/msr/glycresoft/ and Source Code: https://github.com/BostonUniversityCBMS/glycresoft.

Supplementary information

Supplementary data are available at Bioinformatics online.

A Knowledge-Based System for Display and Prediction of O-Glycosylation Network Behaviour in Response to Enzyme Knockouts

O-linked glycosylation is an important post-translational modification of mucin-type protein, changes to which are important biomarkers of cancer. For this study of the enzymes of O-glycosylation, we developed a shorthand notation for representing GalNAc-linked oligosaccharides, a method for their graphical interpretation, and a pattern-matching algorithm that generates networks of enzyme-catalysed reactions. Software for generating glycans from the enzyme activities is presented, and is also available online. The degree distributions of the resulting enzyme-reaction networks were found to be Poisson in nature. Simple graph-theoretic measures were used to characterise the resulting reaction networks. From a study of in-silico single-enzyme knockouts of each of 25 enzymes known to be involved in mucin O-glycan biosynthesis, six of them, β-1,4-galactosyltransferase (β4Gal-T4), four glycosyltransferases and one sulfotransferase, play the dominant role in determining O-glycan heterogeneity. In the absence of β4Gal-T4, all Lewis X, sialyl-Lewis X, Lewis Y and Sda/Cad glycoforms were eliminated, in contrast to knockouts of the N-acetylglucosaminyltransferases, which did not affect the relative abundances of O-glycans expressing these epitopes. A set of 244 experimentally determined mucin-type O-glycans obtained from the literature was used to validate the method, which was able to predict up to 98% of the most common structures obtained from human and engineered CHO cell glycoforms.

Combined biochemical and cytological analysis of membrane trafficking using lectins

We have tested the application of high-mannose-binding lectins as analytical reagents to identify N-glycans in the early secretory pathway of HeLa cells during subcellular fractionation and cytochemistry. Post-endoplasmic reticulum (ER) pre-Golgi intermediates were separated from the ER on Nycodenz-sucrose gradients, and the glycan composition of each gradient fraction was profiled using lectin blotting. The fractions containing the post-ER pre-Golgi intermediates are found to contain a subset of N-linked α-mannose glycans that bind the lectins Galanthus nivalis agglutinin (GNA), Pisum sativum agglutinin (PSA), and Lens culinaris agglutinin (LCA) but not lectins binding Golgi-modified glycans. Cytochemical analysis demonstrates that high-mannose-containing glycoproteins are predominantly localized to the ER and the early secretory pathway. Indirect immunofluorescence microscopy revealed that GNA colocalizes with the ER marker protein disulfide isomerase (PDI) and the COPI coat protein β-COP. In situ competition with concanavalin A (ConA), another high-mannose specific lectin, and subsequent GNA lectin histochemistry refined the localization of N-glyans containing nonreducing mannosyl groups, accentuating the GNA vesicular staining. Using GNA and treatments that perturb ER-Golgi transport, we demonstrate that lectins can be used to detect changes in membrane trafficking pathways histochemically. Overall, we find that conjugated plant lectins are effective tools for combinatory biochemical and cytological analysis of membrane trafficking of glycoproteins.

Mass spectrometric analysis of neutral and anionic N-glycans from a Dictyostelium discoideum model for human congenital disorder of glycosylation CDG IL

The HL241 mutant strain of the cellular slime mold Dictyostelium discoideum is a potential model for human congenital disorder of glycosylation type IL (ALG9-CDG) and has been previously predicted to possess a lower degree of modification of its N-glycans with anionic moieties than the parental wild-type. In this study, we first showed that this strain has a premature stop codon in its alg9 mannosyltransferase gene compatible with the occurrence of truncated N-glycans. These were subject to an optimized analytical workflow, considering that the mass spectrometry of acidic glycans often presents challenges due to neutral loss and suppression effects. Therefore, the protein-bound N-glycans were first fractionated, after serial enzymatic release, by solid phase extraction. Then primarily single glycan species were isolated by mixed hydrophilic-interaction/anion-exchange or reversed-phase HPLC and analyzed using chemical and enzymatic treatments and MS/MS. We show that protein-linked N-glycans of the mutant are of reduced size as compared to those of wild-type AX3, but still contain core α1,3-fucose, intersecting N-acetylglucosamine, bisecting N-acetylglucosamine, methylphosphate, phosphate, and sulfate residues. We observe that a single N-glycan can carry up to four of these six possible modifications. Due to the improved analytical procedures, we reveal fuller details regarding the N-glycomic potential of this fascinating model organism.

Positive selection for gains of N-linked glycosylation sites in hemagglutinin during evolution of H3N2 human influenza A virus

The number of N-linked glycosylation sites in the globular head of hemagglutinin (HA) has increased during evolution of H3N2 human influenza A virus. Here natural selection operating on the gains of N-linked glycosylation sites was examined by using the single-site analysis and the single-substitution analysis. In the single-site analysis, positive selection was not inferred at the amino acid sites where the substitutions generating N-linked glycosylation sites were observed, but was detected at antigenic sites. In contrast, in the single-substitution analysis, positive selection was detected for the amino acid substitutions generating N-linked glycosylation sites. The single-site analysis and the single-substitution analysis appeared to be suitable for detecting recurrent and episodic natural selection, respectively. The gains of N-linked glycosylation sites were likely to be positively selected for the function of shielding antigenic sites from immune responses. At the antigenic sites, positive selection appeared to have operated not only on the radical substitution but also on the conservative substitution in terms of the charge of amino acids, suggesting that the antigenic drift is not a by-product of the evolution of receptor binding avidity in HA of human H3N2 virus.

Compensatory evolution of net-charge in influenza A virus hemagglutinin

The propagation of influenza A virus depends on the balance between the activities of hemagglutinin (HA) for binding to host cells and neuraminidase (NA) for releasing from infected cells (HA-NA balance). Since the host cell membrane and the sialic acid receptor are negatively charged, the amino acid substitutions increasing (charge+) and decreasing (charge−) the positive charge of HA subunit 1 (HA1) enhance and reduce, respectively, the binding avidity and affinity. The positive charge of HA1 in human influenza A virus bearing subtype H3N2 (A/H3N2 virus) was observed to have increased during evolution, but the evolutionary mechanism for this observation was unclear because this may disrupt the HA-NA balance. Here we show, from the phylogenetic analysis of HA for human A/H3N2 and A/H1N1 viruses, that the relative frequencies of charge+ and charge− substitutions were elevated on the branches where the number of N-glycosylation sites (NGS) increased and decreased, respectively, compared to those where the number of NGS did not change. On the latter branches, the net-charge of HA1 appeared to have been largely maintained to preserve its structure and function. Since the charge+ and charge− substitutions in HA1 have opposite effects to the gain and loss of NGS on the binding and release of the virus, the net-charge of HA1 may have evolved to compensate for the effect of the gain and loss of NGS, probably through changing the avidity. Apparently, the relative frequency of charge− substitutions in HA1 of A/H3N2 virus was elevated after the introduction of oseltamivir, and that of charge+ substitutions in HA1 of A/H1N1 virus was elevated after the spread of oseltamivir resistance. These observations may also be explained by the compensatory effect of the net-charge in HA1 on the NA activity for keeping the HA-NA balance.