- Open Access
A highly divergent Encephalomyocarditis virus isolated from nonhuman primates in Singapore
Virology Journalvolume 10, Article number: 248 (2013)
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.
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.
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.
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.
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 , baboons , chimpanzees , lemurs , rhesus macques , and orang utans [3, 6]. EMCV infection was also implicated in the deaths of various animals in Audubon Park Zoo, New Orleans, in 1985 ; and in elephants in a Florida zoo in 1997 . 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 ; 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 . 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 . EMCV has also been recognized as a porcine pathogen, with EMCV infections in European pigs associated with sudden deaths and reproductive failure [12–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 .
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 . 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, 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 .
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 post-mortem 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.
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 reports using the same primer sets [26–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 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 , 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 amino acid (aa) identities of the entire ORF, VP1 capsid and 3Dpol gene regions for Sing-M105-02 and Sing-M100-02 isolates compared with other fully sequenced Cardiovirus strains. In the conserved 3Dpol region, the extent of sequence identities between the Sing-M100-02 and Sing-M105-02 isolates and other EMCV strains ranged from 76.4% to 80.9% (nt level) and 79.4% to 90.5% (aa level). The lowest sequence relationships were observed in the variable VP1 capsid region, with percent identities ranging from 61.2% to 73.6% (nt level) and 66.1% to 87.0% (aa level). 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 BC-loop 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.
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 , 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 was initiated at the Singapore Zoological Gardens as a control measure for the management of the myocarditis outbreak . 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.
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–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–41]. Furthermore, analysis using the VP1 capsid could lead to the discovery of a novel variant such as the discovery reported for Hepatitis A . 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–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 . 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–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 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–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.
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.
Materials and methods
Post-mortem heart and lung tissues were obtained from the two juvenile orang utans. In addition, a total of 122 serum samples from various zoo animals were collected for sero-epidemiological study. These consisted of 6 agile gibbons, 1 cape hunting dog, 1 capuchin monkey, 10 capybaras, 13 chimpanzees, 3 colobus monkeys, 1 coloque monkey, 21 gibbons, 3 howler monkeys, 2 Japanese monkeys, 2 komodo dragons, 19 macaques, 2 mandrills, 14 orang utans, 1 patas monkey, 2 pig-tailed macaques, 1 raccoon, 2 lemurs, 3 siamangs, 1 slow loris, 5 spider monkeys and 4 white nose guenons. Two rats were also captured in a baited trap.
Virus isolation and propagation
Heart and lung tissues were collected from the post-mortem 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% CO2 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., . Virus isolates and EMCV(Aust) were also inoculated into Mardin Darby Canine Kidney (MDCK) cells (ATCC-CCL34) and CPE observed over a few days.
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 antiserum-virus mixtures were next inoculated onto Vero cells seeded in a 96 well formant and examined daily for signs or absence of CPE.
Virus-infected tissue culture fluid was clarified by low-speed 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 . 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.
Virus-infected cells were scraped off from a 25 cm2 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 . 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 MgSO4, 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 neighbor-joining 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  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.
Jones P, Cordonnier N, Mahamba C, Burt FJ, Rakotovao F, Swanepoel R, André C, Dauger S, Bakkali Kassimi L: Encephalomyocarditis virus mortality in semi-wild bonobos (Pan panicus). J Med Primatol 2011, 40: 157-163. 10.1111/j.1600-0684.2010.00464.x
Hubbard GB, Soike KF, Butler TM, Carey KD, Davis H, Butcher WI, Gauntt CJ: An Encephalomyocarditis virus epizootic in a baboon colony. Lab Anim Sci 1992, 42: 233-239.
Reddacliff LA, Kirkland PD, Hartley WJ, Reece RL: Encephalomyocarditis virus infections in an Australian zoo. J Zool Wild Med 1997, 28: 153-157.
Canelli E, Luppi A, Lavazza A, Lelli D, Sozzi E, Martin AM, Gelmetti D, Pascotto E, Sandri C, Magnone W, Cordioli P: Encephalomyocarditis virus infection in an Italian zoo. Virol J 2010, 7: 64. 10.1186/1743-422X-7-64
Masek-Hammerman K, Miller AD, Lin KC, MacKey J, Weissenbock H, Gierbolini L, Burgos A, Perez H, Mansfield KG: Epizootic myocarditis associated with Encephalomyocarditis virus in a group of rhesus macaques (macaca mulatta). Vet Path 2012, 49: 386-392. 10.1177/0300985811409254
Citino SB, Gaskin HJH, Wickham DJ: Fatal Encephalomyocarditis virus infection in Sumatran orangutan ( Pongo pygmaeus abelii ). J Zoo Wildl Med 1998, 19: 214-218.
Wells SK, Gutter AE, Soike KF, Baskin GB: Encephalomyocarditis virus: epizootic in a zoological collection. J Zool Wild Med 1989, 20: 291-296.
Simpson CF, Lewis AL, Gaskin JM: Encephalomyocarditis virus infection of captive elephants. J Am Vet Med Assoc 1977, 171: 902-905.
Grobler DG, Raath JP, Braack LE, Keet DF, Gerdes GH, Barnard BJ, Kriek NP, Jardine J, Swanepoel R: An outbreak of Encephalomyocarditis-virus infection in free-ranging African elephants in the Kruger National Park. Onderstepoort J Vet Res 1995, 62: 97-108.
van Sandwyk JH, Bennett NC, Swanepoel R, Bastos AD: Retrospective genetic characterisation of Encephalomyocarditis viruses from African elephant and swine recovers two distinct lineages in South Africa. Vet Microbiol 2013, 162: 23-31. 10.1016/j.vetmic.2012.08.008
Krylova RI, Dzhikidze EK: Encephalomyocarditis in monkeys. Bull Exp Biol Med 2005, 139: 355-359. 10.1007/s10517-005-0293-3
Joo HS: Encephalomyocarditis virus. In Diseases of Swine. Edited by: Leman AD, Straw BE, Mengeling WL, Mengeling S, D’Allaire S, Taylor DJ. London: Wolfe Publishing; 1992:257-262.
Koenen F, Vanderhallen H, Papadopoulos O, Billinis CP, Paschaleri-Papadopoulou E, Brocchi E, De Simone F, Carra E, Knowles NJ: Comparison of the pathogenic, antigenic and molecular characteristics of two Encephalomyocarditis virus (EMCV) isolates from Belgium and Greece. Res Vet Sci 1997, 62: 239-244. 10.1016/S0034-5288(97)90197-1
Paschaleri-Papadopoulou E, Axiotis I, Laspidis C: Encephalomyocarditis of swine in Greece. Vet Rec 1990, 124: 364-365.
Sidoli L, Barigazzi G, Foni E, Marcato PS, Barbieri G: Encephalomyocarditis (EMC) due to Cardiovirus in Po Valley swines: preliminary observations, clinical aspects, virus isolation, characterization and experimental transmission. Selozoine Vet 1989, 30: 249-260.
Maurice H, Nielen M, Brocchi E, Nowotny N, Kassimi LB, Billinis C, Loukaides P, O’Hara RS, Koenen F: The occurrence of Encephalomyocarditis virus (EMCV) in European pigs from 1990 to 2001. Epidemiol Infect 2005, 133: 547-557. 10.1017/S0950268804003668
Park NY, Ri CY, Chung CY, Kee HY, Bae SY: Pathological findings on Encephalomyocarditis virus infections of swine in Korea. Korea J Vet Res 1992, 32: 99-109.
Hu DG, Chan CH, Shieh WY, Chang CS Research report – animal industry research institute, Taiwan Sugar Corporation. Epidemiological studies of swine Encephalomyocarditis virus infection 1993, 149-154. 1992/93
Miyagi J, Tsuhako K, Kinjo T, Iwamasa T, Kamada Y, Kinju T, Koyanagi Y: Coxsackievirus B4 myocarditis in an orangutan. Vet Path 1999, 36: 452-456. 10.1354/vp.36-5-452
Gajdusek DC: Encephalomyocarditis virus infection in childhood. Pedriatrics 1955, 16: 902-907.
Dick GWA, Best AM, Haddow AJ, Smithburn KC: Mengo encephalomyelitis, a hitherto unknown virus affecting man. Lancet 1948, 252: 286-289. 10.1016/S0140-6736(48)90652-7
Oberste MS, Gotuzzo E, Blair P, Nix WA, Ksiazek TG, Comer JA, Rollin P, Goldsmith CS, Olson J, Kochel TJ: Human febrile illness caused by Encephalomyocarditis virus infection, Peru. Emerg Infect Dis 2009, 15: 640-646. 10.3201/eid1504.081428
Knowles NJ, Hovi T, Hyypiä T, King AMQ, Lindberg AM, Pallansch MA, Palmenberg AC, Simmonds P, Skern T, Stanway G, Yamashita T, Zell R: Picornaviridae. In Virus taxonomy: classification and nomenclature of viruses: ninth report of the international committee on taxonomy of viruses. Edited by: King AMQ, Adams MJ, Carstens EB, Lefkowitz EJ. San Diego: Elsevier; 2012:855-880.
Oberste MS, Maher , Flemister MR, Marchetti G, Kilpatrick DR, Pallansch MA: Comparison of classic and molecular approaches for the identification of untypeable enteroviruses. J Clin Microbiol 2000, 38: 1170-1174.
Rotbart AR: PCR amplification of enteroviruses. In PCR protocols: a guide to methods and amplifications. Edited by: Innis MA, Gelfand DH, Sninsky JJ, White TJ. San Diego: Academic Press; 1990:372-377.
Koenen F, Vanderhallen H, Dickinson ND, Knowles NJ: Phylogenetic analysis of European Encephalomyocarditis viruses: comparison of two genomic regions. Arch Virol 1999, 144: 893-903. 10.1007/s007050050554
Vanderhallen H, Koenen F: Rapid diagnosis of Encephalomyocarditis virus infections in pigs using a reverse transcription-polymerase chain reaction. J Virol Met 1997, 66: 83-89. 10.1016/S0166-0934(97)02214-3
Vanderhallen H, Koenen F: Identification of Encephalomyocarditis virus in clinical samples by reverse transcription-PCR followed by genetic typing using sequence analysis. J Clin Microbiol 1998, 36: 3463-3467.
Philipps A, Dauber M, Groth M, Schirrmeier H, Platzer M, Krumbholz A, Wutzler P, Zell R: Isolation and molecular characterization of a second serotype of the Encephalomyocarditis virus. Vet Microbiol 2012, 161: 49-57. 10.1016/j.vetmic.2012.07.006
Luo M, Vriend G, Kamer G, Minor I, Arnold E, Rossmann MG, Boege U, Scraba DG, Duke GM, Palmenberg AC: The atomic structure of Mengo virus at 3.0 A resolution. Science 1987, 235: 182-191. 10.1126/science.3026048
Boege U, Kobasa D, Onodera S, Parks GD, Palmenberg AC, Scraba DG: Characterization of Mengo virus neutralization epitopes. Virology 1991, 181: 1-13. 10.1016/0042-6822(91)90464-M
Denis P, Liebig HD, Nowotny N, Billinis C, Papadopoulos O, O’Hara RS, Knowles NJ, Koenen F: Genetic Encephalomyocarditis virus (EMCV) isolates. Vet Microbiol 2006, 113: 1-12. 10.1016/j.vetmic.2005.10.032
Zimmerman JJ: Encephalomyocarditis. In CRC handbook series in zoonoses. Section B: viral zoonoses. Edited by: Beran GW, Steele JH. Boca Raton: CRC Press; 1994:423-436.
Brown BA, Oberste MS, Alexander JP Jr, Kennett ML, Pallansch MA: Molecular epidemiology and evolution of enterovirus 71 strains isolated from 1970 to 1998. J Virol 1999, 73: 9969-9975.
Haydon DT, Bastos AD, Knowles NJ, Samuel AR: Evidence for positive selection in foot-and-mouth disease virus capsid genes from field isolates. Genetics 2001, 157: 7-15.
McMinn P, Lindsay K, Perera D, Chan HM, Chan KP, Cardosa MJ: Phylogenetic analysis of enterovirus 71 strains isolated during linked epidemics in Malaysia, Singapore, and Western Australia. J Virol 2001, 75: 7732-7738. 10.1128/JVI.75.16.7732-7738.2001
Mulders MN, Salminen M, Kalkkinen N, Hovi T: Molecular epidemiology of coxsackievirus B4 and disclosure of the correct VP1/2A(pro) cleavage site: evidence for high genomic diversity and long-term endemicity of distinct genotypes. J Gen Virol 2000, 81: 803-812.
Norder H, Bjerregaard L, Magnius LO: Homotypic echoviruses share aminoterminal VP1 sequence homology applicable for typing. J Med Virol 2001, 63: 35-44. 10.1002/1096-9071(200101)63:1<35::AID-JMV1005>3.0.CO;2-Q
Oberste MS, Maher K, Kilpatrick DR, Pallansch MA: Molecular evolution of the human enteroviruses: correlation of serotype with VP1 sequence and application to Picornavirus classification. J Virol 1999, 73: 1941-1948.
Oberste MS, Schnurr D, Maher K, al-Busaidy S, Pallansch MA: Molecular identification of new picornaviruses and characterization of a proposed enterovirus 73 serotype. J Gen Virol 2001, 82: 409-416.
Santti J, Harvala H, Kinnunen L, Hyypia T: Molecular epidemiology and evolution of coxsackievirus A9. J Gen Virol 2000, 81: 1361-1372.
Costa-Mattioli M, Cristina J, Romero H, Perez-Bercof R, Casane D, Colina R, Garcia L, Vega I, Glikman G, Romanowsky V, Castello A, Nicand E, Gassin M, Billaudel S, Ferre V: Molecular evolution of hepatitis A virus: a new classification based on the complete VP1 protein. J Virol 2002, 76: 9516-9525. 10.1128/JVI.76.18.9516-9525.2002
Koonin EW: The phylogeny of RNA-dependent RNA polymerases of positive-strand RNA viruses. J Gen Virol 1991, 72: 2197-2206. 10.1099/0022-1317-72-9-2197
Rodrigo MJ, Dopazo J: Evolutionary analysis of the picornavirus family. J Mol Evol 1995, 40: 362-371. 10.1007/BF00164022
Palmberg AC, Kirby EM, Janda MR, Drake NL, Duke GM, Potratz KF, Collett MS: The nucleotide and deduced amino acid sequences of the Encephalomyocarditis viral polyprotein coding region. Nucleic Acids Res 1984, 12: 2969-2985. 10.1093/nar/12.6.2969
Villaverde A, Martinez-Salas E, Domingo E: 3D gene of foot-and-mouth disease virus: Conservation by convergence of average sequences. J Mol Biol 1988, 204: 771-776. 10.1016/0022-2836(88)90367-1
Domingo E, Holland JJ: RNA virus mutations and fitness for survival. Ann Rev Microbiol 1997, 51: 151-178. 10.1146/annurev.micro.51.1.151
Domingo E: Viruses at the edge of adaptation. Virology 2000, 270: 251-253. 10.1006/viro.2000.0320
Lauring AS, Andino R: Quasispecies theory and the behavior of RNA viruses. PLoS Pathog 2010, 6: e1001005. 10.1371/journal.ppat.1001005
Page RDM: TREEVIEW: an application to display phylogenetic trees on personal computers. Comp Appl Biosci 1996, 12: 357-358.
We thank Dr Cheng Wen-Haur from the Singapore Zoological Gardens, and Dr Paulo Martelli (now the head veterinarian of the Hong Kong Ocean Park), for their assistance in the coordination of sera collection from the zoo animals. We also thank Dr Peter Kirkland, from the MacArthur Institute, Australia, for his gifts of polyclonal antibodies to EMCV, as well as the strain of EMCV(Aust). This project was funded by the Ministry of Defence, Agri-Food and Veterinary Authority of Singapore, and National Environmental Agency, Singapore.
The authors declare that they have no competing interests.
DS-YY, JEL, JCWL, EASL and HCT performed the laboratory experiments. SO performed the clinical assessment. CJF performed the necropsy and sampling. CJF, Y-NL and M-LS performed the virus isolation. B-HT, M-LN, E-EO, K-PC, SO and CJF conceived the study. B-HT and DS-YY were involved in the data analysis and manuscript writing. All authors read and approved the final manuscript.