The global antigenic diversity of swine influenza A viruses
Abstract
Swine influenza presents a substantial disease burden for pig populations worldwide and poses a potential pandemic threat to humans. There is considerable diversity in both H1 and H3 influenza viruses circulating in swine due to the frequent introductions of viruses from humans and birds coupled with geographic segregation of global swine populations. Much of this diversity is characterized genetically but the antigenic diversity of these viruses is poorly understood. Critically, the antigenic diversity shapes the risk profile of swine influenza viruses in terms of their epizootic and pandemic potential. Here, using the most comprehensive set of swine influenza virus antigenic data compiled to date, we quantify the antigenic diversity of swine influenza viruses on a multi-continental scale. The substantial antigenic diversity of recently circulating viruses in different parts of the world adds complexity to the risk profiles for the movement of swine and the potential for swine-derived infections in humans.
https://doi.org/10.7554/eLife.12217.001eLife digest
Influenza viruses, commonly called flu, infect millions of people and animals every year and occasionally causes pandemics in humans. The immune system can neutralise flu viruses by recognising the proteins on the virus surface, generically referred to as antigens. These antigens change as flu viruses evolve to escape detection by the immune system. These changes tend to be relatively small such that exposure to one flu virus generates immunity that is still effective against other related flu viruses. However, over time, the accumulation of these small changes can result in larger differences such that prior infections no longer provide protection against the new virus.
Influenza A viruses infect a wide variety of birds and mammals. Viruses can also transmit from one species to another, which may result in the introduction of viruses with antigens that are new to the recipient species and which have the potential to cause substantial outbreaks. Pig flu viruses have long been considered to be a potential risk for human pandemic viruses and were the source of the 2009 pandemic H1N1 virus. Importantly, humans often transmit flu viruses to pigs. Understanding the dynamics and consequences of this two-way transmission is important for designing effective strategies to detect and respond to new strains of flu.
Influenza A viruses of the H1 and H3 subtypes circulate widely in pigs. However, it was poorly understood how closely related swine and human viruses circulating in different regions were to one another and how much the antigens varied between the different viruses.
Lewis, Russell et al. have now analysed the antigenic variation of hundreds of H1 and H3 viruses from pigs on multiple continents. The antigenic diversity of recent swine flu viruses resembles the diversity of H1 and H3 viruses observed in humans over the last 40 years. A key factor driving the diversity of the H1 and H3 viruses in pigs is the frequent introduction of human viruses to pigs. In contrast, only one flu virus from a bird had contributed to the observed antigenic diversity in pigs in a substantial way.
Once in pigs, human-derived flu viruses continue to evolve their antigens. This results in a tremendous diversity of flu viruses that can be transmitted to other pigs and also to humans. These flu viruses could pose a serious risk to public health because they are no longer similar to the current human flu strains. These findings have important implications not only for developing flu vaccines for pigs but also for informing the development of more-effective surveillance and disease-control strategies to prevent the spread of new flu variants.
https://doi.org/10.7554/eLife.12217.002Introduction
Swine have been hypothesized to be a mixing vessel for influenza viruses and were the source of the 2009 H1N1 human pandemic virus (Garten et al., 2009; Smith et al., 2009). The directionality and relative threat of transmission is generally assumed to be from swine into humans with attention only recently turning to introductions in the opposite direction – from humans into swine – and the role that such anthroponotic introductions play in influenza virus evolution (Nelson et al., 2012; Terebuh et al., 2010). Understanding the dynamics of influenza viruses at the interface between humans and swine is key for designing optimal surveillance and control strategies.
Currently there are three dominant subtypes of influenza A viruses circulating enzootically in swine: H1N1, H1N2, and H3N2. The evolutionary history of these viruses reflects multiple introductions of influenza viruses into swine from other species (Nelson et al., 2012; Vincent et al., 2014; Grøntvedt et al., 2013; Cappuccio et al., 2011; Ottis et al., 1982; Pereda, 2010; Njabo et al., 2012; Karasin et al., 2006; Hofshagen, 2009; Holyoake et al., 2011; Forgie et al., 2011; Pensaert et al., 1981). Much of the work on characterizing swine influenza viruses and their ancestry has been carried out in the U.S and in Europe (Vincent et al., 2014; Marozin et al., 2002; Nelson et al., 2014; Simon et al., 2014; Brown et al., 1993; Brown et al., 1997; Brown et al., 1998; Webby et al., 2000; Webby et al., 2004; Vincent et al., 2009a; Vincent et al., 2009b; Kyriakis et al., 2011; Watson et al., 2015) with evidence for numerous introductions of human seasonal H1 and H3 influenza viruses into swine since 1990. This pattern continued with the 2009 human pandemic H1N1 virus (H1N1pdm09) with at least 49 separate human-to-swine transmission events, some of which appear to have become established in pig populations (Nelson et al., 2012).
In addition to viruses from humans, avian influenza viruses also occasionally infect swine. Notably, in 1979, epidemics of an avian H1N1 influenza virus lineage were reported in Belgian swine leading to the establishment of an ‘avian-like’ H1N1 virus lineage in Europe. This Eurasian avian-like 1C lineage (Zell et al., 2013) has continued to circulate in swine until the present day, and been introduced into other geographic areas through international movement of swine (Pensaert et al., 1981; Zhu et al., 2013; Nelson et al., 2015c). Unsurprisingly, genetic diversity and multi-species origins of influenza viruses in swine are not restricted solely to pig populations in the U.S. and Europe. China has the largest swine population of any country and shows similar patterns of frequent introductions of both H1 and H3 influenza viruses from both avian and human hosts as well as introductions of swine lineages from other geographic regions (Chen et al., 2014; Vijaykrishna et al., 2011; Zhu et al., 2011).
Once established in swine, influenza viruses of human or avian origin pose a substantial threat to swine populations’ health and may spread to other geographically segregated populations (Nelson et al., 2015c). Similarly, these viruses also pose a threat for introduction or re-introduction into the human population (Nelson et al., 2014), with associated outbreak and pandemic risks. This risk is exemplified by the re-introduction of the H3N2v (variant) virus with surface glycoprotein genes of human origin to over 300 people in 2011–12 from an H3N2 virus that has been endemic in swine in the U.S. since the late 1990s (Centers for Disease Control and Prevention, 2011; Bowman et al., 2014).
The influenza virus hemagglutinin (HA) surface glycoprotein is the primary target of the humoral immune response, undergoes antigenic drift over time, and is the major antigenic constituent in influenza vaccines used in both humans and livestock. The genetic heterogeneity among the HAs of influenza viruses circulating in swine in certain geographic regions has been relatively well studied and characterized (Smith et al., 2009; Nelson et al., 2012; Cappuccio et al., 2011; Nelson et al., 2015c; Vijaykrishna et al., 2011; Kitikoon et al., 2013; Anderson et al., 2015; Takemae et al., 2008; 2013; Perera et al., 2013; Pereda et al., 2011). However, the relationships between genetic diversity and antigenic diversity among contemporary swine influenza viruses and their antigenic similarity to human seasonal influenza strains are unknown, largely owing to a lack of antigenic data. Such antigenic data are critical in assessing the phenotype of the virus, any potential cross-immunity with viruses circulating in humans and in swine, and the risk of re-emergence from these hosts. These analyses are also critical in the design of antigenically well-matched vaccines.
In order to help fill this knowledge gap, we compiled the largest and most geographically comprehensive antigenic dataset of swine influenza viruses to date. We then used a Bayesian framework (Bedford et al., 2014) to quantify the antigenic diversity of influenza viruses circulating in humans and swine thus generating a framework for exploring the antigenic component of risk of introduction of influenza A viruses from swine into humans and to swine in other geographic regions based on the antigenic characteristics of swine influenza viruses circulating in different parts of the world.
Results
Dimensionality of influenza virus antigenic evolution in swine
We used hemagglutination inhibition (HI) assays to antigenically characterize a geographically diverse collection of H1 viruses from swine and humans from 1930–2013 (n=194) and H3 viruses from swine and humans from 1968–2013 (n=379) and genetically sequenced all antigenically characterized viruses. Of these, 101 H1 and 73 H3 swine influenza viruses were newly characterized for this study.
In order to accurately quantify the antigenic variation of influenza viruses from swine we first determined the appropriate number of dimensions for accurately representing their antigenic relationships. Previous work on quantifying the antigenic evolution of seasonal influenza viruses showed that the evolution of H3 viruses in humans, as measured in HI assays, could be accurately visualized in two dimensions (2D) using antigenic cartography (Smith et al., 2004). In contrast, previous work on swine and equine H3 influenza viruses found that three dimensions (3D) were necessary to accurately visualize antigenic evolution using antigenic cartography (Lewis et al., 2011; 2014; de Jong et al., 2007; Lorusso et al., 2011). In this study we extensively tested the dimensionality of H1 and H3 influenza virus antigenic maps and found that the antigenic variation of both virus subtypes was most accurately represented in three dimensions.
Antigenic diversity of H1 influenza viruses in swine
To investigate the combined antigenic and genetic evolution of our study viruses we implemented a Bayesian multi-dimensional scaling (BMDS) method for inferring antigenic and genetic relationships (Bedford et al., 2014). As of 2013, three major H1 lineages circulated in swine: 1. the classical 1A lineage – descended from the ancestors of the 1918 human influenza pandemic, was first detected in swine in the 1930’s, and gave rise to the human H1N1pdm09 viruses and swine H1N1pdm09-like sub-lineage viruses; 2. the human seasonal-like 1B lineage – resulting from multiple introductions from human seasonal H1 viruses; 3. the Eurasian avian-like 1C lineage – arising from an introduction from wild birds into swine and first detected in Europe in the 1980’s (Figure 1A, Zell et al., 2013).
Our combined antigenic and genetic analyses of H1 viruses in swine showed an inconsistent relationship between antigenic clustering and genetic lineage through time with occasional emergence of new antigenic variants. In all geographic regions for which we had antigenic data we found unique patterns of antigenic evolution within individual lineages.
In Asia, we observed co-circulation of classical 1A, human seasonal-like 1B, and Eurasian avian-like 1C lineage viruses. However all three lineages were not detected in any one country. In Hong Kong – where samples were derived at slaughter from swine sourced from China – classical swine 1A lineage including the H1N1pdm09 sub-lineage, and Eurasian avian-like 1C lineage viruses were detected (Figure 1A, light red). For the viruses from Hong Kong, the mean pairwise antigenic distance (MPD) between the H1N1pdm09 sub-lineage viruses and the two classical 1A lineage clades was 5.0 and 6.0 AU away (Figure 1—figure supplement 1: Figure 1—figure supplement 2, Supplementary file 2). The Classical swine 1A lineage viruses collected in Japan (Figure 1A, dark red) differed from the Hong Kong H1 viruses with a MPD of 5.6 AU. Within the Hong Kong Eurasian avian-like 1C lineage viruses there were two distinct phylogenetic groups – group 1 and group 2 (as previously reported in Yang et al. (2016) and Figure 1—figure supplement 1, Supplementary file 2). The Group 1 and Group 2 viruses were 6.5 AU and 4.4 AU (MPD) away from the H1N1pdm09 sub-lineage viruses. The MPD between all Asian and European Eurasian avian-like 1C lineage strains was 5.2 AU. All of these distances are sufficiently large that it would be surprising if prior infection or un-adjuvanted vaccines would provide any protection among these virus groups.
In Europe, the Eurasian avian-like 1C lineage co-circulated with the H1N1pdm09 sub-lineage of classical 1A viruses and the human seasonal-like 1B lineage. The human seasonal-like 1B lineage viruses were first detected in European swine in the mid-1990’s but have a putative human seasonal ancestor from the 1980’s. The human seasonal-like 1B viruses differed substantially from the Eurasian avian-like swine 1C viruses and the H1N1pdm09 sub-lineage viruses (MPD: 6.7 AU and 10.4 AU respectively). The H1N1pdm2009 sub-lineage viruses also differed substantially from the human seasonal-like 1B lineage viruses (MPD:10.6 AU).
In Canada, only the classical swine 1A lineages viruses and H1N1pdm09 sub-lineage viruses were detected. Most of the viruses within the classical swine 1A lineage were antigenically similar to each other but differed from H1N1pdm09 sublineages viruses with a MPD of 6.5 AU. However, the alpha sub-lineage of classical 1A lineage viruses (Figure 1B green, Supplementary file 2) antigenically diverged in Canadian swine to form a separate antigenic group with a MPD of 7.0 AU from the H1N1pdm09 strains and 4.2 AU from the other Canadian classical swine 1A lineage viruses.
In the USA, the classical swine 1A lineage viruses co-circulated with the H1N1pdm09 sub-lineage and the human seasonal-like 1B lineage viruses. In the USA the human seasonal-like 1B lineage viruses evolved two sub-lineages – Delta 1 and Delta 2 (Figure 1—figure supplement 1, Supplementary file 2) – likely as the result of two separate but nearly contemporaneous introductions into swine first detected in the early 2000’s. The Delta 1 and 2 sub-lineage viruses antigenically differed from the H1N1pdm2009 lineage viruses with an MPD of 7.8 AU. Representatives of alpha sub-lineage of the classical swine 1A lineage in the USA were 7.0 AU (MPD) from the H1N1pdm2009 sub-lineage viruses. The beta and gamma sub-lineages of the classical swine 1A lineage were 6.6 and 4.7 AU (MPD) respectively from the H1N1pdm2009 viruses.
To assess the risk of virus introduction associated with moving live swine between continents we compared the MPDs of swine influenza H1 viruses circulating in different parts of the world. The human seasonal-like 1B viruses circulating in European swine were on average 11.7 AU (MPD) from the human seasonal-like 1B viruses in USA swine. The Eurasian avian-like 1C lineage viruses were 6.8 AU (MPD) from the human seasonal-like 1B lineage viruses in the USA.
Antigenic diversity of H3 influenza viruses in swine
Just as for the H1 viruses, we used BMDS (Bedford et al., 2014) to investigate the antigenic and genetic evolution of H3 viruses (Figure 1C,D). In European swine there were four H3 introductions included in our antigenic characterization (Figure 1—figure supplement 3: Figure 1—figure supplement 4). Of these, the European 3A lineage (A/Port Chalmers/1973-like) introduction was the most commonly detected lineage, showing evidence for three sequential antigenic clusters and evolving on a markedly different path from the putative ancestor in humans (Figure 1C,D [blue]). A separate introduction from humans into swine of an antigenically distinct 1970’s-like human seasonal virus also circulated in European swine in the 1980s and 1990s (Figure 1—figure supplement 1 [Europe 2 Introduction]). The earliest strain in the European 3A lineage was 6.5 AU from the A/Port Chalmers/1973 reference strain (Figure 1C large grey) and the most recent strain was 9.8 AU (Supplementary file 2). The distance from our representative of recent human seasonal H3 viruses (A/Victoria/361/2011) to the H3 European virus in swine and was 14.9 AU (MPD).
In the USA, we found evidence for three H3 virus introduction events from humans into swine in the mid-1990’s, one of which resulted in the establishment of the 3B lineage (Figure 1C,D) which became the most frequently detected lineage in swine through 2013, and gave rise to the H3N2v viruses which infected >300 humans in the US from 2011–2012 (Centers for Disease Control and Prevention, 2012). We only antigenically characterizes a single representative of the H3N2v viruses – A/Indiana/8/2011 – and found that it differed from A/Victoria/361/2011 by 6.4 AU (Figure 1 pink). The MPD between the swine H3 viruses circulating in the USA and in Europe was 12.6 AU.
In Asia we found evidence for multiple introductions of H3 viruses and subsequent antigenic drift from the six representative viruses in our dataset including an introduction likely mediated by translocation of a virus circulating in USA swine to Asian swine (Nelson et al., 2015d). Two separate introductions from humans were detected in Japan that resulted in circulation in swine beyond the time of circulation of the ancestral strain in humans. The introductions were ~5.2 AU from A/Victoria/361/2011 and 10.9 AU (MPD) from the H3 strains circulating in European swine. The most recent seasonal human H3 introduction detected in Asia was from Vietnam and this virus was 3.8 AU away from A/Victoria/361/2011 (Supplementary file 2).
Rates of antigenic drift in swine
To investigate the rates of swine influenza virus antigenic drift away from the most recent common ancestor for each genetic lineage in each geographic region we used TimeSlicer available in SPREAD (Bielejec et al., 2011). For H1 viruses in swine, the classical 1A lineage evolved at a mean rate of 0.15 AU per year from the earliest antigenically characterized representative in the dataset (A/swine/Iowa/15/1930) (Figure 2A, Supplementary file 3). The human seasonal-like 1B viruses introduced into European swine evolved at a mean rate of 0.17 AU per year (Figure 2B, Supplementary file 3). The introductions of human seasonal H1 viruses into USA swine that gave rise to the Delta-1 and Delta-2 lineages showed slightly different rates of antigenic drift – 0.63 for Delta-1 and 0.85 for Delta-2 – both significantly faster than the human seasonal introduction in European swine (Figure 2B, Supplementary file 3) but they were only observed over a relatively short time period. The Eurasian avian-like swine 1C lineage evolved at a mean rate of 0.15 AU per year (Figure 2C, Supplementary file 3). With the exception of the human seasonal-like 1B lineage viruses in the USA, the overall antigenic diversity of H1 viruses in swine remained remarkably constant over the study period (Figure 2A,B,C). The antigenic diversity of the Delta-1 and Delta-2 sublineages in the USA began to expand rapidly around 2009 for reasons that remain to be elucidated (Figure 2B).
For H3 viruses, we found that the 3B lineages viruses in USA swine evolved antigenically at a faster rate than the H3 viruses introduced into European swine (0.49 AU per year vs 0.28 AU per year [Figure 2D,E, Supplementary file 3]). Additionally, the antigenic diversity of H3 viruses in USA swine, particularly since 2010, has increased at a faster rate than that of H3 viruses in Europeans swine. Thus sustained circulation of antigenic drift variants within US swine has substantially increased the observed antigenic diversity in this lineage relative to the observed antigenic diversity in other geographic regions.
Additional antigenic diversity can be inferred from genetic sequence data
In addition to the BMDS analyses of viruses for which we had both antigenic and genetic data, we also performed large-scale phylogenetic analyses of publicly available swine influenza virus sequence. Through these analyses we inferred potentially important gaps in the global antigenic data that suggest our antigenic dataset underrepresents the true global antigenic diversity of swine influenza viruses.
For H1 viruses that circulated in swine as of 2013, we find that the major genetic lineages have been captured in this study (Figure 1—figure supplement 5: Figure 1—figure supplement 6). However, new genetic evidence recently became available to demonstrate additional human seasonal introductions into swine in previously under-surveilled regions such as Central and South America (Nelson et al., 2015a; Dibárbora et al., 2013; Nelson et al., 2015b) (Figure 1—figure supplement 5). Since genetic lineage does not always predict antigenic phenotype, we cannot make any inferences about their antigenicity.
Even in cases where we have antigenically-characterized clade representatives, the long branch lengths associated with some introductions could conceal antigenic drift away from the characterized clade representatives. The markedly different evolutionary patterns seen in the H3 lineages characterized in swine and the inferred antigenic profile of the viruses introduced to swine, suggests that the true diversity of H3 viruses circulating in swine is likely to be substantially greater than the diversity found in our BMDS analyses. In particular, we identified 19 additional examples of introductions of seasonal H3 viruses from humans in 10 countries (Figure 1—figure supplement 5 – black tips: Figure 1—figure supplement 7). Based on the likely time of these introductions from humans into swine we estimate that the introduced viruses represent at least 12 distinct antigenic variants based on the antigenic variation of their human seasonal influenza ancestors. Similar to the H1 genetic lineages described above, these additional H3 lineages also require antigenic testing to interpret the impact of these introductions on the global antigenic diversity of IAV in swine.
Risk posed by swine influenza viruses to the humans
In addition to the challenge the observed antigenic diversity poses to swine populations, the antigenic diversity of swine influenza viruses also creates substantial risks for human populations. The outbreak and pandemic potential of swine influenza viruses in humans is at least partially determined by human population immunity against swine viruses. This component of risk can be inferred from the antigenic distance of swine influenza viruses to seasonal viruses or vaccine strains in humans to which humans are likely to have immunity. Importantly, individuals are unlikely to have immunity to viruses that only circulated before they were born (Fonville et al., 2014).
Since most of the current swine influenza viruses are the product of human seasonal influenza virus introductions into swine, we anticipate at least some cross-protective immunity in the human population that could potentially interfere with re-introduction of these viruses. For example, the H1N1pdm09 viruses circulating in both humans and swine are very similar antigenically and likely induce at least some cross-immunity in both hosts (Figure 1A,B). However, for the H1 1C, H3 3A, and H3 3B human seasonal lineages in swine, the risk of re-introduction into the human population increases with the proportion of the human population that were born after the human precursor virus circulated in humans and is amplified by antigenic evolution of these viruses in swine. Thus, in terms of antigenic similarity and likely prior exposure, earlier introduced lineages of human H1 and H3 viruses – particularly those with precursors antigenically similar to the H3 1968 pandemic strain – pose the greatest current risk to humans because of the low or negligible predicted levels of cross immunity in individuals born since the 1970s. In comparison, the H3N2v viruses, which ultimately derive from a human seasonal virus introduction in the 1990s, have predominantly infected people born after the mid-1990s.
Discussion
In this study we characterized the antigenic diversity H1 and H3 viruses circulating in both swine and humans on a multi-continental scale. The observed diversity was largely driven by frequent introductions of variants from humans into swine and subsequent antigenic evolution within swine combined with long-term geographic segregation of multiple virus lineages. Much of the antigenic evolution of human-origin viruses in swine was via long and divergent evolutionary paths from those observed in humans. These divergent paths create threats for the re-introduction of viruses into humans as well as the possibility for antigenically novel strains in humans to be introduced into swine. Interestingly, the introduction of viruses from birds to swine has made only one major contribution to the antigenic diversity of influenza viruses in swine. This could be at least partially due to the antigenic similarity of the avian-like and the classical swine 1A lineages found here, although an exhaustive test of antigenic similarity between avian and swine H1 and H3 viruses could not be performed as part of this study and to the best of our knowledge has not been reported elsewhere. Successful introduction of avian influenza viruses into swine will also likely be influenced by other factors including mutations conferring adaption for mammalian transmission (Herfst et al., 2012; Imai et al., 2012) and other physiological and ecological factors.
In addition to multiple introduction events and subsequent antigenic drift, the patterns of antigenic diversity and heterogeneity among geographic regions are also influenced by the relative geographic isolation of swine populations resulting in region-specific patterns of virus evolution and circulation. The observed patterns of antigenic heterogeneity among geographic regions are also likely influenced by livestock production system differences, particularly regional differences in the movement of swine within and between animal holdings. USA producers commonly move swine among otherwise geographically isolated areas at different ages after weaning, whereas in the EU there is negligible intra- or inter-country movement of animals during a single production cycle. Such livestock movement patterns have been shown to influence the evolution of influenza viruses in swine in the USA, with the Midwest serving as an ecological sink for swine influenza virus with the source of genetic diversity located in the south east and south central states for some HA lineages (Nelson et al., 2011).
Vaccination is used extensively as a means to control influenza in swine in the USA and sporadically globally. Control strategies vary by region with some countries having no swine influenza vaccine usage, while others produce autologous (herd specific) vaccine for individual producers, which may be used consecutively with commercially-manufactured multi-strain vaccine when a new outbreak strain requires an update to vaccine composition. Most autologous and commercial products are multivalent due to the diversity of circulating strains. However, there is no formalized system for matching vaccine strains with circulating strains, nor validated protocols for standardization and effective vaccine use. The antigenic diversity quantified here is sufficiently large that it is highly unlikely that one strain per subtype would be efficacious globally, or even within a given region. Multiple within-subtype strains would likely be required to produce a vaccine that afforded adequate protection against most circulating variants using current inactivated vaccine production processes and application protocols.
Interestingly, H3 viruses circulating in European swine showed lower antigenic diversity than other regions. This limited antigenic diversity, combined with the use of oil adjuvanted vaccines, could explain why the prototype A/Port Chalmers/1973 (H3N2)-based vaccine induced protection against contemporary H3N2 swine influenza viruses isolated up to 35 years later (de Jong et al., 2007; Van Reeth and Ma, 2013; De Vleeschauwer et al., 2015). However, movement of swine among geographic regions could lead to the introduction of novel variants into a new sub-population of swine despite current import-export quarantine procedures between some countries. Therefore a more structured approach to vaccine strain selection in swine, where epidemiological, antigenic and genetic data are considered at a global level by the World Organization for Animal Health (OIE), similar to the process undertaken for creating equine influenza vaccine recommendations, could provide an evidence-based rationale for swine vaccine strain selection and improved vaccine efficacy in swine.
Age-related susceptibility within the human population could be one of the limiting factors in onward spread of recent human infections with swine viruses. However there remains a need for focused surveillance in areas with high swine population density and situations where humans and swine have opportunities for close contact to better assess the incidence of cross-species transmission and risk of onward transmission within the human population.
Neutralizing immunity to influenza viruses in both humans and swine is largely directed towards the influenza HA and the antigenic dissimilarity among influenza HAs within and between these two hosts is undoubtedly important for influenza epidemiology. However, the potential for a virus to emerge and spread widely in human or swine populations is likely to be a multi-genic trait comprising factors beyond just antigenicity (Trock et al., 2012). These factors include receptor binding (Matos-Patrón et al., 2015; Elderfield et al., 2014), adaptations to the matrix, non-structural, and neuraminidase proteins (Elderfield et al., 2014) as well as other factors differing between human and swine influenza viruses yet to be identified. There is also a need to better understand the competition dynamics of different influenza virus variants and subtypes in both human and swine populations and how these dynamics affect the potential for invasion by novel viruses.
Quantifying the public health risk of circulating swine influenza viruses in terms of immunological protection could be improved by assessing human population immunity to a selection of the antigenically-diverse swine viruses characterized in this study. A globally coordinated and systematic pipeline of antigenic and genetic analyses of relative antigenic distance between human and swine strains for subsequent testing of highly divergent swine strains against human sera (age-stratified) would add valuable information to the current WHO-led emergent pandemic risk assessment and inform the design of strategic vaccine stockpiles for human populations most at risk of infection with swine influenza virus, whether by age, geographic location or other factors.
Materials and methods
Viruses
Using hemagglutination inhibition (HI) assays, we characterized the antigenic properties of swine influenza A H1 and H3 virus strains circulating 1) in Europe, as part of the third program of the European Surveillance Network for Influenza in Pigs (ESNIP3) (www.esnip3.com), and 2) in Asia in collaboration with laboratories in Hong Kong, Japan, Vietnam and Thailand. We also expanded previously published datasets for North America. To this HI assay dataset we added previously published swine and human HI assay data for H1 and H3 influenza viruses (Lewis et al., 2014; de Jong et al., 2007; Lorusso et al., 2011; Koel et al., 2013; Nfon et al., 2011), resulting in an antigenic dataset consisting of 194 swine and human H1 influenza viruses and 379 swine and human H3N2 influenza viruses.
Viruses were propagated in Madin-Darby canine kidney (MDCK) cells or embryonated fowls’ eggs. Harvested cell culture supernatant or allantoic fluid was clarified by centrifugation. For antisera production virus was concentrated by ultracentrifugation over a 20% w/v sucrose cushion. Virus pellets were resuspended overnight at 4°C in sterile phosphate buffered saline at pH 7.4 and stored at - 70°C for use in immunizing naïve pigs.
Swine antisera production
Request a detailed protocolSwine antisera to influenza A viruses were generated by the United States Department of Agriculture National Animal Disease Centre, Ames, Iowa as previously described (Lorusso et al., 2011). Additional sera were generated by the French Agency for Food, Environmental and Occupational Health & Safety, Ploufragan, France, and by the Technical University of Denmark, Copenhagen, Denmark. To generate these additional antisera each animal was first inoculated intranasally with live influenza virus (7 to 8 log10 EID50 in a volume of 3 ml; 1.5 ml per nostril). Three weeks later equal volumes of the same inoculating virus and Freund’s complete adjuvant (total volume 3 ml) or inoculating virus and Montanide ISA206 (Seppic, Givaudan-Lavirotte, France) (total volume 2 ml) were administered by intramuscular/intradermal injection The animals were killed and bled two weeks after the last immunization. For sera raised in the USA, pigs were cared for in compliance with the Institutional Animal Care and Use Committee of the National Animal Disease Centre and for sera raised in the EU, all animal work was done in accordance with the local rules and procedures with ethical permissions granted following the codes of practice for performing scientific studies using animals.
Virus antigenic characterization
Request a detailed protocolHI assays using swine antisera were performed to compare the antigenic properties of swine and human influenza H1 and H3 viruses. Prior to HI testing, sera used for testing most H1 and all H3 viruses were treated with receptor-destroying enzyme (Sigma-Aldrich, MO, USA) and sera used for testing H1 viruses from the USA were treated with kaolin (Fisher Scientific, Pittsburg, PA, USA). All sera were then heat inactivated at 56°C for 30 min and adsorbed with 50% turkey red blood cells (RBC) to remove nonspecific inhibitors of hemagglutination. HI assays were performed by testing reference antisera raised to swine influenza viruses against selected H1 and H3 swine and human influenza viruses according to standard techniques. Serial 2-fold dilutions starting at 1:10 were tested for their ability to inhibit the agglutination of 0.5% turkey RBC with four hemagglutinating units of swine and human H1 and H3 viruses.
Antigenic cartography
Request a detailed protocolTo combine HI data generated in different laboratories we used a standard reference panel consisting of viruses and swine antisera raised to selected reference strains and shared among laboratories (Lewis et al., 2014; Lorusso et al., 2011). Our initial quantitative analyses of the antigenic properties of swine and human influenza H1 and H3 viruses used antigenic cartography methods as previously described for human (H3) and swine influenza A (H3) and (H1) viruses (Smith et al., 2004; Lewis et al., 2014; de Jong et al., 2007; Lorusso et al., 2011; Nfon et al., 2011).
To identify the most appropriate dimensionality of the antigenic maps for each virus subtype, we made first constructed antigenic maps of H1 and H3 viruses in 1, 2, 3, 4, and 5 dimensions and compared the HI titer differences to the resultant antigenic map distances for each dimension. Increasing from 1 to 2, and from 2 to 3 dimensions, resulted in significant improvement in correlation between HI data and map distances indicating that 3 dimensional (3D) maps provide a more accurate representation of the underlying HI data (Supplementary file 1). Further increasing map dimensionality resulted in smaller and statistically insignificant improvements in the quality of the fit of the map to the underlying HI data. Based on these results, we concluded that 3D maps were sufficient for capturing the antigenic variation of both H1 and H3 viruses in our combined swine and human datasets. Next we removed a random 10%, 20%, 30%, 40%, and 50% of the HI data and re-made the antigenic maps to test the robustness of the antigenic maps and the resolution at which they could visualize antigenic distances among viruses. The precision of point positioning in the 3D antigenic maps was 0.68 antigenic units (AU) for H1 and 0.85 AU for H3. One AU is equivalent to a two-fold difference in HI assay titer. Three AUs are considered sufficient antigenic difference to warrant an update of the human influenza vaccine.
Integrated analysis of antigenic and genetic evolution
Request a detailed protocolFor all isolates with both HA sequence data and HI measurements available, we implemented a Bayesian multidimensional scaling (BMDS) cartographic model to jointly infer antigenic and phylogenetic relationships of the viruses as described by Bedford et al. (2014)
In brief, time-resolved phylogenies were estimated for a total of 194 H1 and 379 H3 sequences using the Bayesian Markov chain Monte Carlo (MCMC) method implemented in BEAST (Drummond et al., 2012). We incorporated the Hasegawa-Kishino-Yano (HKY) model of nucleotide substitution model (Hasegawa et al., 1985) with codon partitioning at all three positions, and a Bayesian Skygrid coalescent model (Gill et al., 2013). A strict molecular clock and an uncorrelated lognormal relaxed clock were used for H3 and H1, respectively, based on linear regression analyses in Path-o-Gen v1.4 (http://tree.bio.ed.ac.uk/software/pathogen/). Ancestral sequence states were reconstructed using the ‘renaissance counting’ method (Lemey et al., 2012). Two independent MCMC chains were run for 200 million states with sampling every 20,000 states and a burn-in of 20 million states, and then combined and further subsampled every 180,000 states for a total of 2000 trees after assessing convergence in Tracer v1.6 (http://tree.bio.ed.ac.uk/software/tracer/). Along with the matching HI measurements for the virus isolates against post-infection ferret and swine antisera, these subsets of 2000 trees were used as an empirical tree distribution for the subsequent BMDS analysis to simultaneously model antigenic locations and evolutionary history.
The BMDS approach implements a Bayesian analog of antigenic cartographic models from Smith et al. (Smith et al., 2004) that arrange virus and serum locations in N-dimensions, such that Euclidean distances between locations are inversely proportional to serological cross-reactivity (Bedford et al., 2014). We generated the antigenic maps by imposing a weakly informative prior on expected antigenic locations such that antigenic distance increases with sampling time along one dimension and incorporated genetic data by modelling changes in antigenic phenotype as a diffusion process along the viral phylogeny, i.e. expected virus locations co-vary with genetic relatedness. Virus avidities and serum potencies were also estimated to account for experimental variation in serum and virus reactivity. MCMC was used to sample virus and serum locations in either two or three antigenic dimensions, as well as virus avidities, serum potencies, antigenic drift rate, MDS precision, virus and serum location precisions, and phylogenetic patterns. MCMC chains were run for 500 million states with sampling every 200,000 states, and checked for convergence by high ESS values in Tracer. Following burn-in of 100 million states, we obtained a total of 2,000 trees from which the maximum clade credibility tree was summarized in TreeAnnotator v1.8.2.
We used TimeSlicer available in SPREAD (Bielejec et al., 2011) to quantify the rates of antigenic evolution and levels of standing diversity over the posterior distribution of trees. This was achieved by ‘slicing’ through the BMDS-derived antigenic map location-tagged phylogenies at year-intervals, imputing the unobserved ancestral antigenic locations for branches that intersect those time points, and summarizing the drift front distance (in AU) from the most recent common ancestor (MRCA) and the mean pairwise distance between locations (in AU) for each lineage at that slice time. The overall antigenic drift rate was calculated by taking the total imputed drift front distance from the lineage MRCA divided by the total time elapsed, measured in AU per year. 95% high posterior density (HPD) estimates were used to measure the uncertainty in these inferences from the posterior sample of trees.
Other analyses
Request a detailed protocolH1 and H3 influenza A hemagglutinin (HA1) sequences representing the swine reference strains as well as viruses from swine populations in all available geographic regions were compiled from the NIAID Influenza Research Database (IRD) through the web site at http://www.fludb.org on 22/9/2014 (Squires et al., 2012) and combined with then unpublished sequences from European swine influenza viruses generated as part of the ESNIP3 project (Watson et al., 2015). We added HA1 human seasonal H1 and H3 influenza virus sequences from IRD that represented the entire time period of the study and all human vaccine strain selections.
Nucleotide alignments of the HA1 domain for both subtypes were generated using default settings in MAFFT with subsequent manual correction. For each alignment the sequences were curated for poorly curated data by inferring maximum likelihood phylogenetic trees in IQ-Tree (Nguyen et al., 2015) having performed the model test procedure to select the best-fit model. The provisional ML trees were then analyzed using Path-O-Gen and outlier sequences in the root-to-tip regression plots were identified and removed. ML phylogenetic trees were then inferred again using IQ-TREE as above and ultrafast bootstrap analyses performed (Minh et al., 2013).
Data availability
-
Data from: The global antigenic diversity of swine influenza A virusesAvailable at Dryad Digital Repository under a CC0 Public Domain Dedication.
References
-
Swine-to-human transmission of influenza A(H3N2) virus at agricultural fairs, Ohio, USA, 2012Emerging Infectious Diseases 20:1472–1480.https://doi.org/10.3201/eid2009.131082
-
Multiple genetic reassortment of avian and human influenza a viruses in European pigs, resulting in the emergence of an H1N2 virus of novel genotypeThe Journal of General Virology 79 (Pt 12):2947–2955.https://doi.org/10.1099/0022-1317-79-12-2947
-
Antigenic and genetic analyses of H1N1 influenza a viruses from European pigsThe Journal of General Virology 78 (Pt 3):553–562.https://doi.org/10.1099/0022-1317-78-3-553
-
Swine influenza outbreaks in England due to a new H1N1 virusThe Veterinary Record 132:461–462.https://doi.org/10.1136/vr.132.18.461
-
Outbreak of swine influenza in Argentina reveals a non-contemporary human H3N2 virus highly transmissible among pigsThe Journal of General Virology 92:2871–2878.https://doi.org/10.1099/vir.0.036590-0
-
Notes from the field: Outbreak of influenza a (H3N2) virus among persons and swine at a county fair--Indiana, July 2012MMWR. Morbidity and Mortality Weekly Report 61:561.
-
Swine-origin influenza a (H3N2) virus infection in two children--Indiana and Pennsylvania, July-August 2011MMWR. Morbidity and Mortality Weekly Report 60:1213–1215.
-
Identification of four genotypes of H3N2 swine influenza virus in pigs from southern ChinaArchives of Virology 159:2705–2709.https://doi.org/10.1007/s00705-014-2040-4
-
Antigenic and genetic evolution of swine influenza a (H3N2) viruses in EuropeJournal of Virology 81:4315–4322.https://doi.org/10.1128/JVI.02458-06
-
Swine influenza: clinical, serological, pathological, and virological cross-sectional studies in nine farms in ArgentinaInfluenza and Other Respiratory Viruses 7 Suppl 4:10–15.https://doi.org/10.1111/irv.12200
-
Bayesian phylogenetics with BEAUti and the BEAST 1.7Molecular Biology and Evolution 29:1969–1973.https://doi.org/10.1093/molbev/mss075
-
Swine outbreak of pandemic influenza a virus on a Canadian research farm supports human-to-swine transmissionClinical Infectious Diseases 52:10–18.https://doi.org/10.1093/cid/ciq030
-
Improving Bayesian population dynamics inference: a coalescent-based model for multiple lociMolecular Biology and Evolution 30:713–724.https://doi.org/10.1093/molbev/mss265
-
Influenza A(H1N1)pdm09 virus infection in Norwegian swine herds 2009/10: the risk of human to swine transmissionPreventive Veterinary Medicine 110:429–434.https://doi.org/10.1016/j.prevetmed.2013.02.016
-
Dating of the human-ape splitting by a molecular clock of mitochondrial DNAJournal of Molecular Evolution 22:160–174.https://doi.org/10.1007/BF02101694
-
Pandemic influenza A (H1N1) v: human to pig transmission in Norway?Euro Surveillance: Bulletin Europeen Sur Les Maladies Transmissibles=European Communicable Disease Bulletin 14:429–433.
-
The first identified case of pandemic H1N1 influenza in pigs in AustraliaAustralian Veterinary Journal 89:427–431.https://doi.org/10.1111/j.1751-0813.2011.00844.x
-
Identification of human H1N2 and human-swine reassortant H1N2 and H1N1 influenza a viruses among pigs in Ontario, Canada (2003 to 2005)Journal of Clinical Microbiology 44:1123–1126.https://doi.org/10.1128/JCM.44.3.1123-1126.2006
-
Genotype patterns of contemporary reassorted H3N2 virus in US swineThe Journal of General Virology 94:1236–1241.https://doi.org/10.1099/vir.0.51839-0
-
Antigenic and genetic evolution of equine influenza a (H3N8) virus from 1968 to 2007Journal of Virology 85:12742–12749.https://doi.org/10.1128/JVI.05319-11
-
DataData from: The global antigenic diversity of swine influenza a virusesDryad Digital Repository.https://doi.org/10.5061/dryad.5ff80
-
Genetic and antigenic characterization of H1 influenza viruses from United States swine from 2008The Journal of General Virology 92:919–930.https://doi.org/10.1099/vir.0.027557-0
-
Antigenic and genetic diversity among swine influenza a H1N1 and H1N2 viruses in EuropeThe Journal of General Virology 83:735–745.https://doi.org/10.1099/0022-1317-83-4-735
-
Ultrafast approximation for phylogenetic bootstrapMolecular Biology and Evolution 30:1188–1195.https://doi.org/10.1093/molbev/mst024
-
Global transmission of influenza viruses from humans to swineThe Journal of General Virology 93:2195–2203.https://doi.org/10.1099/vir.0.044974-0
-
Influenza a viruses of human origin in swine, BrazilEmerging Infectious Diseases 21:1339–1347.https://doi.org/10.3201/eid2108.141891
-
Global migration of influenza a viruses in swineNature Communications 6:6696.https://doi.org/10.1038/ncomms7696
-
Characterization of H1N1 swine influenza viruses circulating in Canadian pigs in 2009Journal of Virology 85:8667–8679.https://doi.org/10.1128/JVI.00801-11
-
IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogeniesMolecular Biology and Evolution 32:268–274.https://doi.org/10.1093/molbev/msu300
-
Pandemic A/H1N1/2009 influenza virus in swine, Cameroon, 2010Veterinary Microbiology 156:189–192.https://doi.org/10.1016/j.vetmic.2011.09.003
-
Evidence for the natural transmission of influenza a virus from wild ducts to swine and its potential importance for manBulletin of the World Health Organization 59:75–78.
-
Pandemic (H1N1) 2009 outbreak on pig farm, ArgentinaEmerging Infectious Diseases 16:304–307.https://doi.org/10.3201/eid1602.091230
-
Swine influenza in Sri LankaEmerging Infectious Diseases 19:481–484.https://doi.org/10.3201/eid1903.120945
-
Influenza research database: an integrated bioinformatics resource for influenza research and surveillanceInfluenza and Other Respiratory Viruses 6:404–416.https://doi.org/10.1111/j.1750-2659.2011.00331.x
-
Antigenic variation of H1N1, H1N2 and H3N2 swine influenza viruses in Japan and VietnamArchives of Virology 158:859–876.https://doi.org/10.1007/s00705-013-1616-8
-
Genetic diversity of swine influenza viruses isolated from pigs during 2000 to 2005 in ThailandInfluenza and Other Respiratory Viruses 2:181–189.https://doi.org/10.1111/j.1750-2659.2008.00062.x
-
Transmission of influenza a viruses between pigs and people, Iowa, 2002-2004Influenza and Other Respiratory Viruses 4:387–396.https://doi.org/10.1111/j.1750-2659.2010.00175.x
-
Development of an influenza virologic risk assessment toolAvian Diseases 56:1058–1061.https://doi.org/10.1637/10204-041412-ResNote.1
-
Swine influenza virus vaccines: to change or not to change-that's the questionCurrent Topics in Microbiology and Immunology 370:173–200.https://doi.org/10.1007/82_2012_266
-
Review of influenza a virus in swine worldwide: a call for increased surveillance and researchZoonoses and Public Health 61:4–17.https://doi.org/10.1111/zph.12049
-
Evolution of swine H3N2 influenza viruses in the United StatesJournal of Virology 74:8243–8251.https://doi.org/10.1128/JVI.74.18.8243-8251.2000
-
Prevalence, genetics, and transmissibility in ferrets of Eurasian avian-like H1N1 swine influenza virusesProceedings of the National Academy of Sciences of the United States of America 113:392–397.https://doi.org/10.1073/pnas.1522643113
-
Genetics, evolution, and the zoonotic capacity of European Swine influenza virusesCurrent Topics in Microbiology and Immunology 370:29–55.https://doi.org/10.1007/82_2012_267
-
History of Swine influenza viruses in AsiaCurrent Topics in Microbiology and Immunology 370:57–68.https://doi.org/10.1007/82_2011_179
Article and author information
Author details
Funding
United States Department of Agriculture-ARS (SCA agreement number 58-3625-2-103F)
- Nicola S Lewis
European Union Seventh Framework Programme (259949)
- Nicola S Lewis
- Paul Kellam
- Gaelle Simon
- Richard J Webby
- Kristien Van Reeth
- Sharon M Brookes
- Ian H Brown
- Lars Larsen
- Simon J Watson
- Amy L Vincent
Medical Research Council Fellowship (MR/K021885/1)
- Judith M Fonville
United States Department of Agriculture-ARS (SCA agreement number 58-3625-4-070)
- Tavis K Anderson
University Research Fellowship from the Royal Society
- Colin A Russell
DEFRA and the UK devolved governments (SV34041)
- Sharon M Brookes
- Ian H Brown
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication. Mention of trade names or commercial products in this article is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the U.S. Department of Agriculture.
Acknowledgements
We acknowledge the co-author contributions of members of the ESNIP3 consortium as follows: Scott Reid (AHPA), Montserrat Agüero Garcia, Laboratorio Central de Veterinaria, Spain; Timm Harder, Friedrich-Loeffler-Institut, Germany; Emanuela Foni, Instituto Zooprofilattico Sperimentale delle Venezie, Italy; Iwona Markowska-Daniel, National Veterinary Research Institute, Poland.
Ethics
Animal experimentation: For sera raised in the U.S., pigs were cared for in compliance with the Institutional Animal Care and Use Committee of the National Animal Disease Centre and for sera raised in the EU, all animal work was done in accordance with the local rules and procedures with ethical permissions granted following the codes of practice for performing scientific studies using animals.
Copyright
© 2016, Lewis et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Epidemiology and Global Health
Artificially sweetened beverages containing noncaloric monosaccharides were suggested as healthier alternatives to sugar-sweetened beverages. Nevertheless, the potential detrimental effects of these noncaloric monosaccharides on blood vessel function remain inadequately understood. We have established a zebrafish model that exhibits significant excessive angiogenesis induced by high glucose, resembling the hyperangiogenic characteristics observed in proliferative diabetic retinopathy (PDR). Utilizing this model, we observed that glucose and noncaloric monosaccharides could induce excessive formation of blood vessels, especially intersegmental vessels (ISVs). The excessively branched vessels were observed to be formed by ectopic activation of quiescent endothelial cells (ECs) into tip cells. Single-cell transcriptomic sequencing analysis of the ECs in the embryos exposed to high glucose revealed an augmented ratio of capillary ECs, proliferating ECs, and a series of upregulated proangiogenic genes. Further analysis and experiments validated that reduced foxo1a mediated the excessive angiogenesis induced by monosaccharides via upregulating the expression of marcksl1a. This study has provided new evidence showing the negative effects of noncaloric monosaccharides on the vascular system and the underlying mechanisms.
-
- Epidemiology and Global Health
- Microbiology and Infectious Disease
Influenza viruses continually evolve new antigenic variants, through mutations in epitopes of their major surface proteins, hemagglutinin (HA) and neuraminidase (NA). Antigenic drift potentiates the reinfection of previously infected individuals, but the contribution of this process to variability in annual epidemics is not well understood. Here, we link influenza A(H3N2) virus evolution to regional epidemic dynamics in the United States during 1997—2019. We integrate phenotypic measures of HA antigenic drift and sequence-based measures of HA and NA fitness to infer antigenic and genetic distances between viruses circulating in successive seasons. We estimate the magnitude, severity, timing, transmission rate, age-specific patterns, and subtype dominance of each regional outbreak and find that genetic distance based on broad sets of epitope sites is the strongest evolutionary predictor of A(H3N2) virus epidemiology. Increased HA and NA epitope distance between seasons correlates with larger, more intense epidemics, higher transmission, greater A(H3N2) subtype dominance, and a greater proportion of cases in adults relative to children, consistent with increased population susceptibility. Based on random forest models, A(H1N1) incidence impacts A(H3N2) epidemics to a greater extent than viral evolution, suggesting that subtype interference is a major driver of influenza A virus infection ynamics, presumably via heterosubtypic cross-immunity.