ABSTRACT
Avian influenza H5N1 virus is a global threat. An emergence of a reassortant virus with a pandemic potential is a major concern. Here we describe a multiplex reverse transcription-PCR assay that is specific for the eight genomic segments of the currently circulating H5N1 viruses to facilitate surveillance for a virus resulting from reassortment between human influenza virus and the H5N1 virus.
The genome of influenza A viruses is composed of eight gene segments, which encode polymerase (PB1, PB2, and PA), nucleocapsid (NP), hemagglutinin (HA), neuraminidase (NA), matrix (M) protein, and nonstructural (NS) protein. Since influenza viruses have segmented genomes, they are able to evolve through reassortment (1). Both the 1957 and 1968 pandemic strains are believed to be reassortants in which HA and PB1 genes, with or without the NA gene, were replaced by the corresponding genomic segments of avian influenza virus strains (5). The reassortments provided the viruses with new HA and NA antigenic profiles and the ability to spread efficiently from human to human. The ongoing outbreak of the highly virulent avian influenza virus H5N1 is concerning because of the possibility of a reassortant virus emerging with pandemic potential (3). Monitoring for a reassortant with mixed genomic segments from human and avian viruses is therefore very important. However, detection of a reassortant requires sequencing of all eight genomic segments, a process which is laborious and time consuming. Therefore, we developed a simple and high-throughput screening method using multiplex reverse transcription-PCR (RT-PCR) specific for eight genomic segments of recent H5N1 viruses to detect the probable reassortant viruses worthy of further characterization.
For the primer design, 100 full-length sequences arbitrarily selected from the Influenza Virus Resource Database (National Center for Biotechnology Information; http://www.ncbi.nlm.nih.gov/genomes/FLU/ ) were used to search for the sequence(s) that was conserved in H5N1 and that differed from those of other subtypes. The viruses used in this study comprised 28 human influenza virus isolates and 103 H5N1 isolates from various animal species (94 isolates from poultry and other birds, 6 isolates from tigers, 2 isolates from leopards, and 1 isolate from a cat). These viruses have been collected since the beginning of the epidemic in the country. In 2003, 2004, 2005, and 2006, 7, 34, 60, and 2 animal H5N1 virus isolates were collected, respectively. The human viruses comprised 6 H5N1, 3 H1N1, 8 H3N2, and 11 influenza B viruses. These viruses were propagated in MDCK cells or embryonated eggs. The H5N1 viruses were originally subtyped by RT-PCR. The subtyping was confirmed by sequencing for all human H5 and some animal H5 viruses (2). Available accession numbers of some of the H5N1 virus isolates used in this study are shown in Table S1 in the supplemental material and also in a report by Puthavathana et al. (4). The human H1, H3, and influenza B viruses were characterized and confirmed with respect to their subtypes by hemagglutination inhibition assay at the WHO Reference Laboratories, Melbourne, Australia, or the Centers for Disease Control and Prevention of the United States via the National Reference Laboratory on Influenza, NIH, Thailand.
RNA was extracted by using a QIAamp viral RNA Mini kit (QIAGEN, Valencia, CA). A 140-μl volume of virus suspension from the infected culture or embryonated egg yielded a volume of 60 μl of extracted RNA. Multiplex RT-PCR was carried out in three reaction tubes with combinations of the following primer sets: NA plus PB1 plus NP, HA plus PA plus PB2, or NS plus M. Sequences of these primers are shown in Table 1. The thermocycling was performed in a GeneAmp PCR system 2400 thermal cycler (Perkin Elmer) using a QIAGEN OneStep RT-PCR kit (QIAGEN, Valencia, CA) with standard reaction conditions (a 50 μl volume containing 1× PCR buffer, 5 to 10 μl of extracted RNA, 400 μM (each) deoxynucleoside triphosphate, 2.5 mM MgCl2, 2 μl enzyme mix, and 2 U RNase inhibitor). The thermocycling profile was 50°C for 30 min, 95°C for 15 min, 35 cycles of 95°C for 45 s, 55°C for 45 s, and 72°C for 45 s, and a final extension of 72°C for 10 min. The PCR products were visualized by ethidium bromide staining after electrophoresis in a 2% agarose gel.
In our multiplex RT-PCR assay, potential reassortants would show the lack of one or more of the eight amplicons in the case of H5N1 viruses and the presence of one or more amplicons in the case of H1N1 and H3N2 viruses. With the designed primer set, we could initially amplify all eight genomic segments from 105 out of 109 H5N1 viruses. Four viruses initially gave negative amplification results for NA. However, these four viruses yielded positive NA band results in separated single-gene amplifications, and multiplex reamplification using 1.5- to 2-fold-higher RNA input gave positive amplification results for all genes. Representative results obtained with the multiplex RT-PCR products are shown in Fig. 1. None of the human H3N2, H1N1, or influenza B virus isolates yielded any amplification products by this method (data not shown). In order to verify the lack of reassortment, all human H5N1 viruses and a subset of animal viruses (62 out of 103) were subjected to a full-genome sequencing protocol as previously described (4). Phylogenetic analyses of the sequence data revealed that all viruses contained their respective genomic segments and that no evidence of reassortment was found (data not shown). In addition to detecting human-avian virus reassortants, our assay may also be able to detect some reassortants that acquire genomic segments from some other unrelated avian viruses. However, the assay was not designed to detect reassortment among avian viruses, which are unlikely to play a role in viral adaptation to the human host. The absence of a band in the amplification of H5N1 viruses is not always due to reassortment. Full-genome sequencing is needed for result confirmation. Because the performance of this assay is directly affected by variation in the sequences of circulating H5N1 strains, the primer sequences should be frequently reevaluated and updated accordingly.
Our multiplex RT-PCR was designed for egg- or MDCK-propagated viruses. Although the assay may be able to detect some of the genes in some clinical specimens, the level of sensitivity is probably not sufficient to reliably detect all the genomic segments simultaneously. This makes our assay less suitable for direct screening using clinical specimens. Nevertheless, the possibility of using a more sensitive version of this assay with clinical specimens should be further explored.
A representative photograph of multiplex RT-PCR products from an H5N1 virus after electrophoresis in 2% agarose and ethidium bromide staining. Each lane shows products from each tube. The first lane is a 100-bp DNA ladder, and the target gene of each amplification product is labeled on the right.
Sequences, concentrations used in the reaction, nucleotide positions, and product sizes of PCR primers for the eight gene segments
ACKNOWLEDGMENTS
This study was supported by research grants from the National Center for Genetic Engineering and Biotechnology and the Thailand Research Fund for Advanced Research Scholar.
We thank Darat Lauhakirti and Sirawat Srichatrapimuk for their collaboration in the preparation of the manuscript.
FOOTNOTES
- Received 19 February 2007.
- Returned for modification 23 February 2007.
- Accepted 14 March 2007.
- Copyright © 2007 American Society for Microbiology