Next Article in Journal
The Cajal Body in Plant-Virus Interactions
Previous Article in Journal
Tick-Borne Encephalitis Virus: An Emerging Ancient Zoonosis?
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 12
Issue 2
10.3390/v12020248
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Substantial Antigenic Drift in the Hemagglutinin Protein of Swine Influenza A Viruses

by
Pia Ryt-Hansen
1,2,*,
Anders Gorm Pedersen
3,
Inge Larsen
2,
Charlotte Sonne Kristensen
4,
Jesper Schak Krog
5,
Silke Wacheck
6 and
Lars Erik Larsen
2
1
National Veterinary Institute, Technical University of Denmark, Kemitorvet Building 204, DK-2800 Kongens Lyngby, Denmark
2
Dpt. of Veterinary and Animal Sciences, University of Copenhagen, Grønnegårdsvej 2, DK-1870 Frederiksberg C, Denmark
3
Department of Health Technology, Section for Bioinformatics, Technical University of Denmark, Kemitorvet Building 208, DK-2800 Kongens Lyngby, Denmark
4
SEGES, Danish Pig Research Centre, Agro Food Park 15, DK-8200 Aarhus, Denmark
5
Statens Serum Institut, Artillerivej 5, DK-2300 Copenhagen S, Denmark
6
Ceva Santé Animale 10 Avenue de la Ballastière, 33500 Libourne, France
*
Author to whom correspondence should be addressed.
Viruses 2020, 12(2), 248; https://0-doi-org.brum.beds.ac.uk/10.3390/v12020248
Submission received: 28 January 2020 / Revised: 18 February 2020 / Accepted: 20 February 2020 / Published: 23 February 2020
(This article belongs to the Section Animal Viruses)
Browse Figures
Review Reports Versions Notes

Abstract

:
The degree of antigenic drift in swine influenza A viruses (swIAV) has historically been regarded as minimal compared to that of human influenza A virus strains. However, as surveillance activities on swIAV have increased, more isolates have been characterized, revealing a high level of genetic and antigenic differences even within the same swIAV lineage. The objective of this study was to investigate the level of genetic drift in one enzootically infected swine herd over one year. Nasal swabs were collected monthly from sows (n = 4) and piglets (n = 40) in the farrowing unit, and from weaners (n = 20) in the nursery. Virus from 1–4 animals were sequenced per month. Analyses of the sequences revealed that the hemagglutinin (HA) gene was the main target for genetic drift with a substitution rate of 7.6 × 10−3 substitutions/site/year and evidence of positive selection. The majority of the mutations occurred in the globular head of the HA protein and in antigenic sites. The phylogenetic tree of the HA sequences displayed a pectinate typology, where only a single lineage persists and forms the ancestor for subsequent lineages. This was most likely caused by repeated selection of a single immune-escape variant, which subsequently became the founder of the next wave of infections.

1. Introduction

Novel influenza A viruses (IAV) can develop through two different mechanisms; genome reassortment (antigenic shift) and gradual accumulation of mutations (antigenic drift). Genome reassortment occurs as a consequence of the segmented genome of IAV, when RNA segments originating from different subtypes/strains are mixed during assembly of progeny virions, leading to the formation of new subtypes/strains that may have novel antigenic properties [1,2]. Antigenic drift is a much slower process where the error prone RNA polymerase causes misincorporation of nucleotides during genome replication [3,4]. Mutations in coding regions of the viral genome are either synonymous or non-synonymous. Non-synonymous mutations occurring in immunogenic epitopes can undergo positive selection driven by host immunity, and may lead to the virus escaping e.g. neutralizing antibodies. The major epitopes of IAV, also termed antigenic sites, are located on the globular head of the hemagglutinin (HA) molecule, which is encoded by the HA1 domain. Several antigenic sites have been identified for both the H1 and H3 subtypes [5,6,7,8,9,10]. For humans, positive selection in these sites has been documented [11,12,13,14], and as little as one mutation in an antigenic site has been shown to affect the vaccine effectiveness [15,16]. Human influenza vaccines are therefore, evaluated twice a year to prevent mismatches between vaccine strain and circulating strains [17].
The rate of antigenic drift of swine IAV (swIAV) has generally been believed to be much lower than that of human IAV, mainly due to the short lifespan of pigs and the acute nature of the infection historically seen in pig herds, which limits the impact of pre-existing immunity [3,18,19,20,21,22]. Consequently, the swine influenza vaccines are updated less frequently [23]. Previous studies on the antigenic drift of swine hemagglutinin of the H1 or H3 subtypes, has mainly focussed on the global or national evolution [19,24,25,26,27,28,29]. For the H1 subtype in swine the nucleotide substitution rates have been estimated to range between 1.9–4.4 × 10−3 per site per year [11,30,31,32,33], whereas the nucleotide substitution rate of the swine H3 subtype has been documented to be as high as 6 × 10−3 per site per year [34]. While the highest rates are comparable to that of human H1 subtypes, the selection pressure expressed as the ratio of nonsynonymous to synonymous mutations has been found to be lower in swine compared to humans [11,22,35,36]. However, a very recent study revealed similar rates of synonymous and non-synonymous nucleotide substitutions within the swine H1 lineages 1B and 1C compared to that of human IAV, when investigating the genetic drift of swIAV sampled between 2003–2015 in Germany [37].
Over the past ten years, the understanding of swIAV circulation in swine herds has changed and it is now recognized that an infection with swIAV is likely to result in an enzootic infected herd [38,39,40,41,42,43,44,45]. This is probably a consequence of increasing herd sizes, which provides a continuous flow of naïve piglets [46,47]. The possible impact that herd level persistence of IAV might have on antigenic drift over time in a specific population of pigs, has, to the best of our knowledge, never been investigated. However, we believe it is highly important to get an increased understanding of antigenic drift occurring within single herds, as it can help explain the high genetic diversity within swIAV lineages documented in large-scale investigations and surveillance programs [25,26,27,29,48,49]. Importantly, if positive selection comparable to that observed in human IAV occurs in swine IAV, the possible effects on herd immunity and vaccination, should be taken into consideration when designing vaccines and evaluating swIAV control programs. We here report the results of a repeated cross-sectional study where we investigated the dynamics and viral evolution within a single sow herd over a one-year period.

2. Materials and Methods

2.1. Herd Description

The herd consisted of 480 sows and stables for 2000 nursery pigs. The farrowing stables were divided into two units and had no sectioning between different age groups, with weekly farrowings. The farrowing unit was cleaned once a year, without the use of disinfectants. At four weeks of age, piglets were weaned into a heated nursery. The nursery stables contained seven separate rooms with separate airflow, and all rooms were cleaned and disinfected between batches. The nursery pigs were housed in the nursery, until they were sold at approx. 30 kilos. Gilts were recruited internally (annual replacement rate of approx. 50%) and were subjected to eight weeks of quarantine from 12 weeks-of-age. Thus, no pigs were introduced from an outside source into the herd during the study period. According to the Danish Specific Pathogen Free program [50], the herd was free from infection with Mycoplasma hyopneumoniae, Actinobacillus pleuropneumoniae serotype, 2, 6, and 12, PRRSv type 1 and 2, Brachyspira hyodysenteriae, Pasteurella multocida, Sarcoptes scabiei var. Suis, and Haematopinus suis. The herd used a low number of nursing sows and minimized cross fostering of piglets. No vaccination against IAV was performed. The herd experienced recurrent respiratory symptoms in both the farrowing and nursery unit, and tested positive for IAV in July 2017, where the subtype “H1avN2sw” was diagnosed by full genome sequencing.

2.2. Study Design

Nasal swabs were collected monthly from November 2017 to October 2018. The nasal swabs were obtained from 20 piglets from four one-week-old litters (five piglets per litter), 20 piglets from four three-week-old litters (five piglets per litter) and 20 pigs from four pens with five-week-old nursery pigs (five pigs per pen). In addition, nasal swabs were collected from each sow of the one-week-old litters, and the parity of the sows was recorded. In total 64 nasal swabs were collected each month corresponding to 768 nasal swabs obtained over the full year. The individual samples and sequences were given an ID ranging from F1–F12 according to which month they were sampled, F1 being the first month (November 2017) and F2 being the second month etc. Moreover, the sample ID also included the age of the pigs, W1, W3, and W5 indicating week 1, 3, and 5, respectively. In addition, sequences obtained previously in July 2017, was included in the genetic analysis and named “W00_W1_01”.

2.3. Sampling

Nasal swabs were collected from the sows and piglets using a large or small rayon dry swab (Medical Wire, UK), respectively. The swabs were inserted into each nostril and turned 360 degrees. Thereafter, the swab was inserted in a tube containing 1ml Sigma Virocult media (Medical Wire, Corsham, UK) and kept at 2–8 °C for maximum 24 hours before being processed.

2.4. Coughing Index

For each litter/pen, included in the sampling, a coughing index was calculated, as described in a previous study [45]. Briefly, the coughing index was calculated by dividing the total number of coughs and sneezes with the number of pigs in the litter/pen multiplied by the time observed (three minutes).

2.5. Pooling and RNA Extraction

All the nasal swabs from pigs were initially pooled prior to extraction. Subsequently, the two most positive pools of each sampling time was identified, and the individual nasal swabs of the pools was subjected to RNA extraction and real-time RT PCR, to identify samples for viral isolation and sequencing. All nasal swabs were mixed using a Vortexer, and 100 µL was extracted for the pool. All individual samples (excluding the samples from sows) were pooled litter-wise, with five nasal swabs in each pool. The pools were mixed and centrifuged. Subsequently, 200 µL suspension was transferred to a new tube containing 400 µL RLT-buffer (QIAGEN, Copenhagen, Denmark) with 2-Mercaptoethanol (Merck, Darmstadt, Germany). Subsequently, the RNA was extracted using the RNeasy mini kit (QIAGEN, Copenhagen, Denmark) automated on the QIAcube (QIAGEN, Copenhagen, Denmark) according to large sample protocol version 2.

2.6. Real-Time RT PCR

In order to determine if a pool was positive for IAV, a previously published real-time RT PCR assay targeting the matrix gene of IAV [51] was adopted. Briefly, the published primers and the OneStep RT-PCR kit (QIAGEN, Copenhagen, Denmark) was used for the PCR mix, which was subsequently run on the Rotor-Gene Q (QIAGEN, Copenhagen, Denmark) using the following program: 50 °C, 30 min; 95 °C, 15 min; cycling 45× (95 °C for 10 s, 60 °C for 20 s, 64 °C for 1 sec, 68 °C for 1 sec, 72 °C for 30 s). All samples were tested in duplicates and the pool was only considered positive if both samples gave a positive result and had a Ct value <36. All positive samples with a Ct value <31 were tested to determine the IAV subtype using the previously published multiplex real-time RT PCR assay [52] with the modifications described in a previous study [45] and run on the Rotor-Gene Q (QIAGEN, Copenhagen, Denmark). The individual samples of two most positive pools of each sampling were also tested. Positive individual samples were selected for viral isolation and sequencing.

2.7. Viral Isolation and NGS

The monthly nasal swab with the lowest Ct value was selected for viral isolation. The nasal swabs were first subjected to sterile filtration using a 0.45 µM Millex-HP Millipore filter (Merck, Darmstadt, Germany) and then grown in Madin-Darby Canine Kidney (MDCK) cells under the conditions as described in a previous study [53]. After incubation, the RNA was extracted from the supernatant of each cell isolate using the same method as described above, however performed manually. Subsequently, the RNA was subjected to PCR amplification of each IAV segment and prepared for sequencing on the Illumina MiSeq platform using the methods described in a previous study [53].

2.8. HA and NA Amplification and Sanger Sequencing

Additional individual nasal swabs, which had Ct values <31, of each sampling were subjected to HA and NA PCR amplification and subsequent Sanger sequencing using the same methods as described in a previous study [53].

2.9. Consensus Sequence Generation

The determination of consensus nucleotide and amino acid sequences based on the Illumina and Sanger sequencing data was done as previously described [53] using the program CLC genomics Workbench version 11.0.1 (Aarhus, Denmark) and CLC main workbench version 8 (Aarhus, Denmark), respectively.

2.10. Characterization of the Herd swIAV Strain

The nucleotide and amino acid consensus sequences of each of the eight gene segments were aligned using the MUSCLE algorithm [54] in CLC main workbench version 8. The sequences of the alignments were compared using the pairwise comparison tool. In addition, the lineage of each gene segment was determined by aligning the respective sequences with contemporary swIAV sequences obtained in the Danish swine IAV surveillance program and subsequently neighbor-joining trees were constructed. Furthermore, HA amino acid sequences were annotated for known antigenic sites (Sa, Sb, Ca1, Ca2, and Cb) [5,6,8,55,56], B-cell epitopes [57,58] and T-cell epitopes [59,60,61,62], which were then manually checked for variation. Potential asparagine-linked glycosylation sites of the HA protein was predicted by the program NetNGlyc 1.0 Server [63].

2.11. Molecular Clock Analysis and Positively Selected Sites

Neighbor joining trees were constructed for each of the eight gene segments using the CLC main workbench version 8 software (Aarhus, Denmark). The eight gene segments obtained from the diagnostic sample of the same herd sampled approx. 4 months earlier (July 2017), were included in all analyses and used as an outgroup. The resulting tree, including information on sampling dates, was subsequently checked for the presence of a temporal signal (i.e., whether nucleotide changes accumulate proportionally to elapsed time) using the program TempEST [64]. Thereafter, the software package BEAST2 version 2.5.2 [65] was used to determine the substitution rate of each of the eight gene segments as previously described [66]. Briefly, the substitution model was specified to be HKY with gamma distributed rates over sites, with a strict clock model, and using tip dates (sampling dates). The chain length was set to 10,000,000 with a log every 1000, and the MCMC was run twice. The program BEAUti [65] was used to set up the analysis with all priors. Summaries of results and checking of MCMC convergence was done using the program Tracer, version 1.7.1 [67].
The program CODEML of the PAML package [68] was used to identify positively selected sites in all 8 genome segments as described in a previous study [66]. Briefly, we did this by comparing the fits of CODEML’s substitution models 1a (M1a) and 2a (M2a) (NSsites = 1 and 2). These substitution models include parameters for the ratio between the rate of non-synonymous substitutions per non-synonymous sites and the rate of synonymous substitutions per synonymous site (the dN/dS ratio, also indicated by ω). A dN/dS ratio above one indicates positive selection (there are more amino acid changing substitutions then expected). M1a includes two categories of codons—some under negative selection (dN/dS ratio < 1) and some codons where mutations are neutral (dN/dS ratio = 1). The model M2a includes three categories of codons—the same two as M1a plus an additional category of codons under positive selection (dN/dS ratio > 1). If M2a fits a dataset significantly better than M1a, then there is evidence of positive selection in some codons (and the identity of these codons is also found during model fitting). The fit of each model was compared using the Akaike Information Criterion (AIC) and likelihood ratio tests [69,70]. In addition, the average dN/dS ratio (global ω ratio) of all HA sequences was also estimated using CODEML (NSsites = 0). We also used CODEML to map codon (and the derived amino acid) mutations to branches in the HA phylogeny. Specifically, we used the M2a model (NSsites = 2) with RateAncestor = 1 to infer the most probable mapping of mutations to branches (mutational mappings are output in the result file “rst” along with ancestral reconstructions). The inferred amino acid changes were subsequently manually plotted on the HA clock-based phylogenetic tree (described below).
The program MrBayes [71] was used for reconstructing clock-based phylogenetic trees using codon-based substitution models, allowing simultaneous estimation of clock rates and detection of positively selected sites for the HA and NA gene segments. This method was described in a previous study [66]. The data analysis was performed using two parallel runs for 3.000.000 generations with a sample frequency of 600. The phylogenetic tree was inferred in a Bayesian framework and with MCMC sampling of posterior probabilities. We used the program Tracer version 1.7.1 to check convergence of MCMC runs [67]. Phylogenies were visualized using FigTree version 1.4.4, and the R-packages “tidyverse”, “treeio”, and “ggtree” in RStudio (version 0.97.551) [72,73,74,75,76,77].

2.12. Statistics

Student’s t-test was used to investigate if the average coughing index was significantly different between the IAV positive and IAV negative litters/pens. A chi-squared test was used to evaluate if 1st parity sows were more likely to be IAV positive compared to older sows, and the same test was also used to test if IAV positive sows were more likely to have an IAV positive litter compared to the IAV negative sows. All calculations were done using the GraphPad Software [78]. A likelihood ratio chi-squared test was used to test if M2a fit the data significantly better than M1a (indicating the presence of positively selected codons) [79]. Statistical significance was considered when the p-value was below 0.05.

3. Results

3.1. Presence of Enzootic IAV

The results of real-time RT PCR targeting the matrix gene revealed that IAV was present at all monthly samplings. Some variations were however observed between months: For example, very few litters/pens tested positive at F11 (September) and the Ct values of the samples did not allow for sequencing (Table 1 and Figure 1). The results of the test of the pooled samples of each month, showed that 60% of the one-week old litters, 69% of the three-week old litters and 60% of the pens with five-week old weaners tested positive over the entire study period. In total, 16 of 48 (33%) sows tested positive for IAV in the nasal swabs over the study period. The majority of IAV positive sows also had a positive litter (14/16). The prevalence of IAV positive litters from IAV positive sows (88%), was significantly higher than the prevalence of IAV positive litters from IAV negative sows (50%) (p = 0.03). For first parity sows, seven of 15 (47%) tested positive for IAV in nasal swabs as opposed to only nine of 33 (27%) of the ≥2nd parity sows. However, this difference was not significant (p value = 0.32).
The percentage of swIAV positive samples at week 1, 3, and 5 were based on pooled samples, whereas samples from sows were tested individually.

3.2. Correlation between IAV and the Coughing Index

For each age group (weeks 1, 3, and 5) the average coughing index (CI) of the IAV positive and negative litters/pens, was calculated (Table 2). Overall, the mean CI was significantly higher in litters/pens that included at least one pig testing positive for IAV in nasal swabs (p value = 0.03). No significant differences were discovered within the individual age groups, but a tendency towards a higher coughing index in the IAV positive litters was most evident at week 1 (p value = 0.07).

3.3. Herd Strain Characterization

In total, ten full genome sequences were obtained from cultured isolates based on individual monthly samples. However, it was not possible to obtain viral isolates from F7 and F11. Full genome sequences from one sample obtained in July 2017 were additionally included. Moreover, 19 HA and NA sequences were obtained from RNA of individual nasal swabs of all samplings except from F11. Thus, in total 30 HA and 29 NA sequences were generated, including 1-4 sequences from each sampling time (except F11).
The herd swIAV isolated throughout the study was of the H1N2 subtype, with a HA gene of Eurasian avian-like origin, and a NA gene of the swine adapted Hong Kong H3N2 origin. All gene segments of the internal gene cassette were of Eurasian avian-like origin. The HA gene segments had a pairwise nucleotide sequence identity ranging between 98.6–100% and, similarly, the NA gene segment had a pairwise nucleotide sequence identity ranging between 98.5–100%. The M, NP, PA, PB1, and PB2 gene segments had a pairwise nucleotide sequence identity ranging between 99.2–100% respectively, whereas the NS gene segments had a pairwise nucleotide sequence identity ranging between 97.5–100%. All sequences are available in NCBI Genbank with the following accession numbers: MN990039-MN990079 (Sanger sequences) and MN990117-MN990192 (NGS including internal genes).

3.4. Phylogenetic Analysis and Substitution Rates

The TempEst analysis revealed that all eight gene-segments, but especially the HA gene (correlation coefficient 0.95), showed association between genetic divergence through time and sampling dates indicating that a phylogenetic molecular clock-analysis (using BEAST and MrBayes) was suitable for the sequences (Table 3). Using BEAST the nucleotide substitution rate for the HA segment was estimated to be 7.6 × 10−3 substitutions/site/year, corresponding to 13 nucleotide substitutions per year for the entire gene (which is 1698 nucleotides long in this dataset). Estimated substitution rates were also high in the NA and NS segments (6.9 × 10−3 and 5.7 × 10−3 substitutions/site/year respectively), while the remaining segments had substitution rates ranging from 1.1 to 2.9 × 10−3 substitutions/site/year.
Interestingly, the phylogenetic tree based on the HA sequences displays the same imbalanced (so-called “comb-like” or pectinate) topology as that which is typical for human influenza trees spanning multiple years. The main feature of this type of topology is the repeated bottlenecks where only a single lineage persists and forms the ancestor for subsequent lineages (Figure 2). This is most likely caused by repeated selection of a single immune-escape variant that becomes the founder of the next wave of infection. The phylogenetic tree based on NA sequences had a somewhat comb-like topology also (Figure 3), however, with fewer bottleneck events (although firm conclusions in this regard is made difficult by greater uncertainty about branching pattern, and the resulting high level of polytomies in the tree).

3.5. Positive Selection

The program CODEML from the PAML package was used to test if positive selection was present in the eight gene segments. The results showed that the M2a model (indicating the presence of positive selection) fitted the HA sequences significantly better, whereas this was not the case for the remaining genes (Table 3). Estimates of dN/dS ratios for individual codons in the HA gene, under the M2a model, strongly indicated the presence of positive selection at position 536 (numbering from first methionine) in the HA2 part of the gene, which encodes the stalk region. Further analysis of the HA amino acid alignment showed that a mutation from tryptophan to arginine at position 536 was present in four pigs at F4, F6 and F8 (F4_W1_13, F6_W1_P16, F8_W1_16, F8_W1_17). Interestingly, this specific position is located in a B-cell epitope identified in the A(H1N1)pdm09 virus and in a T-cell epitope identified among human seasonal H1N1 (Table 4). The average dN/dS ratio of the HA sequences was estimated to be 0.19, while the value for codon 536 was estimated at 1.6, supporting that this position was under positive selection. Analysis of dN/dS ratios for individual codons using MrBayes identified several additional sites having a dN/dS ratio below 1, but significantly higher (p < 0.05) than the average dN/dS ratio (0.19) (Table 4). The majority of these mutations were present in the HA1 part of the HA gene and included seven mutations in known antigenic sites and several of the positions were in other known B- or T-cell epitopes. Further investigation of HA sequences identified seven mutations (D125N, P137L, K155R, V216D, E222K, H283Y, and I404F), that all showed a clear temporal pattern, where the given mutation became established at one time point and remained in all the following sequences until the end of the study.

4. Discussion

This study documented IAV persistence within a swine herd over a one-year period, supporting the increasing number of studies that have shown that IAV may persists within the herd [38,39,40,41,42,43,44,45]. Furthermore, the study confirms previous studies documenting a lack of seasonality of swIAV infections [28,80]. IAV was abundantly present in the one-week old litters, which in addition was the age group that showed the highest average viral load. These findings support the results of a pervious study performed by us, which identified IAV in nasal swabs of piglets from three days of age [45]. In general, a high percentage of litters were positive for IAV in the farrowing unit, probably as a consequence of only two farrowing stables being available, meaning that new-born piglets were housed side-by-side with piglets ready for weaning at 4 weeks of age. In turn, this provided an optimal environment for IAV transmission, as new naïve individuals were readily available for infection. The relatively high percentage of sows (33.3%) that tested positive for IAV in nasal swabs in this study, suggested that the sows had an important role in the transmission dynamics. In addition, a significantly higher number of IAV positive sows also had an IAV positive litter, which suggested a transmission from sow-to-piglet or piglet-to-sow. These findings are in accordance with previous work [43,44,45,81], although more studies are needed to firmly determine the directionality of transmission between sows and piglets. The high number of sows found positive for IAV in this study and the fact that they were positive for approx. 1.5–2 weeks after being introduced into the farrowing unit, emphasize the importance of stimulating sow immunity, especially before entering an environment where IAV is circulating. If the sows are inadequately immunized when entering the farrowing unit, there is a high risk of IAV infection occurring a few days before birth, which potentially could lead to birth complications and a lower milking yield, which in turn will result in compromised animal welfare and production economy. First parity sows were overrepresented among the IAV positive sows, even though this correlation was not statistically significant. In the herd investigated here, the internally recruited gilts were kept in a quarantine stable for eight weeks and were thereby not exposed to the herd strain several weeks prior to re-introduction into the sow-herd. This could explain why the gilts seemingly were more prone to IAV, compared to the older sows continuously exposed to the herd strain. This further highlights that proper gilt immunization is important, through either natural exposure or vaccination. In summary, the results of the study suggest that the constant production of naïve piglets and the continuous introduction of gilts play a key role in the herd-level persistence of swIAV.
The results of the study also highlighted the importance of herd management in the control of viral diseases. The same MCREBEL principles, which are widely used in controlling PRRSv infection [82], could also be a helpful tool in the prevention of IAV transmission within the herd. A clear sectioning and all in/all out management of weekly batches in the farrowing unit would most likely limit the IAV transmission significantly in the present herd.
The overall mean coughing index (CI) was significantly higher in IAV positive litters/pens. This supports our previous findings [45] which found a significant correlation between the CI and the pen/litter testing positive for IAV in nasal swabs. In turn, CI could be a helpful tool in identifying IAV positive litters/pens for diagnostics. Furthermore, the correlation to clinical symptoms of respiratory disease underline that IAV has an impact on the health and welfare of infected pigs.
Throughout the study, all eight genomic segments of the circulating H1avN2sw strain were highly similar (98.5–100%) to the initial consensus sequence, supporting that only a single IAV variant was circulating in the herd during the study. This is consistent with the information that the herd operated as a closed herd with no import of pigs, as the gilts were internally recruited. In turn, this provided the optimal scenario to study the viral drift within a single IAV strain, as the risk of reassortment events was limited. A correlation between sampling time and genetic diversity was found for all eight segments of the IAV strain, with the strongest temporal signal being present in the HA gene. The phylogenetic tree based on HA amino acid sequences had the distinct comb-like topology also known from human HA sequences, with a main trunk that represents the pathway of advantageous (mostly immune escape) mutations, which have been selected over time. Conversely, the shorter side branches represents the isolates that died out, because they were not able to avoid the host immunity [16,83,84,85]. Likelihood-based analysis using CODEML, confirmed the presence of positive selection in the HA gene. Further inspection of the HA amino acid alignment showed how amino acid mutations at specific positions arose at a given time point and remained until the end of the study, consistent with these amino acid changes being advantageous. All the mutations manually identified corresponded to positions in which the dN/dS ratio was estimated to be higher than the average (using Bayesian analysis). Moreover, the CODEML analysis defined the same mutations to be ancestral amino acid changes (presented in Figure 2), thereby supporting the presence of repeated bottlenecks. However, the dN/dS ratios remained below 1, probably due to a low number of sequences obtained at each sampling. Several of these positions were located in previously known antigenic sites (Sa, Sb, Ca1, Ca2, and Cb) or in other known B or T-cell epitopes, suggesting that positive selection occurred in immunogenic important sites. The single positively selected site (dN/dS ratio > 1) identified by the CODEML and Bayesian analysis, was located in the HA2 subunit of the HA gene, which encodes the stalk region. Interestingly, this position was included in a B-cell epitope proposed by a previous study investigating the A(H1N1)pdm09 subtype. In this study the reaction of the peptide encoding the epitope against a panel of swine sera raised against a panel of H1 and H3 subtypes were tested and showed a reaction, indicating this epitope to be highly immune reactive in pigs [58]. Moreover, the position was also included in a T-cell epitope defined in the human seasonal H1N1 subtype [59]. The results of the two studies indicate that the positive selection identified in position 536, could have an impact on the immunity, and thereby be important for immune escape variants of swIAV. However, this specific amino acid change did not show the same evolutionary bottleneck pattern as the other sites with elevated dN/dS ratios, as the change appeared at several independent sampling.
The nucleotide substitution rate calculated for the HA gene in this study exceeded the nucleotide substitution rate of 1.9 to 4.37 × 10−3 substitutions per site per year for swine H1Nx subtypes reported in previous studies [11,30,31,32,33]. The nucleotide substitution rate was also markedly higher than that reported for human H1N1 [11,86,87]. However, it is a well-known phenomenon, that when sequentially viral samples are collected over a short period of time, the clock rate is often estimated to be much higher compared to values estimated from samples collected over a longer period of time [88,89,90]. Therefore, it is difficult to compare the within-herd substitution rate calculated in this study, with the overall substitution rates calculated previously based on samples collected over several years and originating from many different herds or people. However, the nucleotide substitution rate of the NP gene corresponds to the results of an earlier study [91]. The average dN/dS (global ω ratio) of the HA sequences of this study was similar to that found in two studies investigating human H1N1 subtypes (ω ratio: 0.18–0.21) [32,92], but lower than that documented in two other studies (ω ratio: 0.24–0.38) [11,12]. The results of our study need to be confirmed by additional within-herd studies of the viral evolution, however, the data obtained in the present study suggested that the within herd evolution of swine IAV are comparable to that of human seasonal IAV. The rate of nucleotide substitution was high, and there was clear evidence for positive selection on the HA gene, especially in epitopes important for the adaptive and cellular immune response. Moreover, the topology of the phylogenetic tree indicated that immune escape variants were selected over time. Consequently, the results further confirmed that the antigenic drift of IAV in swine is comparable to seasonal human IAV. The extensive antigenic drift with generation of escape variants could possibly lead to acute clinical swIAV outbreaks even in herds without an external IAV introduction. In addition, the intensive viral drift could also have a negative impact on vaccine efficacy, however, this need to be studied in more detail. Nevertheless, it is recommendable for all herds to get their IAV strain/s sequenced. Then, the herd veterinarian will have an additional tool in explaining disease developments and it will contribute to our general understanding of the viral evolution of IAV in swine. From a human health point of view, it is also important to understand the viral evolution of IAV in swine herds, as they represents a reservoir for generation of future human pandemics, as demonstrated by the 2009 human pandemic [93,94]. The risk of generating new human pandemics will probably increase as an increasing number of herds become persistently infected, keeping a constant evolving reservoir of IAV circulating in swine.

5. Conclusions

In conclusion, the present study confirmed other recent studies [40,44,46,81,95] that found that swIAV infections should be regarded as enzootic infections with long term within-herd persistence. Our results also revealed that this change in epidemiology potentially affected the viral evolution, measured as increased diversity and selection of escape mutants. Finally, persistent circulation of swIAV in swine herds increase the likelihood for generation of re-assortments between human and swine IAV strains.

Author Contributions

P.R.-H. conceived the study, performed the laboratory analysis of the samples, interpreted the results and drafted the manuscript. A.G.P. helped conducting the bioinformatics analysis of the study, aided in the interpretation of the results, and corrected the manuscript. I.L. helped interpreting the results and corrected the manuscript. C.S.K. helped interpreting the results and corrected the manuscript. J.S.K. helped interpreting the results and corrected the manuscript. S.W. corrected the manuscript. L.E.L. planned the overall study, helped with the interpretation of the results, corrected, and approved the final draft of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by IDT Biologika GmbH and SEGES pig research center (grant #23352).

Acknowledgments

The authors would like to thank the participating herd and herd-veterinarian.

Conflicts of Interest

The project was partially funded by IDT Biologika GmbH and SEGES pig research center. The funders had no role in the design of the study; in the collection, analyses of data or in the decision to publish the results. However, both funders participated in the elaboration of the manuscript.

References

  1. Webster, R.G.; Laver, W.G.; Air, G.M.; Schild, G.C. Molecular mechanisms of variation in influenza viruses. Nature 1982, 296, 115–121. [Google Scholar] [CrossRef] [PubMed]
  2. Kim, H.; Webster, R.G.; Webby, R.J. Influenza Virus: Dealing with a Drifting and Shifting Pathogen. Viral Immunol. 2018, 31, 174–183. [Google Scholar] [CrossRef] [PubMed]
  3. Yoon, S.-W.; Webby, R.J.; Webster, R.G. Evolution and Ecology of Influenza A Viruses. Curr. Top. Microbiol. Immunol. 2014, 385, 359–375. [Google Scholar] [PubMed] [Green Version]
  4. Parvin, J.D.; Moscona, A.; Pan, W.T.; Leider, J.M.; Palesel, P. Measurement of the Mutation Rates of Animal Viruses: Influenza A Virus and Poliovirus Type 1. J. Virol. 1986, 59, 377–383. [Google Scholar] [CrossRef] [Green Version]
  5. Rudneva, I.; Ignatieva, A.; Timofeeva, T.; Shilov, A.; Kushch, A.; Masalova, O.; Klimova, R.; Bovin, N.; Mochalova, L.; Kaverin, N. Escape mutants of pandemic influenza A/H1N1 2009 virus: Variations in antigenic specificity and receptor affinity of the hemagglutinin. Virus Res. 2012, 166, 61–67. [Google Scholar] [CrossRef]
  6. Matsuzaki, Y.; Sugawara, K.; Nakauchi, M.; Takahashi, Y.; Onodera, T.; Tsunetsugu-Yokota, Y.; Matsumura, T.; Ato, M.; Kobayashi, K.; Shimotai, Y.; et al. Epitope Mapping of the Hemagglutinin Molecule of A/(H1N1)pdm09 Influenza Virus by Using Monoclonal Antibody Escape Mutants. J. Virol. 2014, 88, 12364–12373. [Google Scholar] [CrossRef] [Green Version]
  7. Gerhard, W.; Yewdell, J.; Frankel, M.E.; Webster, R. Antigenic structure of influenza virus haemagglutinin defined by hybridoma antibodies. Nature 1981, 290, 713–717. [Google Scholar] [CrossRef]
  8. Caton, A.J.; Brownlee, G.G.; Yewdell, J.W.; Gerhard, W. The antigenic structure of the influenza virus A/PR/8/34 hemagglutinin (H1 subtype). Cell 1982, 31, 417–427. [Google Scholar] [CrossRef]
  9. Luoh, S. Hemaglutinin mutations related to antigenic variation in H1 swine influenza- viruses. J. Virol. 1992, 66, 1066–1073. [Google Scholar] [CrossRef] [Green Version]
  10. Skehel, J.J.; Wiley, D.C. Receptor Binding and Membrane Fusion in Virus Entry: The Influenza Hemagglutinin. Annu. Rev. Biochem. 2000, 69, 531–569. [Google Scholar] [CrossRef]
  11. Furuse, Y.; Shimabukuro, K.; Odagiri, T.; Sawayama, R.; Okada, T.; Khandaker, I.; Suzuki, A.; Oshitani, H. Comparison of selection pressures on the HA gene of pandemic (2009) and seasonal human and swine influenza A H1 subtype viruses. Virology 2010, 405, 314–321. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Li, W.; Shi, W.; Qiao, H.; Ho, S.Y.; Luo, A.; Zhang, Y.; Zhu, C. Positive selection on hemagglutinin and neuraminidase genes of H1N1 influenza viruses. Virol. J. 2011, 8, 183. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Shen, J.; Ma, J.; Wang, Q. Evolutionary Trends of A(H1N1) Influenza Virus Hemagglutinin Since 1918. PLoS ONE 2009, 4, e7789. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Suzuki, Y. Positive selection operates continuously on hemagglutinin during evolution of H3N2 human influenza A virus. Gene 2008, 427, 111–116. [Google Scholar] [CrossRef] [PubMed]
  15. Petrie, J.G.; Parkhouse, K.; Ohmit, S.E.; Malosh, R.E.; Monto, A.S.; Hensley, S.E. Antibodies Against the Current Influenza A(H1N1) Vaccine Strain Do Not Protect Some Individuals From Infection With Contemporary Circulating Influenza A(H1N1) Virus Strains. J. Infect. Dis. 2016, 214, 1947–1951. [Google Scholar] [CrossRef] [PubMed]
  16. Petrova, V.N.; Russell, C.A. The evolution of seasonal influenza viruses. Nat. Rev. Microbiol. 2018, 16, 47–60. [Google Scholar] [CrossRef]
  17. Boni, M.F. Vaccination and antigenic drift in influenza. Vaccine 2008, 26, C8–C14. [Google Scholar] [CrossRef] [Green Version]
  18. Van Reeth, K.; Labarque, G.; De Clercq, S.; Pensaert, M. Efficacy of vaccination of pigs with different H1N1 swine influenza viruses using a recent challenge strain and different parameters of protection. Vaccine 2001, 19, 4479–4486. [Google Scholar] [CrossRef]
  19. Abente, E.J.; Rajao, D.S.; Santos, J.; Kaplan, B.S.; Nicholson, T.L.; Brockmeier, S.L.; Gauger, P.C.; Perez, D.R.; Vincent, A.L. Comparison of adjuvanted-whole inactivated virus and live-attenuated virus vaccines against challenge with contemporary, antigenically distinct swine H3N2 influenza A viruses. J. Virol. 2018. [Google Scholar] [CrossRef] [Green Version]
  20. McCauley, J.W.; Alexander, D.J.; Brown, I.H.; Olsen, C.W.; Ludwig, S.; Hannoun, C.; Strong, I.; Scholtissek, C.; Hinshaw, V.S.; Harris, P.A. Antigenic and genetic analyses of H1N1 influenza A viruses from European pigs. J. Gen. Virol. 1997, 78, 553–562. [Google Scholar]
  21. Campitelli, L.; Donatelli, I.; Foni, E.; Castrucci, M.R.; Fabiani, C.; Kawaoka, Y.; Krauss, S.; Webster, R.G. Continued Evolution of H1N1 and H3N2 Influenza Viruses in Pigs in Italy. Virology 1997, 232, 310–318. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. de Jong, J.C.; Smith, D.J.; Lapedes, A.S.; Donatelli, I.; Campitelli, L.; Barigazzi, G.; Van Reeth, K.; Jones, T.C.; Rimmelzwaan, G.F.; Osterhaus, A.D.M.E.; et al. Antigenic and Genetic Evolution of Swine Influenza A (H3N2) Viruses in Europe. J. Virol. 2007, 81, 4315–4322. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. O’Neill, K.C.; Shen, H.G.; Lin, K.; Hemann, M.; Beach, N.M.; Meng, X.J.; Halbur, P.G.; Opriessnig, T. Studies on Porcine Circovirus Type 2 Vaccination of 5-Day-Old Piglets. Clin. Vaccine Immunol. 2011, 18, 1865–1871. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Trebbien, R.; Bragstad, K.; Larsen, L.; Nielsen, J.; Bøtner, A.; Heegaard, P.M.; Fomsgaard, A.; Viuff, B.; Hjulsager, C. Genetic and biological characterisation of an avian-like H1N2 swine influenza virus generated by reassortment of circulating avian-like H1N1 and H3N2 subtypes in Denmark. Virol. J. 2013, 10, 290. [Google Scholar] [CrossRef] [PubMed]
  25. Watson, S.J.; Langat, P.; Reid, S.M.; Lam, T.T.-Y.; Cotten, M.; Kelly, M.; Van Reeth, K.; Qiu, Y.; Simon, G.; Bonin, E.; et al. Molecular Epidemiology and Evolution of Influenza Viruses Circulating within European Swine between 2009 and 2013. J. Virol. 2015, 89, 9920–9931. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Kuntz-Simon, G.; Madec, F. Genetic and Antigenic Evolution of Swine Influenza Viruses in Europe and Evaluation of Their Zoonotic Potential. Zoonoses Public Health 2009, 56, 310–325. [Google Scholar] [CrossRef]
  27. Rajao, D.S.; Anderson, T.K.; Kitikoon, P.; Stratton, J.; Lewis, N.S.; Vincent, A.L. Antigenic and genetic evolution of contemporary swine H1 influenza viruses in the United States. Virology 2018, 518, 45–54. [Google Scholar] [CrossRef]
  28. Simon, G.; Larsen, L.E.; Dürrwald, R.; Foni, E.; Harder, T.; Van Reeth, K.; Markowska-Daniel, I.; Reid, S.M.; Dan, A.; Maldonado, J.; et al. European Surveillance Network for Influenza in Pigs: Surveillance Programs, Diagnostic Tools and Swine Influenza Virus Subtypes Identified in 14 European Countries from 2010 to 2013. PLoS ONE 2014, 9, e115815. [Google Scholar] [CrossRef]
  29. Lewis, N.S.; Russell, C.A.; Langat, P.; Anderson, T.K.; Berger, K.; Bielejec, F.; Burke, D.F.; Dudas, G.; Fonville, J.M.; Fouchier, R.A.; et al. The global antigenic diversity of swine influenza A viruses. Elife 2016, 5. [Google Scholar] [CrossRef] [Green Version]
  30. Moreno, A.; Gabanelli, E.; Sozzi, E.; Lelli, D.; Chiapponi, C.; Ciccozzi, M.; Zehender, G.; Cordioli, P. Different evolutionary trends of swine H1N2 influenza viruses in Italy compared to European viruses. Vet. Res. 2013, 44, 112. [Google Scholar] [CrossRef] [Green Version]
  31. Aymard, M.; Cameron, K.; Hay, A.; Marozin, S.; Barigazzi, G.; Foni, E.; Bennett, M.; Gregory, V.; Lin, Y.; Valette, M. Antigenic and genetic diversity among swine influenza A H1N1 and H1N2 viruses in Europe. J. Gen. Virol. 2002, 83, 735–745. [Google Scholar]
  32. Lam, T.-Y.; Hon, C.-C.; Wang, Z.; Hui, R.K.-H.; Zeng, F.; Leung, F.C.-C. Evolutionary analyses of European H1N2 swine influenza A virus by placing timestamps on the multiple reassortment events. Virus Res. 2008, 131, 271–278. [Google Scholar] [CrossRef] [PubMed]
  33. Chastagner, A.; Hervé, S.; Bonin, E.; Quéguiner, S.; Hirchaud, E.; Henritzi, D.; Béven, V.; Gorin, S.; Barbier, N.; Blanchard, Y.; et al. Spatiotemporal Distribution and Evolution of the A/H1N1 2009 Pandemic Influenza Virus in Pigs in France from 2009 to 2017: Identification of a Potential Swine-Specific Lineage. J. Virol. 2018, 92. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Lewis, N.S.; Anderson, T.K.; Kitikoon, P.; Skepner, E.; Burke, D.F.; Vincent, A.L. Substitutions near the Hemagglutinin Receptor-Binding Site Determine the Antigenic Evolution of Influenza A H3N2 Viruses in U.S. Swine. J. Virol. 2014, 88, 4752–4763. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Sugita, S.; Yoshioka, Y.; Itamura, S.; Kanegae, Y.; Oguchi, K.; Gojobori, T.; Nerome, K.; Oya, A. Molecular evolution of hemagglutinin genes of H1N1 swine and human influenza A viruses. J. Mol. Evol. 1991. [Google Scholar] [CrossRef] [PubMed]
  36. Xu, Z.; Zhou, R.; Jin, M.; Chen, H. Selection pressure on the hemagglutinin gene of Influenza A (H1N1) virus: adaptation to human and swine hosts in Asia. Acta Virol. 2010, 54, 113–118. [Google Scholar] [CrossRef]
  37. Zell, R.; Groth, M.; Krumbholz, A.; Lange, J.; Philipps, A.; Dürrwald, R. Displacement of the Gent/1999 human-like swine H1N2 influenza A virus lineage by novel H1N2 reassortants in Germany. Arch. Virol. 2020. [Google Scholar] [CrossRef]
  38. Reynolds, J.J.H.; Torremorell, M.; Craft, M.E. Mathematical Modeling of Influenza A Virus Dynamics within Swine Farms and the Effects of Vaccination. PLoS ONE 2014, 9, e106177. [Google Scholar] [CrossRef]
  39. Chamba Pardo, F.O.; Alba-Casals, A.; Nerem, J.; Morrison, R.B.; Puig, P.; Torremorell, M. Influenza Herd-Level Prevalence and Seasonality in Breed-to-Wean Pig Farms in the Midwestern United States. Front. Vet. Sci. 2017, 4. [Google Scholar] [CrossRef] [Green Version]
  40. Allerson, M.W.; Davies, P.R.; Gramer, M.R.; Torremorell, M. Infection Dynamics of Pandemic 2009 H1N1 Influenza Virus in a Two-Site Swine Herd. Transbound. Emerg. Dis. 2014, 61, 490–499. [Google Scholar] [CrossRef]
  41. Simon-Grifé, M.; Martín-Valls, G.E.; Vilar, M.J.; Busquets, N.; Mora-Salvatierra, M.; Bestebroer, T.M.; Fouchier, R.A.; Martín, M.; Mateu, E.; Casal, J. Swine influenza virus infection dynamics in two pig farms; results of a longitudinal assessment. Vet. Res. 2012, 43, 24. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Cador, C.; Rose, N.; Willem, L.; Andraud, M. Maternally Derived Immunity Extends Swine Influenza A Virus Persistence within Farrow-to-Finish Pig Farms: Insights from a Stochastic Event-Driven Metapopulation Model. PLoS One 2016, 11, e0163672. [Google Scholar] [CrossRef] [PubMed]
  43. Cador, C.; Andraud, M.; Willem, L.; Rose, N. Control of endemic swine flu persistence in farrow-to-finish pig farms: a stochastic metapopulation modeling assessment. Vet. Res. 2017, 48, 58. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Rose, N.; Hervé, S.; Eveno, E.; Barbier, N.; Eono, F.; Dorenlor, V.; Andraud, M.; Camsusou, C.; Madec, F.; Simon, G. Dynamics of influenza A virus infections in permanently infected pig farms: evidence of recurrent infections, circulation of several swine influenza viruses and reassortment events. Vet. Res. 2013, 44, 72. [Google Scholar] [CrossRef] [Green Version]
  45. Ryt-Hansen, P.; Larsen, I.; Kristensen, C.S.; Krog, J.S.; Wacheck, S.; Larsen, L.E. Longitudinal field studies reveal early infection and persistence of influenza A virus in piglets despite the presence of maternally derived antibodies. Vet. Res. 2019, 50, 36. [Google Scholar] [CrossRef] [Green Version]
  46. Pitzer, V.E.; Aguas, R.; Riley, S.; Loeffen, W.L.A.; Wood, J.L.N.; Grenfell, B.T. High turnover drives prolonged persistence of influenza in managed pig herds. J. R. Soc. Interface 2016, 13, 20160138. [Google Scholar] [CrossRef]
  47. Brown, I.H. The epidemiology and evolution of influenza viruses in pigs. Vet. Microbiol. 2000, 74, 29–46. [Google Scholar] [CrossRef]
  48. Bolton, M.J.; Abente, E.J.; Venkatesh, D.; Stratton, J.A.; Zeller, M.; Anderson, T.K.; Lewis, N.S.; Vincent, A.L. Antigenic evolution of H3N2 influenza A viruses in swine in the United States from 2012 to 2016. Influenza Other Respi. Viruses 2019, 13, 83–90. [Google Scholar] [CrossRef]
  49. Ryt-Hansen, P.; Hjulsager, C.K.; Larsen, L.E. Overvågning af Influenza A virus i svin. Available online: https://www.vet.dtu.dk/raadgivning/aarsrapporter/aarsrapport-svineinfluenza (accessed on 1 July 2019).
  50. Svineproduktion, S. SPF Sundhedsstyringen a part of Landbrug & Fødevarer. Available online: http://spfsus.dk/en (accessed on 1 October 2017).
  51. Nagy, A.; Vostinakova, V.; Pirchanova, Z.; Cernikova, L.; Dirbakova, Z.; Mojzis, M.; Jirincova, H.; Havlickova, M.; Dan, A.; Ursu, K.; et al. Development and evaluation of a one-step real-time RT-PCR assay for universal detection of influenza A viruses from avian and mammal species. Arch. Virol. 2010, 155, 665–673. [Google Scholar] [CrossRef]
  52. Goecke, N.B.; Krog, J.S.; Hjulsager, C.K.; Skovgaard, K.; Harder, T.C.; Breum, S.Ø.; Larsen, L.E. Subtyping of Swine Influenza Viruses Using a High-Throughput Real-Time PCR Platform. Front. Cell. Infect. Microbiol. 2018, 8. [Google Scholar] [CrossRef]
  53. Ryt-Hansen, P.; Larsen, I.; Kristensen, C.S.; Krog, J.S.; Larsen, L.E. Limited impact of influenza A virus vaccination of piglets in an enzootic infected sow herd. Res. Vet. Sci. 2019, 127, 47–56. [Google Scholar] [CrossRef] [PubMed]
  54. Edgar, R.C. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Yang, H.; Qiao, C.; Tang, X.; Chen, Y.; Xin, X.; Chen, H. Human Infection from Avian-like Influenza A (H1N1) Viruses in Pigs, China. Emerg. Infect. Dis. 2012, 18, 1144–1146. [Google Scholar] [CrossRef] [PubMed]
  56. Manicassamy, B.; Medina, R.A.; Hai, R.; Tsibane, T.; Stertz, S.; Nistal-Villán, E.; Palese, P.; Basler, C.F.; García-Sastre, A. Protection of Mice against Lethal Challenge with 2009 H1N1 Influenza A Virus by 1918-Like and Classical Swine H1N1 Based Vaccines. PLoS Pathog. 2010, 6, e1000745. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Deem, M.W. The epitope regions of H1-subtype influenza A, with application to vaccine efficacy. Protein Eng 2009, 22, 543–546. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Wang, Z.; Huang, B.; Thomas, M.; Sreenivasan, C.C.; Sheng, Z.; Yu, J.; Hause, B.M.; Wang, D.; Francis, D.H.; Kaushik, R.S.; et al. Detailed mapping of the linear B Cell epitopes of the hemagglutinin (HA) protein of swine influenza virus. Virology 2018. [Google Scholar] [CrossRef]
  59. Babon, J.A.B.; Cruz, J.; Orphin, L.; Pazoles, P.; Co, M.D.T.; Ennis, F.A.; Terajima, M. Genome-wide screening of human T-cell epitopes in influenza A virus reveals a broad spectrum of CD4+ T-cell responses to internal proteins, hemagglutinins, and neuraminidases. Hum. Immunol. 2009, 70, 711–721. [Google Scholar] [CrossRef] [Green Version]
  60. H Gutiérrez, A. In Vivo Validation of Predicted and Conserved T Cell Epitopes in a Swine Influenza Model. PLoS ONE 2019, 11, e0159237. [Google Scholar] [CrossRef] [Green Version]
  61. Pedersen, L.E.; Breum, S.Ø.; Riber, U.; Larsen, L.E.; Jungersen, G. Identification of swine influenza virus epitopes and analysis of multiple specificities expressed by cytotoxic T cell subsets. Virol. J. 2014, 11, 163. [Google Scholar] [CrossRef]
  62. Baratelli, M.; Pedersen, L.E.; Trebbien, R.; Larsen, L.E.; Jungersen, G.; Blanco, E.; Nielsen, J.; Montoya, M. Identification of cross-reacting T-cell epitopes in structural and non-structural proteins of swine and pandemic H1N1 influenza A virus strains in pigs. J. Gen. Virol. 2017, 98, 895–899. [Google Scholar] [CrossRef]
  63. Bioinformatics, D. NetNGlyc 1.0 Server. 2017. Available online: http://www.cbs.dtu.dk/services/NetNGlyc/ (accessed on 23 July 2019).
  64. Rambaut, A.; Lam, T.T.; Max Carvalho, L.; Pybus, O.G. Exploring the temporal structure of heterochronous sequences using TempEst (formerly Path-O-Gen). Virus Evol. 2016, 2, vew007. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Bouckaert, R.; Heled, J.; Kühnert, D.; Vaughan, T.; Wu, C.H.; Xie, D.; Suchard, M.A.; Rambaut, A.; Drummond, A.J. BEAST 2: A software platform for Bayesian evolutionary analysis. PLoS Comput. Biol. 2014, 10, e1003537. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Ryt-Hansen, P.; Pedersen, A.G.; Larsen, I.; Krog, J.S.; Kristensen, C.S.; Larsen, L.E. Acute Influenza A virus outbreak in an enzootic infected sow herd: Impact on viral dynamics, genetic and antigenic variability and effect of maternally derived antibodies and vaccination. PLoS ONE 2019. [Google Scholar] [CrossRef] [PubMed]
  67. Rambaut, A.; Drummond, A.J.; Xie, D.; Baele, G.; Suchard, M.A. Posterior Summarization in Bayesian Phylogenetics Using Tracer 1.7. Syst. Biol. 2018, 67, 901–904. [Google Scholar] [CrossRef] [Green Version]
  68. Yang, Z. PAML 4: Phylogenetic Analysis by Maximum Likelihood. Mol. Biol. Evol. 2007, 24, 1586–1591. [Google Scholar] [CrossRef] [Green Version]
  69. Yang, Z. Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution. Mol. Biol. Evol. 1998, 15, 568–573. [Google Scholar] [CrossRef]
  70. Burnham, K. Model Selection and Multimodel Inference; Burnham, K.P., Anderson, D.R., Eds.; Springer: New York, NY, USA, 2004; ISBN 978-0-387-95364-9. [Google Scholar]
  71. Ronquist, F.; Teslenko, M.; van der Mark, P.; Ayres, D.L.; Darling, A.; Höhna, S.; Larget, B.; Liu, L.; Suchard, M.A.; Huelsenbeck, J.P. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst. Biol. 2012, 61, 539–542. [Google Scholar] [CrossRef] [Green Version]
  72. Rambaut, A. FigTree. Available online: http://tree.bio.ed.ac.uk/software/figtree/ (accessed on 6 June 2019).
  73. R Core Team (2019). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2019. Available online: https://www.R-project.org/http:/ (accessed on 14 February 2020).
  74. Yu, G.; Smith, D.K.; Zhu, H.; Guan, Y.; Lam, T.T.Y. ggtree: an r package for visualization and annotation of phylogenetic trees with their covariates and other associated data. Methods Ecol. Evol. 2017. [Google Scholar] [CrossRef]
  75. Yu, G.; Lam, T.T.Y.; Zhu, H.; Guan, Y. Two methods for mapping and visualizing associated data on phylogeny using GGTree. Mol. Biol. Evol. 2018. [Google Scholar] [CrossRef]
  76. Wang, L.; Lam, T.T.; Xu, S.; Dai, Z.; Zhou, L.; Feng, T.; Guo, P.; Dunn, C.W.; Jones, B.R.; Bradley, T.; et al. treeio: An R package for phylogenetic tree input and output with richly annotated and associated data. Mol. Biol. Evol. 2020, 37, 599–603. [Google Scholar] [CrossRef]
  77. Wickham, H.; Averick, M.; Bryan, J.; Chang, W.; McGowan, L.D.; François, R.; Grolemund, G.; Hayes, A.; Henry, L.; Hester, J.; et al. Welcome to the tidyverse. J. Open Source Softw. 1686, 4. [Google Scholar] [CrossRef]
  78. GraphPad software. Available online: https://www.graphpad.com/quickcalcs/ (accessed on 6 June 2018).
  79. Pedersen, A.G. Chi-squared table. Available online: http://www.cbs.dtu.dk/dtucourse/cookbooks/gorm/phd.phylo.prob/chi-squared.html (accessed on 1 August 2019).
  80. Harder, T.C.; grosse Beilage, E.; Lange, E.; Meiners, C.; Dohring, S.; Pesch, S.; Noe, T.; Grund, C.; Beer, M.; Starick, E. Expanded Cocirculation of Stable Subtypes, Emerging Lineages, and New Sporadic Reassortants of Porcine Influenza Viruses in Swine Populations in Northwest Germany. J. Virol. 2013, 87, 10460–10476. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Diaz, A.; Marthaler, D.; Corzo, C.; Muñoz-Zanzi, C.; Sreevatsan, S.; Culhane, M.; Torremorell, M. Multiple Genome Constellations of Similar and Distinct Influenza A Viruses Co-Circulate in Pigs During Epidemic Events. Sci. Rep. 2017, 7, 11886. [Google Scholar] [CrossRef] [PubMed]
  82. Subcommittee, A. REPORT Control of porcine reproductive and respiratory syndrome (PRRS) virus. Swine Health Prod. 2014, 4. [Google Scholar]
  83. Nelson, M.I.; Holmes, E.C. The evolution of epidemic influenza. Nat. Rev. Genet. 2007, 8, 196–205. [Google Scholar] [CrossRef] [PubMed]
  84. Fitch, W.M.; Leiter, J.M.; Li, X.Q.; Palese, P. Positive Darwinian evolution in human influenza A viruses. Proc. Natl. Acad. Sci. 1991, 88, 4270–4274. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. McHardy, A.C.; Adams, B. The Role of Genomics in Tracking the Evolution of Influenza A Virus. PLoS Pathog. 2009, 5, e1000566. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. Wilson, I.A.; Cox, N.J. Structural Basis of Immune Recognition of Influenza Virus Hemagglutinin. Annu. Rev. Immunol. 1990, 8, 737–787. [Google Scholar] [CrossRef]
  87. Rejmanek, D.; Hosseini, P.R.; Mazet, J.A.K.; Daszak, P.; Goldstein, T. Evolutionary Dynamics and Global Diversity of Influenza A Virus. J. Virol. 2015, 89, 10993–11001. [Google Scholar]
  88. Aiewsakun, P.; Katzourakis, A. Time-Dependent Rate Phenomenon in Viruses. J. Virol. 2016, 90, 7184–7195. [Google Scholar]
  89. Meyer, A.G.; Spielman, S.J.; Bedford, T.; Wilke, C.O. Time dependence of evolutionary metrics during the 2009 pandemic influenza virus outbreak. Virus Evol. 2015, 1, vev006. [Google Scholar] [CrossRef]
  90. Ho, S.Y.W.; Lanfear, R.; Bromham, L.; Phillips, M.J.; Soubrier, J.; Rodrigo, A.G.; Cooper, A. Time-dependent rates of molecular evolution. Mol. Ecol. 2011, 20, 3087–3101. [Google Scholar] [CrossRef] [PubMed]
  91. Jenkins, G.M.; Rambaut, A.; Pybus, O.G.; Holmes, E.C. Rates of Molecular Evolution in RNA Viruses: A Quantitative Phylogenetic Analysis. J. Mol. Evol. 2002, 54, 156–165. [Google Scholar] [CrossRef] [PubMed]
  92. Bragstad, K.; Nielsen, L.P.; Fomsgaard, A. The evolution of human influenza A viruses from 1999 to 2006 - a complete genome study. Virol. J. 2008, 5, 40. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  93. Garten, R.J.; Davis, C.T.; Russell, C.A.; Shu, B.; Lindstrom, S.; Balish, A.; Sessions, W.M.; Xu, X.; Skepner, E.; Deyde, V.; et al. Antigenic and Genetic Characteristics of Swine-Origin 2009 A(H1N1) Influenza Viruses Circulating in Humans. Science 2009, 325, 197–201. [Google Scholar] [CrossRef] [Green Version]
  94. Mena, I.; Nelson, M.I.; Quezada-Monroy, F.; Dutta, J.; Cortes-Fernández, R.; Lara-Puente, J.H.; Castro-Peralta, F.; Cunha, L.F.; Trovão, N.S.; Lozano-Dubernard, B.; et al. Origins of the 2009 H1N1 influenza pandemic in swine in Mexico. Elife 2016, 5. [Google Scholar] [CrossRef]
  95. Diaz, A.; Marthaler, D.; Culhane, M.; Sreevatsan, S.; Alkhamis, M.; Torremorell, M. Complete Genome Sequencing of Influenza A Viruses within Swine Farrow-to-Wean Farms Reveals the Emergence, Persistence, and Subsidence of Diverse Viral Genotypes. J. Virol. 2017, 91, e00745-17. [Google Scholar] [CrossRef] [Green Version]
Viruses 12 00248 g001 550
Figure 1. The average Ct value of the influenza A viruses (IAV) positive litters/pens and sows at each sampling time.
Figure 1. The average Ct value of the influenza A viruses (IAV) positive litters/pens and sows at each sampling time.
Viruses 12 00248 g001
Viruses 12 00248 g002 550
Figure 2. Bayesian strict molecular clock tree of the hemagglutinin (HA) sequences. The x-axis represents time in months from 2017–2019. Branch labels represent posterior clade probabilities. Ancestral amino acid changes are indicated below the branches on which they most probably occurred, with DTLC numbering. F00_W1_01 sampled approx. four months before the actual study was used as an outgroup.
Figure 2. Bayesian strict molecular clock tree of the hemagglutinin (HA) sequences. The x-axis represents time in months from 2017–2019. Branch labels represent posterior clade probabilities. Ancestral amino acid changes are indicated below the branches on which they most probably occurred, with DTLC numbering. F00_W1_01 sampled approx. four months before the actual study was used as an outgroup.
Viruses 12 00248 g002
Viruses 12 00248 g003 550
Figure 3. Bayesian strict molecular clock tree of the NA sequences. The x-axis represents time in months from 2017–2019. Branch labels represent posterior clade probabilities. F00_W1_01 sampled approx. four months before the actual study was used as an outgroup.
Figure 3. Bayesian strict molecular clock tree of the NA sequences. The x-axis represents time in months from 2017–2019. Branch labels represent posterior clade probabilities. F00_W1_01 sampled approx. four months before the actual study was used as an outgroup.
Viruses 12 00248 g003
Table
Table 1. Prevalence of swine influenza A viruses (swIAV) positive litters/pens and sows at each sampling time.
Table 1. Prevalence of swine influenza A viruses (swIAV) positive litters/pens and sows at each sampling time.
Week 1Week 3Week 5Sows
Nov75% (3/4)100% (4/4)75% (3/4)50% (2/4)
Dec75% (3/4)75% (3/4)25% (1/4)50% (2/4)
Jan100% (4/4)75% (3/4)100% (4/4)0% (0/4)
Feb75% (3/4)50% (2/4)100% (4/4)75% (3/4)
Mar75% (3/4)75% (3/4)50% (2/4)100% (4/4)
Apr50% (2/4)75% (3/4)75% (3/4)25% (1/4)
May100% (4/4)50% (2/4)75% (3/4)25% (1/4)
Jun25% (1/4)25% (1/4)25% (1/4)25% (1/4)
Jul25% (1/4)100% (4/4)75% (3/4)25% (1/4)
Aug50% (2/4)100% (4/4)50% (2/4)0% (0/4)
Sep0% (0/4)50% (2/4)50% (2/4)0% (0/4)
Oct75% (3/4)50% (2/4)25% (1/4)25% (1/4)
Total60.5% (29/48)70.8% (34/48)60.5% (29/48)33.3% (16/48)
Table
Table 2. Average coughing index (CI) in IAV positive and negative litters/pens.
Table 2. Average coughing index (CI) in IAV positive and negative litters/pens.
Week 1Week 3Week 5Total
IAV Positive0.12 (SD = 0.14)0.36 (SD = 0.27)0.09 (SD = 0.08)0.32 (SD = 0.23)
IAV Negative0.05 (SD = 0.08)0.30 (SD = 027)0.06 (SD = 0.04)0.12 (SD = 0.17)
p-value0.070.480.120.03
Table
Table 3. The best fitting substitution model, temporal correlation coefficient and nucleotide substitution rate of each gene segment.
Table 3. The best fitting substitution model, temporal correlation coefficient and nucleotide substitution rate of each gene segment.
M1aM2aSignificant DifferenceCorrelation CoefficientSubstitution Rate
AICAkaike WeightAICAkaike Weightp-Value
HA5494.300.11435490.200.8857<0.050.957.6 × 10−3
NA4203.240.94784209.060.0521>0.050.816.9 × 10−3
M2866.580.88102870.580.1189>0.050.822.5 × 10−3
NS2531.020.88102535.020.1189>0.050.685.7 × 10−3
NP4310.140.88104314.140.1189>0.050.771.1 × 10−3
PB16595.060.88106599.060.1189>0.050.862.47 × 10−3
PB26459.840.88106463.840.1189>0.050.942.94 × 10−3
PA6170.360.88106174.360.1189>0.050.602.13 × 10−3
M1a = model 1, which describes neutral and negative selection. M2a = model 2, which describes neutral, negative and positive selection. AIC = Akaike information criterion. The correlation coefficient of the TempEst analysis. The substitutions rate gives the results of the BEAST analysis expressed as nucleotide substitution rate per site per year.
Table
Table 4. Evidence of positive selected amino acid changes and their location in antigenic sites.
Table 4. Evidence of positive selected amino acid changes and their location in antigenic sites.
Position
(from DTLC)		MutationdN/dS RatioPR+Sequences Showing the MutationAntigenic Site/Epitope/Glycosylation Site
1D→E0.79820.27211/3 F8T-cell
43S→N0.82550.28373/3 F8B-cell
125D→N0.78160.26512/3 F6, 1/1 F7 and all F9-F12Sa
137P→L0.87910.30651/4 F2, 1/3 F3, 2/4 F4 and all F5–F12Ca1
155K→R0.76720.25882/3 F6, 1/1 F7 and all F9–F12Sa
156N→D0.79910.2725All F1-F12Ca1
170G→R0.87840.30621/4 F2 and 1/3 F3Ca1
193E→G0.87460.3046All F1–F12Sb
210F→Y0.82150.28201/2 F1B-cell epitope
216V→D0.84070.29032/3 F6 and all F7–F12B-cell epitope
222E→K0.77170.26072/3 F6 and all F7–F12Ca2
276T→N0.82820.28491/2 F1Glyco
283H→Y0.84020.28981/4 F4 and all F5–F12-
288S→N0.82550.28381/4 F4, 1/1 F5 and
1/3 F6		-
404I→F0.81020.2773All F10–F12T-cell epitope
442V→I0.80360.27451/3 F3T-cell epitope
536W→R3.156
1.5882		0.879
0.6380		1/4 F4, 1/3 F6 and
2/3 F8		B-cell and T-cell epitope
Column 1 indicates the position with an elevated dN/dS ratio identified by the Bayesian analysis (normal text) and in both the Bayesian analysis and CODEML (bold text). Column 2 indicates the amino acid change occurring over time. Column 3 and 4 gives the dN/dS ratio and the probability (PR+) of the position being positive selected with the results of the Bayesian analysis in normal text and the CODEML results in bold text. Column 5 presents the sequences wherein the given mutation was identified; herein F1–F12 indicates the sequences where the given mutation was identified. Column 6 specifies if the mutation was located in an antigenic site (Sa, Sb, Ca1, Ca2, and Cb) [5,6,8,55,56], glycosylation site (Glyco) [63] or a B-cell [57,58] or T-cell epitope [59,60,61,62].

Share and Cite

MDPI and ACS Style

Ryt-Hansen, P.; Pedersen, A.G.; Larsen, I.; Kristensen, C.S.; Krog, J.S.; Wacheck, S.; Larsen, L.E. Substantial Antigenic Drift in the Hemagglutinin Protein of Swine Influenza A Viruses. Viruses 2020, 12, 248. https://0-doi-org.brum.beds.ac.uk/10.3390/v12020248

AMA Style

Ryt-Hansen P, Pedersen AG, Larsen I, Kristensen CS, Krog JS, Wacheck S, Larsen LE. Substantial Antigenic Drift in the Hemagglutinin Protein of Swine Influenza A Viruses. Viruses. 2020; 12(2):248. https://0-doi-org.brum.beds.ac.uk/10.3390/v12020248

Chicago/Turabian Style

Ryt-Hansen, Pia, Anders Gorm Pedersen, Inge Larsen, Charlotte Sonne Kristensen, Jesper Schak Krog, Silke Wacheck, and Lars Erik Larsen. 2020. "Substantial Antigenic Drift in the Hemagglutinin Protein of Swine Influenza A Viruses" Viruses 12, no. 2: 248. https://0-doi-org.brum.beds.ac.uk/10.3390/v12020248

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop