Darwinian selection for sites of Asn-linked glycosylation in phylogenetically disparate eukaryotes and viruses

A quest for cytosolic sequons and their functions

Evolution shapes protein sequences for their functions. Here, we studied the moonlighting functions of the N-linked sequon NXS/T, where X is not P, in human nucleocytosolic proteins. By comparing membrane and secreted proteins in which sequons are well known for N-glycosylation, we discovered that cyto-sequons can participate in nucleic acid binding, particularly in zinc finger proteins. Our global studies further discovered that sequon occurrence is largely proportional to protein length. The contribution of sequons to protein functions, including both N-glycosylation and nucleic acid binding, can be regulated through their density as well as the biased usage between NXS and NXT. In proteins where other PTMs or structural features are rich, such as phosphorylation, transmembrane ɑ-helices, and disulfide bridges, sequon occurrence is scarce. The information acquired here should help understand the relationship between protein sequence and function and assist future protein design and engineering.

Defective N-glycosylation in tumor-infiltrating CD8+ T cells impairs IFN-γ-mediated effector function

T cell-mediated antitumor immunity is modulated, in part, by N-glycosylation. However, the interplay between N-glycosylation and the loss of effector function in exhausted T cells has not yet been fully investigated. Here, we delineated the impact of N-glycosylation on the exhaustion of tumor-infiltrating lymphocytes in a murine colon adenocarcinoma model, focusing on the IFN-γ-mediated immune response. We found that exhausted CD8+ T cells downregulated the oligosaccharyltransferase complex, which is indispensable for N-glycan transfer. Concordant N-glycosylation deficiency in tumor-infiltrating lymphocytes leads to loss of antitumor immunity. Complementing the oligosaccharyltransferase complex restored IFN-γ production and alleviated CD8+ T cell exhaustion, resulting in reduced tumor growth. Thus, aberrant glycosylation induced in the tumor microenvironment incapacitates effector CD8+ T cells. Our findings provide insights into CD8+ T cell exhaustion by incorporating N-glycosylation to understand the characteristic loss of IFN-γ, opening new opportunities to amend the glycosylation status in cancer immunotherapies.

Cell Entry and Unusual Replication of SARS-CoV-2

BACKGROUND

SARS-CoV-2 is the causative virus for the CoVID-19 pandemic that has frequently mutated to continue to infect and resist available vaccines. Emerging new variants of the virus have complicated notions of immunity conferred by vaccines versus immunity that results from infection. While we continue to progress from epidemic to endemic as a result of this collective immunity, the pandemic remains a morbid and mortal problem.

OBJECTIVE

The SARS-CoV-2 virus has a very complex manner of replication. The spike protein, one of the four structural proteins of the encapsulated virus, is central to the ability of the virus to penetrate cells to replicate. The objective of this review is to summarize these complex features of viral replication.

METHODS

A review of the recent literature was performed on the biology of SARS-CoV-2 infection from published work from PubMed and works reported to preprint servers, e.g., bioRxiv and medRxiv.

RESULTS AND CONCLUSION

The complex molecular and cellular biology involved in SARS-CoV-2 replication and the origination of >30 proteins from a single open reading frame (ORF) have been summarized, as well as the structural biology of spike protein, a critical factor in the cellular entry of the virus, which is a necessary feature for it to replicate and cause disease.

Evolutionary Dynamics of Indels in SARS-CoV-2 Spike Glycoprotein

SARS-CoV-2, responsible for the current COVID-19 pandemic that claimed over 5.0 million lives, belongs to a class of enveloped viruses that undergo quick evolutionary adjustments under selection pressure. Numerous variants have emerged in SARS-CoV-2, posing a serious challenge to the global vaccination effort and COVID-19 management. The evolutionary dynamics of this virus are only beginning to be explored. In this work, we have analysed 1.79 million spike glycoprotein sequences of SARS-CoV-2 and found that the virus is fine-tuning the spike with numerous amino acid insertions and deletions (indels). Indels seem to have a selective advantage as the proportions of sequences with indels steadily increased over time, currently at over 89%, with similar trends across countries/variants. There were as many as 420 unique indel positions and 447 unique combinations of indels. Despite their high frequency, indels resulted in only minimal alteration of N-glycosylation sites, including both gain and loss. As indels and point mutations are positively correlated and sequences with indels have significantly more point mutations, they have implications in the evolutionary dynamics of the SARS-CoV-2 spike glycoprotein.

Multiple N-glycans cooperate in balancing misfolded BRI1 secretion and ER retention

N-glycans play a protective or monitoring role according to the folding state of associated protein or the distance from structural defects. Asparagine-linked (Asn/N-) glycosylation is one of the most prevalent and complex protein modifications and the associated N-glycans play crucial roles on protein folding and secretion. The studies have shown that many glycoproteins hold multiple N-glycans, yet little is known about the redundancy of N-glycans on a protein. In this study, we used BRI1 to decipher the roles of N-glycans on protein secretion and function. We found that all 14 potential N-glycosylation sites on BRI1 were occupied with oligosaccharides. The elimination of single N-glycan had no obvious effect on BRI1 secretion or function except N¹⁵⁴-glycan, which resulted in the retention of BRI1 in the endoplasmic reticulum (ER), similar to the loss of multiple highly conserved N-glycans. To misfolded bri1, the absence of N-glycans next to local structural defects enhanced the ER retention and the artificial addition of N-glycan could help the misfolded bri1-GFPs exiting from the ER, indicating that the N-glycans might serve as steric hindrance to protect the structure defects from ER recognition. We also found that the retention of misfolded bri1-9 by lectins and chaperones in the ER relied on the presence of multiple N-glycans distal to the local defects. Our findings revealed that the N-glycans might play a protective or monitoring role according to the folding state of associated protein or the distance from structural defects.

Re-annotation of the Theileria parva genome refines 53% of the proteome and uncovers essential components of N-glycosylation, a conserved pathway in many organisms

BACKGROUND

The apicomplexan parasite Theileria parva causes a livestock disease called East coast fever (ECF), with millions of animals at risk in sub-Saharan East and Southern Africa, the geographic distribution of T. parva. Over a million bovines die each year of ECF, with a tremendous economic burden to pastoralists in endemic countries. Comprehensive, accurate parasite genome annotation can facilitate the discovery of novel chemotherapeutic targets for disease treatment, as well as elucidate the biology of the parasite. However, genome annotation remains a significant challenge because of limitations in the quality and quantity of the data being used to inform the location and function of protein-coding genes and, when RNA data are used, the underlying biological complexity of the processes involved in gene expression. Here, we apply our recently published RNAseq dataset derived from the schizont life-cycle stage of T. parva to update structural and functional gene annotations across the entire nuclear genome.

RESULTS

The re-annotation effort lead to evidence-supported updates in over half of all protein-coding sequence (CDS) predictions, including exon changes, gene merges and gene splitting, an increase in average CDS length of approximately 50 base pairs, and the identification of 128 new genes. Among the new genes identified were those involved in N-glycosylation, a process previously thought not to exist in this organism and a potentially new chemotherapeutic target pathway for treating ECF. Alternatively-spliced genes were identified, and antisense and multi-gene family transcription were extensively characterized.

CONCLUSIONS

The process of re-annotation led to novel insights into the organization and expression profiles of protein-coding sequences in this parasite, and uncovered a minimal N-glycosylation pathway that changes our current understanding of the evolution of this post-translational modification in apicomplexan parasites.

Asparagine-linked glycosylation is not directly coupled to protein translocation across the endoplasmic reticulum in Saccharomyces cerevisiae

Mammalian cells express two oligosaccharyltransferase complexes, STT3A and STT3B, that have distinct roles in N-linked glycosylation. The STT3A complex interacts directly with the protein translocation channel to mediate glycosylation of proteins using an N-terminal-to-C-terminal scanning mechanism. N-linked glycosylation of proteins in budding yeast has been assumed to be a cotranslational reaction. We have compared glycosylation of several glycoproteins in yeast and mammalian cells. Prosaposin, a cysteine-rich protein that contains STT3A-dependent glycosylation sites, is poorly glycosylated in yeast cells and STT3A-deficient human cells. In contrast, a protein with extreme C-terminal glycosylation sites was efficiently glycosylated in yeast by a posttranslocational mechanism. Posttranslocational glycosylation was also observed for carboxypeptidase Y-derived reporter proteins that contain closely spaced acceptor sites. A comparison of two recent protein structures indicates that the yeast OST is unable to interact with the yeast heptameric Sec complex via an evolutionarily conserved interface due to occupation of the OST binding site by the Sec63 protein. The efficiency of glycosylation in yeast is not enhanced for proteins that are translocated by the Sec61 or Ssh1 translocation channels instead of the Sec complex. We conclude that N-linked glycosylation and protein translocation are not directly coupled in yeast cells.

N-glycosylation Triggers a Dual Selection Pressure in Eukaryotic Secretory Proteins

Nearly one third of the eukaryotic proteome traverses the secretory pathway and most of these proteins are N-glycosylated in the lumen of the endoplasmic reticulum. N-glycans fulfill multiple structural and biological functions, and are crucial for productive folding of many glycoproteins. N-glycosylation involves the attachment of an oligosaccharide to selected asparagine residues in the sequence N-X-S/T (X ≠ P), a motif known as an N-glycosylation’sequon’. Mutations that create novel sequons can cause disease due to the destabilizing effect of a bulky N-glycan. Thus, an analogous process must have occurred during evolution, whenever ancestrally cytosolic proteins were recruited to the secretory pathway. Here, we show that during evolution N-glycosylation triggered a dual selection pressure on secretory pathway proteins: while sequons were positively selected in solvent exposed regions, they were almost completely eliminated from buried sites. This process is one of the sharpest evolutionary signatures of secretory pathway proteins, and was therefore critical for the evolution of an efficient secretory pathway.

Asparagine-Linked Glycans of Cryptosporidium parvum Contain a Single Long Arm, Are Barely Processed in the Endoplasmic Reticulum (ER) or Golgi, and Show a Strong Bias for Sites with Threonine

Cryptosporidium parvum causes severe diarrhea in infants in developing countries and in immunosuppressed persons, including those with AIDS. We are interested in the Asn-linked glycans (N-glycans) of C. parvum, because (1) the N-glycan precursor is predicted to contain five mannose and two glucose residues on a single long arm versus nine mannose and three glucose residues on the three-armed structure common in host N-glycans, (2) C. parvum is a rare eukaryote that lacks the machinery for N-glycan-dependent quality control of protein folding in the lumen of the Endoplasmic Reticulum (ER), and (3) ER and Golgi mannosidases, as well as glycosyltransferases that build complex N-glycans, are absent from the predicted proteome. The C. parvum N-glycans reported here, which were determined using a combination of collision-induced dissociation and electronic excitation dissociation, contain a single, unprocessed mannose arm ± terminal glucose on the trimannosyl chitobiose core. Upon nanoUPLC-MS/MS separation and analysis of the C. parvum tryptic peptides, the total ion and extracted oxonium ion chromatograms delineated 32 peptides with occupied N-glycan sites; these were derived from 16 glycoproteins. Although the number of potential N-glycan sites with Thr (NxT) is only about twice that with Ser (NxS), almost 90% of the occupied N-glycan sites contain NxT. The two most abundant C. parvum proteins modified with N-glycans were an immunodominant antigen on the surface of sporozoites (gp900) and the possible oocyst wall protein 1 (POWP1). Seven other glycoproteins with N-glycans were unique to C. parvum; five shared common ancestry with other apicomplexans; two glycoproteins shared common ancestry with many organisms. In summary, C. parvum N-glycans are remarkable for the absence of ER and Golgi modification and for the strong bias toward occupancy of N-glycan motifs containing Thr.

Integrated Omics and Computational Glycobiology Reveal Structural Basis for Influenza A Virus Glycan Microheterogeneity and Host Interactions

Despite sustained biomedical research effort, influenza A virus remains an imminent threat to the world population and a major healthcare burden. The challenge in developing vaccines against influenza is the ability of the virus to mutate rapidly in response to selective immune pressure. Hemagglutinin is the predominant surface glycoprotein and the primary determinant of antigenicity, virulence and zoonotic potential. Mutations leading to changes in the number of HA glycosylation sites are often reported. Such genetic sequencing studies predict at best the disruption or creation of sequons for N-linked glycosylation; they do not reflect actual phenotypic changes in HA structure. Therefore, combined analysis of glycan micro and macro-heterogeneity and bioassays will better define the relationships among glycosylation, viral bioactivity and evolution. We present a study that integrates proteomics, glycomics and glycoproteomics of HA before and after adaptation to innate immune system pressure. We combined this information with glycan array and immune lectin binding data to correlate the phenotypic changes with biological activity. Underprocessed glycoforms predominated at the glycosylation sites found to be involved in viral evolution in response to selection pressures and interactions with innate immune-lectins. To understand the structural basis for site-specific glycan microheterogeneity at these sites, we performed structural modeling and molecular dynamics simulations. We observed that the presence of immature, high-mannose type glycans at a particular site correlated with reduced accessibility to glycan remodeling enzymes. Further, the high mannose glycans at sites implicated in immune lectin recognition were predicted to be capable of forming trimeric interactions with the immune-lectin surfactant protein-D.

Mutational and functional analysis of N-linked glycosylation of envelope fusion protein F of Helicoverpa armigera nucleopolyhedrovirus

The envelope fusion (F) protein of baculoviruses is a heavily N-glycosylated protein that plays a significant role in the virus infection cycle. N-Linked glycosylation of virus envelope glycoprotein is important for virus envelope glycoprotein folding and its function in general. There are six predicted N-glycosylation sites in the F (HaF) protein of Helicoverpa armigera nucleopolyhedrovirus (HearNPV). The N-glycosylation site located in the F(2) subunit (N104) of HaF has been identified and functionally characterized previously (Long et al., 2007). In this study, the other five potential N-glycosylation sites located in the HaF1 subunit, namely, N293, N361, N526, N571 and N595, were analysed extensively to examine their N-glycosylation and relative importance to the function of HaF. The results showed that four of these five potential glycosylation sites in the F(1) subunit, N293, N361, N526 and N571, were N-glycosylated in F proteins of mature HearNPV budded viruses (BVs) but that N595 was not. In general, the conserved site N526 was critical to the functioning of HaF, as absence of N-glycosylation of N526 reduced the efficiency of HaF folding and trafficking, consequently decreased fusogenicity and modified the subcellular localization of HaF proteins, and thus impaired virus production and infectivity. The absence of N-glycosylation at other individual sites was found to have different effects on the fusogenicity and subcelluar distribution of HaF proteins in HzAM1 cells. In summary, N-glycosylation plays comprehensive roles in HaF function and virus infectivity, which is further discussed.

Anti-Retroviral Lectins Have Modest Effects on Adherence of Trichomonas vaginalis to Epithelial Cells In Vitro and on Recovery of Tritrichomonas foetus in a Mouse Vaginal Model

Trichomonas vaginalis causes vaginitis and increases the risk of HIV transmission by heterosexual sex, while Tritrichomonas foetus causes premature abortion in cattle. Our goals were to determine the effects, if any, of anti-retroviral lectins, which are designed to prevent heterosexual transmission of HIV, on adherence of Trichomonas to ectocervical cells and on Tritrichomonas infections in a mouse model. We show that Trichomonas Asn-linked glycans (N-glycans), like those of HIV, bind the mannose-binding lectin (MBL) that is part of the innate immune system. N-glycans of Trichomonas and Tritrichomonas bind anti-retroviral lectins (cyanovirin-N and griffithsin) and the 2G12 monoclonal antibody, each of which binds HIV N-glycans. Binding of cyanovirin-N appears to be independent of susceptibility to metronidazole, the major drug used to treat Trichomonas. Anti-retroviral lectins, MBL, and galectin-1 cause Trichomonas to self-aggregate and precipitate. The anti-retroviral lectins also increase adherence of ricin-resistant mutants, which are less adherent than parent cells, to ectocervical cell monolayers and to organotypic EpiVaginal tissue cells. Topical application of either anti-retroviral lectins or yeast N-glycans decreases by 40 to 70% the recovery of Tritrichomonas from the mouse vagina. These results, which are explained by a few simple models, suggest that the anti-retroviral lectins have a modest potential for preventing or treating human infections with Trichomonas.

A sweet code for glycoprotein folding

Glycoprotein synthesis is initiated in the endoplasmic reticulum (ER) lumen upon transfer of a glycan (Glc3Man9GlcNAc2) from a lipid derivative to Asn residues (N-glycosylation). N-Glycan-dependent quality control of glycoprotein folding in the ER prevents exit to Golgi of folding intermediates, irreparably misfolded glycoproteins and incompletely assembled multimeric complexes. It also enhances folding efficiency by preventing aggregation and facilitating formation of proper disulfide bonds. The control mechanism essentially involves four components, resident lectin-chaperones (calnexin and calreticulin) that recognize monoglucosylated polymannose protein-linked glycans, lectin-associated oxidoreductase acting on monoglucosylated glycoproteins (ERp57), a glucosyltransferase that creates monoglucosylated epitopes in protein-linked glycans (UGGT) and a glucosidase (GII) that removes the glucose units added by UGGT. This last enzyme is the only mechanism component sensing glycoprotein conformations as it creates monoglucosylated glycans exclusively in not properly folded glycoproteins or in not completely assembled multimeric glycoprotein complexes. Glycoproteins that fail to properly fold are eventually driven to proteasomal degradation in the cytosol following the ER-associated degradation pathway, in which the extent of N-glycan demannosylation by ER mannosidases play a relevant role in the identification of irreparably misfolded glycoproteins.

Effects of N-glycan precursor length diversity on quality control of protein folding and on protein glycosylation

Asparagine-linked glycans (N-glycans) of medically important protists have much to tell us about the evolution of N-glycosylation and of N-glycan-dependent quality control (N-glycan QC) of protein folding in the endoplasmic reticulum. While host N-glycans are built upon a dolichol-pyrophosphate-linked precursor with 14 sugars (Glc3Man9GlcNAc2), protist N-glycan precursors vary from Glc3Man9GlcNAc2 (Acanthamoeba) to Man9GlcNAc2 (Trypanosoma) to Glc3Man5GlcNAc2 (Toxoplasma) to Man5GlcNAc2 (Entamoeba, Trichomonas, and Eimeria) to GlcNAc2 (Plasmodium and Giardia) to zero (Theileria). As related organisms have differing N-glycan lengths (e.g. Toxoplasma, Eimeria, Plasmodium, and Theileria), the present N-glycan variation is based upon secondary loss of Alg genes, which encode enzymes that add sugars to the N-glycan precursor. An N-glycan precursor with Man5GlcNAc2 is necessary but not sufficient for N-glycan QC, which is predicted by the presence of the UDP-glucose:glucosyltransferase (UGGT) plus calreticulin and/or calnexin. As many parasites lack glucose in their N-glycan precursor, UGGT product may be identified by inhibition of glucosidase II. The presence of an armless calnexin in Toxoplasma suggests secondary loss of N-glycan QC from coccidia. Positive selection for N-glycan sites occurs in secreted proteins of organisms with N-glycan QC and is based upon an increased likelihood of threonine but not serine in the +2 position versus asparagine. In contrast, there appears to be selection against N-glycan length in Plasmodium and N-glycan site density in Toxoplasma. Finally, there is suggestive evidence for N-glycan-dependent ERAD in Trichomonas, which glycosylates and degrades the exogenous reporter mutant carboxypeptidase Y (CPY*).

The intrinsic and extrinsic effects of N-linked glycans on glycoproteostasis

Proteins that traffic through the eukaryotic secretory pathway are commonly modified with N-linked carbohydrates. These bulky amphipathic modifications at asparagines intrinsically enhance solubility and folding energetics through carbohydrate-protein interactions. N-linked glycans can also extrinsically enhance glycoprotein folding by using the glycoprotein homeostasis or 'glycoproteostasis' network, which comprises numerous glycan binding and/or modification enzymes or proteins that synthesize, transfer, sculpt and use N-linked glycans to direct folding and trafficking versus degradation and trafficking of nascent N-glycoproteins through the cellular secretory pathway. If protein maturation is perturbed by misfolding, aggregation or both, stress pathways are often activated that result in transcriptional remodeling of the secretory pathway in an attempt to alleviate the insult (or insults). The inability to achieve glycoproteostasis is linked to several pathologies, including amyloidoses, cystic fibrosis and lysosomal storage diseases. Recent progress on genetic and pharmacologic adaptation of the glycoproteostasis network provides hope that drugs of this mechanistic class can be developed for these maladies in the near future.

A predictive fitness model for influenza

The seasonal human influenza A/H3N2 virus undergoes rapid evolution, which produces significant year-to-year sequence turnover in the population of circulating strains. Adaptive mutations respond to human immune challenge and occur primarily in antigenic epitopes, the antibody-binding domains of the viral surface protein haemagglutinin. Here we develop a fitness model for haemagglutinin that predicts the evolution of the viral population from one year to the next. Two factors are shown to determine the fitness of a strain: adaptive epitope changes and deleterious mutations outside the epitopes. We infer both fitness components for the strains circulating in a given year, using population-genetic data of all previous strains. From fitness and frequency of each strain, we predict the frequency of its descendent strains in the following year. This fitness model maps the adaptive history of influenza A and suggests a principled method for vaccine selection. Our results call for a more comprehensive epidemiology of influenza and other fast-evolving pathogens that integrates antigenic phenotypes with other viral functions coupled by genetic linkage.

Antigenic analyses of highly pathogenic avian influenza a viruses

In response to the ongoing threat to animal and human health posed by HPAI endemic in poultry, Asia (H5N1) and North America (H7N3) have revived efforts to reduce pandemic risk by disease control at the source and improved pandemic vaccines. Discovery of conserved neutralization epitopes in the HA, which mediate broad protection within and across HA subtypes have changed the paradigm of "broadly reactive" or "universal" vaccine design. Development of such vaccines would benefit from comparative antigenic analysis of viruses with increasing divergence within (and between) HA subtypes. A review of recent work to define the antigenic properties of HPAI viruses revealed data generated through an array of experimental approaches. This information has supported diagnostics and vaccine development for animal and human health. Further harmonization of analytical methods is needed to determine the antigenic relationships among multiple lineages of rapidly evolving HPAI viruses.

Glycosylation of closely spaced acceptor sites in human glycoproteins

Asparagine-linked glycosylation of proteins by the oligosaccharyltransferase (OST) occurs when acceptor sites or sequons (N-x≠P-T/S) on nascent polypeptides enter the lumen of the rough endoplasmic reticulum. Metazoan organisms assemble two isoforms of the OST that have different catalytic subunits (STT3A or STT3B) and partially non-overlapping cellular roles. Potential glycosylation sites move past the STT3A complex, which is associated with the translocation channel, at the protein synthesis elongation rate. Here, we investigated whether close spacing between acceptor sites in a nascent protein promotes site skipping by the STT3A complex. Biosynthetic analysis of four human glycoproteins revealed that closely spaced sites are efficiently glycosylated by an STT3B-independent process unless the sequons contain non-optimal sequence features, including extreme close spacing between sequons (e.g. NxTNxT) or the presence of paired NxS sequons (e.g. NxSANxS). Many, but not all, glycosylation sites that are skipped by the STT3A complex can be glycosylated by the STT3B complex. Analysis of a murine glycoprotein database revealed that closely spaced sequons are surprisingly common, and are enriched for paired NxT sites when the gap between sequons is less than three residues.

Comparative glycomics analysis of influenza Hemagglutinin (H5N1) produced in vaccine relevant cell platforms

Hemagglutinin (HA) is the major antigen in influenza vaccines, and glycosylation is known to influence its antigenicity. Embryonated hen eggs are traditionally used for influenza vaccine production, but vaccines produced in mammalian and insect cells were recently licensed. This raises the concern that vaccines produced with different cell systems might not be equivalent due to differences in their glycosylation patterns. Thus, we developed an analytical method to monitor vaccine glycosylation through a combination of nanoLC/MS(E) and quantitative MALDI-TOF MS permethylation profiling. We then used this method to examine glycosylation of HAs from two different influenza H5N1 strains produced in five different platforms, including hen eggs, three different insect cell lines (High Five, expresSF+ and glycoengineered expresSF+), and a human cell line (HEK293). Our results demonstrated that (1) sequon utilization is not necessarily equivalent in different cell types, (2) there are quantitative and qualitative differences in the overall N-glycosylation patterns and structures produced by different cell types, (3) ∼20% of the N-glycans on the HAs produced by High Five cells are core α1,3-fucosylated structures, which may be allergenic in humans, and (4) our method can be used to monitor differences in glycosylation during the cellular glycoengineering stages of vaccine development.

Evaluation of the humoral and cellular immune responses elicited by the live attenuated and inactivated influenza vaccines and their roles in heterologous protection in ferrets

The humoral and cellular immune responses elicited by the trivalent live attenuated influenza vaccine (LAIV) and the trivalent inactivated influenza vaccine (TIV) were evaluated in the ferret model, using newly developed ferret immunological reagents and assays. In contrast to the TIV, which only induced immune responses in primed animals, LAIV induced strong influenza virus-specific serum antibody and T-cell responses in both naive and influenza-seropositive animals. The LAIV offered significant protection against a heterologous H1N1 virus challenge infection in the upper respiratory tract. Influenza virus-specific immunoglobulin A (IgA) and immunoglobulin G (IgG) antibody-secreting cells (ASCs) and influenza virus-specific CD4(+) and CD8(+) T cells were detected in the circulation and local paratracheal draining lymph nodes. The frequency of the influenza-specific ASCs in the local lymph nodes appeared to correlate with the degree of protection in the upper respiratory tract. The protection conferred by the LAIV could be attributed not only to the antibody response but also to the cell-mediated and local mucosal immune responses, particularly in naive ferrets. These findings may explain why the LAIV is immunologically superior and offers immediate protection after a single dose in children.

Extreme C-terminal sites are posttranslocationally glycosylated by the STT3B isoform of the OST

Glycosylation in the C-terminal 50–55 residues of proteins is mediated posttranslocationally by the STT3B isoform of oligosaccharyltransferase, with a preference for NXT sites.

Metazoan organisms assemble two isoforms of the oligosaccharyltransferase (OST) that have different catalytic subunits (STT3A or STT3B) and partially nonoverlapping roles in asparagine-linked glycosylation. The STT3A isoform of the OST is primarily responsible for co-translational glycosylation of the nascent polypeptide as it enters the lumen of the endoplasmic reticulum. The C-terminal 65–75 residues of a glycoprotein will not contact the translocation channel–associated STT3A isoform of the OST complex before chain termination. Biosynthetic pulse labeling of five human glycoproteins showed that extreme C-terminal glycosylation sites were modified by an STT3B-dependent posttranslocational mechanism. The boundary for STT3B-dependent glycosylation of C-terminal sites was determined to fall between 50 and 55 residues from the C terminus of a protein. C-terminal NXT sites were glycosylated more rapidly and efficiently than C-terminal NXS sites. Bioinformatics analysis of glycopeptide databases from metazoan organisms revealed a lower density of C-terminal acceptor sites in glycoproteins because of reduced positive selection of NXT sites and negative selection of NXS sites.

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.

Evidence for N-glycan shielding of antigenic sites during evolution of human influenza A virus hemagglutinin

After the emergence of influenza A viruses in the human population, the number of N-glycosylation sites (NGS) in the globular head region of hemagglutinin (HA) has increased continuously for several decades. It has been speculated that the addition of NGS to the globular head region of HA has conferred selective advantages to the virus by preventing the binding of antibodies (Ab) to antigenic sites (AS). Here, the effect of N-glycosylation on the binding of Ab to AS in human influenza A virus subtype H3N2 (A/H3N2) was examined by inferring natural selection at AS and other sites (NAS) that are located close to and distantly from the NGS in the three-dimensional structure of HA through a comparison of the rates of synonymous (d(S)) and nonsynonymous (d(N)) substitutions. When positions 63, 122, 126, 133, 144, and 246 in the globular head region of HA were non-NGS, the d(N)/d(S) was >1 and positive selection was detected at the AS located near these positions. However, the d(N)/d(S) value decreased and the evidence of positive selection disappeared when these positions became NGS. In contrast, d(N)/d(S) at the AS distantly located from the positions mentioned above and at the NAS of any location were generally <1 and did not decrease when these positions changed from non-NGS to NGS. These results suggest that the attachment of N-glycans to the NGS in the globular head region of HA prevented the binding of Ab to AS in the evolutionary history of human A/H3N2 virus.

Dynamics of glycoprotein charge in the evolutionary history of human influenza

Background

Influenza viruses show a significant capacity to evade host immunity; this is manifest both as large occasional jumps in the antigenic phenotype of viral surface molecules and in gradual antigenic changes leading to annual influenza epidemics in humans. Recent mouse studies show that avidity for host cells can play an important role in polyclonal antibody escape, and further that electrostatic charge of the hemagglutinin glycoprotein can contribute to such avidity.

Methodology/Principal Findings

We test the role of glycoprotein charge on sequence data from the three major subtypes of influenza A in humans, using a simple method of calculating net glycoprotein charge. Of all subtypes, H3N2 in humans shows a striking pattern of increasing positive charge since its introduction in 1968. Notably, this trend applies to both hemagglutinin and neuraminidase glycoproteins. In the late 1980s hemagglutinin charge reached a plateau, while neuraminidase charge started to decline. We identify key groups of amino acid sites involved in this charge trend.

Conclusions/Significance

To our knowledge these are the first indications that, for human H3N2, net glycoprotein charge covaries strongly with antigenic drift on a global scale. Further work is needed to elucidate how such charge interacts with other immune escape mechanisms, such as glycosylation, and we discuss important questions arising for future study.

Do N-glycoproteins have preference for specific sequons?

Protein N-glycosylation requires the presence of asparagine (N) in the consensus tri-peptide NXS/T (where X is any amino acid, S is serine and T is threonine). Several factors affect the glycosylation potential of NXS/T sequons and one such factor is the type of amino acid at position X. While proline was shown to negatively affect N-glycosylation, the nature of other amino acids at this position is not clear. Using Markov chain analysis of tri-peptide NXS/T from viral, archaeal and eukaryotic proteins as well as experimentally confirmed N-glycosylated sequons from eukaryotic proteins, we show here that the occurrence of most sequon types differ significantly from the expected probability. Sequon types with F, G, I, S, T and V amino acids are consistently preferred while those with P and charged amino acids are under-represented in all four groups. Further, proteins contained far fewer number of possible sequon types (maximum 20 types for NXS or NXT taken separately) for any given number of sequons, which may be explained based on random sampling. Consistent with the present finding, majority of the over-represented sequons found in two important viral envelope glycoproteins (hemagglutinin of influenza A H3N2 and glycoprotein120 of HIV-1) are indeed preferred sequon types, which may provide a selective advantage. Accordingly, although there seems to be some preference for sequons, this preference may not be unique to N-glycosylation.

The antiretroviral lectin cyanovirin-N targets well-known and novel targets on the surface of Entamoeba histolytica trophozoites

Entamoeba histolytica, the protist that causes amebic dysentery and liver abscess, has a truncated Asn-linked glycan (N-glycan) precursor composed of seven sugars (Man(5)GlcNAc(2)). Here, we show that glycoproteins with unmodified N-glycans are aggregated and capped on the surface of E. histolytica trophozoites by the antiretroviral lectin cyanovirin-N and then replenished from large intracellular pools. Cyanovirin-N cocaps the Gal/GalNAc adherence lectin, as well as glycoproteins containing O-phosphodiester-linked glycans recognized by an anti-proteophosphoglycan monoclonal antibody. Cyanovirin-N inhibits phagocytosis by E. histolytica trophozoites of mucin-coated beads, a surrogate assay for amebic virulence. For technical reasons, we used the plant lectin concanavalin A rather than cyanovirin-N to enrich secreted and membrane proteins for mass spectrometric identification. E. histolytica glycoproteins with occupied N-glycan sites include Gal/GalNAc lectins, proteases, and 17 previously hypothetical proteins. The latter glycoproteins, as well as 50 previously hypothetical proteins enriched by concanavalin A, may be vaccine targets as they are abundant and unique. In summary, the antiretroviral lectin cyanovirin-N binds to well-known and novel targets on the surface of E. histolytica that are rapidly replenished from large intracellular pools.

Subtle evolutionary changes in the distribution of N-glycosylation sequons in the HIV-1 envelope glycoprotein 120

Many viruses are known to undergo rapid evolutionary changes under selective pressures. The HIV-1 envelope glycoprotein 120 (gp120) shows extreme selection for NXS/T sequons, the potential sites of N-glycosylation. Although the average number of sequons in gp120 appears to be relatively stable in the recent past, even slight changes in the distribution of sequons may potentially play crucial roles in protein interaction and viral infection. This study tracked the prevalence and distribution of NXS/T sequons in gp120 over a period of 29 years (from 1981 to 2009). The gp120 showed location specific distribution of sequons with higher density in the outer domain of the molecule. The NXT sequon density decreased in the outer domain (despite the increase in the sequon specific amino acid threonine), but increased in the inner domain. By contrast, the NXS sequon density increased specifically in the outer domain. Related changes were also seen in the distribution probabilities of sequons in two domains. The results indicate that the gp120, chiefly in subtype B, is redistributing NXS/T sequons within the molecule with specific selection for NXS sequons. The subtle evolution of sequons in gp120 may have implications in viral resistance and infection.

A deeply divergent phosphoglucomutase (PGM) of Giardia lamblia has both PGM and phosphomannomutase activities

Giardia lamblia, which is an important parasitic cause of diarrhea, uses activated forms of glucose to make glycogen and activated forms of mannose to make glycophosphosphoinositol anchors. A necessary step for glucose activation is isomerization of glucose-6-phosphate to glucose-1-phosphate by a phosphoglucomutase (PGM). Similarly, a phosphomannomutase (PMM) converts mannose-6-phosphate to mannose-1-phosphate. While whole genome sequences of Giardia predict two PGM candidates, no PMM candidate is present. The hypothesis tested here is that at least one of the two Giardia PGM candidates has both PGM and PMM activity, as has been described for bacterial PGM orthologs. Nondenaturing gels showed that Giardia has two proteins with PGM activity, one of which also has PMM activity. Phylogenetic analyses showed that one of the two Giardia PGM candidates (Gl-PGM1) shares recent common ancestry with other eukaryotic PGMs, while the other Giardia PGM candidate (Gl-PGM2) is deeply divergent. Both Gl-PGM1 and Gl-PGM2 rescue a Saccharomyces cerevisiae pgm1Delta/pgm2Delta double deletion strain, while only Gl-PGM2 rescues a temperature-sensitive PMM mutant of S. cerevisiae (sec53-ts). Recombinant Gl-PGM1 has PGM activity only, whereas Gl-PGM2 has both PGM and PMM activities. We conclude that Gl-PGM1 behaves as a conventional eukaryotic PGM, while Gl-PGM2 is a novel eukaryotic PGM that also has PMM activity.

Evolutionary Pattern of N-Glycosylation Sequon Numbers in Eukaryotic ABC Protein Superfamilies

Many proteins contain a large number of NXS/T sequences (where X is any amino acid except proline) which are the potential sites of asparagine (N) linked glycosylation. However, the patterns of occurrence of these N-glycosylation sequons in related proteins or groups of proteins and their underlying causes have largely been unexplored. We computed the actual and probabilistic occurrence of NXS/T sequons in ABC protein superfamilies from eight diverse eukaryotic organisms. The ABC proteins contained significantly higher NXS/T sequon numbers compared to respective genome-wide average, but the sequon density was significantly lower owing to the increase in protein size and decrease in sequon specific amino acids. However, mammalian ABC proteins have significantly higher sequon density, and both serine and threonine containing sequons (NXS and NXT) have been positively selected—against the recent findings of only threonine specific Darwinian selection of sequons in proteins. The occurrence of sequons was positively correlated with the frequency of sequon specific amino acids and negatively correlated with proline and the NPS/T sequences. Further, the NPS/T sequences were significantly higher than expected in plant ABC proteins which have the lowest number of NXS/T sequons. Accordingly, compared to overall proteins, N-glycosylation sequons in ABC protein superfamilies have a distinct pattern of occurrence, and the results are discussed in an evolutionary perspective.

Suggestive evidence for Darwinian Selection against asparagine-linked glycans of Plasmodium falciparum and Toxoplasma gondii

We are interested in asparagine-linked glycans (N-glycans) of Plasmodium falciparum and Toxoplasma gondii, because their N-glycan structures have been controversial and because we hypothesize that there might be selection against N-glycans in nucleus-encoded proteins that must pass through the endoplasmic reticulum (ER) prior to threading into the apicoplast. In support of our hypothesis, we observed the following. First, in protists with apicoplasts, there is extensive secondary loss of Alg enzymes that make lipid-linked precursors to N-glycans. Theileria makes no N-glycans, and Plasmodium makes a severely truncated N-glycan precursor composed of one or two GlcNAc residues. Second, secreted proteins of Toxoplasma, which uses its own 10-sugar precursor (Glc(3)Man(5)GlcNAc(2)) and the host 14-sugar precursor (Glc(3)Man(9)GlcNAc(2)) to make N-glycans, have very few sites for N glycosylation, and there is additional selection against N-glycan sites in its apicoplast-targeted proteins. Third, while the GlcNAc-binding Griffonia simplicifolia lectin II labels ER, rhoptries, and surface of plasmodia, there is no apicoplast labeling. Similarly, the antiretroviral lectin cyanovirin-N, which binds to N-glycans of Toxoplasma, labels ER and rhoptries, but there is no apicoplast labeling. We conclude that possible selection against N-glycans in protists with apicoplasts occurs by eliminating N-glycans (Theileria), reducing their length (Plasmodium), or reducing the number of N-glycan sites (Toxoplasma). In addition, occupation of N-glycan sites is markedly reduced in apicoplast proteins versus some secretory proteins in both Plasmodium and Toxoplasma.