Molecular subtyping of European swine influenza viruses and scaling to high-throughput analysis
Virology Journal volume 15, Article number: 7 (2018)
Swine influenza is a respiratory infection of pigs that may have a significant economic impact in affected herds and pose a threat to the human population since swine influenza A viruses (swIAVs) are zoonotic pathogens. Due to the increasing genetic diversity of swIAVs and because novel reassortants or variants may become enzootic or have zoonotic implications, surveillance is strongly encouraged. Therefore, diagnostic tests and advanced technologies able to identify the circulating strains rapidly are critically important.
Several reverse transcription real-time PCR assays (RT-qPCRs) were developed to subtype European swIAVs in clinical samples previously identified as containing IAV genome. The RT-qPCRs aimed to discriminate HA genes of four H1 genetic lineages (H1av, H1hu, H1huΔ146–147, H1pdm) and one H3 lineage, and NA genes of two N1 lineages (N1, N1pdm) and one N2 lineage. After individual validation, each RT-qPCR was adapted to high-throughput analyses in parallel to the amplification of the IAV M gene (target for IAV detection) and the β-actin gene (as an internal control), in order to test the ten target genes simultaneously on a large number of clinical samples, using low volumes of reagents and RNA extracts.
The RT-qPCRs dedicated to IAV molecular subtyping enabled the identification of swIAVs from the four viral subtypes that are known to be enzootic in European pigs, i.e. H1avN1, H1huN2, H3N2 and H1N1pdm. They also made it possible to discriminate a new antigenic variant (H1huN2Δ146–147) among H1huN2 viruses, as well as reassortant viruses, such as H1huN1 or H1avN2 for example, and virus mixtures. These PCR techniques exhibited a gain in sensitivity as compared to end-point RT-PCRs, enabling the characterization of biological samples with low genetic loads, with considerable time saving. Adaptation to high-throughput analyses appeared effective, both in terms of specificity and sensitivity. This new development opens novel perspectives in diagnostic capacities that could be very useful for swIAV surveillance and large-scale epidemiological studies.
Swine influenza (SI) is a highly contagious viral respiratory infection of pigs that has become enzootic in areas densely populated with this species . It is responsible for significant economic losses in affected herds due to morbidity, stunted growth and secondary infections. Swine influenza A viruses (swIAVs) are also of public health concern due to their zoonotic potential. Although human infections remain subclinical or commonly produce only mild symptoms, swIAVs may be responsible for fatal cases and/or contribute, e.g. by reassortment, to pandemics, as illustrated by the pandemic that occurred in 2009 [2, 3].
swIAVs are enveloped viruses belonging to the genus Influenzavirus A, family Orthomyxoviridae. Their genome consists of eight segments of single-stranded RNA of negative polarity. They are typed according to their major antigenic determinants, i.e. the surface glycoproteins hemagglutinin (HA) and neuraminidase (NA) . SwIAVs can evolve by antigenic drift (nucleotide point mutations, substitutions, deletions or insertions in HA and NA) and/or by antigenic shift (HA and/or NA reassortments) [5, 6]. Three subtypes of swIAVs, H1N1, H3N2 and H1N2, are simultaneously circulating among pigs worldwide, whereas lineages may vary within each subtype depending on the region, i.e. North America, Europe and Asia .
In Europe, the predominant H1N1 viruses are entirely of avian origin, introduced from wild ducks to pigs in 1979 and referred as “avian-like” swine H1N1 (H1avN1) [8, 9], with HA genes belonging to the clade 1C defined recently . H3N2 viruses became widespread in European pigs in the mid-1980s. They are “human-like reassortant swine H3N2” viruses that combine the HA and NA genes from a descendant of the 1968 Hong Kong pandemic H3N2 virus and the internal genes of the H1avN1 virus . The dominant H1N2 viruses are “human-like reassortant swine H1N2” (H1huN2) viruses that emerged in 1994. They retained the genotype of European H3N2 swIAVs, except for the HA gene that they acquired from a human H1N1 virus from the 1980s  and which belongs to clade 1B . Thus, these three swIAV lineages share common internal genes, but have clearly distinguishable HA- and NA- encoded genes . Since 2009, the H1N1 pandemic (H1N1pdm) virus has been transmitted from humans to pigs and has become a fourth enzootic swIAV lineage in European pigs . The H1N1pdm virus originates from a reassortment between a North American triple reassortant swIAV and a Eurasian H1avN1 swIAV . The H1pdm gene (clade 1A ) and the N1pdm gene, can be distinguished genetically from H1 and N1 genes encountered in other European swIAVs . Due to co-circulation and co-infections, reassortant viruses originating from these four European enzootic lineages are also sporadically but regularly detected in several European countries, i.e. H1avN2, H1huN1, H1pdmN2, and some of them have even become enzootic in certain regions [8, 14, 15]. Moreover, a new variant resulting from an antigenic drift in the HA gene of the enzootic H1huN2 virus has been identified on several occasions over the last few years in Italy [16, 17] and France, where it was named H1huN2Δ146–147 .
Due to the increasing genetic diversity of swIAVs in Europe and because novel reassortants or variants may become enzootic in swine or have zoonotic implications, surveillance is strongly encouraged . Therefore, diagnostic tests and advanced technologies able to identify the circulating strains rapidly are critically important. The development and harmonization of diagnostic tools for detection (M gene) and subtyping (HA and NA genes) of European swIAVs have been initiated in some European countries, especially through the concerted action “European surveillance network for influenza in pigs (ESNIP) 3” [8, 19]. One issue, following a multi-center study using a panel of reference strains from several countries, was the selection of the best-performing primers/probe sets and amplification protocols that could be used in reverse transcription real-time PCR assays (RT-qPCRs) for the detection of HA and NA genes from the different genetic lineages . In this paper, we describe RT-qPCRs aimed specifically at identifying HA and NA genes from European swIAVs, taking into account local strains isolated recently in France. Following an original approach, we also report their adaptation to high-throughput analyses, in parallel to the amplification of the IAV M gene and the β-actin gene, which provide a comprehensive tool for detection and subtyping of European swIAVs in large biological sample sets.
Clinical samples and virus strains
Nasal swabs (MW950Sent2mL Virocult®, Kitvia, Labarthe-Inard, France) were collected from herds in France by veterinary practitioners or ANSES personnel, in the context of virological diagnosis passive surveillance programs or specific epidemiological investigations. They were taken from pigs during outbreaks of acute respiratory disease. The swabs were mixed vigorously and supernatants were stored at −70 °C until analysis. Some lung samples were also obtained from necropsy of fatal cases. Total viral RNAs were extracted from 200 μL of nasal swab supernatants or 20–30 mg of lung homogenate using the NucleoSpin RNA or Nucleospin 8 RNA kits (Macherey-Nagel, Hoerdt, France) or the RNeasy Mini Kit (Qiagen, Courtaboeuf, France). IAV genome detection was performed by M gene RT-qPCR using one of the two ready-to-use commercial kits previously validated by the French National Reference Laboratory (NRL), i.e. the LSI VetMAX™ Swine Influenza A-A/H1N1/2009- included kit (Life technologies, Carlsbad, CA, USA) or the ADIAVET™ SIV REALTIME kit (Bio-X Diagnostics, Rochefort, Belgique). Both steps, i.e. viral RNA purification and IAV genome detection, followed the instructions provided by RT-qPCR manufacturers .
Virus isolation was attempted from positive clinical samples in Madin Darby canine kidney (MDCK) cells, according to standard procedures . Reference strains, previously selected as representative of the main viral subtypes and lineages encountered in European pigs over the last few years, were also retrieved from the collection of the French NRL or kindly provided by the European surveillance network for influenza in pigs (ESNIP) [8, 23]. A/California/04/2009 (H1N1pdm) was obtained from the French National Reference Center for Influenza, Institut Pasteur, Paris.
Initially, swIAV-positive samples and/or virus isolates were subjected to molecular subtyping using conventional RT-PCR assays, as previously described [19, 23]. One multiplex RT-PCR enabled the specific detection of HA genes of the H1av (clade 1C ), H1hu (clade 1B.1 ) and H3 lineages, while another enabled the amplification of NA genes of the N1 and N2 lineages. Since 2010, H1 and N1 genes from H1N1pdm have been screened using ready-to-use commercial kits previously validated by the French NRL for Swine Influenza , according to the manufacturer’s instructions. The H1pdm gene (clade 1A.3.3.2 ) was detected using the real-time RT-PCR LSI VetMAX™ Swine Influenza A-A/H1N1/2009-H1 detection kit (Life technologies, Carlsbad, CA, USA) or the ADIAVET™ A/H1N1(2009) REALTIME kit (Bio-X Diagnostics, Rochefort, Belgique). The N1pdm gene was detected using the real-time RT-PCR LSI VetMAX™ Swine Influenza A-A/H1N1/2009-N1 detection kit (Life technologies, Carlsbad, CA, USA).
Some virus strains were further propagated in 9-day-old specific pathogen-free (SPF) embryonated chicken eggs and subjected to antigenic characterization using hemagglutination inhibition (HI) tests, according to standard procedures . Allantoic fluids were tested with hyperimmune sera produced in SPF pigs against reference strains representative of the different subtypes: A/Sw/Cotes d’Armor/0388/2009 (H1avN1), A/Sw/Flandres/1/1998 (H3N2), A/sw/Scotland/410440/94 (H1huN2), A/Sw/Cotes d’Armor/113/2006 (H1huN2), A/Sw/Cotes d’Armor/0070/2010 (H1huN1), A/Sw/Cotes d’Armor/0186/2010 (H1avN2), A/Sw/Sarthe/0255/2010 (H1N1pdm) and A/Sw/France/22–130212/2013 (H1huN2Δ146–147) [23, 24].
The subtype and lineage of some selected swIAVs and some positive field samples were determined and/or confirmed by sequencing the HA and NA gene segments. The viral genome was reverse transcribed using universal primers  and SuperScript II Reverse Transcriptase (Life technologies, Carlsbad, CA, USA) following the manufacturer’s instructions. Full length HA and NA genes were amplified using Platinum Taq DNA Polymerase High Fidelity (Life technologies, Carlsbad, CA, USA) and in-house designed primers (sequences available on request). Amplified products were then separated in agarose gel and purified with the QIAquick Gel Extraction Kit (Qiagen, Hilden, Germany). Both strands of the amplicons were sequenced with the same primers used for the amplification, as well as with additional internal primers (sequences available on request). Sequencing reactions were performed using the BigDye Terminator v3.1 Cycle Sequencing Kit (Life technologies, Carlsbad, CA, USA) on an automatic DNA sequencer ABI 3130 Genetic Analyzer (Life technologies, Carlsbad, CA, USA). Manual editing of sequences and assembly refinement were done with Vector NTI Advance 11.0 software (Life technologies, Carlsbad, CA, USA). The Influenza Research Database (IRD) was screened with BLASTN2 to identify sequences closely related to the HA and NA sequences.
Primers and probes for HA and NA RT-qPCRs
Sets of primers (forward and reverse) and the probe were designed to specifically amplify the different HA and NA genes from European swIAVs, on the basis of alignments of nucleotide sequences of viruses isolated since 2000 and retrieved from the IRD, GISAID (Global Initiative on Sharing All Influenza Data) or Genbank databases. Melting temperatures and/or basic properties of oligonucleotides were approximated using OligoCalc  or Vector NTI Advance® sequence analysis and design software (Thermo Fischer Scientific, Waltham, MA, USA) (Table 1). Whereas ready-to-use commercial kits have previously been validated for H1pdm and N1pdm detection , new sets of primers and the probe specific to these two genes were designed in our study in order to include H1pdm and N1pdm RT-qPCRs in the high-throughput PCR array. Thus, the H1av, H1hu, H1pdm and H3 sets specifically amplify HA genes from the European H1avN1, H1huN2, H1N1pdm and H3N2 lineages, respectively. The N1 set matches with NA genes from both the H1avN1 and the H1N1pdm lineages whereas the N1pdm set specifically amplifies the NA gene from the H1N1pdm lineage. The N2 primers and probe can detect NA genes from both the H1huN2 and the H3N2 lineages. In addition, an H1huΔ146–147 set was manually designed to identify the novel antigenic variant that have emerged among H1huN2 viruses following antigenic drift . Alignments comprised full length HA sequences of H1huN2 and H1huN2Δ146–147, as well as various H1huN1 reassortant viruses. Thus, the H1huΔ146–147 set design took into account a deletion of 6 nucleotides and a mutation located in the receptor binding site of the H1huN2Δ146–147 variant (manuscript in preparation). Finally, M and β-actin sets were included, since in-house RT-qPCRs for swIAV and reference gene detection were used, respectively, as controls in high-throughput PCR assays. All RT-qPCRs used standard TaqMan DNA probes except those for H1pdm, N1pdm and H1huΔ146–147 which used highly specific DNA probes with conjugated minor groove binder (MGB) groups at the 3‘-end. The fluorophore covalently attached to the 5’-end of the oligonucleotide probe was 6-carboxyfluorescein (FAM), hexachlorofluorescein (HEX) or cyanine5 (Cy5) (Table 1). All standard TaqMan probes used black hole quencher (BHQ1) at the 3′-end. Primers and labelled standard TaqMan probes were purchased from Sigma-Aldrich (Saint-Louis, MO, USA). MGB-labelled probes were from Life technologies (Carlsbad, CA, USA).
Real-time RT-PCR assays for European swIAV subtyping
Similarly to virus detection by M gene RT-qPCR, the HA/NA subtyping methods were assessed on RNA extracted from amplified viruses or clinical samples using the NucleoSpin RNA kit (Macherey-Nagel, Hoerdt, France) or the RNeasy Mini Kit (Qiagen, Courtaboeuf, France) following the manufacturer’s instructions. Real-time RT-PCR assays were run either on Chromo4 (Bio-Rad Laboratories, Hercules, CA, USA) or MX3005P (Stratagene, Agilent Technologies, La Jolla, CA, USA). Each RT-qPCR mixture contained 5 μL of RNA extract and 20 μL of master mix containing 1X GoTaq Probe qPCR Master Mix, dUTP (2X) (Promega, Fitchburg, WI, USA), 2X GoScript RT Mix for 1-step RT-qPCR (50X) (Promega, Fitchburg, WI, USA) as well as primers and probe in different concentrations depending on the target gene. Thus, the master mix contained 800 nM of each primer and 100 nM of probe for H1av, H3, N1 and N2, 800 nM of each primer and 140 nM of probe for H1hu, 400 nM of each primer and 500 nM of probe for H1huΔ146–147, and 800 nM of each primer and 250 nM of probe for H1pdm and N1pdm. All RT-qPCRs were run as singleplex assays, except N1 and N2 that were run as a duplex assay. The cycling conditions used for H1av, H1hu, H3 and N1/N2 were: 15 min reverse transcription at 45 °C, 2 min of denaturation at 95 °C and 42 cycles of 15 s at 95 °C and 1 min at 56 °C. The program used for H1huΔ146–147, H1pdm and N1pdm (with MGB probe) was 15 min reverse transcription at 45 °C, 2 min of denaturation at 95 °C, and 40 cycles of 15 s at 95 °C and 1 min at 60 °C.
Evaluation of real-time RT-PCR subtyping performance
Several panels incorporating virus strains and/or clinical samples were constituted to evaluate the performance of the RT-qPCRs developed for the amplification of the H1av, H1hu, H1huΔ146–147, H3, N1 and N2 genes.
A first panel (panel 1) of 42 swIAV isolates was set up to assess the analytical specificity of the RT-qPCRs designed for H1av, H1hu, H3, N1 and N2 gene amplifications. Among these virus strains, 32 belong to one of the four European enzootic lineages, i.e. 12 H1avN1, 10 H1huN2, 8 H1N1pdm and 2 H3N2, whereas 10 others are isolates originating from reassortments between enzootic viruses, i.e. 5 H1huN1 and 5 H1avN2 (Table 2). RNA extracts were diluted in RNAse free-water to adjust the quantification cycle (Cq) value to 24–27 when tested with one of the M gene RT-qPCRs described above, in order to reduce the risk of cross-reactions between H1av and H1hu.
Panel 2 was assembled to test the specificity of the H1huΔ146–147 gene RT-qPCR, i.e. its ability to differentiate the H1 gene of the antigenic variant H1huN2Δ146–147 virus among HA genes previously identified as belonging to the H1hu lineage. This panel consisted of RNA extracts obtained from 46 swIAV strains previously identified as H1huN2 viruses (44/46) or H1huN1 viruses (2/46) (Additional file 1).
The analytical sensitivity of the assays was evaluated by testing 10-fold serial dilutions (from 10−2 to 10−9) of RNA extracted from reference strains (Table 3). The range of linearity, the coefficient of linear regression (R2) and the efficacy were then calculated to evaluate the performance (Table 4). The repeatability of each RT-qPCR was estimated by testing, in 5–6 independent assays, the RNA extract dilution immediately below the last dilution in which the target gene was detected by the assay, and calculating the inter-assay coefficient of variation (CV) (Table 5).
Panels 3 and 4 were set up to evaluate the diagnostic abilities of the various RT-qPCRs, i.e. their potential to identify HA and NA genes directly from clinical samples previously shown to contain swIAV genome (M gene RT-qPCR positive). Panel 3 included 60 RNA extracts obtained from 59 nasal swab supernatants and one lung tissue homogenate (Table 6). The swIAV genome was previously subtyped in 37 nasal swabs and the lung sample (11 samples were positive for H1avN1, 4 for H1avN2, 15 for H1huN2, 2 for H1huN1, 1 for H3N2, 4 for H1N1pdm, and 1 for both H1av and H1hu genes mixed with N2), whereas in 22 nasal swab supernatants the IAV genome was only partially subtyped using conventional RT-PCRs (18 samples tested positive for H1av but not for N1 or N2; 3 samples tested negative for any HA but were found positive for N2, and 1 sample was found to be negative for HA and NA). Panel 4 was specifically assembled to assess the diagnostic ability of the H1huΔ146–147 RT-qPCR. It was composed of 78 RNA extracts, i.e. 56 obtained from nasal swab supernatants and 22 from swIAV strains, that were previously found either positive for H1hu but not sequenced, or positive for N2 but in which HA was not identified, or in which neither HA nor NA were identified (see Additional file 2).
High-throughput real-time RT-PCR procedure
High-throughput real-time RT-PCR amplifications (M, β-actin, H1av, H1hu, H1huΔ146–147, H1pdm, H3, N1, N1pdm and N2) were conducted on a LightCycler®1536 real-time PCR system (Roche, Meylan, France). A Bravo automated liquid handling platform equipped with a chiller and a PlateLoc thermal microplate sealer (Agilent Technologies, La Jolla, CA, USA) was used as a microplate dispenser. Each RT-qPCR mixture contained 1 μL of RNA extract and 1 μL of master mix containing 1X RealTime ready DNA Probes Master (Roche, Meylan, France), 2X GoScript RT Mix for 1-step RT-qPCR (50X) (Promega, Fitchburg, WI, USA), 800 nM of each concerned primer, with the exception of the M gene forward primer that was included at a final concentration of 400 nM, and 250 nM of standard or MGB-labelled TaqMan probe. The one-step real-time RT-PCR program involved a 30 min reverse transcription of RNA at 45 °C, followed by a 2 min denaturation step at 95 °C, and lastly 40 cycles of 0 s (i.e. a pulse) at 95 °C and 1 min at 58 °C.
Evaluation of high-throughput RT-qPCR performance
Panel 5 was established to evaluate the analytical specificity of the real-time RT-PCRs (including new in-house H1pdm and N1pdm RT-qPCRs) when using the high-throughput RT-qPCR procedure (see Additional file 3). This panel included of RNA extracts obtained from 88 swIAVs and 30 nasal swab supernatants. The isolates as well as the viruses present in clinical samples were first subtyped using RT-qPCRs run under low-throughput conditions. Serial 10-fold dilutions (10−1 to 10−8) of RNA obtained from 8 reference strains were prepared in water to evaluate the analytical sensitivity of each RT-qPCR when run on the LightCycler®1536 (Table 7). The M gene copy number present in each dilution point was measured using a duplex M/β-actin gene RT-qPCR as described separately . Each dilution point was tested in duplicate.
Analytical performances of the real-time RT-PCRs for H1av, H1hu, H1huΔ146–147, H3, N1 and N2 gene amplification
The analytical performances of the novel RT-qPCRs developed for European swIAV subtyping in addition to the H1pdm and N1pdm commercial kits previously validated were first evaluated in low-throughput analyses.
Analyses of 42 swIAV isolates from panel 1 showed perfect concordance between the results from H1av, H1hu, H3, N1 and N2 RT-qPCRs and those expected from conventional, i.e. end-point RT-PCRs and/or sequencing data, indicating 100% analytical specificity for the assays (Table 2). Despite non negligible viral loads, i.e. Cq-values around 25 in M gene RT-qPCR, none of these assays showed any undesired amplification of genes from influenza viruses other than those expected. According to plan, the N1 RT-qPCR was able to detect NA genes from the H1avN1 and H1N1pdm lineages, while the N2 RT-qPCR amplified NA genes of viruses from the H1huN2 and H3N2 lineages.
Among the 46 swIAV strains included in panel 2 and previously identified as viruses that bear an H1hu gene, 25/44 H1huN2 strains were detected by the H1huΔ146–147 RT-qPCR indicating they were variant viruses, in accordance with HA gene sequencing (see Additional file 1). All other H1huN2 strains did not exhibit the 6 nt-deletion encountered in the variant and were not detected by specifically designed RT-qPCR H1huΔ146–147, showing that there were no false-positive results. By contrast, the two H1huN1 viruses were not detected even though they also showed the 146–147 amino acid deletion, due to a mutation in the MGB probe target region. Indeed, sequencing revealed that these viruses retained a thymine residue (T) in probe position 12, instead of having mutated into a cytosine residue (C) as compared to the parental H1hu gene. These reassortant viruses, 100% identical in their HA genes, were isolated from the same farm 2.5 months apart.
Real-time RT-PCRs developed for HA or NA subtyping of previously detected swIAVs exhibited lower sensitivity levels than the M gene assays used for the detection, except those amplifying the H1huΔ146–147 and the H3 genes (Table 3). Thus, assays conducted on 10-fold dilutions of reference strains indicated that the detection limits based on M gene RT-qPCR Cq-values would be >35 for H1huΔ146–147 and H3, around 33 for H1av, around 27 for N1, between 24 and 28 for N2 depending on the virus subtype, and around 24 for H1hu. Based on calculation of the ranges of linearity, it appeared that all genes were correctly detected from the first (10−1) strain dilution, i.e. at Cq-value <13.5, except the H3 gene for which the lowest limit was estimated at the second (10−2) dilution, at a Cq-value equivalent to 16.7 (Table 4). Between the limits of the range of linearity, the slopes of the standard curves varied from −3.3483 to −3.266. The coefficients of linear regression (R2) varied from 0.9957 to 0.9998 and the efficiencies were calculated between 96.81% and 104.16%. Repeatability and intermediate precision, assessed by calculating the inter-assay CVs, gave satisfactory results (Table 5). The inter-assay CVs were all <10%, ranging from 1.20 for the H1av RT-qPCR to 6.51 for the N2 RT-qPCR run on the H3N2 strain.
Diagnostic abilities of the RT-qPCR assays targeting the H1av, H1hu, H1huΔ146–147, H3, N1 and N2 genes of European swIAVs
Nasal swabs from animals naturally infected with H1avN1, H1huN2 or H3N2 swIAVs, as well as with H1avN2 or H1huN1 reassortant viruses (panel 3), were detected positive for both HA and NA corresponding genes in 100% of cases tested for both genes (Table 6). Additional samples tested either only in the simplex HA assays or in the N1/N2 multiplex assay were also qualified positive according to the expected HA or NA gene. Inclusion of 4 samples positive for H1N1pdm confirmed N1 amplification without any H1av gene amplification. Sample 130103–2, which was previously known to contain both the H1av and H1hu genes together with the N2 gene, was also proved here to contain the N1 gene, enabling us to confirm the hypothesis of a mixture of H1avN1 and H1huN2 viruses in this sample. Analyses of 18 samples previously partially subtyped as “H1avN?” confirmed the H1av subtype in 9/18 cases. Others did not give a positive signal, probably due to limiting genome amounts (7/11 were samples with M gene RT-qPCR Cq-values >32, while the H1av RT-qPCR detection limit was estimated around M gene Cq-value = 33). Conversely, the N1 gene was identified in 10/18 H1avN? samples, even in samples where M gene Cq-value >27, which was estimated as the detection limit. The N2 gene was not amplified in any of these samples. Panel 3 also comprised 3 H?N2 samples. The N2 gene was detected in all of them, whereas the H1hu gene was identified in 2/3. The last one (sample 130277–3) was shown to have an H1hu gene after propagation in cell culture (data not shown). Finally, a sample (130133–5) that was left fully unsubtyped using conventional RT-PCRs was here classified as containing the H1av, H1hu and N2 genes, confirming the ability to detect virus mixtures using these RT-qPCRs.
The diagnostic ability of the H1huΔ146–147 RT-qPCR was evaluated by testing 78 samples of mainly H1N2 or H?N2 subtypes included in panel 4. The H1huΔ146–147 gene was identified in 3/15 H1huN2 isolates (not sequenced previously), in 16/46 nasal swab supernatants in which the H1huN2 virus was identified, as well as in 1/6 H1huN? samples (Additional file 2). Although the specificity of the H1huΔ146–147 was previously demonstrated, these results were further validated by sequencing the HA genes of the three isolates that tested positive for the H1huΔ146–147 gene as well as the HA gene of strain A/Sw/France/29–150034/15 that tested negative. Sequencing confirmed that the three positive strains contained the 6 nt deletion by contrast to the other one that contained a H1hu gene without this deletion (Genbank accession numbers KJ128334, KY241154, KY241155 and KY241115). Thus, 20/78 viruses or clinical samples previously identified as H1hu-positive were rapidly identified as having an H1huΔ146–147 gene.
Analytical performances of the HA and NA RT-qPCRs in high-throughput analyses
Real-time RT-PCRs for H1av, H1hu, H1huΔ146–147, H3, N1 and N2 gene amplifications, as well as RT-qPCRs designed to amplify the H1pdm and N1pdm genes specifically, were evaluated in high-throughput analyses. Real-time RT-PCRs targeting the IAV M gene and the porcine β-actin genes were included in order to check the amount and quality of the genetic material and/or demonstrate any potential PCR inhibition. These last RT-qPCRs were previously validated in low-throughput analyses [21, 27]. Thus, 10 RT-qPCR assays were run together on a LightCycler®1536 to analyze the 118 samples from panel 5, i.e. 83 swIAVs, 5 virus mixtures and 30 nasal swab supernatants, all previously analyzed in low-throughput PCR assays (see Additional file 3). As expected, all samples were detected positive for the M gene. RNA extracts from virus strains exhibited M gene Cq-values ranging from <5 to 16.69, whereas extracts from clinical samples exhibited higher Cq-values, which mostly ranged between 14 and 30. Only 2/118 samples (1.7%) (both RNA extracts from MDCK cell culture supernatants) were found negative for the β-actin gene.
All HA and NA RT-qPCRs exhibited 100% specificity, as no undesirable cross-reaction or false-positive results were obtained, irrespective of the nature of the sample, i.e. RNA extracts from virus stock containing high amounts of genetic material, or RNA extracts from nasal swab supernatants. As expected, the H1av RT-qPCR amplified H1 genes originating from the H1avN1 lineage, i.e. H1 genes from H1avN1 enzootic viruses and from H1avN2 reassortants. The H1hu RT-PCR amplified H1 genes from H1huN2 viruses, H1huN2Δ146–147 variants and H1huN1 reassortants. By contrast, the H1huΔ146–147 RT-PCR only detected the H1huN2Δ146–147 antigenic variants. In accordance with the results obtained in low-throughput analyses, it did not detect H1huN1 reassortants bearing H1hu genes with one or two amino acid deletions close to the receptor-binding site (RBS) (called H1huN1Δ146–147 and H1huN1Δ147). The H1pdm RT-qPCR amplified HA genes of viruses from the H1N1pdm lineage but also any other H1 gene, and the H3 RT-qPCR amplified HA from H3N2 viruses only. The N1 RT-qPCR amplified N1 genes from both the H1avN1 and the H1N1pdm lineages, whereas the N1pdm RT-qPCR amplified the N1 genes from H1N1pdm viruses specifically. The N2 RT-qPCR amplified NA genes from all HxN2 viruses. Thus, all samples were fully subtyped according to expected results, including the virus mixtures, except one of them in which the N2 gene was not amplified in parallel to the H1av, H1hu and N1 genes, as expected. As observed for the M gene, the Cq-values obtained for the HA and NA RT-qPCRs were lower on RNAs extracted from virus isolates as compared to clinical samples. In several cases, the Cq-value was too low to be precisely determined (< 5). A sample without any viral RNA but rather consisting of water was added in each PCR array as a negative control and no unspecific detection was observed (data not shown).
In order to go further in evaluating the sensitivity of the methods when run as high-throughput analyses, each RT-qPCR (except the one targeting the β-actin gene) was tested against 10-fold serial dilutions (10−1 to 10−8) of reference strains (Table 7). Depending on the virus stock and its initial genomic load, the M gene was still detected to the last dilution or to the 10−6 or the 10−7 dilution. In all cases, it was detected to the 10−5 dilution. The H1av, H1huΔ146–147 and H1pdm genes were detected in samples exhibiting a corresponding M gene Cq-value until 32, approximately, while the H3 gene detection limit was found to exceed an M gene Cq-value >35. The detection limit for the H1hu gene varied from 25 to 33 depending on the virus (parental virus, antigenic variant, reassortant or reassortant with deletion), that of N1 from 29 to 31, and that of N2 from 30 to 33. The N1pdm was detected until an M gene Cq-value of 31. The ranges of linearity are given in Table 8. The slopes calculated over these ranges were all (except one) comprised between −3.2 and −3.7, which means the efficacies of the RT-PCRs varied from 85% to 105%. The N2 RT-qPCR run on the reference H3N2 strain was the only one found to be outside these limits, showing an efficacy of 80.51%.
Robustness of the RT-qPCR subtyping tool in low- and high-throughput analyses
The full RT-qPCR subtyping tool began to be used routinely by the French NRL for Swine Influenza from January 2014, as part of the analytical workflow for swIAV surveillance. Thus, nasal swab supernatants previously selected as containing M gene from a swIAV were subjected to H1av, H1hu and H3 monoplex assays, to the N1/N2 duplex assay, as well as to commercial H1pdm and N1pdm RT-PCRs. Samples that were found to be positive for H1hu were subjected to the H1huΔ146–147 RT-PCR in a second step. Considering the detection limits evaluated above, subtyping was undertaken on samples with M gene Cq-values <35 only. Samples exhibiting M gene Cq-values >35 were considered to be not typable. Looking at the proportions of HA and NA genes that were successfully subtyped in these samples from January 2014 to August 2016, it appeared that both genes were identified in almost all samples with M gene Cq-value <25 (Table 9). For 25 ≤ Cq-value <30, the proportion of unidentified HA and/or NA genes increased slightly. In these samples, the N1/N2 multiplex appeared less sensitive than HA simplex assays, as NA genes were less frequently amplified. In samples with Cq-value >30, the proportion of viruses fully or partially subtyped fell to about 50%.
Finally, the diagnostic ability of the high-throughput subtyping RT-qPCR assay was tested on RNAs extracted from 919 nasal swabs collected during a longitudinal study conducted in three pig herds located in Brittany and previously selected as M gene-positive samples . These herds were known to be affected by H1avN1 and H1huN2 viruses only. Among them, HA and NA genes were both identified in 697 samples, whereas only the HA or NA gene was detected in 77 others, depending on the apparent amount of virus genome in the samples (Table 9). Almost all samples with an M gene Cq-value <30 were fully subtyped (Table 9). In samples with an M gene Cq-value >30, the NA gene was more frequently amplified than the HA gene (Table 9). Fourteen samples from one herd were shown to contain virus mixtures and/or reassortant viruses (data not shown).
Eight real-time RT-PCRs were evaluated for European swIAV molecular subtyping in clinical samples demonstrated to contain IAV genome. They aimed to discriminate HA genes for four H1 genetic lineages (H1av, H1hu, H1huΔ146–147, H1pdm) and one H3 lineage, as well as NA genes of two N1 lineages (N1, N1pdm) and one N2 lineage. Altogether, the RT-qPCRs enabled the identification of swIAVs from the four viral subtypes known to be enzootic in European pigs, i.e. H1avN1, H1huN2, H3N2 and H1N1pdm. They also made it possible to quickly identify a new antigenic variant (H1huN2Δ146–147) among H1huN2 viruses, as well as reassortant viruses, i.e. H1huN1 or H1avN2, and virus mixtures. All assays were optimized to take into account the genetic diversity encountered among swIAVs isolated in France in recent years, whereas keeping in mind usefulness on an European level. They exhibited a gain in sensitivity as compared to conventional RT-PCRs, allowing the characterization of biological samples with low genetic loads, with considerable time saving. Reagents and amplification procedures were harmonized to run all assays in parallel, using a unique amplification cycle profile irrespective of the equipment.
Evaluation in low-throughput analyses showed perfect specificities for H1av, H1hu and H3 RT-qPCRs conducted in parallel as simplex assays. While multiplexing may facilitate the diagnostic procedure, preliminary studies showed that individual RT-qPCR sensitivities were somewhat affected when run in a triplex assay (data not shown). Based on M gene Cq-values, RT-qPCRs appeared slightly less sensitive than assays involving MGB-labelled probes, such as previously validated commercial kits or other in-house methods aimed at amplifying the M, H1pdm or N1pdm genes . Nevertheless, the detection levels were satisfactory as HA subtyping was nearly 100% successful in clinical samples up to M gene Cq-value = 30. The H1hu RT-qPCR was the least efficient, probably due to higher genetic variability among H1hu genes, as compared to the H1av and H3 genes . The H1huΔ146–147 RT-qPCR, run in a second step on H1huNy-positive samples also demonstrated excellent specificity and sensitivity, leading to rapid discrimination of novel antigenic variants among H1huN2 viruses. While H1huN1 reassortants are rare events in most European countries, they are sporadically detected in France [8, 29] and some of them may also have an H1huΔ146–147 gene. In this study, two H1huN1 reassortants exhibiting 2 amino acid deletions at positions 146–147 of the RBS, thus bearing H1 genes genetically and antigenically closer to the H1huΔ146–147 variant than to the parental H1hu, were not detected by the H1huΔ146–147 RT-PCR due to a mismatch within the “H1hu_var” MGB probe. As a result, classification of the H1hu gene from H1huN1 reassortants as a “Δ146–147 variant” (or not) would be better confirmed by HA sequencing.
N1 and N2 RT-qPCRs, adapted from protocols previously evaluated at the European level, confirmed the specificities expected from in silico analyses . While designed to amplify N2 genes from both European H1huN2 and H3N2 lineages, the N2 RT-PCR proved to detect N2 genes from the novel H1huN2Δ146–147 variants. By contrast to HA RT-qPCRs, they were run in a duplex assay because preliminary studies indicated no reduction in sensitivity as compared to corresponding simplex assays (data not shown). Nevertheless, it should be noted that in routine diagnosis, the NA gene was less frequently subtyped than the HA gene in samples with an M gene Cq-value >30, leading to the hypothesis that this difference could be related, at least partially, to multiplexing.
Scaling M and β-actin RT-qPCRs as well as HA (H1av, H1hu, H1huΔ146–147, H1pdm, H3) and NA (N1, N1pdm, N2) subtyping assays to higher throughput, with concurrent miniaturization of individual reactions, was successful. The LightCycler®1536 system designed for automated high-throughput laboratory workflows proved in this context to provide very high well-to-well homogeneity, good inter-assay reproducibility and low inter-plate variability. The 10 RT-qPCRs, run together onto one microplate using the same amplification procedure, showed very good efficiency. Each RT-qPCR retained its specificity. They all exhibited comparable and acceptable sensitivities. When run on clinical samples with an M gene Cq-value <30, the techniques were able to fully subtype more than 90% of detected swIAV genomes, showing an increased proportion of characterized NA genes in samples with Cq-values of 25–30 as compared to equivalent samples using classical thermocyclers. By contrast, when run on clinical samples with an M gene Cq-value >30, H1av, H1hu, N1 and N2 RT-PCRs appeared slightly less efficient than in low-throughput analyses. This could be related to the low volume of RNA extract included into the RT-qPCR mixture when run on LightCycler®1536 system as compared to other thermocyclers (1 μL instead of 5 μL), but this would need to be further investigated and confirmed for other subtyping RT-qPCRs.
The emergence of the pandemic A/H1N1 virus of swine origin in 2009 highlighted the need for global surveillance of influenza A viruses in pigs. Its subsequent introduction into the pig population, its co-circulation with other enzootic swIAVs, and genomic reassortment events have led to an increase in the genetic diversity of swIAVs. Altogether, these RT-qPCR assays provide a rapid and simple genotyping method to identify viruses that infect the European pig population, including a novel H1huN2 variant identified in France, sporadic HA/NA reassortants and virus mixtures, as a first characterization step before full genome sequencing and/or antigenic subtyping. Screening of individual samples against the 10 target genes in a high-throughput scenario opens novel perspectives in diagnostic abilities, which will be very useful for swIAV surveillance and large-scale epidemiological studies.
Van Reeth K, Brown IH, Olsen CW. Influenza Virus. In: Zimmerman JJ, KLA, Ramirez A, Schwartz KJ, Stevenson GW, editors. Diseases of Swine. Tenth Edition. Wiley; 2012. p. 557-71.
Dawood FS, Jain S, Finelli L, Shaw MW, Lindstrom S, Garten RJ, Gubareva LV, Xu X, Bridges CB, Uyeki TM. Emergence of a novel swine-origin influenza a (H1N1) virus in humans. New Eng J Med. 2009;360(25):2605–15.
Myers KP, Olsen CW, Gray GC. Cases of swine influenza in humans: a review of the literature. Clin Infect Dis. 2007;44(8):1084–8.
Cheung TK, Poon LL. Biology of influenza a virus. Ann N Y Acad Sci. 2007;1102:1–25.
Nelson MI, Holmes EC. The evolution of epidemic influenza. Nat Rev Genet. 2007;8(3):196–205.
Webster RG, Bean WJ, Gorman OT, Chambers TM, Kawaoka Y. Evolution and ecology of influenza a viruses. Microbiol Rev. 1992;56(1):152–79.
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(6–7):310–25.
Simon G, Larsen LE, Durrwald R, Foni E, Harder T, Van Reeth K, Markowska-Daniel I, Reid SM, 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(12):e115815.
Watson SJ, Langat P, Reid SM, Lam TT, 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(19):9920–31.
Anderson TK, Macken CA, Lewis NS, Scheuermann RH, Van Reeth K, Brown IH, Swenson SL, Simon G, Saito T, Berhane Y, et al. A phylogeny-based global nomenclature system and automated annotation tool for H1 hemagglutinin genes from swine influenza a viruses. mSphere. 2016;1(6):e00275-16.
Castrucci MR, Campitelli L, Ruggieri A, Barigazzi G, Sidoli L, Daniels R, Oxford JS, Donatelli I. Antigenic and sequence analysis of H3 influenza virus haemagglutinins from pigs in Italy. J Gen Virol. 1994;75(2):371–9.
Brown IH, Harris PA, McCauley JW, Alexander DJ. Multiple genetic reassortment of avian and human influenza a viruses in European pigs, resulting in the emergence of an H1N2 virus of novel genotype. J Gen Virol. 1998;79(12):2947–55.
Smith GJ, Vijaykrishna D, Bahl J, Lycett SJ, Worobey M, Pybus OG, Ma SK, Cheung CL, Raghwani J, Bhatt S, et al. Origins and evolutionary genomics of the 2009 swine-origin H1N1 influenza a epidemic. Nature. 2009;459(7250):1122–5.
Trebbien R, Bragstad K, Larsen LE, Nielsen J, Botner A, Heegaard PM, Fomsgaard A, Viuff B, Hjulsager CK. 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. Virology J. 2013;10:290.
Lange J, Groth M, Schlegel M, Krumbholz A, Wieczorek K, Ulrich R, Koppen S, Schulz K, Appl D, Selbitz H-J, et al. Reassortants of the pandemic (H1N1) 2009 virus and establishment of a novel porcine H1N2 influenza virus, lineage in Germany. Vet Microbiol. 2013;167(3–4):345–56.
Moreno A, Chiapponi C, Boniotti MB, Sozzi E, Foni E, Barbieri I, Zanoni MG, Faccini S, Lelli D, Cordioli P. Genomic characterization of H1N2 swine influenza viruses in Italy. Vet Microbiol. 2012;156(3):265–76.
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.
Garin E, Hervé S, Rose N, Locatelli C, Ngwa-Mbot D, Wendling S, Bournez L, Calavas D, Simon G. French network for the surveillance of influenza a viruses in pigs (Résavip) - results of the surveillance carried out in 2015. Bull Epid Santé Anim Alim. 2016;75:16–9.
Chiapponi C, Moreno A, Barbieri I, Merenda M, Foni E. Multiplex RT-PCR assay for differentiating European swine influenza virus subtypes H1N1, H1N2 and H3N2. J Virol Methods. 2012;184(1–2):117–20.
Henritzi D, Zhao N, Starick E, Simon G, Krog JS, Larsen LE, Reid SM, Brown IH, Chiapponi C, Foni E, et al. Rapid detection and subtyping of European swine influenza viruses in porcine clinical samples by haemagglutinin- and neuraminidase-specific tetra- and triplex real-time RT-PCRs. Influenza Other Resp Viruses. 2016;10(6):504–17.
Pol F, Quéguiner S, Gorin S, Deblanc C, Simon G. Validation of commercial real-time RT-PCR kits for detection of influenza a viruses in porcine samples and differentiation of pandemic (H1N1) 2009 virus in pigs. J Virol Methods. 2011;171:241–7.
OIE. Influenza A virus of swine. In: Manual of diagnostic tests and vaccines for terrestrial animals; 2015. Chapter 2.8.7. http://www.oie.int/fileadmin/Home/fr/Health_standards/tahm/2.08.07_INF_A_SWINE.pdf.
Kyriakis CS, Brown IH, Foni E, Kuntz-Simon G, Maldonado J, Madec F, Essen SC, Chiapponi C, Van Reeth K. Virological surveillance and preliminary antigenic characterization of influenza viruses in pigs in five European countries from 2006 to 2008. Zoonoses Public Health. 2011;58:93–101.
Lewis NS, Russell CA, Langat P, Anderson TK, Berger K, Bielejec F, Burke DF, Dudas G, Fonville JM, Fouchier RA, et al. The global antigenic diversity of swine influenza a viruses. elife. 2016;5:e12217.
Hoffmann E, Stech J, Guan Y, Webster RG, Perez DR. Universal primer set for the full-length amplification of all influenza a viruses. Arch Virol. 2001;146(12):2275–89.
Kibbe WA. OligoCalc: an online oligonucleotide properties calculator. Nucleic Acids Res. 2007;35:W43–6.
Deblanc C, Delgado-Ortega M, Gorin S, Berri M, Paboeuf F, Berthon P, Herrler G, Meurens F, Simon G. Mycoplasma hyopneumoniae does not affect the interferon-related anti-viral response but predisposes the pig to a higher inflammation following swine influenza virus infection. J Gen Virol. 2016;97:2501–15.
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(1):72.
Simon G, Hervé S, Rose N. Epidemiosurveillance of swine influenza in France from 2005 to 2012: programs, viruses and associated epidemiological data. Bull Epid Santé Anim Alim. 2013;56:17–22.
Weingartl HM, Berhane Y, Hisanaga T, Neufeld J, Kehler H, Emburry-Hyatt C, Hooper-McGreevy K, Kasloff S, Dalman B, Bystrom J, et al. Genetic and pathobiologic characterization of pandemic H1N1 2009 influenza viruses from a naturally infected swine herd. J Virol. 2010;84:2245–56.
Duvigneau JC, Hartl RT, Groiss S, Gemeiner M. Quantitative simultaneous multiplex real-time PCR for the detection of porcine cytokines. J Immunol Methods. 2005;306(1–2):16–27.
The authors would like to thank the members of the French national network for Influenza A surveillance in pigs (Résavip) and, more broadly, all practicing veterinarians who send swine samples for influenza A virus monitoring.
This study was partially funded by the European FP7 project ESNIP3 (#259949) and the CoVetLab project “Development and validation of molecular tools for sub-typing swine influenza viruses” (#CF0010).
Availability of data and materials
All relevant information is provided in this current manuscript.
Ethics approval and consent to participate
Consent for publication
None of the authors of this paper have a financial or personal relationship with other people or organizations that could inappropriately influence or bias the content of the paper. The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Specificity of the real-time RT-PCR developed to identify H1huΔ146–147 antigenic variants among H1huNY swIAVs (panel 2). (DOCX 18 kb)
Identification of new antigenic variants among H1huNY swIAVs by real-time RT-PCR targeting the H1huΔ146–147 gene. (DOCX 22 kb)
Specificity of real-time RT-PCRs for detection and subtyping of swIAVs when run simultaneously as simplex assays on LightCycler®1536. (DOCX 61 kb)
About this article
Cite this article
Bonin, E., Quéguiner, S., Woudstra, C. et al. Molecular subtyping of European swine influenza viruses and scaling to high-throughput analysis. Virol J 15, 7 (2018). https://doi.org/10.1186/s12985-018-0920-z
- High-throughput real-time RT-PCR