Volume 20, Issues 3–4, June–August 2006, Pages 245–249

A multiplex RT-PCR for detection of type A influenza virus and differentiation of avian H5, H7, and H9 hemagglutinin subtypes

  • Zhixun Xiea, Corresponding author contact information, E-mail the corresponding author,
  • Yao-shan Panga,
  • Jiabo Liua,
  • Xianwen Denga,
  • Xiaofei Tanga,
  • Jianhua Suna,
  • Mazhar I. Khanb, 1, Corresponding author contact information, E-mail the corresponding author
  • a Guangxi Veterinary Research Institute, 51 You Ai North Road, Nanning, Guangxi 530001, People's Republic of China
  • b Department of Pathobiology and Veterinary Science, University of Connecticut, 61 North Eagleville Road, Storrs, CT 06269-3089, USA

Abstract

A multiplex reverse transcriptase-polymerase chain reaction (mRT-PCR) was developed and optimized for the detection of type A influenza virus; the assay simultaneously differentiates avian H5, H7 and H9 hemagglutinin subtypes. Four sets of specific oligonucleotide primers were used in this test for type A influenza virus, H5, H7 and H9 heamagglutinin subtypes. The mRT-PCR DNA products were visualized by gel electrophoresis and consisted of fragments of 860 bp for H5, 634 bp for H7, 488 bp for H9 hemagglutinin subtypes, and 244 bp for type A influenza virus. The common set primers for type A influenza virus were able to amplify a 244 bp DNA band for any of the other subtypes of AIV. The mRT-PCR assay developed in this study was found to be sensitive and specific. Detection limit for PCR-amplified DNA products was 100 pg for the subtypes H5, H7, and H9 and 10 pg for type A influenza virus in all subtypes. No specific amplification bands of the same sizes (860, 634 and 488 bp) could be amplified for RNA of other influenza hemagglutinin subtypes, nor specific amplification bands of type A influenza (244 bp) for other viral or bacterial pathogens.

Keywords

  • Multiplex;
  • Polymerase chain reaction;
  • Avian influenza virus;
  • H5 subtype;
  • H7 subtype;
  • H9 subtype;
  • Hemagglutinin

1. Introduction

Influenza is a zoonotic disease, infecting a wide variety of warm-blooded animals, including birds and mammals. Influenza viruses are classified into types A, B, and C. Influenza A viruses are responsible for major disease problem in birds, as well as in humans [1], [2], [3], [4], [5], [6], [7], [8], [9] and [10]. Infections among domestic or confined birds have been associated with a variety of disease syndromes ranging from sub-clinical to mild upper respiratory disease, to loss of egg production, to acute generalized fatal disease. In domestic avian species, influenza viruses have caused considerable economic losses [1] and [2].

Influenza A are enveloped, negative-sense RNA viruses. Influenza A viruses are further classified into subtypes on the basis of the antigenic properties of their two surface glycoproteins, hemagglutinin (HA) and neuraminidase (NA). To date, 15 HA and 9 NA subtypes have been identified. All influenza A virus subtypes have been found in aquatic and domestic birds, but only a few subtypes have been recovered from mammals and humans.

Among 15 HA subtypes, only H5 and H7 are highly virulent in poultry [1]. Historically, highly pathogenic avian influenza viruses (AIV) of poultry only belonged to the H5 and H7 hemagglutinin (HA) subtypes. Therefore, because there is a greater risk for these subtypes to become highly pathogenic, it is important to identify them specifically in surveillance programs [10] and [16]. Recently, H9 subtypes have been seen to cause infections in poultry [3] and [6]. Separate from the H5N1 influenza virus, another subtype of influenza virus, H9N2, has become panzootic in the last decade and has been isolated from different types of terrestrial poultry worldwide [3], [6] and [9].

Diagnosis of influenza A virus infection is routinely done by the isolation and identification of the virus. Serotyping is required to differentiate the subtypes of the AI viruses and is laborious and time-consuming. Furthermore, other tests required to determine the HA cleavage site sequence must be done to determine its potential virulence [12] and [13]. Although single band PCR has been used to detect and differentiate subtypes, it only recognizes one specific subtype at a time [14]. Standard RT-PCR has been previously applied to the detection of avian influenza virus [10], [15] and [16] and each of the 15 HA subtypes [14]. In addition, real time-RT-PCR assays for influenza virus have been developed for the detection of influenza virus types A and B [17] and differentiation of two subtypes H5, and H7 [11]. Equipment costs and specific technical training requirements limit usefulness of these assays as routine laboratory tests. In our study, we have developed a specific and sensitive multiplex RT-PCR that can simultaneously detect and differentiate the three most important subtypes of avian influenza viruses.

2. Materials and methods

2.1. Avian pathogens and culture conditions

The avian pathogens used in this study are listed in Table 1. All avian influenza virus (AIV) subtypes, Newcastle disease virus (NDV), and infectious bronchitis virus (IBV) were propagated in the allantoic cavity of 10-day-old specific-pathogen-free (SPF) embryonated chicken eggs, whereas infectious laryngotracheitis virus (ILT) was propagated on the chorioallantoic membrane in 10-day-old SPF embryonated chicken eggs as described [18]; the allantoic fluids from embryonated eggs infected with AI, NDV, ILT and IBV were harvested after 36 h of incubation at 37 °C [18]. Mycoplasma gallisepticum (MG) was propagated in Frey's broth and incubated at 37 °C as described [19].

Table 1. Avian pathogens used in multiplex PCR

Avian pathogenSubtypeSourceResults of mRT-PCR
Type AH5H7H9
Duck/HK/717/79-d1H1N3Uni. Of HK+
Duck/HK/717/79-d7H1N3Uni. of HK+
Human/NJ/8/76H1N1Uni. of HK+
Duck/HK/77/76 d77/3H2N3Uni. of HK+
Duck/HK/77/76H2N3Uni. of HK+
Duck/HK/526/79/2BH3N6Uni. of HK+
Duck/HK/526/79/2B-1H3N6Uni. of HK+
Duck/HK/668/79H4N5Uni. of HK+
Duck/HK/668/79-1H4N5Uni. of HK+
Duck/HK/313/78H5N3Uni. of HK++
Duck/HK/313/78-1H5N3Uni. of HK++
Duck/Guangxi/1/04H5N1GVRI++
Duck/Guangxi/2/04H5N1GVRI++
Duck/Guangxi/3/04H5N1GVRI++
Chicken/Guangxi/1/04H5N1GVRI++
Chicken/Guangxi/2/04H5N1GVRI++
Goose/Guangxi/1/04H5N1GVRI++
Goose/Guangxi/2/04H5N1GVRI++
Duck/HK/531/79-1H6N8Uni. of HK+
Duck/HK/531/79H6N8Uni. of HK+
Duck/HK/47/76H7N2Uni. of HK++
Duck/HK/47/76-1H7N2Uni. of HK++
Turkey/ont/6118/68H8N4Uni. of HK+
Turkey/ont/6118/68-1H8N4Uni. of HK+
Duck/HK/147/77H9N6Uni. of HK++
Duck/HK/147/77-1H9N6Uni. of HK++
Duck/Guangxi/1/00H9N2GVRI++
Duck/Guangxi/2/00H9N2GVRI++
Duck/Guangxi/3/00H9N2GVRI++
Duck/HK/876/80H10N3Uni. of HK+
Duck/HK/876/80-1H10N3Uni. of HK+
Duck/HK/661/79H11N3Uni. of HK+
Duck/HK/661/79-1H11N3Uni. of HK+
Duck/HK/862/80H12N5Uni. of HK+
Duck/HK/862/80-1H12N5Uni. of HK+
Gull/MD/704/77H13N5Uni. of HK+
Gull/MD/704/77-1H13N5Uni. of HK+
NDVF68-E9CIVDC, Beijing
ILTVFieldCIVDC, Beijing
IBVM41CIVDC, Beijing
MGS6UCDavis, Calif
LiverSPF chicken
LungSPF chicken
Small intestineSPF chicken

Uni. of HK, University of Hong Kong, China; GVRI, Guangxi Veterinary Research Institute; CIVDC, China Institute of Veterinary Drug Control; UCDAVIS, University of California, Davis.

Full-size table

2.2. Extraction of RNA and DNA

Allantoic fluids from embryonated eggs infected with AI, NDV, IBV and ILT were first clarified by centrifugation at 500g for 15 min. The supernatants were transferred to new tubes and then centrifuged at 45,000g for 30 min at 4 °C. The supernatants were discarded and the pellets treated with 25 μl Rnase-free TE buffer (10 mM Tris–HCl and 1 mM ethylendiaminetetraacetic acid (EDTA)). The RNA extraction from AI, NDV, and IBV was carried out according to the Trizol LS manufacturer's protocol (Trizol, Invitrogen, Carlsbad, CA, USA). DNA from ILT and MG was extracted using the phenol:chloroform:isoamyl alcohol (25:24:1 v/v) (Amersham Life Science, Cleaveland, OH, USA) method as described by Pang et al. [20]. The concentrations of the RNA or DNA were determined by spectrophotometry using the UV2501PC (Shimadzu Corporation, Tokyo, Japan) and stored at −20 °C. An aliquot of 250 μl each was used for RNA extraction from the AI virus isolates listed in Table 1 according to the Trizol LS manufacturer's protocol (Trizol, Invitrogen, Carlsbad, CA, USA). Four-week-old SPF chicken (Jinan City, Shandong, China) was euthanized and 30 mg tissue samples from lung, liver, and small intestine were minced in TE buffer and RNA extraction was carried out according to the Trizol LS manufacturer's protocol for tissue samples (Trizol, Invitrogen, Carlsbad, CA, USA).

2.3. Primers designs and selection

Four sets of primers that specifically amplify type A influenza virus and simultaneously detect and differentiate H5, H7 and H9 subtypes of hemagglutinins are listed in Table 2. The primers for type A avian influenza were designed after reviewing the matrix gene sequences of over 10 different subtypes of AIV from chickens, ducks, geese, swine and humans; whereas primers for H5 subtypes were designed from hemagglutinin genes of avian origin. Primers for H7 and H9 subtypes were used from published data [14]. All four sets of oligonucleotide primers were synthesized at the Takara Shuzo Co., Ltd (Dalian, Shandong, China). The primers were aliquoted to a final concentration of 100 pmol/μl and stored at −20 °C.

Table 2. Multiplex RT-PCR primers

Primer's namePrimer's oligonucleotide sequenceaProduct (bp)
Type A influenza virus
XZ145-25′-CTTCTAACCGAGGTCGAAAC-3′
XZ1465′-AGGGCATTTTGGACAAAKCGTCTA-3′244

Subtype H5
XZ H5-15′-ACACATGCYCARGACATACT-3′
XZ H5-55′-CAGGAACGYTCWCCTGAKTCT-3′860

Subtype H7
XZ H7-15′-GGGATACAAAATGAAYACTC-3′
XZ H7-25′-CCATABARYYTRGTCTGYTC-3′634

Subtype H9
XZ H9-15′-CTYCACACAGARCACAATGG-3′
XZ H9-25′-GTCACACTTGTTGTTGTRTC-3′488
a

Codes for mixed bases position: Y, C/T; R, A/G; W, A/T; B, G/C/T; K, G/T.

Full-size table

2.4. Optimization of Reverse transcription and multiplex PCR reaction

The mRT-PCR consists of a two-step procedure, which includes reverse transcription (RT) and PCR amplification. An RT-PCR kit (Takara Shuzo Co., Ltd, Dalian, Shandong, China) was used for the reverse transcription reaction. RT was performed in 20 μl volumes, in which the reaction mixture contained 2 μl RNA in different concentrations and 5 mM MgCl2, in PCR buffer Buffer (500 mM KCl, 100 mM Tris–HCl, pH 8.3), 1 mM of each dinucleoside triphosphate (dNTP), 2 units RNase inhibitor, 0.25 units avian myeloblastosis virus (AMV) reverse transcriptase, 1.25 pmol upstream primers of XZ145-2, XZ H5-1, XZ H7-1 and XZ H9-1. DEPC treated water was added to bring the final volume to 20 μl. RT was performed in a thermal cycler (Model 9600, Perkin Elmer Cetus, Norwalk, CT) for one cycle at 42 °C for 25 min, 99 °C for 3 min and 4 °C for 5 min.

The multiplex PCR was performed in a 50 μl volume using a PCR kit (Takara Shuzo Co., Ltd, Dalian, Shandong, China). The reaction contained 5 mM MgCl2, 1×PCR Buffer, 10 mM of each dNTP, 0.5 pmol of each down stream primer XZ146, XZ H5-5, XZ H7-2 and XZ H9-2, and 1 unit TaKaRa LA Taq™ (Takara Shuzo Co., Ltd, Dalian, Shandong, China). This mixture was added to the RT reaction tubes. Sterile deionized water was added to the mixture to bring the total volume to 50 μl. The mPCR was carried out in the same thermal cycler used for RT. The cycling protocol consisted of an initial denaturing at 94 °C for 5 min, then 35 cycles that each consisted of denaturing at 94 °C for 45 s, annealing at 55 °C for 45 s, and extension at 72 °C for 105 s. The sample was then heated at 72 °C for 10 min for a final extension. A negative control did not contain template cDNA and consisted of PCR master mix, all four sets of primers and deionized water.

2.5. Detection of amplified nucleic acid products

Agarose gel electrophoresis was used to detect mRT-PCR nucleic acid products. A volume of 10 μl of amplified PCR nucleic acid products was subjected to electrophoresis at 80 V in horizontal gels containing 1% agarose with Tris–borate buffer (45 mM Tris–borate, 1 mM EDTA) as described [21]. The Gels were stained with ethidium bromide (0.5 mg/ml or 0.5 μg/ml−l), and exposed to UV light to visualize the amplified nucleic acid products, and photographed using a Bio-vision Post-electrophoresis instrument (Vilber Lourmat, Paris, France).

2.6. Multiplex-RT-PCR sensitivity and specificity

Specificity of the mRT-PCR was determined by examining the ability of the test to detect type A influenza viruses and simultaneously differentiate H5, H7, and H9 subtypes. The mRT-PCR was tested using other avian pathogens that produce similar clinical signs or that can be present in mixed infections with AI subtypes. These pathogens and AIV subtypes are listed in Table 1. Random samples from the amplified mRT-PCR DNA bands were isolated from the gel and purified using a DNA Glass-milk purification kit (BioVed, Beijing, China). The purified DNA products were sent out to Takara Shuzo Co., Ltd, Dalian, Shandong, China for DNA sequencing. DNA sequences of the amplified products were than analyzed using the DNAStar software to confirm the amplified DNA sequence in the products. To determine the ability of the multiplex PCR assay to detect type A influenza virus and differentiate three subtypes in the same reaction, we used a mixture of RNA concentrations ranging from 200 ng to 10 fg RNA in various combinations of all subtypes of avian influenza virus. Sensitivity of the mRT-PCR was determined by making 10-fold serial dilutions of a mixture containing 10 ng of the RNA templates of each of the three viruses of AIV.

3. Results

Multiplex RT-PCR was developed to detect group A avian influenza viruses and simultaneously differentiate three hemagglutinin subtypes H5, H7, and H9 in a single reaction through 35 cycles of PCR. The mRT-PCR products were 244 bp for type A avian influenza viruses (group specific), 860 bp for H5, 634 bp for H7, 488 bp for H9. These products were visualized by electrophoresis (Fig. 1 and Fig. 2). This mRT-PCR was found to be a specific assay for type A avian influenza and hemagglutinin subtypes H5, H7, and H9, with no amplification of nucleic acid from NDV-F48, ILTV, IBV-M41, MG-S6 and SPF chicken tissues. The mRTPCR specifically detected the H5, H7, and H9 hemagglutinin subtypes along with their group specific type A specification. Thirty-seven avian influenza type A, including nine H5 subtypes, two H7 subtypes, and five H9 subtypes (listed in Table 1) were detected by the mPCR with the amplification of 244 bp DNA products. Limit of detection by visualization of PCR-amplified DNA products was 100 pg for the subtypes hemagglutinins, and 10 pg for type A avian influenza viruses (Fig. 2). Throughout the development of mRT-PCR, various modifications were made to the annealing temperature, extension time, cycle quantity, and primer concentrations in order to obtain the optimal conditions for this mRT-PCR. No spurious PCR amplification reactions among all influenza subtypes and other pathogens were noticed with various amounts of DNA and RNA mixtures. All the negative controls including RNA samples from the SPF chickens tissues were negative. DNAStar software analysis indicated that the mRT-PCR amplified DNA products were similar to the group A, subtypes H5, H7, and H9 genes sequences of avian influenza viruses.

Fig. 1. Agarose gel elecrophoresis of multiplex RT-PCR amplified products from purified DNA and RNAs of known avian influenza subtypes and other avian pathogens. Lane 1, molecular size marker; lane 2, PCR reagent buffer as a negative control; lane 3, H5N1 (Duck/Guangxi/2/04), H7N2 (Duck/HK/47/76-1) and H9N2 (Duck/Guangxi/3/00) subtypes of AIV; lane 4, H5N3 (Duck/HK/313/78); lane 5, H5N1 (Duck/Guangxi/2/04); lane 6, H7N2 (Duck/HK/47/76); lane 7, H9N2 (Duck/Guangxi/2/00); lane 8, H9N6 (Duck/HK/147/77); lane 9, NDVF68-E9; lane 10, IBV-M41; lane11, ILTV; lane 12, MG-S6.

View thumbnail images

Fig. 2. Sensitivity of multiplex RT-PCR. Agarose gel electrophoresis of multiplex RT-PCR amplified products from purified RNAs from avian influenza subtypes H5N1 (Duck/Guangxi/2/04), H7N2 (Duck/HK/47/76-1) and H9N2 (Duck/Guangxi/3/00) subtype of AIV. Lane 1, molecular size marker; lane 2, PCR reagent buffer as a negative control; lane 3, 10 ng; lane 4, 1 ng; lane 5, 100 pg; lane 6, 10 pg; lane 7, 1 pg; lane 8, 100 fg; lane 9, 10 fg; lane10, 1 fg.

View thumbnail images

4. Discussion and conclusions

The mRT-PCR developed here was able to detect type A influenza virus hemagglutin subtypes from H1 to H13 as tested (Table 1) and simultaneously detected and differentiated the very important hemagglutinin subtypes H5, H7 and H9 in one single reaction (Fig. 1 and Fig. 2). These hemagglutin subtypes recently have been a cause of human infections [2], [4], [5], [6] and [7]. Therefore, a mRT-PCR which can rapidly identify type A influenza as well as subtypes H5, H7 and H9, will be very important for the control of disease transmission in poultry and in humans. Use of this assay will also help reduce the economic losses in poultry associated with an AIV outbreak. This mRT-PCR is sensitive, specific, cost effective and it may be useful in diagnosis, screening and surveillance of poultry, including live bird market population. This mRT-PCR has the added benefits of being time saving, and using fewer reagents. Studies will be carried out to further test the specificity and sensitivity of this mRT-PCR on avian influenza subtypes from the various diagnostic and research laboratories as well as on clinical samples originated in the USA and South East China.

Acknowledgements

The Guangxi Science and Technology Bureau and Guangxi Aquaculture and Animal Husbandry Bureau supported this work. We also thank Dr K.F. Shortridge, Department of Microbiology at the University of Hong Kong for providing the avian influenza subtypes. Authors thank Drs Sandy Bushmich and Herbert Van Kruiningen for their constructive suggestions and critical review of the manuscript.

References

    • [1]
    • D.J. Alexander
    • A review of avian influenza in different bird species

    • Vet Microbiol, 74 (2000), pp. 3–13

    • [2]
    • D. Normile
    • North Korea collaborates to fight bird flu

    • Science, 308 (2005), p. 175

    • [3]
    • K.R. Cameron, V. Gregory, J. Banks, I.H. Brown, D.J. Alexander, A.J. Hay et al.
    • H9N2 subtype influenza viruses in poultry in Pakistan are closely related to the H9N2 viruses responsible for human infection in Hong Kong

    • Virology, 278 (2000), pp. 36–41

    • [4]
    • R.A. Fouchier, P.M. Schneeberger, F.W. Rozendaal, J.M. Broekman, S.A. Kemink, V. Munster et al.
    • Avian influenza A virus (H7N7) associated with human conjunctivitis and a fatal case of acute respiratory distress syndrome

    • Proc Natl Acad Sci USA, 101 (2004), pp. 1356–1361

    • [5]
    • D. Normile
    • Avian flu. First human case in Cambodia highlights surveillance shortcomings

    • Science, 307 (2005), p. 1027

    • [6]
    • Y. Guo, J. Li, X. Chen, M. Wang, Y. Zhou, X.H. Li et al.
    • Discovery of men infected by avian influenza A (H9N2) virus

    • Chin J Exp Clin Virol, 13 (1999), pp. 105–108

    • [7]
    • T.H. Tran, T.L. Nguyen, T.D. Nguyen, T.S. Luong, P.M. Pham, V.C. Nguyen et al.
    • Avian influenza A (H5N1) in 10 patients in Vietnam

    • N Engl J Med, 350 (2004), pp. 1179–1188

    • [8]
    • M. Koopmans, B. Wilbrink, M. Conyn et al.
    • Transmission of H7N7 avian influenza A virus to human being during a large outbreak in commercial poultry farm in the Netherlands

    • Lancet, 363 (2004), pp. 587–593

    • [9]
    • K.S. Li, K.M. Xu, J.S. Peiris, L.L. Poon, K.Z. Yu, K.Y. Yuen et al.
    • Characterization of H9 subtype influenza viruses from the ducks of Southern China: a candidate for the next influenza pandemic in humans?

    • J Virol, 77 (2003), pp. 6988–6994

    • [10]
    • R.G. Webster, Y. Kawaoka
    • Avian influenza

    • Rev Poult Biol, 1 (1987), pp. 212–246

    • [11]
    • E. Spackman, D.A. Senne, T.J. Myers, L.L. Bulaga, L.P. Garber, M.L. Perdue et al.
    • Development of a real-time reverse transcriptase PCR assay for type a influenza virus and the avian H5 and H7 hemagglutinin subtypes

    • J Clin Microbiol, 40 (2002), pp. 3256–3260

    • [12]
    • D.A. Senne, B. Panigrahy, Y. Kawaoka, J.E. Pearson, J. Suss, M. Lipkind et al.
    • Survey of the hemagglutinin (HA) cleavage site sequence of H5 and H7 avian influenza viruses: amino acid sequence at the HA cleavage site as a marker of pathogenicity potential

    • Avian Dis, 40 (1996), pp. 425–437

    • [13]
    • T. Horimoto, Y. Kawaoka
    • Direct reverse transcriptase PCR to determine virulence potential of influenza A viruses in birds

    • J Clin Microbiol, 33 (1995), pp. 748–751

    • [14]
    • S. Lee, P. Chang, J. Shien, M. Cheng, H. Shieh
    • Identification and subtyping of avian influenza viruses by reverse transcription-PCR

    • J Virol Methods, 97 (2001), pp. 13–22

    • [15]
    • M. Munch, L. Nielsen, K. Handberg, P. Jorgensen
    • Detection and sub typing (H5 and H7) of avian type A influenza virus by reverse transcription-PCR and PCR-ELISA

    • Arch Virol, 146 (2001), pp. 87–97

    • [16]
    • E. Starick, A. Romer-Oberdorfer, O. Werner
    • Type- and subtype-specific RT-PCR assays for avian influenza A viruses (AIV)

    • J Vet Med B Infect Dis Vet Public Health, 47 (2000), pp. 295–301

    • [17]
    • L. Van Elden, M. Nijhuis, P. Schipper, R. Schuurman, A. Van Loon
    • Simultaneous detection of influenza viruses A and B using real-time quantitative PCR

    • J Clin Microbiol, 39 (2001), pp. 196–200

    • [18]
    • D. Swayne, D.A. Senne, C. Beard
    • Avian influenza

    • D.E. Swayne, J.R. Glisson, M.W. Jackson, J.E. Pearson, W.M. Reed (Eds.), A laboratory manual for the isolation and identification of avian pathogens, American association of avian pathologists (4th ed), Kennett Square, Pennsylvania (1998), pp. 150–155

    • [19]
    • M.I. Frey, R.P. Hanson, D.P. Anderson
    • A medium for the isolation of avian mycoplasma

    • Am J Vet Res, 29 (1968), pp. 2163–2171

    • [20]
    • Pang, Yao-shan, H. Wang, T. Girshick, Z. Xie, M.I. Khan
    • Development and application of a multiplex polymerase chain reaction for avian respiratory agents

    • Avian Dis, 46 (2002), pp. 691–699

    • [21]
    • J. Sambrook, E.F. Fritsch, T. Maniatis
    • Molecular cloning

    • A laboratory manual (2nd ed)Cold Spring Harbor Laboratory Press, New York (1989) p. 6.6–7 [B.23; E.8–9]

Corresponding author contact information
Corresponding authors. Tel.: +1 860 486 0228; fax: +1 860 486 2794.
Corresponding author contact information
Tel.: +86 771 3120371.
1

Visiting Scholar at the Guangxi Veterinary Research Institute, Nanning, Guangxi, People's Republic of China.