A highly divergent Encephalomyocarditis virus isolated from nonhuman primates in Singapore

Background In 2001 and 2002, fatal myocarditis resulted in the sudden deaths of four, two adult and two juvenile, orang utans out of a cohort of 26 in the Singapore Zoological Gardens. Methods Of the four orang utans that underwent post-mortem examination, virus isolation was performed from the tissue homogenates of the heart and lung obtained from the two juvenile orang utans in Vero cell cultures. The tissue culture fluid was examined using electron microscopy. Reverse transcription and polymerase chain reaction with Encephalomyocarditis virus (EMCV)-specific primers targeting the gene regions of VP3/VP1 and 3D polymerase (3Dpol) confirmed the virus genus and species. The two EMCV isolates were sequenced and phylogenetic analyses of the virus genes performed. Serological testing on other animal species in the Singapore Zoological Gardens was also conducted. Results Electron microscopy of the two EMCV isolates, designated Sing-M100-02 and Sing-M105-02, revealed spherical viral particles of about 20 to 30 nm, consistent with the size and morphology of members belonging to the family Picornaviridae. In addition, infected-Vero cells showed positive immunoflorescence staining with antiserum to EMCV. Sequencing of the viral genome showed that the two EMCV isolates were 99.9% identical at the nucleotide level, indicating a similar source of origin. When compared with existing EMCV sequences in the VP1 and 3Dpol gene regions, the nucleotide divergence were at a maximum of 38.8% and 23.6% respectively, while the amino acid divergence were at a maximum of 33.9% and 11.3% respectively. Phylogenetic analyses of VP1 and 3Dpol genes further grouped the Sing-M100-02 and Sing-M105-02 isolates to themselves, away from existing EMCV lineages. This strongly suggested that Sing-M100-02 and Sing-M105-02 isolates are highly divergent variants of EMCV. Apart from the two deceased orang utans, a serological survey conducted among other zoo animals showed that a number of other animal species had neutralizing antibodies to Sing-M105-02 isolate, indicating that the EMCV variant has a relatively wide host range. Conclusions The etiological agent responsible for the fatal myocarditis cases among two of the four orang utans in the Singapore Zoological Gardens was a highly divergent variant of EMCV. This is the first report of an EMCV infection in Singapore and South East Asia.


Introduction
Between July 2001 and January 2002, four (3 Bornean, Pongo pygmaeus and 1 Sumatran, Pongo abelii) orang utans, two (male and female) adults and two (male and female) juveniles, from a cohort of 26 in the Singapore Zoological Gardens died suddenly. Post-mortem examination of the four orang utans revealed pathology characterized by multifocal myocarditis, pulmonary congestion and edema, hydropericardium, hydrothorax and Ascites.
A literature review revealed several reports describing Encephalomyocarditis virus (EMCV) as the etiological agent of viral myocarditis in captive primates. These primates included semi-wild bonobos [1], baboons [2], chimpanzees [3], lemurs [4], rhesus macques [5], and orang utans [3,6]. EMCV infection was also implicated in the deaths of various animals in Audubon Park Zoo, New Orleans, in 1985 [7]; and in elephants in a Florida zoo in 1997 [8]. Outside the United States, sporadic outbreaks of EMCV involving a variety of zoo animals had occurred in Taronga Zoo in Australia from 1987 to 1995 [3]; in free-living elephants in Kruger National Park, South Africa between 1993 and 1994 [9,10]; and in an Italian zoo affecting 15 different primates between 2006 and 2008 [4]. In addition, EMCV infections were reported in Russia from monkeys bred from the Sukhumi Breeding Center in 1974; and in the Adler Breeding Center, since 2001 [11]. EMCV has also been recognized as a porcine pathogen, with EMCV infections in European pigs associated with sudden deaths and reproductive failure [12][13][14][15][16]. In Asia, EMCV was isolated from pigs in South Korea and implicated as the cause of reproductive failure in pigs in Taiwan [17,18]. Apart from EMCV, there has been one other report describing Coxsackie virus B4 as the etiological cause of fatal myocarditis in a female orang utan at the Okinawan Zoo in 1999 [19].
Although the potential for EMCV to cross the species barrier has been demonstrated as seen from the various zoo outbreaks described above, human cases have fortunately been rare. Sporadic human EMCV infections and disease have been documented by virus isolation from different specimen types such as serum, stool samples, cerebrospinal fluid and throat washings [20,21]. A recent study describing the etiology of acute febrile disease in locations across South America concluded that there is evidence supporting a role for EMCV in human infection and febrile illness [22]. However, the extent of the effect of EMCV on human health is still largely unknown because the disease is so infrequent in humans.
We hypothesized that the etiological agent behind the deaths of the orang utans in the Singapore Zoological Gardens was either a Coxsackie virus B, or more likely, an EMCV. Although both viruses are members of the same family, Picornaviridae [23], they belong to different genera as Coxsackie virus B is an Enterovirus, and EMCV, a Cardiovirus. The Picornaviridae family is one of the largest and most diverse families of RNA viruses, and includes etiological agents that are responsible for a wide variety of human and animal diseases [23].
The picornavirions are small, non-enveloped and spherical, with a diameter of about 20 to 30 nm. Picornaviruses have a single-stranded, positive sense RNA genome that is between 7.2 and 9 kb in length and contain a large opening reading frame (ORF). The ORF encodes for a polyprotein that comprises both non-structural and structural elements divided into three primary precursor molecules, namely P1, P2 and P3, encoding for 11 distinct proteins. The structural proteins VP4, VP2, VP3 and VP1 make up the viral capsid and are encoded in the P1 region towards the 5′-end of the genome. Non-structural proteins are derived from the P2 and P3 regions and are encoded towards the 3′-end of the genome, the largest of which is the RNA-dependent RNA polymerase (3Dpol). In addition, cardioviruses code for an L (Leader) protein at the N-terminus of their polyproteins. The genomic RNA also contains a highly structured 5′-UTR (untranslated region) that includes an internal ribosome entry site (IRES) from which viral protein translation is initiated in a cap-independent manner. The shorter 3′-UTR terminates with a heterogeneous poly(A) tail that is known to be involved in the binding process of the viral replicase complex.
In this paper, we described the cell culture isolation and identification of two viral isolates obtained postmortem from the two juvenile Bornean orang utans. The viral genomes of the two isolates were >95% sequenced, with the complete genome sequence lacking about 200 nucleotides (nt) from the 5′ end. Phylogenetic and sequence analyses suggested that the newly isolated viruses are highly divergent variants of EMCV and possibly a new serotype of the virus. This is the first report of EMCV infection in Singapore and South East Asia and its molecular characterization.

Virus isolation and preliminary diagnosis
Virus isolates, designated Sing-M105-02 and Sing-M100-02, were successfully isolated in African Green Monkey kidney (Vero) cells from the heart tissues of the first juvenile orang utan, and lung tissues of the second autopsied juvenile orang utan, respectively. Cytopathic effect (CPE) was observed in Vero cells inoculated with the heart homogenate, at 2 days post-infection (pi). Infected-Vero cells were distorted from its normal shape, became granular and finally rounded-up (results not shown). Similar observations of CPE were made with the inoculation of lung homogenate at 3 days pi. Uninfected Vero cells did not show any signs of CPE. Both virus isolates were also inoculated into Mardin Darby Canine Kidney (MDCK) cells (data not shown). At 1 to 2 days pi, CPE in virus infected-MDCK cells was characterized by rounded and wrinkled appearances.
The viruses in infected-Vero tissue culture fluid were negatively stained and examined with electron microscopy. The electron micrographs revealed numerous small round virus particles with smooth appearances in the range of 20 to 30 nm for both Sing-M105-02 and Sing-M100-02 isolates ( Figure 1). The morphology and size observed are consistent with virus members belonging to the family of Picornaviridae [23].
A second passage of the virus isolates in Vero cells resulted in rapid and advanced CPE a day after inoculation. Sing-M105-02 ( Figure 2a) and Sing-M100-02 (Figure 2b) virus-infected Vero cells reacted to polyclonal antiserum raised to EMCV with specific green and yellow fluorescence staining in the cytoplasm. The perinuclear regions were stained more intensely. These results were consistent with observations made for EMCV(Aust)-infected Vero cells (Figure 2c), indicating that the Sing-M105-02 and Sing-M100-02 viruses could be EMCV. Cytoplasmic areas appeared red with counter stain Evans Blue, and no green fluorescence was observed in uninfected Vero cells (Figure 2d).

Virus identification and full-genome molecular characterization
Viral RNA, extracted from the virus-infected tissue culture fluid, was subjected to RT-PCR analysis using one set of primers that were specific for enterovirus detection [24,25], while another two sets of primers specific for the detection of EMCV were also used [13,26] (Table 1). No products were detected from the enterovirus RT-PCR, while approximately 855 bp and 286 bp products were observed with the VP3/VP1 and 3Dpol primer sets respectively for both Sing-M105-02 and Sing-M100-02 viruses (results not shown). The sizes were consistent with

Sing-M105-02 virus
Sing-M100-02 virus 100 nm  reports using the same primer sets [26][27][28]. No results were observed on RT-PCR performed on sample processed from uninfected Vero cells or when water was used to replace the viral RNAs. All the PCR assays were repeated at low stringency so as to detect potentially related viruses. The RT-PCR results indicated that Sing-M105-02 and Sing-M100-02 isolates belong to Cardiovirus genus and to EMCV species. We managed to sequence more than 95% of the full virus genome (i.e.~7.6 kb) of both Sing-M105-02 and Sing-M100-02 isolates by primer walking ( Table 2). This included the entire ORF and 3′-UTR of the virus. Despite  EMCV-P2 GGTGAGAGCAAGCCTCGCAAAGACAG 7370 Forward [13] numerous attempts to amplify out the 5′-UTR region, we were unable to obtain about 200 nt of virus sequence from the 5′ end of the genome, including the poly(C) tract. This difficulty in sequencing the 5′-UTR region was also reported for another EMCV strain [29], possibly due to the high GC content and complex secondary structure in this region. The Sing-M105-02 and Sing-M100-02 isolates showed a high sequence identity of 99.9% in the entire ORF (Tables 3 and 4), indicating that both viruses represent separate isolates of a single virus strain, probably originating from a similar source.

Sequence comparison with existing EMCV strains
To investigate the genetic relationship of Sing-M105-02 and Sing-M100-02 isolates with existing cardioviruses (EMCV and Theiloviruses), we performed sequence comparison for the entire ORF, VP1 capsid and 3Dpol gene regions. Table 3 reflects the percent nt and deduced  Table 4 compares the relationship of the Sing-M105-02 and Sing-M100-02 isolates against selected fully sequenced EMCV strains in the different coding regions across the entire ORF (namely L, P1, P2 and P3), and also the 3′-UTR region. We note that nt sequence identities were all below 80% in the P1, P2 and P3 regions between our isolates and existing EMCV strains, once again demonstrating the high divergence of the Sing-M105-02 and Sing-M100-02 isolates. The aa sequence alignment of the P1 structural region was next examined and we observed significant differences in the BC-loop and loop I regions in the VP1 capsid of the Sing-M105-02 and Sing-M100-02 isolates with other EMCV strains (see boxed region in Figure 3). The BCloop and loop I regions of the VP1 capsid protein are two putative neutralizing antigenic sites proposed in the EMCV Mengo virus strain [30,31]. Taken together, our sequence data indicated that the Sing-M105-02 and Sing-M100-02 isolates have a highly divergent nature compared to any known strains of EMCV.

Phylogenetic analysis
We included fully sequenced EMCV and Theilovirus strains in our phylogenetic trees to demonstrate the divergence of the 2 species in the genus Cardiovirus, as well as to emphasize the diversity and clustering of various EMCV strains (Figure 4). Our phylogenetic analyses showed branching of the EMCV strains into four main lineages, A, B, C and D, at the nt level in the VP1 and 3Dpol regions, as well as for the entire ORF ( Figure 4). Lineages A, B and C concurred with similar clusterings of EMCV strains as described previously [10,26,32]. However, the Sing-M105-02 and Sing-M100-02 viruses clustered separately as a distinct group by themselves in lineage D (Figure 4). A recent report of a second serotype of EMCV, isolate RD1338 [29], also clustered by itself in lineage E, highly distinct from the Sing-M105-02 and Sing-M100-02 viruses (lineage D) and other EMCV in lineages A, B and C. Our phylogenetic studies further demonstrated that Sing-M105-02 and Sing-M100 viruses are highly divergent strains of EMCV.

Serological survey of the zoo animals
In order to investigate whether other zoo animals had been exposed to the EMCV variant, a serological survey was conducted. Sera were collected from a variety of zoo animals, in particular animals residing near the orang utans, and titrated for neutralizing antibodies to both Sing-M105-02 virus and EMCV(Aust). A selection of capybaras (80%), chimpanzees (31%) and orang utans (50%) exhibited neutralizing antibody titers to Sing-M105-02 virus, strongly indicating that they had been exposed to the virus (Table 5). In total, 22 of 122 zoological animals (18%) exhibited neutralizing antibody titers to the Sing-M105-02 isolate. All of the chimpanzees, except one, and an orang utan also had high neutralizing antibody titer to EMCV(Aust) (results not shown). Our serological survey indicated that there may have been two strains of EMCV circulating among the zoo animals. One of the viruses could be an EMCV, closely related to EMCV (Aust), although no attempts at virus isolation were made from these exposed animals. The second virus could be a variant of the EMCV, similar to Sing-M105-02 and Sing-M100-02.
As previous reports have implicated rodents in the initiation of EMCV outbreaks, a rodent control program The proposed EMCV lineage is also indicated. The virus strains names are shown followed by the GenBank accession number in brackets. NA indicates sequence data not available.
was initiated at the Singapore Zoological Gardens as a control measure for the management of the myocarditis outbreak [33]. However, only two rats were caught, and these were bled. The two rats were sampled and demonstrated high titers of neutralizing antibodies to Sing-M105-02 virus ( Table 5); one of the rats had neutralizing antibody to both Sing-M105-02 virus and EMCV(Aust) (results not shown). Our data provided evidence that the rats were exposed to either the Sing-M105-02 virus, or both the Sing-M105-02 virus and another EMCV.

Discussion
The virus outbreak in the Singapore Zoological Gardens took place over an estimated period of 8 months, and occurred in orang utans of various ages, as well as sporadically among other zoo animals. This is the first report of an EMCV infection in Singapore. There were altogether four deaths in the orang utan cohort, caused by acute myocarditis. Other animals, two zebras, one spider monkey, one guanaco, and two capybaras also died of similar clinical signs. One adult orang utan displayed respiratory problems but recovered with intensive supportive therapy. From the two juvenile orang utans that died, viruses were isolated successfully from the heart (designated Sing-M105-02) and lung (designated Sing-M100-02) homogenates in Vero cells. Virus isolation was not performed on the two autopsied adult orang utans. The viral agents isolated in Vero cells were examined under transmission electron microscopy and the viruses showed similar morphology to the picornaviruses (Figure 1). In addition, virus-infected Vero cells showed positive immunofluorescence to polyclonal antiserum raised to EMCV (Figure 2).   Our preliminary characterization of the Sing-M105-02 and Sing-M100-02 viruses indicated that the etiological agent in the Singapore zoo outbreak was an EMCV. EMCV-specific primers targeting the VP3/VP1 capsid and 3Dpol genes were able to PCR-amplify products from the viral isolates. These primer sets had been tested on a wide variety of EMCV isolates from different geographic locations, different host species, and at different times [26][27][28]. Structural proteins such as the VP1 capsid, have been shown to be suitable for phylogenetic analyses within different genera and species of picornaviruses [24,[33][34][35][36][37][38][39][40][41]. Furthermore, analysis using the VP1 capsid could lead to the discovery of a novel variant such as the discovery reported for Hepatitis A [42]. Analysis at the VP1 capsid region revealed a maximum of 38.8% nt and 33.9% aa divergence between the Sing-M105-02 and Sing-M100-02 isolates and existing EMCV strains (Table 3). In addition, aa alignment of Sing-M105-02 and Sing-M100-02 against other EMCV strains showed considerable variation in two putative neutralization antigenic sites in VP1 (Figure 3), further supporting our hypothesis that these viral isolates can represent a new divergent group of EMCV and possibly a new serotype in the EMCV species.
The picornaviral polymerases or 3Dpol genes, on the other hand, are well conserved during evolution, and have been used as a common marker for comparisons between genera in phylogenetic studies [43][44][45]. They can also be used to determine the degree of variation between the newly isolated viruses with existing EMCV strains. For example in Foot-and-Mouth-Disease viruses, high sequence divergence found in the 3Dpol regions is an indication of the virus strains diverging from the same isolate [46]. Sequence comparison for Sing-M105-02 and Sing-M100-02 isolates with existing EMCV strains showed a maximum divergence of 24.6% and 11.3% at the nt and aa levels respectively ( Table 3), suggesting that they are indeed highly divergent variants of EMCV.
In addition, the phylogenetic analyses based on complete virus gene sequences (as compared to partial gene sequences) afforded us increased confidence that the Sing-M105-02 and Sing-M100-02 isolates constituted a divergent group of EMCV variants, as they clustered distinctly away from the other lineages of EMCV strains (Figure 4). This high degree of divergence within the EMCV species is perhaps not surprisingly as viruses with RNA genomes are generally reported to have a high rate of mutation [47][48][49].
Our serological survey of the animal population further suggested that there were at least two virus strains circulating in the Singapore Zoo; the first, an EMCV (Aust)-like virus, and the second, a variant of EMCV. High seropositive rates for the EMCV variant was found in capybaras, chimpanzees, orang utans and two rats caught in baited traps. In addition, neutralizing antibodies to EMCV(Aust) were also found in all except one of the chimpanzee, one orang utan and one of the rats, indicating prior exposure to single or dual infections with the viruses. Although no attempt was made to recover the virus from these animals, we postulated that the EMCV(Aust)-like virus could have been the parental   (18) virus from which Sing-M105-02 and Sing-M100-02 viruses gradually evolved. Rats have been considered to be the primary host reservoir and disseminators of EMCV, and they have been implicated in other zoological outbreaks of EMCV [2][3][4]7,9]. The rats could also have initiated the sporadic infections, by spreading both the parental and its variant virus in the zoo. However, this hypothesis remains to be investigated.

Conclusions
The etiological agent responsible for the outbreak for 2 of the 4 cases of fatal myocarditis in orang utans at the Singapore Zoological Gardens was an EMCV. High divergence in nt sequences of the VP1 capsid and 3Dpol genes as compared to other known EMCV strains, indicated that the outbreak was caused by a highly divergent variant of EMCV and possibly a new serotype of the virus. This is the first report of an EMCV infection in Singapore and South East Asia.

Virus isolation and propagation
Heart and lung tissues were collected from the postmortem analyses of two juvenile orang utans. The tissues were homogenised in Hank's balanced salt solution and the cellular debris cleared by low-speed centrifugation at about 3,000 to 4,000 x g. The homogenates were filtered with a 45 μm pore size and the clarified homogenate inoculated onto monolayers of African Green Monkey kidney (Vero) cells (ATCC-CCL81). The infection was carried out at 37°C and infected cells were observed for cytopathic cell effect (CPE) under a light microscope. The viruses isolated from the heart and lung tissues were designated as Sing-M105-02 and Sing-M100-02, respectively. Sing-M105-02, Sing-M100-02 and EMCV(Aust) were propagated in Vero cells in minimal Eagle's medium (MEM) (Invitrogen) supplemented with 5% fetal calf serum (FCS) (Hyclone, Logan) at 37°C, 5% CO 2 for 3 days. The culture supernatant was then harvested, centrifuged at 4,000 x g and stored at -80°C.
EMCV(Aust) strain was imported into Singapore from Dr Peter Kirkland, MacArthur Institute, Australia. The EMCV(Aust) strain was isolated from a rodent in an Australian zoo as reported in Reddacliff et al., [3]. Virus isolates and EMCV(Aust) were also inoculated into Mardin Darby Canine Kidney (MDCK) cells (ATCC-CCL34) and CPE observed over a few days.

Virus neutralization
The antiserum was serially diluted in two-folds and mixed with a known titer of either Sing-M105-02 or EMCV (Aust) in equal volume for 1 hr at 37°C. The antiserumvirus mixtures were next inoculated onto Vero cells seeded in a 96 well formant and examined daily for signs or absence of CPE.

Electron microscopy
Virus-infected tissue culture fluid was clarified by lowspeed centrifugation at 2,000 to 3,000 x g for 10 min.
The virus was pelleted from the supernatant by ultracentrifugation at 159,000 x g for 1.5 h at 4°C and resuspended in PBS [20]. Drops of virus suspension were placed on carbon-formvar coated grids. Negative staining was performed with 2% phosphotungstic acid, pH 7.4, for 1 min. Grids were examined with a JEOL-1010 electron microscope.

Immunofluorescence assay
Virus-infected cells were scraped off from a 25 cm 2 tissue culture flask at 24 h pi and washed in phosphate-buffered saline (PBS), pH 7.4, by low speed centrifugation at 2,000 to 3,000 x g for 10 min. The cell pellet was resuspended in 1 ml of PBS and 5 μl of cell suspension was distributed into each well of a Teflon-coated slide and air-dried. Cells were fixed in cold acetone for 10 min and virus antigens detected by immunofluorescence using polyclonal anti-EMCV raised in pig as the primary antibody (a gift from Dr Peter Kirkland, MacArthur Institute, Australia). The secondary antibody was goat anti-porcine IgG, conjugated to fluorescein isothiocynate (Chemicon). Evans Blue (Sigma Diagnostics) was used as counter-stain and the slide examined with a Olympus DX51 microscope.

Primers for virus identification and sequencing
For the initial determination of the virus genus and species, we used a total of 4 sets of PCR primers. Two published primer sets, RotbartFor/RotbartRev and EV1/EV2, both for the detection of enterovirus [24,25], and another two primer sets, P9/P10 (targeting the capsid VP3/VP1 gene) and P1/P2 (targeting the 3Dpol gene), for the detection of EMCV were selected [13]. The sequences and references of these primers sets are described in Table 1. The primers used for subsequent full genome assembly and sequencing were designed by primer walking, and are shown in Table 2.
RNA extraction, cDNA synthesis and PCR Viral RNA was extracted from 100 μl of virus-infected Vero cell culture supernatant using the RNeasy Mini Kit (Qiagen) as per the manufacturer's instructions. RNA was eluted in 50 μl of RNase-free water. For initial virus identification, the first strand cDNA was synthesized with gene-specific primers (2 μM) ( Table 1) using the Superscript III First Strand Synthesis System for RT-PCR (Invitrogen, USA) as per the manufacturer's protocol. For virus genome amplification and sequencing, the first strand cDNA was synthesized using using either random hexamers (50 ng/ μl) or gene-specific primers (2 μM) ( Table 2) using the Superscript III First Strand Synthesis System for RT-PCR (Invitrogen, USA) as per the manufacturer's protocol. For PCR, 2 μl of each viral cDNA sample was amplified in a 50 μl reaction mixture (33.5 μl RNase-free water, 5 μl 10x High Fidelity (HiFi) PCR Buffer, 5 μl 2 μM dNTPs, 2 μl 50 mM MgSO 4 , 1 μl 10 μM forward primer, 1 μl 10 μM reverse primer, 0.5 μl 50 U/μl HiFi Platinum Taq). Thermal cycling conditions included initial denaturation at 95°C for 5 min; followed by 45 cycles of 95°C for 30 s, 45°C for 30 s, 68°C for 2 min; and a final extension step of 68°C for 10 min. Touchdown PCR was also performed in several instances to improve PCR amplification specificity.
To obtain the 3′ end sequence of the virus, the GeneRacer kit and core module primers (Invitrogen, USA) were employed as per the manufacturer's instructions. We were unable to obtain about 200 nt of the virus 5′ UTR sequence despite numerous attempts using the GeneRacer kit, other reverse transcriptases, DNA polymerases, PCR additives (such as DMSO and betaine), and various PCR thermocycling conditions.

Sequencing and phylogenetic analysis
PCR products were resolved on 1-2% agarose gels (depending on amplicon size), and the subsequent PCR clean-up and sequencing of the amplicons were performed by AITbiotech (Singapore) using BigDye chemistry (Life Technologies, USA). In brief, the genomes of EMCV isolates Sing-M100-02 and Sing-M105-02 were determined from a series of PCR amplicons that contained overlapping cDNA fragments. Each gene region was sequenced at least two times. Raw sequences were first assembled and aligned with the reference EMCV Mengo M virus genome (GenBank accession number L22089) using the SeqMan Pro software found in the Lasergene Version 8 (DNASTAR, USA).
Final nt sequences were aligned against other strains of EMCV and picornaviruses using the Clustal W method as implemented in the MegAlign program (DNASTAR, USA). Percent identity and phylogenetic trees were constructed for the entire ORF as well as for individual genes coding for the various proteins at the nt level and neighborjoining trees were constructed after alignment by the same program. The consistency and robustness of tree topology were bootstrapped using 1000 trials and a seed of 111. The TreeView X 1.6.6 program [50] was used to visualize trees.

Nucleotide sequence accession numbers
The partial sequence data for the VP3/VP1 genes have been deposited in GenBank under the accession numbers AF525466 for Sing-M105-02 and AY162280 for Sing-M100-02. For the partial 3Dpol genes, the nt sequences were deposited under the accession numbers AF510055 for Sing-M105-02, and AY162279 for Sing-M100-02. The full virus genome sequences (excluding about 200 nt from the 5′ end) have been deposited in GenBank under accession numbers KC310737 for Sing-M100-02, and KC310738 for Sing-M105-02.