Genetic Diversity of Type A Influenza Viruses Found in Swine Herds in Northwestern Poland from 2017 to 2019. The One Health Perspective

1. Introduction

Influenza A virus (IAV) belongs to the Orthomyxoviridae family and is the main causative agent of human and animal flu. The severity and geographical spread of the disease depend primarily on the origin and genetic evolution of the virus. Seasonal strains circulating among certain hosts typically cause waves of incidence. When some slight changes are introduced in viral surface glycoproteins due to accumulation of mutations, such strains may be a reason for local epidemics, whereas major changes, or even replacement of these proteins, can potentially lead to a pandemic [1].

The difference in the peptide sequences of two surface glycoproteins, hemagglutinin (HA) and neuraminidase (NA), is the basis of IAV classification into subtypes. Today, 18 types of HA and 11 types of NA are recognized. The vast majority of subtypes circulate among wild waterfowl and constitute a natural reservoir of IAV [2]. In the course of viral evolution several “spill-over” events followed by virus adaptation introduced the IAV into susceptible species like marine mammals, bats, horses, ferrets, pigs and humans, where certain adapted strains continue circulating and undergo selective environmental pressure. In the case of human infections, zoonosis is usually restricted to infected people and rarely develops into human-to-human transmitted disease, although it is possible in certain conditions and when it happens, it might be dangerous for the human population. [3].

To become a transmissible disease, novel IAV must replicate efficiently in a new host, and for this reason, they must undergo some adaptive changes that will allow them to cross the host species barrier [4]. The genome architecture of IAV and the process of viral assembly facilitate these events[5]. Eight segments of negative-polarity RNA encode viral proteins that have a significant impact on the specificity of the host, the efficiency of viral RNA replication, and evading the immune response of the host [6]. RNA-dependent RNA polymerase, which is responsible for the replication of the genetic material of the virus, is an enzyme without proofreading properties that leads to the formation of point mutations in each replication cycle [7]. It also occurs in segments coding for major viral antigens, causing the accumulation of changes. This so-called antigenic drift may result in peptide sequence alteration, which translates directly to structural changes in the surface proteins. If the changes are sufficiently large to avoid specific antibodies raised against these proteins, it may lead to an epidemic. The mutation frequency is so high that in a few years, we can observe changes in the major virus variants circulating worldwide [8]. Such differences necessitate continuous monitoring of the variability of the virus to design the most effective vaccine each year [9].

When various strains of the virus circulate in the same area, one host cell may be infected with two different viral strains. After possible reasortation, completely new set of genome segments originating from two different viruses would be incorporated into a viral particle giving rise to a whole new virus. This interchange in genetic segments is known as an antigenic shift. When such reassortment occurs between strains of different species specificity, there is a chance for a new dangerous strain with a pandemic potential to emerge. It would become widespread in a naive population with no previous contact with this strain, and due to other adaptations or changes (e.g., in viral polymerase subunits and receptor site recognition), it could efficiently replicate in new hosts, becoming a threat to the whole population [10].

One of the factors responsible for host specificity of IAV is the ability of HA to recognize sialic acids, which are specific receptors on the membrane of the host epithelial cells. These receptors play a crucial role in binding of viral particles and are linked to a sugar galactose by either SA-α-2,3-Gal or SA-α-2,6-Gal bond which is host specific [11]. Because most of the receptor bonds in the respiratory tract of birds are -2,3 bonds, whereas in the upper airways of humans, we find predominantly -2,6 bonds, strains circulating among humans predominantly recognize -2,6 bonds, whereas avian strains need -2,3 bonds for efficient binding. Absence of alpha-2,3 bonds in the upper respiratory tract of the human seems to be the main reason for limitation of the transmission of "avian" virus between humans [12]. For efficient "bird to human" virus transmission, "avian" strains must adapt and acquire the ability to recognize alpha-2,6 bonds [13]. Genetic reassortment during co-infection with different strains could be the basis for such adaptations. The presence of sialic acids linked by both alpha-2,3 and alpha-2,6 bonds has been observed in the airways of pigs. Thus, pigs susceptible to infection with both "avian" and "human" strains are believed to play the role of a closed subtype-mixing environment. This allows co-infection and creates conditions for rearrangement of the genetic segments of the virus [14].

The epidemiology of swine influenza A virus in Europe is complex and dynamic. Although only four major swIAV subtypes, H1N1, H3N2, H1N2, and H1N1pdm09, are the most prevalent, numerous constellations of internal genes with genes coding for HA and NA are possible. This abundance of genotypes is caused by different origins of RNA segments as a result of spill over or reassortment events in the past [15]. Until the end of 70’s a “classical swine” H1N1 from Eurasian lineage was the dominant strain of swIAV in European countries. In 1979, this “classical” strain was entirely replaced by the introduction of an antigenically distinct H1N1 strain of avian origin (H1avN1) that seemed to have a selective advantage over its predecessor [16]. The circulating H3N2 strain emerged after genomic reassortment, which took place in the early 1980s, between H1avN and human-like swine H3N2 viruses. The most recent reassortment event in 1994 in Great Britain gave rise to “human seasonal” H1N2 (H1huN2) strain, derived from human seasonal H1N1, swine H1avN1 and “human like” swine H3N2 strains [16,17]. H1N1pdm09 of pandemic origin has been circulating in the swine population since 2009 [18]. All these strains circulate among swine populations in Europe, and their different combinations of swine, human seasonal, human pandemic, and avian genes contribute to a plethora of reassortants. Some of these supplanted or even replaced endemic strains found in Denmark, Germany, or Sweden (H1avN2dk, H1pdm09N2, H3huN2, H1avN2) [19,20,21]. Similar events have been reported in Italy, Belgium, and France [22,23,24]. The high genetic variability of circulating strains and the growing number of passive and active surveillance programs have created a need for a uniform nomenclature for the classification of HA lineages and swIAV genotypes. In 2015 Watson et al. (2015) [25] proposed a complex genotype nomenclature for 23 combinations of genes (A-W) based on 290 swIAVs collected in 14 European countries between 2009 and 2013. Phylogenetic analysis was based on nine major swIAV lineages that were circulating worldwide. With additional genotypes identified by Henritzi et al. (2020) [26], this genotyping system is well-established and widely used. A new classification for H1 lineages was developed by Anderson et al. (2016) [27], who divided H1 into three main lineages: Classical swine lineage-1A (with H1pdm09), Human seasonal lineage-1B, and Eurasian avian lineage-1C. These two classification systems proved to be useful for organizing swIAV diversity.

Since the influenza pandemic in 2009, numerous cases of swine-to-human zoonotic influenza infections have been reported as well as a growing number of swine infections with human strains of the virus [28,29,30,31,32,33,34,35,36]. These findings highlight the importance of research in the field of influenza interhost relationships. In this context, surveillance and genome sequence analysis will prove useful for understanding the molecular evolution dynamics of influenza. The data we present here are part of a project aimed at assessing influenza A virus subtype diversity in Poland. In this study, we analyzed influenza A virus subtype diversity among swine herds in northwestern Poland along with phylogeny based on HA, NA, and PB2 genome segments. We discuss the necessity of large-scale surveillance of the common human-animal interface and possible ways to raise awareness of the threat posed by the emergence of a new viral strain with pandemic potential.

2. Materials and Methods

2.1. Collection of samples from swine and isolation of viral RNA.

Samples were collected from farms in the northwestern region of Poland from pigs of various ages exhibiting influenza-like symptoms. Nasal swabs were collected with FLOQSwabs^® (COPAN Diagnostics Inc.) and placed in 3 ml Universal Transport Medium (UTM™, COPAN Diagnostics Inc.). The specimens were transported to the laboratory in a portable refrigerator at 4°C – 8°C. In the laboratory, the swabs were vortexed for 10 s and centrifuged for 5 min at 1000 rpm. Aliquots of 1 ml were used immediately for viral RNA isolation or stored at 4°C – 6°C for no longer than 48 hours. The rest of the medium was stored at - 70°C for further use as an inoculum for virus isolation from the cell culture. Viral RNA isolation was performed using high-load volume funnel columns NucleoSpin^® RNA Virus F Kit (Macherey-Nagel GmbH) according to the manufacturer’s instructions.

2.2. Confirmation of influenza A virus infection in collected samples.

All samples were initially tested for the presence of influenza A type genetic material using the virotype^® Influenza A RT-PCR Kit (QIAGEN). Viral RNA isolated from swabs, positive control, negative control and internal control of β-actin mRNA, were subjected to Real Time PCR protocol according to manufacturer’s instructions on Light Cycler^® 480 System (ROCHE).

2.3. Virus isolation from Madin-Darby Canine Kidney (MDCK) cell culture.

Monolayers of MDCK cells (MERCK) were cultured in T-25 flasks in D-MEM growth media (Corning, [+] 4.5 g/L glucose, L-glutamine, [-] sodium pyruvate) supplemented with antibiotic/antimycotic solution (Invitrogen) and 10% fetal bovine serum (FBS). When the cell layer was 90% confluent, the growth medium was decanted and the cells were washed three times with 6 ml of pure D-MEM containing 2 μg/ml of TPCK-trypsin. One milliliter of each specimen (swabs in UTM) was used to inoculate T-25 flasks. The control flask was inoculated with 1 ml of pure UTM™. After one hour of incubation at 37°C and 5% CO_2, the inoculum was replaced with 6 ml of viral growth medium (D-MEM growth media with 2 μg/ml of TPCK-trypsin and without FBS). Each flask was inspected daily and cultured in 37°C and 5% CO₂ until a cytopathic effect (CPE) was observed, usually after 3-7 days. The influenza virus was harvested by collecting the supernatant supplemented with 0.5% glycerol for stabilization. If no CPE was detected after 10 days, blind passage was performed on a new T-25 flask 2 times before the sample was considered unable to isolate the virus from the specimen.

2.4. Nucleic acid extraction from viral isolates and amplification of the whole influenza A virus genome segments.

RNA from the viral isolates was extracted using the RNeasy^® Mini Kit (QIAGEN) according to the manufacturer’s instructions. Extracted RNA was used for whole genome amplification using the SuperScript™ III One-Step RT-PCR System with Platinum™ Taq High-Fidelity DNA Polymerase (Thermo Fisher Scientific). Briefly: 10 μl of extracted RNA from each isolate was mixed with 25 μl of master mix (containing 0.4 mM dNTPs, 2.4 mM MgSO4), 1 μl of SuperScript™ III RT/Platinium™ Taq High-Fidelity Enzyme Mix, 3 μl of 5μM mix of sense and antisense primers (sense: 5´-CTGGATACGCCAGCRAAAGCAGG-3´; antisense: 5´-GACCTGATGCGGAGTAGAAACAAGG-3´) and 11 μl of PCR grade water. The reactions were placed in a pre-heated ProFlex PCR System thermal cycler (Applied Biosystems) and amplified as follows: cDNA synthesis and denaturation: 1 cycle at 45°C for 30 min, 95°C for 1 min; PCR amplification, 5 cycles at 95°C for 15 s, 47°C for 30 s, 68°C for 3 min; 23 cycles at 95°C for 15 sec, 57°C for 30 sec, 68°C for 3 min; final extension 1 cycle 68°C for 7 min. The amplified genome was analyzed on a 1% agarose gel stained with SimplySafe dye (EURx).

2.5. Preparation of sequencing library.

PCR products were purified after enzymatic reaction on magnetic beads using Agencourt AMPure XP (Beckman Coulter) according to the manufacturer’s instructions, utilizing 96 deep well microtiter plates. The concentration of DNA in the samples after purification was measured using a Quantus™ Fluorometer (Promega) and the samples were diluted to a concentration of 0.2 ng/μl. Diluted samples were used to create a sequencing library using the Nextera XT DNA Library Prep Kit in 96 well plate. Briefly, genomic DNA in samples was tagmented; after addition of indexes, all libraries (each sample is a library with unique index combination) were amplified, amplified libraries were purified on Agencourt AMPure XP magnetic beads, and libraries were checked for right library size distribution on Agilent 4150 TapeStation System (Agilent Technologies). All libraries were normalized using a bead-based normalization method, pooled, and diluted to a loading concentration suitable for the MiniSeq System (Illumina).

2.6. Whole-genome sequencing using MiniSeq platform (Illumina), followed by bioinformatical and phylogenetic analysis of gathered data.

Whole-genome sequencing was performed on the MiniSeq platform from Illumina, and raw read data were basecalled and demultiplexed by the sequencing platform. Individual samples, represented as FASTQ files, each containing paired reads of 150 nucleotides in length, underwent bioinformatic analysis through two avenues: an automated process via the INSaFLU platform (https://doi.org/10.1186/s13073-018-0555-0), and an in-house pipeline. Qualitative verification was performed using FastQC (v.0.11.5), which served as a preliminary step for both methods. The automated process yielded the output sequences and a subtype-specific influenza virus report. The in-house pipeline employed tools integrated within Geneious Prime (v.2022) (Biomatters; http://www.geneious.com). For instance, BBDuk (v.38.84) was utilized for trimming, whereas normalization and error correction were performed using BBNorm (v.38.84). De novo assembly was accomplished through the assembly utility embedded within Geneious Prime, using a medium sensitivity option. The obtained contigs were subsequently compared using MAFFT (v.7), and, as necessary (such as a discrepancy at the nucleotide position or an incomplete contig sequence), reads were mapped to the consensus using Minimap2 (v.2.17), thereby establishing a final output sequence. The full sequences corresponding to the segments of the influenza virus were compared, establishing that there were five genetic variants present in all sequenced samples. Preliminary HA and NA subtype classifications were performed using multiplex PCR assays, as described in detail by Chiapponi et al. (2021) [22]. H1 clade classification was performed as described in Anderson et al. (2016) [27] using the swine H1 influenza classification tool available on http://www.fludb.org, accessed on 22 February 2022. For phylogenetic inference, the HA, NA, and PB2 sequences were extracted and aligned using MAFFT (v.7), following which maximum likelihood (ML) trees were calculated using RAxML (v.8.2). The correct substitution model and its parameters for inferred evolutionary history was chosen based on Bayesian Information Criterion and corrected Akaike Information Criterion using MEGAX. For phylogenetic analysis, the GTR GAMMA model was utilized, featuring rapid bootstrapping and a search for the best-scoring ML tree, validated by 1000 bootstrap replicates. For the segment-specific collection utilized for analyses, additional sequences were included, specifically relevant swIAVs representing major European lineages, human reference strains associated with “human-like” swIAVs and recent Polish isolates. Sequences were retrieved from Genebank InfluenzaVirus Resource Database (http://www.ncbi.nlthe m.nih.gov/genomes/FLU/Database/nph-select.cgi?go=database), and EpiFlu™ Database (http://www.gisiaid.org) accessed on 4 of May 2023. The obtained data were used for genotyping according to the method described by Watson et al. (2015) [25]. Sequences of all eight genome segments were used for phylogenetic inference. The origin of each analyzed segment was named by clustering with the reference strains. A list of reference strains together with accession numbers can be found in Supplementary Table S1.

3. Results

A total of 376 samples were collected from 11 farms. All farms were located in the northwestern region of Poland. Sampling was conducted between May 2017 and August 2019. RNA was isolated from all samples for subsequent confirmation of the presence of influenza A virus using a virotype® Influenza A RT-PCR Kit. We were able to confirm influenza A virus infection in 112 of 376 (29.8%) samples collected from 5 of the 11 facilities (45.5%). The infection prevalence in the farms ranged from 42.8% to 55%, with the exception of one farm, where the prevalence level was 80%. To isolate the virus, all of the influenza-positive samples served as an inoculum for MDCK cells, which were then cultured to check whether CPE would occur. We observed a strong CPE in 92 cultures. In the next four cultures, CPE was observed after the first blind passage, and 16 cultures showed no signs of CPE. The virus was harvested from all 96 CPE-positive samples. Viral RNA was extracted from all obtained isolates and used for whole genome amplification of influenza A virus. We successfully amplified all eight segments of the viral genome from 93 isolates. The PCR products were then used to create a set of sequencing libraries, and 93 samples were sequenced using a MiniSeq NGS sequencing system from Illumina. After raw data analysis, we assembled 93 sets of influenza A virus full genome segment sequences. Using the INSaFLU algorithm, we assigned a subtype to all 93 isolates. Thirty-four isolates were recognized as the A/H1N1 subtype and 59 as the A/H1N2 subtype. Isolates were further identified and classified as described by Chiapponi et al. (2021) [22]. Briefly, the HA and NA genes from each of the sequenced isolates were assigned to one of the HA lineage or NA subtype, based on PCR analysis using linage and subtype specific primers. Multiplex subtyping revealed that the HA gene from 59 of the isolates belonged to the H1-1C lineage, and all 59 corresponding NA genes were N2 subtypes. The remaining 34 isolates were classified as H1-1A and N1 subtypes, which matched in full with the initial INSaFLU algorithm subtyping. Due to the high variability of HA segments between circulating strains, we decided to use a phylogeny-based global nomenclature system and an automated annotation tool for H1 hemagglutinin genes from swine influenza A viruses available online in the Influenza Research Database. The H1 clade classification showed that 59 isolates classified as HA-1C by multiplex PCR belonged to the 1.C.2 clade (63.4%) was related to the Eurasian avian lineage, and the hemagglutinin from the rest of the isolates could be assigned to clade 1.A.3.3.2. (36.6 %) of H1N1pdm09 origin. A summary of the sample collection process together with clade and lineage classifications for HA and NA subtypes are described in Table 1.

After multiple alignments of all segments, we found five different strains (I-V) representing all 93 isolates (three among H1N2 isolates and two in the group of H1N1 isolates), which we chose for further analysis. The names of the strains corresponding to the analyzed genetic variants and accession numbers to the full genome sequences submitted to the GISAID EpiFlu database (http://www.gisiaid.org) are listed in Table 2.

The phylogeny based on HA, NA, and PB2 segments revealed the relationship of the strains analyzed in this study to relevant swIAVs representing major lineages circulating in Europe and recent Polish isolates. As shown in Figure 1, three H1N2 strains clustered with swIAVs, which harbor the HA gene derived from the Eurasian avian lineage, while two H1N1 strains were closely related to human pandemic strains and swIAVs with the HA gene of pandemic origin. NA gene phylogeny (Figure 2) shows a clear division into H1N1 strains with A/California/07/2009 related clusters and H1N2 strains grouped with the A/swine/Gent/1/1984-like strains. When analyzing Figure 3, we can see that the PB2 segment from all analyzed strains originated from the H1N1pdm09 line. In all phylogenetic trees, the analyzed strains were located in the vicinity of recent Polish isolates classified into the corresponding genotypes.

To study the possible reassortment events, we applied the genotype classification system developed by Watson et al. (2015) [25]. The origin of internal gene cassette (IGC) was determined and results showed that ICG’s from all strains ale closely related to A/California/07/2009, therefore, H1N1 strains from this study can be assigned to a P genotype and all H1N2 strains are classified as a genotype T (Table 3). Samples represented by strains A/swine/Poland/GIV-S21/2018(H1N1) and A/swine/Poland/GV-S23/2019(H1N1) were collected from the same farm 9 months apart and showed evidence of circulation of the swIAV of the same genotype but genetically different. The three H1N2 genetic variants, A/swine/Poland/GI-S01/2019(H1N2), A/swine/Poland/GII-S02/2017(H1N2), and A/swine/Poland/GIII-S38/2019(H1N2), represent different sampling locations and collection times

4. Discussion

It is estimated that seasonal flu epidemics cause 2–5 million cases and 250,000–500,000 deaths worldwide annually. The emergence of a new strain of influenza A virus in 2009 resulted in a pandemic that lasted over two years and was responsible for more than 60 million cases in the United States of America alone [8]. These figures show that newly emerging strains of influenza A virus pose a serious threat to public health. It is sufficient to mention human cases of highly pathogenic H5N1 and H7N9 avian influenza in Asia and Europe or swine origin H3N2v outbreaks in North America [37,38,39,40]. The lack of resistance to such strains in the human population is caused by continuous genetic variation in the virus, which is based on phenomena such as reassortment, antigenic shift, or antigenic drift [41]. Prophylactic vaccination against currently circulating human strains shows high efficiency, but due to the variability of surface antigens, hemagglutinin (HA), and neuraminidase (NA), the treatment must be repeated every year with a new vaccine, the composition of which is determined on the basis of antigenic drift. Furthermore, in the case of the emergence of an entirely new strain resulting from reassortment, seasonal vaccines will most likely prove ineffective [42].

Reassortment of genetic segments during coinfection of one organism with strains from different hosts could serve as a potential source of new, dangerous strains of influenza A virus and may lead to overcoming species barriers and adaptation of the virus to new hosts. The majority of avian influenza viruses originating from natural reservoirs cannot cause serious diseases in humans. However, some zoonotic strains that acquired the possibility of infection in the human population caused serious disease and, in some cases, fatalities [43,44,45,46].

During the pandemic, the new A (H1N1) pdm09 strain was also detected in animals, alarming farmers, and organizations involved in public and animal health, pushing them to identify the origin of the pandemic strain [47]. It was found that it was a reasortant combined with four different strains of the virus. It contained a combination of genes from "human" influenza virus, "swine" influenza virus from North America, Asian "swine" influenza virus of avian origin and the "avian" influenza virus from North America. This reassortment of influenza virus strains has never been recorded previously [48]. In the short time the new pandemic strain was responsible for majority of influenza cases worldwide. Before the beginning of the post-pandemic period, the seasonal H1N1 circulating among humans was almost entirely displaced. Simultaneously, A(H1N1)pdm09 infected swine and spread to pig populations in Europe [49,50]. There has been evidence of specific strain circulation between humans and swine, as a high degree of homology was observed between the sequences of isolates collected during local epidemics from both hosts [51,52,53]. Infection of pigs with "human" virus poses a serious risk of the emergence of a new virus strain capable of efficient transmission between humans [3]. Most likely, this scenario led to an outbreak of H3N2v in the U.S. in 2011 and 2012. The newly emerged strain was a swine reasortant that acquired gene segments from the swine H3N2 triple reasortant virus, circulating in the U.S in 1998, and an M gene from pandemic H1N1 [54]. Previous contact with pigs was the main etiology of the infections, but some cases of human-to-human transmission have also been reported, probably due to the presence of the M gene from the pandemic strain. [55]. A similar event occurred independently in Europe in 2012. A triple reasortant A/swine/Spain/28778/2012 was isolated in Spain, and similar to the North American strain, it contained external glycoprotein genes from H3N2, in this case from the Gent/84-like strain, an almost complete IGC from the Eurasian avian lineage, and the M gene from the H1N1pdm09-like strain. This shows how the independent evolution of the virus on two continents can lead to the emergence of a strain with similar epidemiological features [25]. Since the end of the 2009 pandemic, cases of human infection with swine-origin influenza virus have been reported in numerous locations in Europe and Asia, some of which are fatal [28,30,36,56,57,58,59]. Conversely, reverse zoonosis cases have been known for a long time. Some of them are well documented as pivotal changes in swine influenza diversity [60]; for example, the introduction of H3N2 “Port Chalmer’s like” virus into swine populations in Europe or the sustained presence of the A(H1N1)pdm09 virus among herds in the US and Europe [51,61,62]. Furthermore, there is some evidence supporting the fact that human-to-swine transmission frequency exceeds the reverse scenario, and the human-swine interface must be perceived as bidirectional to fully understand complex influenza A ecology and evolution [36].

The initiative to assess the epidemiological situation among swine influenza viruses in Europe gained momentum after the H1N1pdm09 lineage became well established in the pig population. Together with evolving sequencing and bioinformatics methods, new insights into the origin of circulating strains were made possible. Moreno et al. (2013) [63] and Trebbien et al. (2013) [19] show different evolutionary trends and reasortation events in Italy and Denmark when analyzing samples from 2000-2012 period. Although focused on a specific country and with a small number of samples, they gave a good perspective on the epidemiological situation of SIV in the pre-pandemic years [19,63]. To date, the largest and most diverse surveillance programs have been carried out by the European Surveillance Network for Influenza in Pigs (ESNIP-3) and Watson et al. (2015) [25], followed by a study described by Henritzi et al. (2020) [26]. The first program analyzed 290 isolates collected from14 European countries between 2009-2013. They identified 23 distinct genotypes (A-W) and proposed a classification system to distinguish between the main lineages and reasortant strains. At the animal level, H1avN1 (genotype A) was the most common strain (29%), followed by H3N2 (genotype B) (13%), H1N1pdm09 (genotype P) (12%), and H1huN2 (genotype C) (9%). Evidence showed a significant level of replacement of the endemic Eurasian avian viruses with H1N1pdm09-like strains and a high number of reasortation events, leading to high genotype diversity among European countries [25]. The latter program was based on over 18,000 samples collected from approximately 2,500 farms in 17 countries between 2015 and 2018. The study showed that 56.6% of swine holdings were infected with viruses belonging to the four major lineages. Twelve distinct HA/NA combinations were added to those reported by Watson et al. (2015) [25]. At the farm level, 39.2% of isolates belonged to the H1avN1av subtype, which is the most dominant subtype. The H1avN2 reasortant had a prevalence of 12.6% with H1huN2 at almost the same level (11.4%), and H1N1pdm09 accounted for 4.2% of cases. H3N2 subtype was more rarely detected (3.2%) and H3 strains also shown the lowest antigenic diversity, especially in comparison with H1pdm subtypes showing complex phylogenetic patterns and high evolution dynamics [26].

Our previous study carried out on virus isolates obtained from pigs in Poland revealed coinfection of animals with several different genetic variants of influenza A virus, some associated with human influenza A strains, most likely A/California/07/2009-like strains [64]. There is evidence that the number of pigs in Poland, which are prone to the risk of mixed infections, may constitute a reservoir of potentially dangerous strains of influenza A virus, resulting from genetic reassortment. Our findings were supported by study described in Czyzewska-Dors et al. (2017) [32], where over 5,900 samples from 145 farms in Poland were analyzed to assess the seroprevalence of swIAV in farrow-to-finish pig herds and in almost 18% of herds antibodies against all four major lineages (H1N1, H1N2, H3N2 and H1N1pdm09) were detected [32]

The data presented here are part of a study in which we analyzed over 370 samples taken from 11 swine herds in northwestern Poland. We were able to confirm the infection in 112 samples from 5 of the 11 herds, and we obtained 93 full genome sequences. All strains belonged to either A/H1N2 (63.4%) or A/H1N1(36.6%) swine influenza subtypes. Molecular subtyping and phylogenetic analysis showed that H1N2 strains found in three herds (27.3%) represented the T genotype, a reasortant strain carrying the H1av gene, N2 gene related to A/swine/Gent/1/1984, and IGC of pandemic origin. All H1N1 strains found in one herd (9.1%) represented the P genotype, where all the genetic segments were derived from H1N1pdm09. Multiple alignment of all segment sequences revealed five distinct genetic variants of H1N2 and H1N1 isolates (three and two, respectively).

The general number of infections in farms (45.5%) and animals (29.8%) was very similar to the latest large-scale European survey conducted by Henritzi et al. (2020) [26] (56.6% and 30.5% respectively) although when we look at particular genotype prevalence, surprisingly we did not find any H1avN1 subtype, dominant in Europe in recent years. In contrast, the incidence of H1avN2 in our study was almost the same (27.3% to 25%), and H1N1pdm09 was slightly lower (9.1% to12.5%). Taking into account the fact, that IGC’s from all strains described in this study are of pandemic origin, the lower incidence of H1N1pdm09 may be the result of the high rate of external glycoproteins segments exchange during reassortment events. The most recent data concerning the Polish swIAV are those described by Vereecke et al. (2023) [65]. Nineteen isolates from 14 Polish farms collected between 2017-2020 were analyzed and classified as genotypes T (58.3%), U (25%), P (16.7%), and A (7.1%). The incidence levels of particular subtypes differ from our findings, but genotypes T and P also constitute the majority of the analyzed strains, and the genotype T/P ratio is comparable (33.4% to 28.6%). We must also consider that these two studies were conducted on a small number of isolates, and without further investigation and country-wide sampling on a much larger scale, we cannot properly assess the genetic diversity of swIAVs circulating among herds in Poland; therefore, there is a need for a broad surveillance program coordinated by veterinary services.

5. Conclusions

The pandemic potential of zoonotic strains cannot be underestimated, public health measures must be implemented to prevent and manage infections in human and animal ecosystems, particularly in farm animals.

There are very few epidemiological studies involving full genetic characterization of strains and detection of co-infections with influenza A virus in Poland, that has been carried out, among livestock population so far. In the interest of public health, we should assess the risk within the human-animal interface, as it may serve as a reservoir of potentially dangerous, new genetic variants of the influenza A virus.