To avoid any probable contamination, each PCR positive sample was separately amplified using specific primers targeting the VP1/VP3 junction region (
19,
20). Touch down PCR was carried out in a 25 μL reaction mixture consisting of 2.5 μL 10x PCR buffer (Fermentas Thermo Scientific, Lithuania), 0.75 μL MgCl
2 50 mM (Fermentas, Thermo Scientific, Lithuania), 0.5 μL 10 mM dNTPs (Fermentas, Thermo Scientific, Lithuania), 0.5 μL of 0.04 μmol forward primer, 5'-GCA GCC CAA AAG AAC TTC AC-3', 0.5 μL 0.04 μmol of reverse primer, 5'-ATT TCA GCA GCT TGG AGT GC-3', 0.5 μL Taq DNA polymerase 100 U (Fermentas, Thermo Scientific, Lithuania), 1 μL DMSO, and 2 μL of the EV-positive cDNA (
18,
19). The amplification reactions were as follows: initial denaturation at 95°C for 5 minutes, followed by 15 cycles comprising denaturation at 95°C for 35 seconds, annealing at 50°C for 16 seconds, and extension at 72°C for 35 seconds, followed by 40 cycles consisting of denaturation at 95°C for 35 seconds, annealing at 48°C for 16 seconds, and extension at 72°C for 35 seconds, and final extension at 72°C for 10 minutes. The PCR product underwent electrophoresis on a 1.5 % agarose gel containing 0.5 μg/mL ethidium bromide (Cina Gen, Iran) using a 50-bp DNA ladder as a molecular marker to detect a product size of 226bp. Diethylpyrocarbonate water and the EV71, which was kindly provided by Dr. Qibin Leng, were used as negative and positive controls, respectively.
The product was purified using a QIAquich gel extraction kit (Qiagen, Inc. GmbH Hilden, Germany) according to the manufacturer’s instructions. It was then sent to a dealer company in Iran, which sent it on to a laboratory in South Korea for sequencing in both directions using a BigDye terminator cycle sequencing kit (Perkin Elmer-Applied Biosystems, Inc. CA), and specific primers on an ABI Prism 3730/3100 DNA analyzer ( Perkin Elmer-Applied Biosystems, Inc.). The sequence data were analyzed using BioEdit software (version 7.0.5). Phylogenic analysis was conducted using the molecular evolutionary genetics analysis software package, version 5.0, and the neighbor-joining method with 1000 bootstrap repetitions (
Figure 1).
The polio virus type 1 was used as an out-group. The tree was elaborated by the neighbor-joining method using molecular evolutionary genetics analysis 4, and the bootstrap (1,000 replicates) supporting each node is indicated. The scale bar corresponds to 0.1 substitutions/site references and GenBank accession numbers for the sequences used in VP1 comparisons were as follows: EV71B4: gi [AF376067.1] (Australia), EV71B3: gi [AF376117.1] (Australia), EV71B2: gi [AF009530.1] (USA), EV71B1: gi [AF135886.1] (USA), EV71A: gi [AB204852.1] (Japan), EV71A: gi [AB204853.1] (Japan), EV71C3: gi [DQ341356.1] (Malaysia), EV71C1: gi [AF376098.1] (Australia), EV71C1: gi [AF376080.1] (Australia), EV71C4: gi [FJ606447.1] (China), EV71C4: gi [FJ606448.1] (China), EV71C4: gi [FJ606450.1] (China), EV71C4: gi [HQ129932.1] (China), EV71C4: gi [AF302996.1] (China), EV71C4: gi [HQ188292.1] (China). The types of EV 71 analyzed in this study are indicated by diamond with their genotypes.
These EV71 sequences have been deposited in the GenBank database under GenBank accession numbers: gi/JX495150/EV71, to gi/JX495164/EV71.
Statistical analysis was carried out using the chi-square test, with statistical package for the social sciences software, version 19.0. P < 0.05 was considered as the level of statistical significance.