- Short report
- Open Access
Real-time reverse transcription PCR-based sequencing-independent pathotyping of Eurasian avian influenza A viruses of subtype H7
Virology Journalvolume 14, Article number: 137 (2017)
Low pathogenic avian influenza viruses (LPAIV) of the subtypes H5 and H7 are known to give rise to highly pathogenic (HP) phenotypes by spontaneous insertional mutations which convert a monobasic trypsin-sensitive endoproteolytical cleavage site (CS) within the hemagglutinin (HA) protein into a polybasic subtilisin-sensitive one. Sporadic outbreaks of notifiable LPAIV H7 infections are continuously recorded in Europe and in Asia, and some lineages showed zoonotic transmission. De novo generation of HPAIV H7 from LPAIV precursors has been reported several times over the past decade. Rapid differentiation between LP and HP H7 virus strains is required as a prerequisite to emplace appropriate control measures. Here, reverse transcription real-time PCR assays (RT-qPCR) were developed and evaluated that allow LP and HP pathotype identification and distinction by probe-assisted detection of the HACS. These new RT-qPCRs allow a sensitive and highly specific pathotype identification of Eurasian subtype H7 AIV in allantoic fluids as well as in diagnostic field samples. RT-qPCR assisted pathotyping presents a rapid and sensitive alternative to pathotyping by animal inoculation or nucleotide sequencing.
Avian influenza viruses (AIV) are members of the family Orthomyxoviridae, specified as influenza virus type-A. These viruses are further classified by the serologically defined subtypes of the predominant viral surface glycoproteins, the hemagglutinin (HA) and neuraminidase (NA) . Their genome is composed of single-stranded, negative-sense RNA and comprises eight genome segments which encode at least ten proteins . All 16 HA and nine NA AIV subtypes can be detected in populations of aquatic wild birds which form the natural reservoir of these viruses .
Based on their pathogenicity in chickens, two phenotypes of AIV are distinguished: highly pathogenic (HP) AIV and AIV of low pathogenicity (LPAIV). In nature, HP phenotypes have been restricted to viruses of subtypes H5 and H7. HPAIV arises from LPAI precursor viruses by spontaneous mutations leading to the insertion of basic amino acids into the cleavage site (CS) of the hemagglutinin protein (HA) which renders the HACS processible to subtilisin-like host proteases that are ubiquitous in all host tissues. Such viruses, therefore, gain competence for fatal systemic infections in avian hosts. LPAIV, in contrast, depends on local provision of trypsin-like proteases at the epithelial surfaces of the respiratory and/or gastrointestinal tracts and per se do not cause severe clinical signs . All LPAIV and HPAIV infections of subtypes H5 and H7 in poultry are notifiable to the World Organization for Animal Health (O.I.E.).  Determination of the type of HACS is of utmost importance for the diagnosis of these infections. This can be achieved biologically by determination of the intravenous pathogenicity index (IVPI) in experimentally inoculated chickens or molecularly by nucleotide sequence analysis of the site encoding the HACS . Since animal experiment facilities or expensive equipment are required for either pathway, solutions for alternative techniques have been sought in the past: These included restriction enzyme cleavage patterns , probe hybridization  and real time RT-PCR (RT-qPCR) approaches . Based on the widespread availability of RT-qPCR technology in diagnostic laboratories and its recent favorable use in pathotyping of HPAIV H5 of the goose/Guangdong (gs/GD) lineage , this study was conducted to develop and validate sequencing-independent RT-qPCRs for pathotyping of Eurasian H7 AI viruses.
Over the past two decades, several incursions into poultry of subtype H7 LPAIV as well as the de novo generation and (in one case) spread of H7 HPAI viruses have been reported from Europe (Table 1). Other H7 LPAIV lineages have arisen in Eastern Asia, and one of them (H7N9/China) showed significant zoonotic propensities in annual waves of poultry-to-human transmission with more than 550 fatal human cases [11, 12]. Recently, the H7N9 lineages has also yielded an HP mutant which is spreading in southern China . Considering the annual presence of LPAIV of subtype H7 in Eurasian wild bird populations  risks of new incursions into poultry in Europe are perpetuating.
Based on the alignments of the HA H7 gene of a comprehensive selection of sequences from LP (n = 60) and HPAIV (n = 21) of Eurasian origin collected over the last decade in sequence databases (GenBank at NCBI; EpiFlu of the Global Initiative on Sharing Avian Influenza Data (GISAID)), a set of six primers was designed (Table 2). The selected primers targeted a short fragment of the HA gene that spans the endoproteolytic CS region [15,16,17]. The primers were designed for the broadest possible reactivity with recent Eurasian H7 sequences.
For validation of the assays, viral RNA from reference H7 LPAIV and HPAIV was used. Moreover, non-H7 influenza subtypes H5 and H9 as well as other avian respiratory viruses (infectious bronchitis virus (IBV), Newcastle disease virus (NDV)) were tested (Table 3). Viral RNA was purified with the QIAamp®Viral RNA Mini Kit (Qiagen, Hilden Germany) according to the manufacturer’s instructions. Primers were first evaluated in conventional RT-PCRs. The PCR reactions were carried out on a CFX96 thermocycler machine (Bio-Rad) using the following temperature profile: 30 min at 50 °C (RT), 2 min at 94 °C (inactivation of reverse transcriptase/activation of Taq polymerase), followed by 42 cycles of 30 s at 94 °C (denaturation), 30 s at 56 °C (annealing), and 30 s at 68 °C (elongation). Twenty-five μL per reaction were prepared using the SuperScript III One-Step RT-PCR system with Platinum Taq DNA polymerase (Invitrogen, Carlsbad, CA): For one reaction, 6.5 μl of RNase-free water, 12.5 μl reaction mix (2×), 1 μl of SuperScript III RT/Platinum Taq, and 5 μl of template RNA were mixed. Pre selected primers were then screened for their specifity using non H7-subtypes. Amplificates of the expected sizes were generated from both LP and HP phenotypes of subtype H7 viruses by conventional RT-PCR and visualized on an 2% agarose gel (Fig. 1).
Having assured the broad but exclusive specificity of the selected primers for Eurasian H7 viruses, matching probes for use in the RT-qPCR assays were developed. Initially, probes were designed with the aim to universally differentiate between LP and HP Eurasian H7 CS sequences. Probes were therefore placed directly across the sequence stretch encoding the CS. Closer inspection of the alignments, and taking into account also the list of HP H7 CS sequences provided by OFFLU , revealed that HP H7 CS sequences of Eurasian origin viruses were highly divergent: Viruses of separate outbreaks and epizootics represented unique CS sequences with little homology to viruses of other outbreaks. Within an outbreak series, however, HP H7 CS sequences proved to be conserved. This situation is opposed to HPAI H5 viruses of the gs/GD lineage which show considerable conservation even across different clades and allowed designing of a universal conserved probe for the HP phenotype of these viruses . In contrast to HP H7, the HA CS of LP H7 viruses of Eurasian origin origin appeared to be fairly conserved . Therefore, two strategies were followed to prove that sequencing-independent pathotyping by RT-qPCRs is principally possible also for Eurasian H7 viruses:
For HP H7, probes were designed that are specific for viruses of distinct outbreaks. Two distinct HP H7 outbreaks were selected: Isolates from a historic epizootic (Italy 1999, H7N1) and from the most recent HP H7 outbreak in Germany (referred to as ‘Emsland’; a region in the Northwest of Germany where a very high density of poultry population is reared) affecting a single holding in 2015 (H7N7) were chosen and specific Taqman probes matching the HA CS consensus sequences of each of these outbreaks designed (Table 2).
For Eurasian LP H7 a universal probe was developed and several universal Taqman probes were synthesized for comparison.
The same PCR conditions as described above for conventional RT-PCR were used for RT-qPCR, however, 2 μL of the RNAse-free water were replaced by 2 μl specific primer-probe mix. The HP mixes were composed of 1,25 pmol probe/μl and 3,75 pmol/μl for each forward and reverse primer.
Specificity was initially confirmed only for the two HP probes which specifically reacted with their homologous sequences but did not cross react with LP H7 or other HP H7 viruses (Table 3). The standard Taqman LP probes, however, did not sharply distinguish between pathotypes and cross reacted with various HP H7 viruses (not shown). Closer inspection of the alignments revealed that a single G/A mutation in the HA CS distinguished between LP and HP pathotypes (Fig. 2). Consequently, an LNA probe was designed placing the critical nucleotide position at the centre of the respective probe. Using this probe at a concentration of 2,5 pmol in the reaction mix finally allowed clear-cut distinction between LP and HP pathotypes by RT-qPCR (Table 3).
The detection limit of the H7 pathotyping RT-qPCRs was determined by testing ten-fold serial dilutions of viral RNAs extracted from representative H7 LPAI and HPAI viruses. Average values of three independent runs were used for comparisons to a generic RT-qPCR for the M gene of these viruses . A standard curve of each assay was generated showing a linear relationship between the log dilution of the viral RNA and the cycle quantification (Cq) value for both the specific and the generic assays (Fig. 3a-c). Considering the universal LP as well as the ‘Emsland’-specific HP probe, no significant difference between the median Cq values of each specific assay and the M RT-qPCR was found indicating that the newly developed and the generic RT-qPCRs have a similar analytical sensitivity. In contrast, the RT-qPCR detecting the historic Italian H7 HP lineage showed slightly higher sensitivity than the generic M RT-qPCR.
Furthermore, we determined the ability of the H7 pathotyping RTqPCRs to detect mixtures of RNAs of LPAIV and HPAIV derived from the Emsland outbreak in Germany, 2015, and compared it to the M gene-specific generic RT-qPCR (Fig. 3d). Different concentrations of LP/HPAIV-mixtures (0, 0.1%, 1%, 10%, 50% and 100% LP) were generated, and HP H7 RNA was added to 100%. Both RNA species were detected by the specific RT-qPCRs in the mixtures, and the respective Cq values reflected the concentration of the RNA species in the mixtures (Fig. 3d). H7 LP RNA was not detected in the sample containing 100% H7 HP RNA, and vice versa, once more confirming the specificity of the pathotyping RT-qPCRs (Fig. 3d). Thus, these PCRs can be used to study the generation and co-circulation of H7 HPAIV from its LPAI precursor viruses.
Assessment of the diagnostic performance characteristics of the established RT-qPCRs was carried out with a collection of H7 AIV isolates (n = 48) and H7-positive field samples (n = 27) collected between 1999 and 2016. Samples were obtained from the virus repository of the German National Reference Laboratory for Avian Influenza at the Friedrich-Loeffler-Institut, Germany, or kindly provided by the OIE Reference Laboratory for Newcastle Disease and Avian Influenza in Italy, ISZVe, Padua, the Central Veterinary Research laboratory at Dubai, United Arab Emirates, the National Centre for Foreign Animal Disease, Winnipeg, Canada and the WHO Collaborating Centre, London, United Kingdom, under the patronage of the global influenza programme (Table 4). Amplificates produced from these viral RNAs by H7-specific RT-qPCR analysis were also further processed for sequence analysis using the H7-specific reverse primer mix (Table 2) for Sanger sequencing: Following agarose gel electrophoresis and amplicon purification using the QIAquick Gel Extraction Kit (Qiagen, Hilden, Germany) they were cycle-sequenced (BigDye Terminator v1.1 Cycle Sequencing Kit, Applied Biosystems, California, United States) and analysed on an ABI PRISM 3130 Genetic Analyzer (Life Technologies, Darmstadt, Germany) . Partial HA sequences of the diagnostic samples are shown in the sequence alignment of Additional file 1: Figure S1; in all cases subtype H7 was confirmed. Pathotypes were assigned as based on the deduced amino acid sequence of the HACS according to the list of published H7 CS sequences (Additional file 1: Figure S1 and Table 4).
In total, 75 samples positive for AIV of subtype H7 were used. Based on nucleotide sequence analysis and/or IVPI, 49 samples were classified as LPAIV and 15 as HPAIV (Table 4, Additional file 1: Figure S1). They were of both historic and recent origin and mainly derived from European locations. Four samples originated from North America, nine from the United Arab Emirates/Dubai and one represented the Chinese LP H7N9 lineage. The samples mainly consisted of egg-derived isolates or native combined oropharyngeal and cloacal swabs obtained from poultry or wild birds. Seven samples were taken from the environment during a recent HPAIV outbreak in a chicken layer holding in Germany (referred to as ‘Emsland’). For the H7 LP RT-qPCR, 48 out of 56 samples were correctly identified as LP (Table 4, Fig. 4), also including the Chinese LP H7N9 reference virus. Three historic LP isolates (Table 4, nos. 1–3) and the two North American LP H7 viruses (Table 4, nos. 71–72) were not detected despite high viral loads. Sequence mismatches affected binding of either probe and/or primers in these cases. In three further samples (Table 4, nos. 26, 37, 38) low virus loads were detected by the generic M RT-qPCR and these were missed by the H7 LP specific RT-qPCR. However, in most samples, the H7 LP specific RT-qPCR proved to be more sensitive as compared to the generic M specific one (Table 4, Fig. 4). Since none of the HP H7 positive samples cross reacted in the H7 LP RT-qPCR, complete specificity was achieved.
A total of 19 samples harbored HP H7 RNA. None of them was detected by the LP specific RT-qPCR (Table 4, Fig. 4). Two isolates originating from the Italian HP H7N1 epizootic of 1999 were detected by the H7 HP ‘Italy’-specific RT-qPCR (Table 4, nos. 57–58); no further viruses were identified by this PCR. This includes another HPAIV H7N1 isolate from Italy originating from 2002 and distinguished from the 1999 viruses by 13 mutations in the primer and probe binding sites (Table 4, no. 56). Thus, the ‘Italy 99’ RT-qPCR proved to be highly lineage-specific. The second H7 HP RT-qPCR aimed at detecting HP AIV related to the most recent outbreak in Germany in 2015. All nine samples classified to harbor HP H7 were identified by this PCR with a high sensitivity (Table 4, nos. 43–51). At similarily high sensitivity four historic European HP H7 viruses (Table 4, nos. 39, 40, 42, 69), but none of the Italian HP viruses, or an isolate (Table 4, no. 41) representing the large HP H7N7 epizootic affecting the Netherlands, Belgium and Germany in 2003, reacted with either of the two HP specific RT-qPCRs. No cross reactivity to any of the LP H7 samples was detected indicating excellent performance values regarding sensitivity and specificity. Due to our results, the threshold distinguishing reliably between positive and negative samples was set at Cq = 38.
Although not all of the LP and HP H7 samples did show a positive signal with the respective RT-qPCR due to mismatches in the probe binding regions, the newly developed set of primers produced a sequenceable amplificate even of those virus strains. Consequently, pathotype confirmation of a H7 positive sample that tested negative by the LP and HP RT-qPCRs is still possible by nucleotide sequence analysis using the amplificate produced by these RT-qPCRs. In this respect, the newly developed RT-qPCRs resemble the one introduced by Slomka et al.  which also spanned the H7 HACS but its probe targeted a highly conserved sequence stretch outside the CS.
The pathotype-specific RT-qPCRs developed here for avian influenza viruses of subtype H7 proved to be a useful, sensitive and highly specific alternative to nucleotide sequence analysis for the characterization of LPAI and HPAI H7 viruses of European origin. Proper detection of HP H7 viruses required knowledge of the HACS of the specific lineage, and specific probes are to be used for each distinct lineage. Thus, initial characterization of an H7 HP virus still depends on nucleotide sequence analysis of its HACS. However, in case of on-going spread of the identified HP H7 lineage a lineage-specific probe can then be used in a pathotyping RT-qPCR for the swift examination and pathotyping of further cases and outbreaks. Furthermore, the LP LNA probe introduced here was universally usuable for Eurasian LP H7 viruses circulating in Europe over the past decade. In conclusion, these here described RT-qPCRs complement a sequencing-independent approach, and allow a high-speed pathotyping helping the authorities to install necessary control measures in time.
Avian influenza virus
Global initiative on sharing Avian influenza data
Highly pathogenic avian influenza virus
Infectious bronchitis virus
Intravenous pathogenicity index
Low pathogenic avian influenza virus
Newcastle disease virus
Reverse transcription real-time PCR
Spackman E. A brief introduction to the avian influenza virus. Methods Mol Biol. 2008;436:1–6.
Cheung TK, Poon LL. Biology of influenza a virus. Ann N Y Acad Sci. 2007;1102:1–25.
Webster RG, et al. Evolution and ecology of influenza a viruses. Microbiol Rev. 1992;56(1):152–79.
Alexander DJ. An overview of the epidemiology of avian influenza. Vaccine. 2007;25(30):5637–44.
(OIE), W.O.f.A.H., Avian Influenza. Chapter 10.4. Article 10.4.1. Available from: http://web.oie.int/eng/normes/mcode/en_chapitre_1.10.4.pdf. 2010.
(OIE/FAO), W.O.f.A.H.F.a.A.O.o.t.U.N., Network of expertise on animal influenza. Influenza a cleavage sites. Paris/Rome: OIE/FAO; 2015. [Accessed 01 Dec 2016] Available from: http://www.offlu.net/fileadmin/home/en/resource-centre/pdf/Influenza_A_Cleavage_Sites.pdf. 2015.
Fereidouni SR, et al. Rapid molecular subtyping by reverse transcription polymerase chain reaction of the neuraminidase gene of avian influenza a viruses. Vet Microbiol. 2009;135(3–4):253–60.
Leijon M, et al. Rapid PCR-based molecular pathotyping of H5 and H7 avian influenza viruses. J Clin Microbiol. 2011;49(11):3860–73.
Hoffmann B, et al. Rapid and highly sensitive pathotyping of avian influenza a H5N1 virus by using real-time reverse transcription-PCR. J Clin Microbiol. 2007;45(2):600–3.
Naguib MM, et al. Novel real-time PCR-based patho- and phylotyping of potentially zoonotic avian influenza a subtype H5 viruses at risk of incursion into Europe in 2017. Euro Surveill. 2017;22(1). doi: http://dx.doi.org/10.2807/1560-7917.ES.2017.22.1.30435.
(WHO), W.H.O., WHO RISK ASSESSMENT of human infections with avian influenza a(H7N9) virus. [Accessed 23 Feb 2015] Available from: http://www.who.int/influenza/human_animal_interface/influenza_h7n9/RiskAssessment_H7N9_23Feb20115.pdf?ua=1. 2015.
Li Q, et al. Epidemiology of human infections with avian influenza a(H7N9) virus in China. N Engl J Med. 2014;370(6):520–32.
Cima G. Another H7N9 epidemic, with possible HPAI isolates. Am J Vet Res. 2017. Available from: https://www.avma.org/News/JAVMANews/Pages/170415l.aspx.
EC, E.c., Annual report for avian influenza in poultry and wild birds in member states of the European Union in 2015; Available from: https://ec.europa.eu/food/sites/food/files/ad_control-measures_ai_surv-rslt_pltry-wld-brds_2015.pdf. 2015.
Steinhauer DA. Role of hemagglutinin cleavage for the pathogenicity of influenza virus. Virology. 1999;258(1):1–20.
Bosch FX, et al. Proteolytic cleavage of influenza virus hemagglutinins: primary structure of the connecting peptide between HA1 and HA2 determines proteolytic cleavability and pathogenicity of avian influenza viruses. Virology. 1981;113(2):725–35.
Gall A, et al. Universal primer set for amplification and sequencing of HA0 cleavage sites of all influenza a viruses. J Clin Microbiol. 2008;46(8):2561–7.
Hoffmann B, et al. New real-time reverse transcriptase polymerase chain reactions facilitate detection and differentiation of novel a/H1N1 influenza virus in porcine and human samples. Berl Munch Tierarztl Wochenschr. 2010;123(7–8):286–92.
Slomka MJ, et al. Validated real time reverse transcriptase PCR methods for the diagnosis and pathotyping of Eurasian H7 avian influenza viruses. Influenza Other Respir Viruses. 2009;3(4):151–64.
OIE, W.O.f.A.H., Update on highly pathogenic avian influenza in animals (type H5 and H7); Available from: http://www.oie.int/animal-health-in-the-world/update-on-avian-influenza/2015/.2017.
The authors thank Aline Maksimov, Diana Wessler and Mahmoud Naguib for excellent support. We are grateful to Christine Luttermann and her team for conducting the Sanger sequencing analyses, and we wish to acknowledge veterinarians and diagnosticians in Germany for submitting diagnostic samples. We are grateful to colleagues from various research institutions as mentioned in the text who made H7 samples and/or isolates available to the FLI virus repository. We kindly acknowledge all researchers who made available sequences, used in this study, to the EpiFlu database of GISAID.
The study was funded by an intramural grant of the Federal Ministry for Food and Agriculture to the Friedrich-Loeffler-Institut.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
Consent for publication
The authors that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Nucleotide sequences encoding the HA endoproteolytic cleavage site of H7N7 low pathogenic avian influenza viruses generated within this study. (PDF 90 kb)
Nucleotide sequences encoding the HA endoproteolytic cleavage site of H7N7 highly pathogenic avian influenza viruses generated within this study. (PDF 57 kb)