Evolutionary changes in influenza B are not primarily governed by antibody selection

Spatial and Temporal Coevolution of N2 Neuraminidase and H1 and H3 Hemagglutinin Genes of Influenza A Virus in United States Swine

The neuraminidase (NA) and hemagglutinin (HA) are essential surface glycoproteins of influenza A virus (IAV). In this study, the evolution of subtype N2 NA paired with H1 and H3 subtype HA in swine was evaluated to understand if genetic diversity of HA and NA were linked. Using time-scaled Bayesian phylodynamic analyses, the relationships of paired swine N2 with H1 or H3 from 2009 to 2018 were evaluated. These data demonstrated increased relative genetic diversity within the major N2 clades circulating in swine in the United States (N2.1998 between 2014-2017 and N2.2002 between 2010-2016). Preferential pairing was observed among specific NA and HA genetic clades. Gene reassortment between cocirculating influenza A strains resulted in novel pairings that persisted. The changes of genetic diversity in the NA gene were quantified using Bayesian phylodynamic analyses and increases in diversity were observed subsequent to novel NA-HA reassortment events. The rate of evolution among NA-N2 clades and HA-H1 and HA-H3 clades were similar. Bayesian phylodynamic analyses demonstrated strong spatial patterns in N2 genetic diversity, but frequent interstate movement of rare N2 clades provided opportunity for reassortment and emergence of new N2-HA pairings. The frequent regional movement of pigs and their influenza viruses is an explanation for the documented patterns of reassortment and subsequent changes in gene diversity. The reassortment and evolution of NA and linked HA evolution may result in antigenic drift of both major surface glycoproteins, reducing vaccine efficacy, with subsequent impact on animal health.

Altered NKp46 Recognition and Elimination of Influenza B Viruses

Every year, millions of people worldwide are infected with influenza, causing enormous health and economic problems. The most common type of influenza is influenza A. It is known that Natural Killer (NK) cells play an important role in controlling influenza A infection, mostly through the recognition of the viral protein hemagglutinin (HA) by the activating receptor, NKp46. In contrast, little is known regarding NK cell recognition of influenza B viruses, even though they are responsible for a third of all pediatric influenza deaths and are therefore included in the seasonal vaccine each year. Here we show that NKp46 also recognizes influenza B viruses. We show that NKp46 binds the HA protein of influenza B in a sialic acid-dependent manner, and identified the glycosylated residue in NKp46, which is critical for this interaction. We discovered that this interaction has a binding affinity approximately seven times lower than NKp46 binding of influenza A’s HA. Finally, we demonstrated, using mice deficient for the mouse orthologue of NKp46, named NCR1, that NKp46 is not important for influenza B elimination. These findings enable us to better understand the interactions between the different influenza viruses and NK cells that are known to be crucial for viral elimination.

Fitness cost of reassortment in human influenza

Reassortment, which is the exchange of genome sequence between viruses co-infecting a host cell, plays an important role in the evolution of segmented viruses. In the human influenza virus, reassortment happens most frequently between co-existing variants within the same lineage. This process breaks genetic linkage and fitness correlations between viral genome segments, but the resulting net effect on viral fitness has remained unclear. In this paper, we determine rate and average selective effect of reassortment processes in the human influenza lineage A/H3N2. For the surface proteins hemagglutinin and neuraminidase, reassortant variants with a mean distance of at least 3 nucleotides to their parent strains get established at a rate of about 10⁻² in units of the neutral point mutation rate. Our inference is based on a new method to map reassortment events from joint genealogies of multiple genome segments, which is tested by extensive simulations. We show that intra-lineage reassortment processes are, on average, under substantial negative selection that increases in strength with increasing sequence distance between the parent strains. The deleterious effects of reassortment manifest themselves in two ways: there are fewer reassortment events than expected from a null model of neutral reassortment, and reassortant strains have fewer descendants than their non-reassortant counterparts. Our results suggest that influenza evolves under ubiquitous epistasis across proteins, which produces fitness barriers against reassortment even between co-circulating strains within one lineage.

The genome of the human influenza virus consists of 8 disjoint RNA polymer segments. These segments can undergo reassortment: when two viruses co-infect a host cell, they can produce viral offspring with a new combination of segments. In this paper, we show that reassortment within a given influenza lineage induces a fitness cost that increases in strength with increasing genetic distance of the parent viruses. Our finding suggests that evolution continuously produces viral proteins whose fitness depends on each other; reassortment reduces fitness by breaking up successful combinations of proteins. Thus, selection across proteins constrains viral evolution within a given lineage, and it may be an important factor in defining a viral species.

Host proteostasis modulates influenza evolution

Predicting and constraining RNA virus evolution require understanding the molecular factors that define the mutational landscape accessible to these pathogens. RNA viruses typically have high mutation rates, resulting in frequent production of protein variants with compromised biophysical properties. Their evolution is necessarily constrained by the consequent challenge to protein folding and function. We hypothesized that host proteostasis mechanisms may be significant determinants of the fitness of viral protein variants, serving as a critical force shaping viral evolution. Here, we test that hypothesis by propagating influenza in host cells displaying chemically-controlled, divergent proteostasis environments. We find that both the nature of selection on the influenza genome and the accessibility of specific mutational trajectories are significantly impacted by host proteostasis. These findings provide new insights into features of host–pathogen interactions that shape viral evolution, and into the potential design of host proteostasis-targeted antiviral therapeutics that are refractory to resistance.

Influenza viruses, commonly called flu, can evade our immune system and develop resistance to treatments by changing frequently. Specifically, mutations in their genome cause influenza proteins to change in ways that can help the virus evade our defences. However, these mutations come at a cost and can prevent the viral proteins from forming functional and stable three-dimensional shapes – a process known as protein folding – thereby hampering the virus’ ability to replicate.

In human cells, proteins called chaperones can help our other proteins fold properly. Influenza viruses do not have their own chaperones and, instead, hijack those of their host. Host chaperones are therefore crucial to the virus’ ability to replicate. However, until now, it was not known if host chaperones can influence how these viruses evolve.

Here, Phillips et al. used mammalian cells to study how host chaperones affect an evolving influenza population. First, cells were engineered to either have normal chaperone levels, elevated chaperone levels, or inactive chaperones. Next, the H3N2 influenza strain was grown in these different conditions for nearly 200 generations and sequenced to determine how the virus evolved in each distinctive host chaperone environment.

Phillips et al. discovered that host chaperones affect the rate at which mutations accumulate in the influenza population, and also the types of mutations in the influenza genome. For instance, when a chaperone called Hsp90 was inactivated, mutations became prevalent in the viral population more slowly than in cells with normal or elevated chaperone levels. Moreover, some specific mutations fared better in cells with high chaperone levels, whilst others worked better in cells with inactivated chaperones.

These results suggest that influenza evolution is affected by host chaperone levels in complex and important ways. Moreover, whether chaperones will promote or hinder the effects of any single mutation is difficult to predict ahead of time. This discovery is significant, as the chaperones available to influenza can vary in different tissues, organisms and infectious conditions, and may therefore influence the virus' ability to change and evolve in a context-specific manner.

The findings are likely to extend to other viruses such as HIV and Ebola, which also hijack host chaperones for the same purpose. More work is now needed to systematically quantify these effects so that we can better predict how specific chaperones will affect the ability of viruses to adapt, especially in pathologically relevant conditions like fever or viral host-switching. In the future, such insights could help shape the design of treatments to which viruses do not evolve resistance.

Modes of transmission of influenza B virus in households

Introduction

While influenza A and B viruses can be transmitted via respiratory droplets, the importance of small droplet nuclei “aerosols” in transmission is controversial.

Methods and Findings

In Hong Kong and Bangkok, in 2008–11, subjects were recruited from outpatient clinics if they had recent onset of acute respiratory illness and none of their household contacts were ill. Following a positive rapid influenza diagnostic test result, subjects were randomly allocated to one of three household-based interventions: hand hygiene, hand hygiene plus face masks, and a control group. Index cases plus their household contacts were followed for 7–10 days to identify secondary infections by reverse transcription polymerase chain reaction (RT-PCR) testing of respiratory specimens. Index cases with RT-PCR-confirmed influenza B were included in the present analyses. We used a mathematical model to make inferences on the modes of transmission, facilitated by apparent differences in clinical presentation of secondary infections resulting from aerosol transmission. We estimated that approximately 37% and 26% of influenza B virus transmission was via the aerosol mode in households in Hong Kong and Bangkok, respectively. In the fitted model, influenza B virus infections were associated with a 56%–72% risk of fever plus cough if infected via aerosol route, and a 23%–31% risk of fever plus cough if infected via the other two modes of transmission.

Conclusions

Aerosol transmission may be an important mode of spread of influenza B virus. The point estimates of aerosol transmission were slightly lower for influenza B virus compared to previously published estimates for influenza A virus in both Hong Kong and Bangkok. Caution should be taken in interpreting these findings because of the multiple assumptions inherent in the model, including that there is limited biological evidence to date supporting a difference in the clinical features of influenza B virus infection by different modes.

Identification of conserved RNA secondary structures at influenza B and C splice sites reveals similarities and differences between influenza A, B, and C

BACKGROUND

Influenza B and C are single-stranded RNA viruses that cause yearly epidemics and infections. Knowledge of RNA secondary structure generated by influenza B and C will be helpful in further understanding the role of RNA structure in the progression of influenza infection.

FINDINGS

All available protein-coding sequences for influenza B and C were analyzed for regions with high potential for functional RNA secondary structure. On the basis of conserved RNA secondary structure with predicted high thermodynamic stability, putative structures were identified that contain splice sites in segment 8 of influenza B and segments 6 and 7 of influenza C. The sequence in segment 6 also contains three unused AUG start codon sites that are sequestered within a hairpin structure.

CONCLUSIONS

When added to previous studies on influenza A, the results suggest that influenza splicing may share common structural strategies for regulation of splicing. In particular, influenza 3' splice sites are predicted to form secondary structures that can switch conformation to regulate splicing. Thus, these RNA structures present attractive targets for therapeutics aimed at targeting one or the other conformation.

Hemagglutinin and neuraminidase genes of influenza B viruses circulating in Riyadh, Saudi Arabia during 2010-2011: evolution and sequence analysis

Influenza viruses are known as continuing threats to human public health every year worldwide. Evolutionary dynamics of influenza B viruses in humans are in a unique progression having two lineages; B/Yam and B/Vic-like viruses, which are circulating simultaneously worldwide. There is a considerable lack of data on influenza B viruses circulating in Saudi Arabia. During the winter-spring season of 2010-2011, 80 nasopharyngeal aspirates were collected from hospitalized patients with flu-like symptoms in Riyadh. Screening of samples by one-step RT-PCR identified three (3.8%) influenza B viruses. Sequencing of hemagglutinin (HA) and neuraminidase (NA) genes was performed to analyze influenza B viruses circulating in Riyadh as compared to the globally circulating strains. Several common and six unique amino acid substitutions were observed for both HA and NA genes of influenza B Saudi strains. Three unique substitutions (T182A, D196N, and K254R) were identified in HA gene of the B/Yam-like Riyadh strains. In NA gene, a unique common substitution (D53G) was found in all Riyadh strains, while two unique substitutions (L38P, G233R) were recognized only in B/Vic-like Riyadh strains. Riyadh strains were also found to contain N-glycosylation site in HA gene of both B/Vic and B/Yam lineages at positions 197-199 (NET) and 196-198 (NNK/DNK), respectively. The significance of these mutations on the antigenicity of both lineages is discussed herein. The unique changes observed in HA and NA genes of influenza B Riyadh strains support strongly the need for continuous surveillance and monitoring of new evolving strains that might pose threat to the Saudi community.

Phylogenetic and evolutionary history of influenza B viruses, which caused a large epidemic in 2011-2012, Taiwan

The annual recurrence of the influenza epidemic is considered to be primarily associated with immune escape due to changes to the virus. In 2011-2012, the influenza B epidemic in Taiwan was unusually large, and influenza B was predominant for a long time. To investigate the genetic dynamics of influenza B viruses during the 2011-2012 epidemic, we analyzed the sequences of 4,386 influenza B viruses collected in Taiwan from 2004 to 2012. The data provided detailed insight into the flux patterns of multiple genotypes. We found that a re-emergent TW08-I virus, which was the major genotype and had co-circulated with the two others, TW08-II and TW08-III, from 2007 to 2009 in Taiwan, successively overtook TW08-II in March and then underwent a lineage switch in July 2011. This lineage switch was followed by the large epidemic in Taiwan. The whole-genome compositions and phylogenetic relationships of the representative viruses of various genotypes were compared to determine the viral evolutionary histories. We demonstrated that the large influenza B epidemic of 2011-2012 was caused by Yamagata lineage TW08-I viruses that were derived from TW04-II viruses in 2004-2005 through genetic drifts without detectable reassortments. The TW08-I viruses isolated in both 2011-2012 and 2007-2009 were antigenically similar, indicating that an influenza B virus have persisted for 5 years in antigenic stasis before causing a large epidemic. The results suggest that in addition to the emergence of new variants with mutations or reassortments, other factors, including the interference of multi-types or lineages of influenza viruses and the accumulation of susceptible hosts, can also affect the scale and time of an influenza B epidemic.

Subtype- and antigenic site-specific differences in biophysical influences on evolution of influenza virus hemagglutinin

Background

Influenza virus undergoes rapid evolution by both antigenic shift and antigenic drift. Antibodies, particularly those binding near the receptor-binding site of hemagglutinin (HA) or the neuraminidase (NA) active site, are thought to be the primary defense against influenza infection, and mutations in antibody binding sites can reduce or eliminate antibody binding. The binding of antibodies to their cognate antigens is governed by such biophysical properties of the interacting surfaces as shape, non-polar and polar surface area, and charge.

Methods

To understand forces shaping evolution of influenza virus, we have examined HA sequences of human influenza A and B viruses, assigning each amino acid values reflecting total accessible surface area, non-polar and polar surface area, and net charge due to the side chain. Changes in each of these values between neighboring sequences were calculated for each residue and mapped onto the crystal structures.

Results

Areas of HA showing the highest frequency of pairwise changes agreed well with previously identified antigenic sites in H3 and H1 HAs, and allowed us to propose more detailed antigenic maps and novel antigenic sites for H1 and influenza B HA. Changes in biophysical properties differed between HAs of different subtypes, and between different antigenic sites of the same HA. For H1, statistically significant differences in several biophysical quantities compared to residues lying outside antigenic sites were seen for some antigenic sites but not others. Influenza B antigenic sites all show statistically significant differences in biophysical quantities for all antigenic sites, whereas no statistically significant differences in biophysical quantities were seen for any antigenic site is seen for H3. In many cases, residues previously shown to be under positive selection at the genetic level also undergo rapid change in biophysical properties.

Conclusions

The biophysical consequences of amino acid changes introduced by antigenic drift vary from subtype to subtype, and between different antigenic sites. This suggests that the significance of antibody binding in selecting new variants may also be variable for different antigenic sites and influenza subtypes.

Surveillance and molecular characterization of human influenza B viruses during 2006-2010 revealed co-circulation of Yamagata-like and Victoria-like strains in eastern India

Acute respiratory illness (ARI) is one of the major health problems in tropical countries of Asia, like India where approximately 0.5 million children in the age group of < 5 years die annually. Previously we have reported the genetic characterization of influenza A (Inf-A) strains circulating in Kolkata, eastern India. This study was initiated to characterize the genetic diversity of the circulating influenza B (Inf-B) viruses. Of 3035 nasal/throat swabs, 494 (16.3%) samples were identified as influenza A/B positive by real time RT-PCR, of which 244 samples were confirmed having Inf-B infection. Comparison of nucleotide (nt) and amino acid (aa) sequences of HA and NA gene of Inf-B viruses revealed co-circulation of B/Yamagata and B/Victoria lineages. Of the 32 randomly selected Inf-B strains from Kolkata, seventeen strains possessed reassorted NA gene. There was a single Histidine to Asparagine substitution in the 131st position which is a part of 120 loop on HA1 region along with a deletion at position 178 in the Kolkata strains belonging to the Yamagata lineage. Amino acid substitution was observed at position 198 on NA gene in the strains B/Kol/542/2006, B/Kol/1373/2008, B/Kol/1880/2008, B/Kol/2044/2008 and in all the representative strains isolated during 2009 with respect to the circulating vaccine strains. This substitution is responsible for reduced sensitivity of neuraminidase inhibitors. The results highlight the importance of monitoring Inf-B viruses for development of antiviral resistance among circulating strains.

Influenza control in the 21st century: Optimizing protection of older adults

Older adults (> or =65 years of age) are particularly vulnerable to influenza illness. This is due to a waning immune system that reduces their ability to respond to infection, which leads to more severe cases of disease. The majority ( approximately 90%) of influenza-related deaths occur in older adults and, in addition, catastrophic disability resulting from influenza-related hospitalization represents a significant burden in this vulnerable population. Current influenza vaccines provide benefits for older adults against influenza; however, vaccine effectiveness is lower than in younger adults. In addition, antigenic drift is also a concern, as it can impact on vaccine effectiveness due to a mismatch between the vaccine virus strain and the circulating virus strain. As such, vaccines that offer higher and broader protection against both homologous and heterologous virus strains are desirable. Approaches currently available in some countries to meet this medical need in older adults may include the use of adjuvanted vaccines. Future strategies under evaluation include the use of high-dose vaccines; novel or enhanced adjuvantation of current vaccines; use of live attenuated vaccines in combination with current vaccines; DNA vaccines; recombinant vaccines; as well as the use of different modes of delivery and alternative antigens. However, to truly evaluate the benefits that these solutions offer, further efficacy and effectiveness studies, and better correlates of protection, including a precise measurement of the T cell responses that are markers for protection, are needed. While it is clear that vaccines with greater immunogenicity are required for older adults, and that adjuvanted vaccines may offer a short-term solution, further research is required to exploit the many other new technologies.

Diversifying selective pressure on influenza B virus hemagglutinin

Influenza B virus hemagglutinin (HA) is a major surface glycoprotein with frequent amino acid substitutions. However, the roles of antibody selection in the amino acid substitutions of HA were still poorly understood. In order to gain insights into this important issue, an analysis was conducted on a total of 271 HA1 sequences of influenza B virus strains isolated during 1940-2007. In this analysis, phylogenetic analysis by maximum likelihood (PAML) package was used to detect the existence of positive selection and to identify positively selected sites on HA1. Strikingly, all the positively selected sites were located in the four major epitopes (120-loop, 150-loop, 160-loop, and 190-helix) of HA identified in previous studies, thus supporting a predominant role of antibody selection in HA evolution. Of particular significance is the involvement of the 120-loop in positive selection, which may become increasingly important in future field isolates. Despite the absence of different subtypes, influenza B virus HA continued to evolve into new sublineages, within which the four major epitopes were targeted selectively in positive selection. Thus, any newly emerging strains need to be placed in the context of their evolutionary history in order to understand and predict their epidemic potential.

Heterogeneous selective pressure acting on influenza B Victoria- and Yamagata-like hemagglutinins

As a consequence of immune pressure, influenza virus hemagglutinin presents some of its amino acids under positive selection. Several authors have reported the existence of influenza A hemagglutinin codons under positive selective pressure (PSP). In this framework, the present work objectives were to demonstrate the presence of PSP and evaluate its effects on Victoria- and Yamagata-like influenza B viruses. Methodology adopted consisted in estimating the acceptance rate of nonsynonymous substitutions (omega = dN/dS) that describe the strength of selective pressure and identifying codons that may be positively selected, applying a set of continuous-time Markov chain codon-substitution models. Two groups of HA1 sequences (140 from Yamagata and 60 from Victoria lineage) were used. All the model maximum-likelihood estimates were obtained using codeml software application (PAML 3.15). The hypothesis of no existence of sites under PSP was rejected for both lineages (p < 0.001), using likelihood ratio tests. These results demonstrate the presence of positive selection acting on hemagglutinin of both Yamagata- and Victoria-like influenza B viruses. Several different sites were identified to be under PSP on Yamagata and Victoria hemagglutinins. Sites found with a posterior probability > 0.95 were codons 197 and 199 in both lineages, codon 75 in the Yamagata lineage, and codon 129 in the Victoria lineage. The detected amino acids are located at or near antigenic sites in influenza A virus H3 hemagglutinin.

Live, attenuated influenza virus (LAIV) vehicles are strong inducers of immunity toward influenza B virus

Historically, vaccines developed toward influenza viruses of the B type using methodologies developed for influenza A viruses as a blueprint have not been equally efficacious or effective. Because most influenza research and public attention concerns influenza A viruses, these shortcomings have not been adequately addressed. In this manuscript, we utilized different influenza vaccine vehicles to compare immunogenicity and protection in mice and ferrets after vaccination against an influenza B virus. We report that plasmid DNA vaccines demonstrate low immunogenicity profiles and poor protection compared to either whole, inactivated influenza virus (IIV) or, live, attenuated influenza virus (LAIV) vaccines. When mixed prime:boost regimens using LAIV and IIV were studied, we observed a boosting effect in mice after priming with LAIV that was not seen when IIV was used as the prime. In ferrets LAIV induced high antibody titers after a single dose and provided a boost in IIV-primed animals. Regimens including LAIV as a prime demonstrated enhanced protection, and adjuvantation was required for efficacy using the IIV preparation. Our results differ from generally accepted influenza A virus vaccine models, and argue that strategies for control of influenza B virus should be considered separately from those for influenza A virus.

The evolutionary dynamics of human influenza B virus

Despite their close phylogenetic relationship, type A and B influenza viruses exhibit major epidemiological differences in humans, with the latter both less common and less often associated with severe disease. However, it is unclear what processes determine the evolutionary dynamics of influenza B virus, and how influenza viruses A and B interact at the evolutionary scale. To address these questions we inferred the phylogenetic history of human influenza B virus using complete genome sequences for which the date (day) of isolation was available. By comparing the phylogenetic patterns of all eight viral segments we determined the occurrence of segment reassortment over a 30-year sampling period. An analysis of rates of nucleotide substitution and selection pressures revealed sporadic occurrences of adaptive evolution, most notably in the viral hemagglutinin and compatible with the action of antigenic drift, yet lower rates of overall and nonsynonymous nucleotide substitution compared to influenza A virus. Overall, these results led us to propose a model in which evolutionary changes within and between the antigenically distinct 'Yam88' and 'Vic87' lineages of influenza B virus are the result of changes in herd immunity, with reassortment continuously generating novel genetic variation. Additionally, we suggest that the interaction with influenza A virus may be central in shaping the evolutionary dynamics of influenza B virus, facilitating the shift of dominance between the Vic87 and the Yam88 lineages.

Crystal structure of unliganded influenza B virus hemagglutinin

Here we report the crystal structure of hemagglutinin (HA) from influenza B/Hong Kong/8/73 (B/HK) virus determined to 2.8 A. At a sequence identity of approximately 25% to influenza A virus HAs, B/HK HA shares a similar overall structure and domain organization. More than two dozen amino acid substitutions on influenza B virus HAs have been identified to cause antigenicity alteration in site-specific mutants, monoclonal antibody escape mutants, or field isolates. Mapping these substitutions on the structure of B/HK HA reveals four major epitopes, the 120 loop, the 150 loop, the 160 loop, and the 190 helix, that are located close in space to form a large, continuous antigenic site. Moreover, a systematic comparison of known HA structures across the entire influenza virus family reveals evolutionarily conserved ionizable residues at all regions along the chain and subunit interfaces. These ionizable residues are likely the structural basis for the pH dependence and sensitivity to ionic strength of influenza HA and hemagglutinin-esterase fusion proteins.

Genetic analysis and evaluation of the reassortment of influenza B viruses isolated in Taiwan during the 2004-2005 and 2006-2007 epidemics

Influenza B viruses were predominant in Taiwan during the 2004-2005 epidemic and both Victoria and Yamagata lineage viruses co-circulated. A reassortant influenza B virus that contained a Victoria lineage hemagglutinin (HA) gene and Yamagata lineage neuraminidase (NA) gene appeared first in 2002 and became predominant during the 2004-2005 epidemic. During the 2006-2007 epidemic, an influenza B outbreak occurred in Taiwan and only Victoria lineage viruses circulated. We characterized the viruses isolated in the 2006-2007 epidemic and found that the HA genes of influenza B viruses from that epidemic were highly similar to those from the 2004-2005 epidemic. We also analyzed the NA genes of isolates from the 2006-2007 epidemic and found that they all belonged to the Yamagata lineage and formed a new genetic subclade. Comparison of isolates from the 2004-2005 and 2006-2007 epidemics revealed four substitutions, N220K, E320D, K343R and E404K in NA genes. Although the HA sequences from the 2006-2007 epidemic were similar to those from the 2004-2005 epidemic, the NA sequences differed, suggesting distinct patterns of evolution of the HA and NA genes from 2004-2007 in Taiwan. This study emphasizes that the evolution of the NA genes may contribute to reemergence of influenza B viruses.

Structural basis for receptor specificity of influenza B virus hemagglutinin

Receptor-binding specificity of HA, the major surface glycoprotein of influenza virus, primarily determines the host ranges that the virus can infect. Influenza type B virus almost exclusively infects humans and contributes to the annual "flu" sickness. Here we report the structures of influenza B virus HA in complex with human and avian receptor analogs, respectively. These structures provide a structural basis for the different receptor-binding properties of influenza A and B virus HA molecules and for the ability of influenza B virus HA to distinguish human and avian receptors. The structure of influenza B virus HA with avian receptor analog also reveals how mutations in the region of residues 194 to 196, which are frequently observed in egg-adapted and naturally occurring variants, directly affect the receptor binding of the resultant virus strains. Furthermore, these structures of influenza B virus HA are compared with known structures of influenza A virus HAs, which suggests the role of the residue at 222 as a key and likely a universal determinant for the different binding modes of human receptor analogs by different HA molecules.

Emergence of influenza A virus variants after prolonged shedding from pheasants

We previously demonstrated the susceptibility of pheasants to infection with influenza A viruses of 15 hemagglutinin (HA) subtypes: 13/23 viruses tested were isolated for >or=14 days, all in the presence of serum-neutralizing antibodies; one virus (H10) was shed for 45 days postinfection. Here we confirmed that 20% of pheasants shed low-pathogenic influenza viruses for prolonged periods. We aimed to determine why the antibody response did not clear the virus in the usual 3 to 10 days, because pheasants serve as a long-term source of influenza viruses in poultry markets. We found evidence of virus replication and histological changes in the large intestine, bursa of Fabricius, and cecal tonsil. The virus isolated 41 days postinfection was antigenically distinct from the parental H10 virus, with corresponding changes in the HA and neuraminidase. Ten amino acid differences were found between the parental H10 and the pheasant H10 virus; four were in potential antigenic sites of the HA molecule. Prolonged shedding of virus by pheasants results from a complex interplay between the diversity of virus variants and the host response. It is often argued that vaccination pressure is a mechanism that contributes to the generation of antigenic-drift variants in poultry. This study provided evidence that drift variants can occur naturally in pheasants after prolonged shedding of virus, thus strengthening our argument for the removal of pheasants from live-bird retail markets.

The evolution of epidemic influenza

Recent developments in complete-genome sequencing, antigenic mapping and epidemiological modelling are greatly improving our knowledge of the evolution of human influenza virus at the epidemiological scale. In particular, recent studies have revealed a more complex relationship between antigenic evolution, natural selection and reassortment than previously realized. Despite these advances, there is much that remains to be understood about the epidemiology of influenza virus, particularly the processes that determine the virus's strong seasonality. We argue that a complete understanding of the evolutionary biology of this important human pathogen will require a genomic view of genetic diversity, including the acquisition of polymorphism data from within individual hosts and from geographical regions, particularly the tropics, which have been poorly surveyed to date.

Comparison of the mutation rates of human influenza A and B viruses

Human influenza A viruses evolve more rapidly than influenza B viruses. To clarify the cause of this difference, we have evaluated the mutation rate of the nonstructural gene as revealed by the genetic diversity observed during the growth of individual plaques in MDCK cells. Six plaques were studied, representing two strains each of type A and B viruses. A total of 813,663 nucleotides were sequenced, giving rates of 2.0 x 10(-6) and 0.6 x 10(-6) mutations per site per infectious cycle, which, when extended to 1 year, agree well with the published annual evolutionary rates.

Molecular characterization of the HA gene of influenza type B viruses

Nucleotide sequences of the HA1 subunit of influenza B viruses isolated in Portugal between 1994 and 2003 influenza winter seasons were analyzed by the Neighbor-Joining algorithm and rates of HA1 evolution estimated by linear regression. From 1994 to 2002, all influenza B viruses studied were of the Yamagata lineage. Strains isolated from 1994 to 1996, 1996 to 1999, and 1999 to 2002 revealed a high similarity with B/Beijing/184/93, B/Yamanashi/166/98, and B/Sichuan/379/99, respectively, and strains isolated during 1994-1995, 1996-1997, and 1998-1999 clustered in more than one branch of the phylogenetic tree. Victoria-related strains reappeared during 2002/2003 and formed only one branch in the phylogenetic tree revealing a closer relationship to B/Shandong/7/97. Evolutionary rates for strains from the Yamagata lineage were estimated as 3.82x10(-3) nucleotides/site/year and 2.62x10(-3) nucleotides/site/year for Victoria-related strains. In order to identify putative influenza B HA1 codons under selective pressure, a codon-substitution model for heterogeneous selective pressure at amino acid sites was used. A percentage of 97.3% of codons under negative selective pressure and 2.7% of codons under positive selective pressure (omega=dN/dS=2.65) were estimated, with posterior probability higher than 0.90. Amino acid sites 75, 197, and 199 were found more likely to be under positive selective pressure.

A single amino acid change in the C-terminal domain of the matrix protein M1 of influenza B virus confers mouse adaptation and virulence

Serial passage of an initially avirulent influenza B virus, B/Memphis/12/97, resulted in the selection of a variant which was lethal in mice. Virulence correlated with improved growth in vivo and prolonged replication. Sequencing of the complete coding regions of the parent and mouse-adapted viruses revealed 8 amino acid differences. Sequencing and characterization of intermediate passages suggested that one change in the C-terminal domain of the M1 protein, an asparagine to a serine at position 221, was responsible for acquisition of virulence and lethality. Site-directed mutagenesis of the M segment of a different virus, B/Yamanashi/166/98, to change this amino acid residue confirmed its importance by conferring improved growth and virulence in mice. This observation suggests a role for the C domain of the M1 protein in growth and virulence in a mammalian host.

Genetic diversity of influenza B virus: the frequent reassortment and cocirculation of the genetically distinct reassortant viruses in a community

To characterize the genetic diversity of influenza B viruses isolated during one influenza season, the antigenic and genetic relationships among 20 strains of influenza B virus isolated in February and March 2001 at one pediatric clinic in Yamagata City, Japan, were investigated. The HA gene and seven other gene segments were phylogenetically divided into three distinct sublineages (Harbin/7/94-, Tokyo/6/98-, and Shiga/T30/98-related lineage) of the Yamagata/16/88-like lineage. The NS genes of the viruses belonging to the Harbin/7/94-related lineage have additional three nucleotides at positions 439-447, and were phylogenetically distinguishable from those of the currently circulating Yamagata/16/88- and Victoria/2/87-like lineages, but were closely related to that of the Yamagata/16/88-like lineage isolated before 1994. Moreover, four strains of influenza B virus isolated in the same community between 2002 and 2003 were further examined. Phylogenetic analysis revealed that a virus of Victoria/2/87-like lineage isolated in 2003 had acquired the NA, NS, M, and PA gene segments from a Shiga/T30/98-like virus, and two strains of Harbin/7/94-related lineage had acquired the various gene segments from Shiga/T30/98-like virus through a reassortment event. These results indicate that genetically distinct multiple viruses can combine to cause an influenza B epidemic in a community and that the frequent reassortment among these viruses plays a role in generating the genetic diversity of influenza B viruses.

Multiple genotypes of influenza B virus circulated between 1979 and 2003

The segmented genome of influenza B virus allows exchange of gene segments between cocirculating strains. Through this process of reassortment, diversity is generated by the mixing of genes between viruses that differ in one or more gene segments. Phylogenetic and evolutionary analyses of all 11 genes of 31 influenza B viruses isolated from 1979 to 2003 were used to study the evolution of whole genomes. All 11 genes diverged into two new lineages prior to 1987. All genes except the NS1 gene were undergoing linear evolution, although the rate of evolution and the degree to which nucleotide changes translated into amino acid changes varied between lineages and by gene. Frequent reassortment generated 14 different genotypes distinct from the gene constellation of viruses circulating prior to 1979. Multiple genotypes cocirculated in some locations, and a sequence of reassortment events over time could not be established. The surprising diversity of the viruses, unrestricted mixing of lineages, and lack of evidence for coevolution of gene segments do not support the hypothesis that the reassortment process is driven by selection for functional differences.

Evolution of surface and nonstructural-1 genes of influenza B viruses isolated in the Province of Québec, Canada, during the 1998-2001 period

After 2 minor winter seasons, influenza B viruses were predominantly isolated in the Province of Quebec, Canada, during the 2000-2001 season representing 74% of laboratory-confirmed influenza viruses. We performed an antigenic study of the hemagglutinin (HA) protein and a molecular characterization of the HA1 region, nonstructural-1 (NS1) and neuraminidase (NA)/NB genes of 20 influenza B strains isolated in the Province of Quebec during the 1998-2001 period. Our isolates were compared to recent vaccine strains (B/Harbin/7/94 in 1998-1999, B/Yamanashi/166/98 in 1999-2000 and 2000-2001, and B/Sichuan/379/99 in 2001-2002). The hemagglutination inhibition (HI) test revealed that all isolates were different from B/Harbin/7/94 and were more related to the 2 other vaccine strains although precise identification was often impossible. Molecular analysis of the HA1 gene revealed that both B/Yamanashi/166/98-like and B/Sichuan/379/99-like isolates co-circulated during the 1998-1999 season whereas isolates from the 2 subsequent years were more related to B/Sichuan/379/99. Most isolates (8/9) of the 2000-2001 season contained a N126D substitution recently associated with altered antigenicity in recent influenza B/Yamagata/16/88-related viruses. Although the HA1 and NS1 protein sequences of viruses isolated during the 1998-1999 season were clearly different from those of the respective vaccine strain (B/Harbin/7/94), the NA protein sequence of those isolates was slightly more related to B/Harbin/7/94 than B/Yamanashi/166/98 suggesting distinct patterns of evolution for these genes. This study confirms the importance of a detailed molecular analysis for understanding the evolution of influenza B viruses.

Molecular characterization of influenza B viruses circulating in northern Italy during the 2001-2002 epidemic season

During the 2001-2002 influenza season, virological surveillance highlighted the predominant circulation of B viruses (86% of isolates) in Italy, in contrast to many other countries in Europe and North America where AH3N2 viruses were isolated most frequently, and in contrast to the infrequent isolation of B viruses in Italy during the previous two years. Associated with this predominance of influenza B was the re-emergence of B/Victoria/2/87-lineage viruses, closely related to B viruses prevalent during the 1980s, which are distinct antigenically and genetically from circulating B/Sichuan/379/99-like viruses of the B/Yamagata/16/88 lineage, which predominated in most parts of the world during the last 10 years. Ninety-four viruses isolated in two regions of northern Italy were characterized, 50 by direct sequencing of haemagglutinin (HA). Viruses of both Victoria and Yamagata lineages co-circulated throughout the 12 weeks of the influenza season. The HAs of the Yamagata-lineage viruses were heterogeneous and comprised two sublineages, represented by B/Sichuan/379/99 and B/Harbin/7/94, whereas the Victoria-lineage viruses were more homogeneous and closely related to B/Hong Kong/330/01, the current prototype vaccine strain. The antigenic and genetic characteristics of the viruses correlated with certain epidemiological features. In particular, the low age (<14 years) of individuals infected with B/Hong Kong/330/01-like viruses is likely to reflect the greater susceptibility of the youngest cohort, due to lack of previous exposure to Victoria-lineage viruses, and is consistent with the conclusion that vaccination with a B/Sichuan/379/99-like virus would give poor protection against infection with B/Hong Kong/330/01-like (Victoria-lineage) viruses.

Origin and evolution of influenza virus hemagglutinin genes

Influenza A, B, and C viruses are the etiological agents of influenza. Hemagglutinin (HA) is the major envelope glycoprotein of influenza A and B viruses, and hemagglutinin-esterase (HE) in influenza C viruses is a protein homologous to HA. Because influenza A virus pandemics in humans appear to occur when new subtypes of HA genes are introduced from aquatic birds that are known to be the natural reservoir of the viruses, an understanding of the origin and evolution of HA genes is of particular importance. We therefore conducted a phylogenetic analysis of HA and HE genes and showed that the influenza A and B virus HA genes diverged much earlier than the divergence between different subtypes of influenza A virus HA genes. The rate of amino acid substitution for A virus HAs from duck, a natural reservoir, was estimated to be 3.19 x 10(-4) per site per year, which was slower than that for human and swine A virus HAs but similar to that for influenza B and C virus HAs (HEs). Using this substitution rate from the duck, we estimated that the divergences between different subtypes of A virus HA genes occurred from several thousand to several hundred years ago. In particular, the earliest divergence time was estimated to be about 2,000 years ago. Also, the A virus HA gene diverged from the B virus HA gene about 4,000 years ago and from the C virus HE gene about 8,000 years ago. These time estimates are much earlier than the previous ones.

Immunity to influenza in the elderly

Influenza is caused by a constantly varying segmented RNA virus that necessitates yearly review of vaccine composition. Humans over the age of 65 years are considered at high risk from influenza; during influenza epidemics the rate of hospitalization in the elderly is very high and up to 90% mortality can occur. Vaccination of the elderly has been shown to be efficacious and cost effective but immunological senescence in the institutionally confined frail elderly is demonstrated by failure to induce herd immunity after vaccination. Reductions in B- and T-cell immunity and in levels of interleukin-2 are age related. Attempts to increase the immunoresponsiveness of the elderly to influenza vaccines have given mixed results. The most convincing evidence is in rodents where dietary caloric restriction has been shown to enhance viral immunity.

Genetic characterization of the hemagglutinin of two strains of influenza B virus co-circulated in Taiwan

Two isolates of influenza B virus were obtained in the spring of 1997. One strain, B/Taiwan/21706/97, was isolated from a patient who had acute tonsillitis. The other, B/Taiwan/3143/97, was isolated from a patient who was diagnosed with meningoencephalitis. This implies that the influenza B viruses not only cause respiratory symptoms but may also cause inflammation of the nervous system. Sequence analysis of the hemagglutinin (HA) gene, HA1 domain, indicated that there were remarkable amino acid changes in the strain B/Taiwan/3143/97 compared to B/Victoria/2/87, B/Yamagata/16/88, and B/Taiwan/7/88. The changes in the positions 116, 200, 238, 242, and 271 were correlated with receptor binding. Furthermore, a potential glycosylation site at position 233 was lost. In total, 30 amino acid changes were noted at positions ranging from 116 to 295. These changes may affect the antigenicity of the virus. Phylogenetic analyses also showed that the B/Taiwan/3143/97 was located in an independent lineage, when compared to the reference strains belonging to B/Victoria/2/87 and B/Yamagata/16/88 lineages. This supports the hypothesis that influenza B viruses with distinct genetic characteristic were co-circulated in Taiwan.

Reassortment and insertion-deletion are strategies for the evolution of influenza B viruses in nature

The evolution of influenza B viruses is poorly understood. Reassortment of influenza B viruses in nature as a means of genetic variation has not been considered to be a major contributor to their evolution. However, the current practice of assigning evolutionary relationships by antigenic analysis of the hemagglutinin of influenza B viruses would fail to detect reassortants. In this study, influenza B viruses isolated within the past 10 years from sites in the United States and China were studied by nucleotide sequencing of the hemagglutinin and neuraminidase genes and construction of phylogenetic trees to assess evolutionary relationships. A group of viruses represented by B/Houston/1/92 possess a hemagglutinin derived from a B/Yamagata/16/88-like strain and a neuraminidase derived from a B/Victoria/2/87-like strain. A second reassortment event between the hemagglutinin of a B/Yamagata/16/88-like virus closely related to the B/Beijing/184/93 strain and the neuraminidase of a B/Victoria/2/87-like strain is represented by a single virus, B/Memphis/3/93. The neuraminidase of the reassortant viruses is most closely related to that of B/Victoria/2/87-like viruses currently circulating in Nanchang, China. A pattern of insertions and deletions in the hemagglutinin and the neuraminidase of different strains of influenza B viruses is observed. Reassortment plays a role in the evolution of influenza B viruses and may necessitate a change in the methods used to assess and identify new influenza viruses.

Genetic variability among group A and group B respiratory syncytial viruses in a children's hospital

Respiratory syncytial (RS) viruses isolated over three epidemic periods in a children's hospital in the United States were analyzed. The viruses (n = 174) were characterized as to major antigenic group (group A or B) by a PCR-based assay. Group A RS viruses were dominant the first 2 years, followed by a year with group B dominance (ratios of group A to group B viruses for epidemic periods, 56/4 for 1993-1994, 42/3 for 1994-1995, and 19/50 for 1995-1996). Genetic variability within the groups was assessed by restriction fragment analysis of PCR products; 79 isolates were also analyzed by nucleotide sequence determination of a variable region of the glycoprotein G gene. Among the group A RS virus isolates, this G-protein variable region had amino acid differences of as great as 38%. The G-protein amino acids of the group A viruses differed by up to 31% from the G-protein amino acids of a prototype (A2) group A virus. Among the group B RS virus G proteins, amino acid differences were as great as 14%. The G-protein amino acids of the group B viruses differed by up to 27% from the G-protein amino acids of a prototype (18537) group B virus. The group A and group B RS viruses demonstrated genetic variability between years and within individual years. Phylogenetic analysis revealed that there were multiple evolutionary lineages among both the group A and group B viruses. Among the recent group B isolates, variability was less than that seen for the group A viruses. However, comparisons to prototype strains revealed that the group B RS viruses may vary more extensively than was observed over the 3 years studied in the present investigation.

Serological and genetic characterisation of bovine respiratory syncytial virus (BRSV) indicates that Danish isolates belong to the intermediate subgroup: no evidence of a selective effect on the variability of G protein nucleotide sequence by prior cell culture adaption and passages in cell culture or calves

Danish isolates of bovine respiratory syncytial virus (BRSV) were characterised by nucleotide sequencing of the G glycoprotein and by their reactivity with a panel of monoclonal antibodies (MAbs). Among the six Danish isolates, the overall sequence divergence ranged between 0 and 3% at the nucleotide level and between 0 and 5% at the amino acid level. Sequence divergences of 7-8%, 8-9% and 2-3% (nucleotide) and 9-11%, 12-16% and 4-6% (amino acid) were obtained in the comparison made between the group of Danish isolates and the previously sequenced 391-2USA, 127UK and 220-69Bel isolates, respectively. Phylogenetic analysis showed that the Danish isolates formed three lineages within a separate branch of the phylogenetic tree. Nevertheless, the Danish isolates were closely related to the 220-69Bel isolate, the prototype of the intermediate antigenic subgroup. The sequencing of the extracellular part of the G gene of additional 11 field BRSV viruses, processed directly from lung samples without prior adaption to cell culture growth, revealed sequence variabilities in the range obtained with the propagated virus. In addition, several passages in cell culture and in calves had no major impact on the nucleotide sequence of the G protein. These findings indicated that the previously established variabilities of the G protein of RS virus isolates were not attributable to mutations induced during the propagation of the virus. The reactivity of the Danish isolates with G protein-specific MAbs were similar to that of the 220-69Bel isolate. Furthermore, the sequence of the immunodominant region was completely conserved among the Danish isolates on one side and the 220-69Bel isolate on the other. When combined, these data strongly suggested that the Danish isolates belong to the intermediate subgroup.

Evolution of the hemagglutinin gene of influenza B virus was driven by both positive and negative selection pressures

Nucleotide sequences of the hemagglutinin (HA) gene of twelve influenza B strains and their deduced amino acid sequences were extracted from the GenBank and compared to those of early isolates. Separate analyses of nonsynonymous and synonymous substitutions for the HA1 and HA2 regions individually indicate that the percentage of nonsynonymous substitutions in the HA1 ranges from 44.4-56.0% but in the HA2 from 0.0-6.5%. The results suggest that positive selection as well as negative selection played a role in the evolution of the influenza B HA gene with the former acting on the HA1 and the latter on the HA2.

The three subunits of the polymerase and the nucleoprotein of influenza B virus are the minimum set of viral proteins required for expression of a model RNA template

The genes encoding the nucleoprotein, PB1, PB2, and PA proteins of the influenza virus strain B/Panamá/45/90 have been cloned under control of the T7 RNA polymerase promoter of plasmid pGEM-3. Transfection of the recombinant plasmids obtained into mammalian cells, which had been infected with a vaccinia virus encoding the T7 RNA polymerase, resulted in expression of the expected influenza B virus polypeptides. Moreover, it is shown that coexpression of the four recombinant core proteins in COS-1 cells reconstituted a functional polymerase capable of expressing a synthetic influenza B virus-like CAT RNA. By using the influenza B virus recombinant plasmids and a set of pGEM-derived plasmids encoding the homologous core proteins of the influenza A virus A/Victoria/3/75 (I. Mena et al. (1994). J. Gen. Virol. 75, 2109-2114), the capabilities of homo- and heterotypic mixtures of the four core proteins to express synthetic type A and B CAT RNAs were analyzed. Both the influenza A and B virus polymerases were active in expressing, albeit with reduced efficiencies, the heterotypic model CAT RNAs. However, none of all possible heterotypic mixtures of the core proteins reconstituted a functional polymerase. In order to fully characterize the recombinant plasmids obtained, the nucleotide sequences of the cloned genes were determined and compared to sequences of other type B virus isolates. The results obtained from these latter analyses are discussed in terms of the conservation and evolution of the influenza B virus core genes.

Influenza B viruses with site-specific mutations introduced into the HA gene

We have succeeded in engineering changes into the genome of influenza B virus. First, model RNAs containing the chloramphenicol acetyltransferase gene flanked by the noncoding sequences of the HA or NS genes of influenza B virus were transfected into cells which were previously infected with an influenza B helper virus. Like those of the influenza A viruses, the termini of influenza B virus genes contain cis-acting signals which are sufficient to direct replication, expression, and packaging of the RNA. Next, a full-length copy of the HA gene from influenza B/Maryland/59 virus was cloned. Following transfection of this RNA, we rescued transfectant influenza B viruses which contain a point mutation introduced into the original cDNA. A series of mutants which bear deletions or changes in the 5' noncoding region of the influenza B/Maryland/59 virus HA gene were constructed. We were able to rescue viruses which contained deletions of 10 or 33 nucleotides at the 5' noncoding region of the HA gene. The viability of these viruses implies that this region of the genome is flexible in sequence and length.

Evolution and ecology of influenza A viruses

In this review we examine the hypothesis that aquatic birds are the primordial source of all influenza viruses in other species and study the ecological features that permit the perpetuation of influenza viruses in aquatic avian species. Phylogenetic analysis of the nucleotide sequence of influenza A virus RNA segments coding for the spike proteins (HA, NA, and M2) and the internal proteins (PB2, PB1, PA, NP, M, and NS) from a wide range of hosts, geographical regions, and influenza A virus subtypes support the following conclusions. (i) Two partly overlapping reservoirs of influenza A viruses exist in migrating waterfowl and shorebirds throughout the world. These species harbor influenza viruses of all the known HA and NA subtypes. (ii) Influenza viruses have evolved into a number of host-specific lineages that are exemplified by the NP gene and include equine Prague/56, recent equine strains, classical swine and human strains, H13 gull strains, and all other avian strains. Other genes show similar patterns, but with extensive evidence of genetic reassortment. Geographical as well as host-specific lineages are evident. (iii) All of the influenza A viruses of mammalian sources originated from the avian gene pool, and it is possible that influenza B viruses also arose from the same source. (iv) The different virus lineages are predominantly host specific, but there are periodic exchanges of influenza virus genes or whole viruses between species, giving rise to pandemics of disease in humans, lower animals, and birds. (v) The influenza viruses currently circulating in humans and pigs in North America originated by transmission of all genes from the avian reservoir prior to the 1918 Spanish influenza pandemic; some of the genes have subsequently been replaced by others from the influenza gene pool in birds. (vi) The influenza virus gene pool in aquatic birds of the world is probably perpetuated by low-level transmission within that species throughout the year. (vii) There is evidence that most new human pandemic strains and variants have originated in southern China. (viii) There is speculation that pigs may serve as the intermediate host in genetic exchange between influenza viruses in avian and humans, but experimental evidence is lacking. (ix) Once the ecological properties of influenza viruses are understood, it may be possible to interdict the introduction of new influenza viruses into humans.

Evolutionary pattern of the H 3 haemagglutinin of equine influenza viruses: multiple evolutionary lineages and frozen replication

The nucleotide and deduced amino acid sequences of the haemagglutinin genes coding for the HA 1 domain of H3N8 equine influenza viruses isolated over wide regions of the world were analyzed in detail to determine their evolutionary relationships. We have constructed a phylogenetic model tree by the neighbour-joining method using nucleotide sequences of 15 haemagglutinin genes, including those of five viruses determined in the present study. This gene tree revealed the existence of two major evolutionary pathways during a twenty five-year period between 1963 to 1988, and each pathway appeared to consist of two distinct lineages of haemagglutinin genes. Furthermore, our analysis of nucleotide sequences showed that two distinct lineages of equine H3N8 viruses were involved in an equine influenza outbreak during the period of December 1971-January 1972 in Japan. The number of nucleotide changes between strains was proportional to the length of time (in years) between their isolation except for three of the HA genes. However, there are three exceptional strains isolated in 1971, 1987, and 1988, respectively. The haemagglutinin gene in these strains showed a small number of nucleotide substitutions after they branched off around 1963, suggesting an example of frozen replication. Although the estimated rate (0.0094/site/year) of synonymous (silent) substitutions of the haemagglutinin gene of equine H3N8 viruses was nearly the same as that of human H 1 and H 3 haemagglutinin genes, the rate of nonsynonymous (amino-acid changing) substitutions of the former equine virus gene was estimated to be 0.00041/site/year--that is about 5 times lower than that estimated for the human H 3 haemagglutinin gene. The present study is the first demonstration that multiple evolutionary lineages of equine H3N8 influenza virus circulated since 1963.

Evolution of the H3 influenza virus hemagglutinin from human and nonhuman hosts

The nucleotide and amino acid sequences of 40 influenza virus hemagglutinin genes of the H3 serotype from mammalian and avian species and 9 genes of the H4 serotype were compared, and their evolutionary relationships were evaluated. From these relationships, the differences in the mutational characteristics of the viral hemagglutinin in different hosts were examined and the RNA sequence changes that occurred during the generation of the progenitor of the 1968 human pandemic strain were examined. Three major lineages were defined: one containing only equine virus isolates; one containing only avian virus isolates; and one containing avian, swine, and human virus isolates. The human pandemic strain of 1968 was derived from an avian virus most similar to those isolated from ducks in Asia, and the transfer of this virus to humans probably occurred in 1965. Since then, the human viruses have diverged from this progenitor, with the accumulation of approximately 7.9 nucleotide and 3.4 amino acid substitutions per year. Reconstruction of the sequence of the hypothetical ancestral strain at the avian-human transition indicated that only 6 amino acids in the mature hemagglutinin molecule were changed during the transition between an avian virus strain and a human pandemic strain. All of these changes are located in regions of the molecule known to affect receptor binding and antigenicity. Unlike the human H3 influenza virus strains, the equine virus isolates have no close relatives in other species and appear to have diverged from the avian viruses much earlier than did the human virus strains. Mutations were estimated to have accumulated in the equine virus lineage at approximately 3.1 nucleotides and 0.8 amino acids per year. Four swine virus isolates in the analysis each appeared to have been introduced into pigs independently, with two derived from human viruses and two from avian viruses. A comparison of the coding and noncoding mutations in the mammalian and avian lineages showed a significantly lower ratio of coding to total nucleotide changes in the avian viruses. Additionally, the avian virus lineages of both the H3 and H4 serotypes, but not the mammalian virus lineages, showed significantly greater conservation of amino acid sequence in the internal branches of the phylogenetic tree than in the terminal branches. The small number of amino acid differences between the avian viruses and the progenitor of the 1968 pandemic strain and the great phenotypic stability of the avian viruses suggest that strains similar to the progenitor strain will continue to circulate in birds and will be available for reintroduction into humans.

In vitro selection of lymphocytic choriomeningitis virus escape mutants by cytotoxic T lymphocytes

Cytotoxic T lymphocyte (CTL)-mediated cytolysis is induced via the interaction of the specific T-cell antigen receptor and the peptidic viral antigen associated with the major histocompatibility complex class I antigen. Here we demonstrate in vitro that lymphocytic choriomeningitis virus (LCMV) can escape the cytotoxic activity of LCMV-specific cloned CTLs by single amino acid changes within the recognized T-cell epitope defined by residues 275-289 of the LCMV glycoprotein [LCMV-GP-(275-289)]. LCMV-infected fibroblasts at a multiplicity of infection of 10(-3) exposed to virus-specific CTL at an effector-to-target cell ratio of 4:1 4 hr after infection was optimal for virus mutant selection. The selections were carried out with three LCMV-GP-(275-289)-specific CTL clones expressing T-cell antigen receptors containing the identical variable gene segments V alpha 4 and V beta 10 but different junctional regions; selection was also possible with LCMV-GP-(275-289)-specific cytotoxic polyclonal T cells. The most common escape mutation was an amino acid change of asparagine (AAT) to aspartic acid (GAT) at position 280; an additional mutation was glycine (GGT) to aspartic acid (GAT) at position 282. The results presented show that relevant point mutations within the T-cell epitope of LCMV-GP-(275-289) occur frequently and that they are selectable in vitro by CTLs.

Rapid antigenic-type replacement and DNA sequence evolution of canine parvovirus

Analysis of canine parvovirus (CPV) isolates with a panel of monoclonal antibodies showed that after 1986, most viruses isolated from dogs in many parts of the United States differed antigenically from the viruses isolated prior to that date. The new antigenic type (designated CPV type 2b) has largely replaced the previous antigenic type (CPV type 2a) among virus isolates from the United States. This represents the second occurrence of a new antigenic type of this DNA virus since its emergence in 1978, as the original CPV type (CPV type 2) had previously been replaced between 1979 and 1981 by the CPV type 2a strain. DNA sequence comparisons showed that CPV types 2b and 2a differed by as few as two nonsynonymous (amino acid-changing) nucleotide substitutions in the VP-1 and VP-2 capsid protein genes. One mutation, resulting in an Asn-Asp difference at residue 426 in the VP-2 sequence, was shown by comparison with a neutralization-escape mutant selected with a non-CPV type 2b-reactive monoclonal antibody to determine the antigenic change. The mutation selected by that monoclonal antibody, a His-Tyr difference in VP-2 amino acid 222, was immediately adjacent to residue 426 in the three-dimensional structure of the CPV capsid. The CPV type 2b isolates are phylogenetically closely related to the CPV type 2a isolates and are probably derived from a common ancestor. Phylogenetic analysis showed a progressive evolution away from the original CPV type. This pattern of viral evolution appears most similar to that seen in some influenza A viruses.

Genetic diversity of the attachment protein of subgroup B respiratory syncytial viruses

Respiratory syncytial (RS) virus causes repeated infections throughout life. Between the two main antigenic subgroups of RS virus, there is antigenic variation in the attachment protein G. The antigenic differences between the subgroups appear to play a role in allowing repeated infections to occur. Antigenic differences also occur within subgroups; however, neither the extent of these differences nor their contributions to repeat infections are known. We report a molecular analysis of the extent of diversity within the subgroup B RS virus attachment protein genes of viruses isolated from children over a 30-year period. Amino acid sequence differences as high as 12% were observed in the ectodomains of the G proteins among the isolates, whereas the cytoplasmic and transmembrane domains were highly conserved. The changes in the G-protein ectodomain were localized to two areas on either side of a highly conserved region surrounding four cysteine residues. Strikingly, single-amino-acid coding changes generated by substitution mutations were not the only means by which change occurred. Changes also occurred by (i) substitutions that changed the available termination codons, resulting in proteins of various lengths, and (ii) a mutation introduced by a single nucleotide deletion and subsequent nucleotide insertion, which caused a shift in the open reading frame of the protein in comparison to the other G genes analyzed. Fifty-one percent of the G-gene nucleotide changes observed among the isolates resulted in amino acid coding changes in the G protein, indicating a selective pressure for change. Maximum-parsimony analysis demonstrated that distinct evolutionary lineages existed. These data show that sequence diversity exists among the G proteins within the subgroup B RS viruses, and this diversity may be important in the immunobiology of the RS viruses.

NS 1 gene sequences from eight dengue-2 viruses and their evolutionary relationships with other dengue-2 viruses

The nucleotide sequences of the NS 1 genes from five Thai and three Sri Lankan dengue-2 viruses were determined by sequencing the viral RNA using synthetic oligonucleotide primers. The results were shown to be similar to four published dengue-2 NS 1 sequences and the classification of these genes was compared with the one obtained for the envelope genes of the same viruses. The classification was similar and showed that the Thai isolates could be divided into two separate groups and that the Sri Lankan isolates were distinct. We found no correlation between disease severity, serological response (1 degree or 2 degrees), or year of isolation and various aspects of NS 1 protein sequence variation; and no particular amino acid changes were correlated with virulence. The sequences were combined with those published and classified elsewhere to provide a comprehensive E/NS 1 gene taxonomy of dengue-2 virus isolates.

Quantitation of relative fitness and great adaptability of clonal populations of RNA viruses

We describe a sensitive, internally controlled method for comparing the genetic adaptability and relative fitness of virus populations in constant or changing host environments. Certain monoclonal antibody-resistant mutants of vesicular stomatitis virus can compete equally during serial passages in mixtures with the parental wild-type clone from which they were derived. These genetically marked "surrogate wild-type" neutral mutants, when mixed with wild-type virus, allow reliable measurement of changes in virus fitness and of virus adaptation to different host environments. Quantitative fitness vector plots demonstrate graphically that even clones of an RNA virus are composed of complex variant populations (quasispecies). Variants of greater fitness (competitive replication ability) were selected within very few passages of virus clones in new host cells or animals. Even clones which were well adapted to BHK21 cells gained further fitness during repeated passages in BHK21 cells.

Positive Darwinian evolution in human influenza A viruses

We earlier suggested that type A human influenza virus genes undergo positive Darwinian selection through immune surveillance. This requires more favorable amino acid replacements fixed in antigenic sites among the surviving lineages than among the extinct lineages. We now show that viral hemagglutinins fix proportionately more amino acid replacements in antigenic sites in the trunk of the evolutionary tree (survivors) than in the branches (nonsurvivors), demonstrating that type A human influenza virus is undergoing positive Darwinian evolution. The hemagglutinin gene is evolving 3 times faster than the nonstructural gene and the average age of the sampled nonsurvivors is only 1.6 years, so that extinction is not only common but rapid.

Molecular evidence for naturally occurring single VP7 gene substitution reassortant between human rotaviruses belonging to two different genogroups

Twenty four stool rotaviruses that comprised 22 distinct electropherotypes were selected for genome analysis from the collection of diarrheal specimens obtained over an eight-year period. These 22 electropherotypes were found in 46% of the total electropherotypes identified during the previous studies and represented 328 (64%) of rotavirus specimens in the collection. When genomic RNAs from these stool rotaviruses were hybridized to the 32P-labeled transcription probes prepared from prototypes representing three human rotavirus genogroups, Wa, DS-1, and AU-1, any one of the isolates showed a high degree of homology only with one of the three probes, which data confirmed and extended our previous observation on the existence of three distinct genogroups among human rotaviruses. Two stool rotaviruses which had an unusual combination of serotype (G1), subgroup (I) and RNA pattern (an identical short electropherotype), however, yielded the hybridization pattern indicative of an intergenogroupic single VP7 gene substitution reassortant. When they were cell culture adapted and analyzed by RNA-RNA hybridization, molecular evidence was obtained indicating that their VP7 gene derived from viruses belonging to the Wa genogroup whereas the remaining 10 genes hybridized with viruses belonging to the DS-1 genogroup. Interestingly, these natural reassortants emerged in the midst of the rotavirus season in which G1 strains predominated.

Molecular evolution of hemagglutinin genes of H1N1 swine and human influenza A viruses

The hemagglutinin (HA) genes of influenza type A (H1N1) viruses isolated from swine were cloned into plasmid vectors and their nucleotide sequences were determined. A phylogenetic tree for the HA genes of swine and human influenza viruses was constructed by the neighbor-joining method. It showed that the divergence between swine and human HA genes might have occurred around 1905. The estimated rates of synonymous (silent) substitutions for swine and human influenza viruses were almost the same. For both viruses, the rate of synonymous substitution was much higher than that of nonsynonymous (amino acid altering) substitution. It is the case even for only the antigenic sites of the HA. This feature is consistent with the neutral theory of molecular evolution. The rate of nonsynonymous substitution for human influenza viruses was three times the rate for swine influenza viruses. In particular, nonsynonymous substitutions at antigenic sites occurred less frequently in swine than in humans. The difference in the rate of nonsynonymous substitution between swine and human influenza viruses can be explained by the different degrees of functional constraint operating on the amino acid sequence of the HA in both hosts.