INTRODUCTION
The surface hemagglutinin (HA) glycoprotein of influenza A virus (IAV) has a pivotal role in initiating the virus life cycle by binding to the virus receptor on target cells. HA binds to sialic acid (SA) residues that occur as terminal monosaccharides in glycoproteins and glycolipids on the cell surface. SA receptors engaged by IAV are linked to galactose (Gal) in an α2-3 (SAα2-3Gal) or α2-6 (SAα2-6Gal) linkage configuration (
1). Located in a small depression on the globular head of HA, the receptor binding site (RBS) is composed of the 130-loop, the 150-loop, the 190-helix, and the 220-loop. A series of conserved residues, including Tyr98, Ser136, Trp153, and His183 (H3 numbering), forms the base of the RBS and is important for SA interaction (
2,
3). Although some HA residues on the RBS are critical for receptor specificity (
4–6), other residues may influence binding by modulating virus-receptor binding avidity (
7,
8).
IAV remains an important pathogen for humans and swine (
9). While influenza vaccines are commercially available, the relative effectiveness of these vaccines is heavily dependent on the antigenic match of vaccine strains to circulating virus strains (
10,
11). Most of the humoral immune response elicited by influenza vaccination or natural exposure is directed against HA to block virus infection (
12,
13). Through a mechanism known as antigenic drift, IAV can circumvent the preexisting antibody response by generating genetic variants during replication with amino acid substitutions in key HA epitopes on the globular head, allowing for the emergence of escape mutant viruses (
14). Antigenic drift is a frequent cause of reduced vaccine effectiveness, especially for H3N2 IAVs (
15).
Defining the molecular basis of antigenic drift has important implications for understanding IAV evolution and has been facilitated by methodological advances, such as antigenic cartography (
16). A recent study identified seven residues (residues 145, 155, 156, 158, 159, 189, and 193) on the globular head of the HA as the major determinants of antigenic drift during the evolution of human H3N2 IAVs (
17). Interestingly, all seven residues were located adjacent to the RBS (
17), and these residues therefore may also influence receptor binding. The importance of this small set of HA residues as major drivers of antigenic evolution has been demonstrated for IAVs circulating in other hosts, including swine H3N2 IAVs (
18).
Antigenic changes in human H3N2 IAVs over the course of time were shown to be frequently caused by a single amino acid substitution in one of the seven residues, with specific substitutions involved in antigenic change more than once (
17). One of those substitutions, N 145 K, was shown in two separate instances to be the major determinant of antigenic drift during the evolution of human H3N2 IAVs (
16,
17). While there is evidence that N 145 K caused large antigenic changes (
7,
16–18) and also altered receptor binding avidity (
7), it remains unclear whether there are functional constraints to substitutions at residue 145 or how alternative amino acid substitutions at this residue may affect receptor binding and antigenic properties.
In the present study, we carried out site-directed mutagenesis to prepare a panel of H3 HA mutant viruses carrying a single amino acid substitution at residue 145 and subsequently evaluated the impact of substitutions at this key HA residue on receptor binding and antibody recognition. Our data indicate that residue 145 displayed remarkable amino acid plasticity in vitro, tolerating multiple amino acid substitutions, many of which have not yet been observed in nature. Mutant viruses carrying substitutions at residue 145 showed no major impairment in virus replication. While all substitution mutants retained binding to SAα2-6Gal glycans compared to the wild type, mutant viruses with substitutions not commonly found in nature displayed diminished receptor binding capacity. Antigenic characterization confirmed the impact of HA residue 145 K in antibody immunodominance. These findings have important implications for understanding virus evolution and aiding the development of novel vaccine design approaches.
DISCUSSION
HA engagement with terminal SA residues on host cells is an essential step in the IAV replication cycle. Not surprisingly, receptor binding specificity is a major host range restriction factor. In a simplistic view, IAVs of avian origin prefer binding to SAα2-3Gal glycans, whereas those that circulate in humans favor binding to SAα2-6Gal receptors (
5,
26–28). During human adaptation of avian-origin IAVs, the receptor specificity switch is accompanied by specific substitutions on the HA (E 190 D/G 225 D in the H1 subtype and Q 226 L/G 228 S in the H2 and H3 subtypes) (
5,
26,
29,
30). While these residues are critical to receptor binding specificity, other residues on the HA have the potential to affect receptor binding functions, including receptor binding avidity (
7,
8,
23,
31). The identification of up to 7 residues (residues 145, 155, 156, 158, 159, 189, and 193) near the RBS as the major determinants of antigenic drift during the evolution of human and swine H3N2 IAVs (
17,
18,
32) has led to the hypothesis that emerging substitutions in these residues must drive antigenic change and immune escape without disrupting receptor binding properties, potentially limiting the flexibility of amino acid residues in these positions (
17,
33). In agreement with this hypothesis, analysis of HA sequences showed the occurrence of only a small subset of amino acid substitutions at these residues that are recycled sporadically (
17,
18,
32,
34,
35).
Here we examined the amino acid plasticity at residue 145, one of the identified key determinants of antigenic drift, by mutating this residue to create a panel of H3 HA mutant viruses carrying a single amino acid substitution, representing every possible amino acid in a swine H3N2 backbone. HA residue 145 showed extraordinary plasticity, with 16 out of 19 substitutions (145 A, 145 C, 145 F, 145 G, 145 H, 145 I, 145 K, 145 L, 145 M, 145 P, 145 Q, 145 R, 145 S, 145 T, 145 V, and 145 Y) being well tolerated
in vitro (with 90 to 100% of sequenced reads possessing the mutated codon at residue 145). Whether these substitutions are tolerated
in vivo is beyond the scope of the present report and remain to be determined. Consistent with the intricate balance between immune escape and receptor binding, three substitutions (145 D, 145 E, and 145 W) were not tolerated and led to the emergence of partial reversion at residue 145 and additional substitutions in HA (T 128 A) and/or NA (T 148 I). HA T 128 A disrupted a potential glycosylation site (
36,
37) and was shown to emerge along with a substitution at residue 145 in human H3N2 IAVs during the 2013 to 2014 season (
34). NA T 148 I has been associated with reduced NA activity, decreased susceptibility to neuraminidase inhibitors, and NA-mediated binding and agglutination (
38,
39). The previously reported role of these substitutions suggests that they may be compensatory.
With the exception of 145 C virus, there was no discernible difference in the abilities of the tested mutant viruses to agglutinate turkey RBCs. Turkey and chicken RBCs carry both SAα2-3Gal and SAα2-6Gal glycans on their cell surface (
21–23), and some of tested mutant viruses exhibited an impaired ability to agglutinate chicken RBCs. Loss of the ability to agglutinate RBCs has plagued the antigenic characterization of recent human H3N2 IAVs (
40). As previously reported for chicken RBCs (
41), an in-depth analysis of major glycan structures present on RBCs from other hosts will help shed light on current issues on antigenic characterization of IAVs.
Using independent assays, we demonstrated that all tested mutant viruses retained some binding to SAα2-6Gal glycans. However, viruses possessing naturally occurring substitutions (145 N [wt], 145 K, and 145 S) showed the strongest receptor binding. With the exception of 145 M virus, all viruses carrying alternative substitutions (145 A, 145 C, 145 F, 145 G, 145 H, 145 I, 145 L, 145 P, 145 Q, 145 R, 145 T, 145 V, and 145 Y) displayed a decreased ability to bind SAα2-6Gal receptors. In contrast to other reports, no compensatory substitutions were identified on the HA or NA by next-generation sequencing (
42,
43). This indicates that the observed changes in receptor binding reported here are directly related to substitutions at residue 145. It is important to emphasize that although viruses carrying alternative substitutions exhibited decreased receptor binding avidity
in vitro, the lower threshold for biologically relevant avidity
in vivo is unknown. Glycan array analysis revealed a broad range of SAα2-6Gal glycan interactions that were modulated by substitutions at residue 145. Interestingly, one of the viruses with the lowest binding to 6-modified fetuin glycans in the panel (145 F virus) displayed extensive binding to SAα2-6Gal glycans, including extended glycans (
25). While the expanded glycan array provides important insights on the binding preference, the question of which receptors are relevant for IAV attachment and transmission
in vivo remains to be elucidated. Furthermore, the density, distribution, and organization of glycans on host tissues are poorly defined.
As observed in human H3N2 IAVs, the emergence of the HA N 145 K substitution in swine H3N2 IAVs led to a change in the antibody immunodominance. It is astonishing that reactivity to sera raised against a virus possessing HA 145 N was indistinguishable among all of the viruses carrying substitutions at residue 145, but reactivity to sera raised against a virus possessing HA 145 K was highly skewed to recognize an epitope bearing HA residue 145 K even in the context of a distantly related H3 HA, demonstrating the dramatic impact that N 145 K mutants can have on population immunity. Antibody immunodominance may be key for understanding antigenic drift and refers to the immunological phenomenon in which the immune system preferentially mounts a response to complex antigens in a dynamic hierarchical order (
14). Immunodominance hierarchy can occur at the level of viruses within multivalent immunogens, proteins within viruses, antigenic sites within proteins, epitopes within antigenic sites, and, as shown in the present report and by others, single amino acid substitutions within epitopes (
7,
14,
16,
17,
44,
45). The impact of substitutions at residue 145 on immunogenicity and antibody immunodominance remains to be fully characterized. While there may be a fitness cost associated with the emergence of such substitutions in nature, the ability to manipulate the amino acid plasticity near the RBS may offer an alternative approach to induce broader protection and, perhaps, potentially elicit receptor mimicry and broad-spectrum neutralizing antibodies (
46).
The complex and continuous antigenic evolution of IAVs remains a major hurdle for vaccine selection and effective vaccination. We provide a better understanding of the functional effects of amino acid substitutions near the RBS that are implicated in antigenic drift and their consequences for receptor binding. The mutation analyses presented in this report represent a significant data set to aid and test the ability of computational approaches to predict binding of glycans.
MATERIALS AND METHODS
Cell culture.
Madin-Darby canine kidney (MDCK) and human embryonic kidney 293T cells were maintained in Dulbecco’s modified Eagle’s medium (DMEM) supplemented with 10% fetal bovine serum (FBS). Cells were propagated at 37˚C in a humidified incubator under a 5% CO2 atmosphere.
Molecular cloning and virus rescue.
The OH/04 wt strain is a prototypic swine-origin virus amenable to genetic manipulation by an established reverse genetics (RG) system (
18,
47,
48). Single amino acid substitutions representing each of the 20 naturally occurring amino acids (
Table 2) were inserted by site-directed mutagenesis at the codon corresponding to HA residue 145 (H3 numbering). To prevent the carryover of OH/04 wt HA plasmid DNA during PCR amplification, the cloned OH/04 wt HA segment was split into two overlapping plasmids, pDP-SD1 and pDP-SD2. The pDP-SD1 plasmid carries the mouse RNA polymerase I terminator (t1) (
20) followed by nucleotides 1 to 522 of the OH/04 wt HA segment. The pDP-SD2 plasmid contains nucleotides 500 to 1762 of the OH/04 wt HA segment followed by the human RNA polymerase I promoter (polI) (
20). To introduce the desired mutations into OH/04 wt HA, a PCR fragment containing t1 followed by the 5′ portion of the HA was generated with forward primer 5′-ACC GGA GTA CTG GTC GAC CTC CGA AGT TGG GGG GGA GCA AAA GCA GG-3′ and the respective reverse primer (
Table 2) using pDP-SD1 as the template. Similarly, a PCR fragment comprising the 3′ portion of HA followed by polI was amplified from pDP-SD2 using reverse primer 5′-ATG CTG ACA ACG TCC CCG GCC CGG CGC TGC T-3′ and the respective forward primer (
Table 2). All PCR products were purified by gel extraction using the QIAquick gel extraction kit (Qiagen, Valencia, CA) and combined to produce RG-ready PCR-based HA segments for individual OH/04 HA aa 145 single mutants by overlapping PCR as previously described (
19). All PCRs were performed with Phusion high-fidelity DNA polymerase (New England Biolabs, Ipswich, MA), and products were confirmed to be free of unwanted mutations by sequencing. Viruses were rescued by PCR-based RG using a coculture of 293T and MDCK cells as previously described (
19,
20). To generate OH/04 HA aa 145 single mutant viruses, the respective RG PCR-based HA segment was paired with the seven plasmids representing the remaining OH/04 wt gene segments. Following transfection, cells were incubated at 35˚C. After 24 h of incubation, medium was replaced with Opti-MEM I (Life Technologies, Carlsbad, CA) containing 1 μg/ml TPCK (tosylsulfonyl phenylalanyl chloromethyl ketone)-trypsin (Worthington Biochemicals, Lakewood, NJ) and 1% antibiotic-antimycotic solution (Sigma-Aldrich, St. Louis, MO). Following virus rescue, virus stocks were amplified in MDCK cells. Virus stocks were titrated by TCID
50, and virus titers were determined by the method of Reed and Muench (
49).
Whole-genome sequencing.
Virus RNA from tissue culture supernatant virus stocks was purified using the RNeasy minikit (Qiagen, Valencia, CA) or the MagNA Pure LC RNA isolation kit (Roche Life Science, Mannheim, Germany). Isolated virus RNA served as the template in a one-step reverse transcriptase PCR (RT-PCR) for multisegment, whole-genome amplification (
50). Amplicon sequence libraries were prepared as previously described (
50) or using the Nextera XT DNA Library Prep kit (Illumina, San Diego, CA) according to the manufacturer’s protocol. Barcoded libraries were multiplexed and sequenced on the high-throughput Illumina MiSeq sequencing platform in a paired-end 150-nucleotide (nt) run format. De novo genome assembly was performed as described previously (
50), and HA- and NA-specific reads were mapped to OH/04 reference sequences using Geneious 10.1.3 (
51).
In vitro growth kinetics.
Confluent monolayers of MDCK cells were inoculated at a multiplicity of infection (MOI) of 0.01 for each virus. After 1 h of incubation, the virus inoculum was removed and cells washed twice with 1× phosphate-buffered saline (PBS). Opti-MEM I (Life Technologies, Carlsbad, CA) containing TPCK-trypsin (Worthington Biochemicals, Lakewood, NJ) and an antibiotic-antimycotic solution (Sigma-Aldrich, St. Louis, MO) was then added to the cells. At the indicated time points, tissue culture supernatant from inoculated cells was collected for virus titer quantification. Virus RNA from tissue culture supernatant was isolated using the MagMAX-96 AI/ND viral RNA isolation kit (Thermo Fisher Scientific, Waltham, MA). Virus titers were determined using a real-time reverse transcriptase PCR (rRT-PCR) assay based on the influenza A matrix gene (
52). The rRT-PCR was performed in a LightCycler 480 real-time PCR instrument (Roche Diagnostics, Rotkreuz, Switzerland) using the LightCycler 480 RNA master hydrolysis probe kit (Roche Life Science, Mannheim, Germany). A standard curve was generated using 10-fold serial dilutions from an OH/04 wt virus stock of known titer to correlate quantitative PCR (qPCR) crossing-point (Cp) values with virus titers, as previously described (
53). Virus titers were expressed as log
10 TCID
50/ml equivalents.
HA assay.
Chicken (Poultry Diagnostic and Research Center, Athens, GA), turkey (Poultry Diagnostic and Research Center, Athens, GA), and horse (Lampire Biologicals, Pipersville, PA) red blood cells (RBCs) were prepared from whole blood preparations using standard techniques (
54). Virus HA assays were carried out using 0.5% (vol/vol in PBS) chicken RBCs, 0.5% (vol/vol in PBS) turkey RBCs, or 1% (vol/vol in PBS) horse RBCs. Briefly, 50 μl of RBC preparations was added to 50 μl of 2-fold serial dilutions of virus stocks and allowed to incubate at room temperature for 45 min. After incubation, agglutination was measured, and data were expressed as the inverse of the highest dilution that allowed full agglutination.
Receptor cell binding assay.
The receptor cell binding assay was performed as previously described (
8,
18,
55). Briefly, 10% (vol/vol in PBS) turkey RBCs were pretreated for 1 h at 37°C with 2-fold serial dilutions of bacterial neuraminidase from
Arthrobacter ureafaciens (New England BioLabs, Ipswich, MA) or
Clostridium perfringens (New England BioLabs, Ipswich, MA). Following neuraminidase treatment, treated RBCs were washed twice with cold PBS and then resuspended to 1% (vol/vol in PBS), and then 50 μl of 1% RBCs treated with the different neuraminidase concentrations was added to 50 μl of each virus (8 HA units [HAU], as determined on untreated RBCs) and allowed to incubate at room temperature for 45 min. After incubation, agglutination was measured, and data were expressed as the maximal concentration of neuraminidase that allowed full agglutination.
Solid-phase assay of receptor binding specificity.
The receptor-binding specificity was determined in a solid-phase direct binding assay using monospecific preparations of horseradish peroxidase (HRP)-conjugated fetuin (Fet-HRP). Monospecific preparations of Fet-HRP were synthesized using α-2,3-sialyltransferase from
Pasteurella multocida (Sigma, St. Louis, MO) for 3-modified fetuin (3-Fet-HRP) or α-2,6-sialyltransferase from
Photobacterium damsela (Sigma, St. Louis, MO) for 6-modified fetuin (6-Fet-HRP), essentially as described previously (
24). Ninety-six-well native fetuin-coated flat-bottom plates (Greiner Bio-One, Monroe, NC) were incubated overnight at 4°C with 128 HAU of each virus in 0.02 M Tris-buffered saline (TBS), pH 7.2 to 7.4. Virus samples were run in duplicate. Plates were washed three times with PBS and blocked with blocking solution (BS) (PBS containing 0.1% neuraminidase-treated bovine serum albumin [BSA-NA]) for 2 h at room temperature. After blocking, plates were washed twice with ice-cold washing solution (WS) (PBS containing 0.02% Tween 80) and incubated with 2-fold serial dilution of 3-Fet-HRP (SAα2-3Gal) or 6-Fet-HRP (SAα2-6Gal) in reaction solution (RS) (PBS containing 0.02% Tween 80, 0.1% BSA-NA, and 2 μM oseltamivir carboxylate) for 1 h at 4°C. After incubation, plates were washed five times with ice-cold WS before adding freshly prepared substrate solution (SS) (0.01% 3,3′,5,5′-tetramethylbenzidine in 0.05 M sodium acetate with 0.03% H
2O
2). Reactions were allowed to proceed at room temperature for 30 min unless otherwise stated. Reactions were stopped with 3% (vol/vol in double-distilled water [ddH
2O]) H
2SO
4. Absorbance readings were obtained at 450 nm using a Victor ×3 multilabel plate reader (PerkinElmer, Waltham, MA).
Glycan array analysis.
Virus stocks were grown in MDCK cells, clarified by low-speed centrifugation, and inactivated by treatment with 0.1% β-propiolactone (BPL) for 1 day at 4°C. Inactivated virus stocks were concentrated as described previously (
56). Concentrated virus stocks were resuspended in PBS with 5% glycerol to 1280 to 5120 HAU, aliquoted, and stored at −80°C. HA titers of concentrated virus stocks were determined by HA assay using 0.5% (vol/vol in PBS) turkey RBCs. Glycan array analysis was performed using an
N-hydroxysuccinimide (NHS) ester-coated glass microarray slide containing six replicates of 128 diverse synthetic sialic acid-containing glycans, including terminal sequences as well as N-linked and O-linked glycans found on mammalian and avian glycoproteins and glycolipids (
25). Whole influenza virus samples were diluted to 256 HAU in PBS containing 3% BSA (PBS-BSA) and incubated on the array surface for 1 h at room temperature in a humidity-controlled chamber. After incubation, slides were washed in PBS and incubated with OH/04 wt swine antiserum (
18) diluted 1:200 in PBS-BSA for 1 h. Slides were washed in PBS and incubated for 1 h in goat anti-pig IgG conjugated to fluorescein isothiocyanate (FITC) (Thermo Fisher Scientific, Waltham, MA) diluted in PBS-BSA (20 μg/ml, final concentration). Slides were washed twice in PBS and in distilled water (dH
2O) and then dried prior to detection. Slide scanning to detect bound virus was conducted using an InnoScan 1100AL (Innopsys, Carbonne, France) fluorescent microarray scanner. Fluorescent signal intensity was measured using Mapix (Innopsys, Carbonne, France), and the mean intensity minus mean background of 4 replicate spots was calculated. A complete list of the glycans present in the array is presented in Table S1 in the supplemental material. The array is comprised of nonsialoside control (1 to 10; gray), SAα2-3Gal (11 to 76; yellow), and SAα2-6Gal (77 to 128; green) glycans. Glycans are grouped by structure type: L, linear; O, O linked; N, N linked; and L
x, sialyl Le
x.
Antisera.
Swine antisera against A/swine/New York/A01104005/2011 (H3N2) (NY/11), A/swine/Iowa/A01480656/2014 (H3N2) (IA/14), and A/turkey/Ohio/313053/2004 (H3N2) (OH/04) were generated in previous studies (
18,
32). For each virus, two pigs were primed and boosted intramuscularly with UV-inactivated whole-virus vaccine prepared with commercial oil-in-water adjuvant (1:5 ratio) (Emulsigen D; MVP Laboratories, Inc., Ralston, NE). The pigs were humanely euthanized for blood collection.
HI assay.
Hemagglutination inhibition (HI) assays were performed as previously described (
54). Prior to HI testing, swine antisera were treated overnight with receptor-destroying enzyme (Denka Seiken, Tokyo, Japan) and heat inactivated at 56°C for 30 min. Serial 2-fold dilutions starting at 1:10 were tested for the ability to inhibit the agglutination of 0.5% turkey RBCs with 4 HAU of each virus. HI titers were recorded as the inverse of the highest dilution that inhibited hemagglutination.
ELISA.
For enzyme-linked immunosorbent assay (ELISA), 96-well flat-bottom plates (Greiner Bio-One, Monroe, NC) were incubated overnight at 4°C with 16 HAU of each virus in 1× coating solution (Seracare, Milford, MA). Virus samples were run in duplicate. Plates were blocked with StartingBlock (PBS) blocking buffer (Thermo Fisher Scientific, Waltham, MA) for 1 h at room temperature. After blocking, plates were washed three times with PBS containing 0.05% Tween 20 (PBS-T). Twofold serial dilutions of swine antisera were added and allowed to incubate for 1 h at room temperature. After incubation, plates were washed three times with PBS-T before adding the secondary goat anti-swine IgG-HRP HRP-conjugated polyclonal antibody (Seracare, Milford, MA). Plates were incubated at room temperature for 1 h. After incubation, plates were washed three times with PBS-T before treatment with freshly prepared SS buffer or 2, 2'-azino-di(3-ethylbenzthiazoline-6-sulfonate) (ABTS) 1-component microwell peroxidase substrate (Seracare, Milford, MA) for 1 h at room temperature. The reaction was stopped by adding 3% (vol/vol in ddH2O) H2SO4 or ABTS peroxidase stop solution (Seracare, Milford, MA), respectively. Absorbance readings were obtained at 405 nm using a Victor ×3 multilabel plate reader (PerkinElmer, Waltham, MA).
Statistical analysis.
All statistical analyses were performed using the GraphPad Prism software, version 7 (GraphPad Software Inc., San Diego, CA). For multiple comparisons, one-way or two-way analysis of variance (ANOVA) followed by a post hoc Tukey test was performed. When indicated, a P value below 0.05 (P < 0.05) was considered significant.
Structure modeling.
A model of the structure of the HA of A/turkey/Ohio/313053/2004 (H3N2) was built by homology modeling using Modeller v9.16 (
57) based upon the crystal structures of multiple H3 HA proteins (Protein Data Bank [PDB] codes
2YP7,
1HA0,
2YP2,
4WE8, and
4WE5). The generated model was subsequently rendered with PyMOL v2.1 (
58).
Computational analysis of HA sequences.
The frequency distribution of amino acid identities at residue 145 (H3 numbering) was computed using the sequence variation analysis tool in the Influenza Research Database (IRD) (
59). HA amino acid sequences from IAVs of the H3 subtype isolated from swine, human, avian, canine, and equine hosts and publicly available in the IRD as of 6 September 2017 were analyzed. For humans, the frequency was calculated from precomputed data in the IRD database. Only amino acids that reached a frequency of at least 1% are labeled in the plot legend.