Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Discovery of a Bovine Enterovirus in Alpaca

  • Shasta D. McClenahan,

    Affiliation Center for Biologics Evaluation and Research, Food and Drug Administration, Bethesda, Maryland, United States of America

  • Gail Scherba,

    Affiliations Department of Pathobiology, University of Illinois, Urbana, Illinois, United States of America, Veterinary Diagnostic Laboratory, University of Illinois, Urbana, Illinois, United States of America

  • Luke Borst,

    Current address: Department of Population Health and Pathobiology, North Carolina State University, Raleigh, North Carolina, United States of America

    Affiliation Veterinary Diagnostic Laboratory, University of Illinois, Urbana, Illinois, United States of America

  • Richard L. Fredrickson,

    Affiliation Veterinary Diagnostic Laboratory, University of Illinois, Urbana, Illinois, United States of America

  • Philip R. Krause,

    Affiliation Center for Biologics Evaluation and Research, Food and Drug Administration, Bethesda, Maryland, United States of America

  • Christine Uhlenhaut

    uhlenhautc@rki.de

    Current address: Highly Pathogenic Viruses, Centre for Biological Threats and Special Pathogens, Robert Koch Institute, Berlin, Germany

    Affiliation Center for Biologics Evaluation and Research, Food and Drug Administration, Bethesda, Maryland, United States of America

Abstract

A cytopathic virus was isolated using Madin-Darby bovine kidney (MDBK) cells from lung tissue of alpaca that died of a severe respiratory infection. To identify the virus, the infected cell culture supernatant was enriched for virus particles and a generic, PCR-based method was used to amplify potential viral sequences. Genomic sequence data of the alpaca isolate was obtained and compared with sequences of known viruses. The new alpaca virus sequence was most similar to recently designated Enterovirus species F, previously bovine enterovirus (BEVs), viruses that are globally prevalent in cattle, although they appear not to cause significant disease. Because bovine enteroviruses have not been previously reported in U.S. alpaca, we suspect that this type of infection is fairly rare, and in this case appeared not to spread beyond the original outbreak. The capsid sequence of the detected virus had greatest homology to Enterovirus F type 1 (indicating that the virus should be considered a member of serotype 1), but the virus had greater homology in 2A protease sequence to type 3, suggesting that it may have been a recombinant. Identifying pathogens that infect a new host species for the first time can be challenging. As the disease in a new host species may be quite different from that in the original or natural host, the pathogen may not be suspected based on the clinical presentation, delaying diagnosis. Although this virus replicated in MDBK cells, existing standard culture and molecular methods could not identify it. In this case, a highly sensitive generic PCR-based pathogen-detection method was used to identify this pathogen.

Introduction

Alpaca (Vicugna pacos, also known as Lama guanicoe pacos) are domesticated members of the New World camelid species (Lamini), which also include guanaco (Lama guanicoe), vicuna (Vicugna vicugna), and llama (Lama glama). The natural habitat for alpaca is at high altitude (3500–5000 m) in South America (Peru, Ecuador, Bolivia, and Chile) where they are kept as livestock in herds and their fiber is used much like wool. Approximately 300,000 animals [1] are in the U.S. Compared to other livestock, e.g., about 96 million cattle [2], their number is still relatively small.

Previously reported viral infections in domestic alpaca include adenovirus, equine viral arteritis virus, rabies, bluetongue virus, foot-and-mouth disease virus, bovine respiratory syncytial virus, influenza A virus, rotavirus, orf virus, bovine papillomavirus, vesicular stomatitis virus, coronavirus, bovine parainfluenza-3 virus, West Nile virus, equine herpesvirus-1 [1], [3], [4] and bovine viral diarrhea virus [5][13]. Bovine enteroviruses (BEV) have not previously been reported to infect alpaca. The bovine enterovirus species previously contained two types, BEV-A and BEV-B [14], [15] although a new classification structure was ratified recently, redesignating these as species Enterovirus E (EV-E) and Enterovirus F (EV-F), respectively [14], [16]. Each of the new BEV species includes multiple serotypes, with EV-E comprising four described serotypes (previously A1–4, renamed E1–E4), and EV-F containing six reported serotypes (previously B1–6, renamed F1–F6).

Recently developed approaches to virus detection have the potential to further expand understanding of viral disease in animals, including alpaca. Many of these approaches are based on non-specific PCR amplification used in conjunction with standard or high-throughput sequencing to identify PCR products.

We utilized such a method [17][19] to investigate an outbreak of a respiratory infection in alpaca, identifying a bovine enterovirus (EV-F), named Enterovirus F, strain IL/Alpaca, after other techniques had failed to detect any pathogen.

Results

Four out of 32 alpaca in an Illinois herd ranging in age from 1.5 to 14 years of age died from an acute respiratory infection (with some evidence of systemic spread in two of the animals) of unknown etiology or origin. The other animals in the herd remained clinically healthy. Necropsy revealed grossly moderate acute diffuse interstitial pneumonia in all four animals and acute renal cortical infarcts in two of the alpaca. Microscopically, marked pulmonary congestion and edema were noted in all lungs, as well as moderate erosive gastritis, acute renal infarcts, mild esophageal erosion and ulceration with suppurative esophagitis in two of the alpaca. Quantitative RT-PCR for bovine viral diarrhea virus 1 and 2 failed to detect viral genomes. A cytopathic virus was isolated on subpassage from pulmonary tissue of one affected animal using MDBK cells. Cytopathic effect (CPE) was not observed in inoculated bovine turbinate, rabbit kidney or uninoculated cells, and therefore these isolation attempts were not pursued. FITC-conjugated fluorescent antibodies against several bovine viruses (adenovirus types 1 and 5; bluetongue; bovine viral diarrhea virus; coronavirus; herpesvirus types 1, 2, and 5; parainfluenza virus 3; parvovirus; reovirus, rotavirus and respiratory syncytial virus) failed to detect a virus in the infected cell cultures. Negative staining electron microscopy (EM) of frozen and thawed infected MDBK cell culture revealed the presence of numerous, uniformly shaped, non-enveloped virus particles approximately 25 to 30 nm in diameter (Figure 1).

thumbnail
Figure 1. Electron photomicrograph of the alpaca virus isolate.

Virus isolated from lung tissue and propagated in MDBK cell culture was imaged by negative staining EM. Virus particles are 25 to 30 nm in diameter. Image was taken at 100, 000×. Size bar is 50 nm.

https://doi.org/10.1371/journal.pone.0068777.g001

In order to identify the cytopathic virus isolated from the alpaca, a generic, degenerate oligonucleotide primer (DOP) PCR-based virus detection assay [17][19] was utilized. Infected and uninfected cell culture supernatants were enriched for viral capsids by nuclease digestion and ultracentrifugation. Extracted nucleic acids were subjected to reverse-transcription, amplified by DOP-PCR, and separated by agarose gel electrophoresis (Figure 2). The gel electrophoresis pattern of these amplified nucleic acids differed between infected and uninfected MDBK cells. Ten bands were excised each from the infected cell lane and from the uninfected cell lane, cloned, and sequenced.

thumbnail
Figure 2. DOP-PCR products obtained from cell culture supernatant with DOP-PCR.

cDNA was amplified by DOP-PCR and the obtained products were analyzed by gel electrophoresis. Marker (M): 100 bp ladder, Invitrogen, Lane 1- Infected MDBK cells, Lane 2- Uninfected MDBK cells, Lane 3- no template control. The arrowheads indicate PCR products that were sequenced and shown to be BEV-related.

https://doi.org/10.1371/journal.pone.0068777.g002

Sequencing of nucleic acid from the infected cell lane revealed 47 distinct products, encompassing regions with homology to approximately 46% of the EV-F genome, with sequences showing greatest homology to serotypes 1 and 3 (Figure 3). One sequence was of cellular origin due to residual MBDK cell DNA. Thirty-six distinct sequences obtained from the DOP-PCR amplicons of the negative control MDBK cells were consistent with amplification of MDBK DNA and with amplification of residual DNA in DOP-PCR reagents that we have observed previously, with no enterovirus-like sequences observed.

thumbnail
Figure 3. Diagram of BEV sequences identified by DOP-PCR, aligned with the EV-F genome.

Sequencing of the DOP-PCR amplicons identified 47 different amplicons. These sequences had greatest homology to serotypes 1–3 of the BEV-B species, and are displayed relative to the EV-F genome. These DOP-PCR products represented 46% of the alpaca-sourced virus genome.

https://doi.org/10.1371/journal.pone.0068777.g003

The remainder of the complete viral genome was identified by specific PCRs and RACE, based on primers designed from the already-obtained sequence and from BEV sequences in GenBank (Table 1). The complete 7433 bp genome for this virus, named Enterovirus F, strain IL/Alpaca, has been deposited in the GenBank database under accession KC748420.

thumbnail
Table 1. Primer sequences use for genomic sequencing and RACE PCRs of the alpaca-sourced virus.

https://doi.org/10.1371/journal.pone.0068777.t001

Based on recently changed nomenclature [14][16], the genome of the novel virus was most closely related to EV-F (previously BEV type B) serotypes 1 and 3, with homology to EV-F complete genome sequences ranging from 75–83%. Homology with EV-E sequences was 67–68% at the genome level. Since the capsid is used for typing picornaviruses, the virus identified in this study has to be considered as type 1.

In order to analyze the virus for potential recombination and to describe it more accurately, we performed more detailed phylogenetic analyses on several proteins of the novel virus, deduced from the translated nucleic acid sequences. Analysis of the full polyprotein and the individual capsid, 2A protease, 3C protease, and polymerase proteins of the alpaca-infecting virus relative to sequences of other representative enteroviruses from bovine EV-E (BEV-A serotypes 1–4) and EV-F (BEV-B serotypes 1–4), and sequences from three unclassified EV-F viruses [16], two from bovine sources (AY724744 and AY724745) [20], and one from a capped langur (JX538037) [21], possum, porcine (PEV), and human (HEV) hosts. These analyses revealed the alpaca virus to be most closely related to EV-F (Figures 4 to 8). Based on analysis of the full polyprotein, the alpaca-sourced virus clusters most closely with the EV-F, with homologies exceeding 85%, highest with serotype 1 viruses (Figure 4), and is more distantly related to the EV-F serotypes 2 and 3, followed by the EV-E species. The more diverse capsid protein (comprising the external surface of the virus) sequence of the alpaca-sourced virus was also most closely related to EV-F, serotype 1 (Figure 5) sharing 81% and 97% identity at the nucleotide and amino acid levels respectively. The amino acid homology with serotype 2 viruses was 86–87%, and 79% with serotype 3, and 78% with a serotype 4 possum isolate. As compared with serotype 1, the capsid sequence also was less similar to the partial capsid sequences of the unclassified EV-F viruses from bovine (AY424745) and capped langur species (JX538037), each with 87% amino acid identity, although the incomplete nature of these sequences makes it impossible to be certain of the degree of relatedness. Because the capsid gene is used for serotyping picornaviruses, this virus is thus considered a type 1. However, based on 2A protease (which cleaves the viral polypeptide into its individual components) sequences, the alpaca-sourced virus groups most closely with serotype 3 with 95% homology, followed by serotype 1 with 89% homology (Figure 6), indicating that the virus had attributes of type 3 and thus could have been a recombinant between types 1 and 3. The less diverse 3C protease (which the virus also uses to cleave the polypeptide into its individual components) of the enteroviruses groups the alpaca-sourced virus most closely with EV-F, serotypes 1 and 3 (Figure 7), with 97% amino acid identity. The gene for the polymerase enzyme (which the virus uses to transcribe its RNA after infection of a cell) is also highly conserved among the enteroviruses and cannot be used to clearly delineate a serotype for the alpaca-sourced virus, which still clusters most closely with the EV-F species (Figure 8), with amino acid identity greater than 97% for serotypes 1–3 and 94% for serotype 4.

thumbnail
Figure 4. Neighbor-joining phylogenetic tree of the deduced amino acid sequences from the complete enterovirus polyprotein.

Enteroviruses representing the bovine enterovirus (BEV) in species EV-E and EV-F, porcine enterovirus (PEV)/Enterovirus G, and human enterovirus (HEV)/Enterovirus A–D groups are included and species and serotypes are indicated. The amino acid sequences were aligned with the Clustal W program, and bootstrap confidence values were determined by 1000 replications. The scale bar represents the number of amino acid substitutions per site.

https://doi.org/10.1371/journal.pone.0068777.g004

thumbnail
Figure 5. Neighbor-joining phylogenetic tree of the deduced amino acid sequences from the capsid gene.

Enteroviruses representing the bovine enterovirus (BEV) in species EV-E and EV-F, porcine enterovirus (PEV)/Enterovirus G, and human enterovirus (HEV)/Enterovirus A–D groups are included and species and serotypes are indicated. The capsid sequences for the capped langur (JX538037) and unclassified AY24745 sequences (denoted with asterisks) are partial, so the correct placement of these sequences in these trees may change as more sequence data become available. The amino acid sequences were aligned with the Clustal W program, and bootstrap confidence values were determined by 1000 replications. The scale bar represents the number of amino acid substitutions per site.

https://doi.org/10.1371/journal.pone.0068777.g005

thumbnail
Figure 6. Neighbor-joining phylogenetic tree of the deduced amino acid sequences from the 2A protease gene.

Enteroviruses representing the bovine enterovirus (BEV) in species EV-E and EV-F, porcine enterovirus (PEV)/Enterovirus G, and human enterovirus (HEV)/Enterovirus A–D groups are included and species and serotypes are indicated. The amino acid sequences were aligned with the Clustal W program, and bootstrap confidence values were determined by 1000 replications. The scale bar represents the number of amino acid substitutions per site.

https://doi.org/10.1371/journal.pone.0068777.g006

thumbnail
Figure 7. Neighbor-joining phylogenetic tree of the deduced amino acid sequences from the 3C protease gene.

Enteroviruses representing the bovine enterovirus (BEV) in species EV-E and EV-F, porcine enterovirus (PEV)/Enterovirus G, and human enterovirus (HEV)/Enterovirus A–D groups are included and species and serotypes are indicated. The amino acid sequences were aligned with the Clustal W program, and bootstrap confidence values were determined by 1000 replications. The scale bar represents the number of amino acid substitutions per site.

https://doi.org/10.1371/journal.pone.0068777.g007

thumbnail
Figure 8. Neighbor-joining phylogenetic tree of the deduced amino acid sequences from the polymerase gene.

Enteroviruses representing the bovine enterovirus (BEV) in species EV-E and EV-F, porcine enterovirus (PEV)/Enterovirus G, and human enterovirus (HEV)/Enterovirus A–D groups are included and species and serotypes are indicated. The amino acid sequences were aligned with the Clustal W program, and bootstrap confidence values were determined by 1000 replications. The scale bar represents the number of amino acid substitutions per site.

https://doi.org/10.1371/journal.pone.0068777.g008

We also compared the 5′ untranslated region (UTR) of the alpaca virus with the bovine enteroviruses and found the greatest homology with EV-F strains, the highest with serotypes 1–3 and unclassified bovine sequence AY24744 at 87–90% homology. The alpaca virus homology with the EV-E 5′ UTR was 75–78%.

Several attempts to perform enterovirus-specific PCRs, using primers developed for the alpaca-source enterovirus and published EV-E and EV-F primers, were made on RNA extracted from paraffin-embedded lung tissues from the two of the alpaca (data not shown). In addition we also performed DOP-PCR on RNA extracted from these embedded tissues. Some non-specific PCR bands were evident, but sequencing of these PCR products revealed no enterovirus sequences. PCRs for the housekeeping genes β–actin and glyceraldehyde 3-phosphate dehydrogenase (GAPDH) were negative for some of the tissue samples, indicating that RNA quality was low in these fixed tissues.

Discussion

In this report, we describe an enterovirus that was isolated on subpassage from pulmonary tissue of an alpaca that died with evidence of respiratory and systemic infection. Using a universal virus detection assay, we identified a significant portion of the genome of this picornavirus with a single PCR. This finding is consistent with the EM data that visualized non-enveloped viral particles of approximately 25–30 nm in diameter. This is the first report of a BEV isolation from alpaca.

All four of the diseased animals had similar clinical symptomatology and had similar pulmonary histology on autopsy. Because EV-F was the only potential pathogen isolated from any of these animals, the alpaca-adapted EV is a potential cause of this syndrome, although these experiments clearly did not fulfill Koch's postulates and limitations in sensitivity of the other tests that were performed do not exclude the potential for other causes. Attempts to identify EV-F by PCR of paraffin-embedded pulmonary tissue samples obtained from these animals failed. This could be due to low copy number of EV-F RNA in the sections of the paraffin blocks that were examined or low stability of the EV-F RNA under the conditions of paraffin block storage.

Enteroviruses comprise one of the nine genera of picornaviruses; all of which include members that infect vertebrates. Picornaviridae members are small, non-enveloped viruses with a single-stranded RNA genome of positive polarity. Members of the Enterovirus genus include human pathogenic poliovirus, coxsackieviruses, enteroviruses, and echoviruses. Other mammalian enteroviruses, including those infecting bovine, simian and porcine species, also have been described [22]. The only picornavirus previously reported to infect alpaca is the foot-and-mouth disease virus (FMDV), which belongs to the Aphthovirus genus. However, FMDV does not usually cause severe disease in alpaca [12], [23].

EV-E and EV-F are globally prevalent infections in cattle, and while virus can be shed in high titers in the feces [24], such infections are usually subclinical and their ability to cause disease in any animal is unclear. Earlier studies described enteroviruses isolated from calves suffering from respiratory disease [25][27]. However, in these studies, respiratory disease could not be reproduced using viral isolates from the infected calves. Subsequent studies in cattle have not been reported.

We hoped to be able to identify sequences that could account for the alpaca infection. While there are insufficient data to determine whether or not the virus adapted to alpaca, the frequent housing of alpaca with cattle without other such reports suggests that these infections are unusual. The alpaca-sourced virus has the interesting characteristic of possessing sequences that are most similar to serotype 1 (including the capsid region that is used to determine picornavirus serotype), but in at least one gene is closest to serotype 3, suggesting that this virus could have arisen by recombination of other EV-F serotypes. It is thus possible that recombination of viruses from two EV-F serotypes led to this unusual infection. The isolate has approximately 80–85% homology in its protein sequence to previously described EV-F strains, which is similar to the degree of homology shared among protein sequences from previously sequenced EV-F strains isolated from cattle, which ranges from 79–99% for EV-F strains, and 50 to 95% when EV-E strains are also considered [15]. The sequence of the alpaca-infecting virus isolate is divergent enough from previously reported strains that it does not provide clear evidence for the basis of its pathogenicity. While it is possible that this virus was transmitted directly from cattle to alpaca, it also is possible that there were one or more intermediate hosts. Besides cattle, EV-F has been reported as an infection of possum and of capped langur [28], [29]. We suspect that the absence of previously reported bovine enterovirus infections in U.S. alpaca is related to the relative isolation of alpaca herds, making it less likely that an alpaca-adapted virus would be further transmitted among alpaca.

Introducing new species (as livestock or as pets) to a habitat potentially increases the risk of an indigenous pathogen causing infections in new species. While most pathogens do not cross the species barrier due to adaptive constraints, those that succeed often cause more severe disease in the new host. Notable human examples of this phenomenon are yellow fever virus, HIV and more recently Nipah virus [30], Hendra virus [31] and SARS virus [32]. Animal examples include the devastating infections of canine distemper virus in raccoons and African lions [33][38]. U.S. alpaca are outside their native South American habitat and are exposed to viruses endemic to the U.S., especially those from U.S. domestic farm animals to which U.S. alpaca herds often have close proximity. Thus, recently described alpaca infections include bovine viral diarrhea virus, equine herpesvirus 1, and bluetongue virus. Newly introduced animals also can potentially carry pathogens that are relatively benign to them, but not to the indigenous fauna. The risk obviously increases if newly introduced and indigenous livestock are kept in close proximity and if their pathogens are able to remain stable in the environment or persist in the host species. Since picornaviruses are non-enveloped viruses, they often are very stable under environmental conditions, increasing the opportunity for infection of different hosts over a prolonged period of time.

While the enterovirus infection described in this report was temporally associated with illness in three other alpaca in the affected herd that may have represented limited spread of the virus, the sparse distribution of alpaca, together with the severe and rapid course of disease likely prevented further dissemination of the virus, as evidenced by the absence of other reports of similar illnesses in the herd or other alpaca in the region. However, even though it appears that this outbreak was controlled, bovine enteroviruses should be added to the list of viruses that can infect alpaca, and that could potentially be associated with severe respiratory and systemic infections in alpaca. Furthermore, considering the relative stability of enteroviruses, the ubiquity of cattle and likely frequent co-location of domestic cattle with alpaca, it is quite plausible that similar outbreaks may occur in the future. Therefore, this alpaca virus infection serves to remind us that viral species are constantly evolving and that the opportunity to infect new hosts may hasten that process.

Materials and Methods

Samples

Samples were obtained from animals within the affected commercial herd that had been submitted after death for a diagnostic necropsy at the University of Illinois Veterinary Diagnostic Laboratory. Grossly affected tissues were harvested at necropsy for routine histopathological examination using 10% neutral buffered formalin-fixed, paraffin-embedded, hematoxylin and eosin-stained sections. Based on the consistent gross necropsy and microscopic findings of acute diffuse interstitial pneumonia in all four alpaca, lung tissue was used for virus isolation.

Ethics statement

The animals used in this study met the definition of “farm animals”, which are not covered by the U.S. Animal Welfare Act (9 CFR 1). Thus, IACUC or ethics committee approval was not required for these studies. The owner of the animals provided permission for these studies.

Cell culture, virus isolation and direct fluorescent antibody detection

MDBK cells were maintained in growth media consisting of Eagle's modified essential medium (MEM; Sigma, St. Louis, MO) supplemented with a 10% horse serum (Sigma), 100 U/ml penicillin, 0.1 mg/ml streptomycin, 0.05 mg/ml gentamicin and 0.0025 mg/ml amphotericin B (Sigma) and kept in a humidified incubator at 37°C with 5% CO2. The following cell lines were used for virus isolation attempts: Madin-Darby bovine kidney (MDBK, ATCC CCL-22), bovine turbinate (gift from Dr. C.L. Kanitz, ADDL, Purdue University, West Lafayette, IN) and rabbit kidney (gift from Dr. C.L. Kanitz). Approximately of 105 cells in one mL was plated for each cell line into separate wells of 24-well tissue culture dishes (Midwest Scientific, St. Louis, MO) containing sterile glass coverslips (12 mm, 31.5 round; Fisher Scientific, Hanover Park, IL). Approximately one gram of lung was homogenized in 10 ml of MEM containing 2× antibiotics using a Stomacher Lab Blender (Fisher Scientific). The homogenate was clarified by centrifugation at 930× g for 20 min at 4°C followed by filtration through a 0.45 µm syringe filter (Minisart NY 25; Fisher Scientific) for inoculation of 200 µl unto each of the individual cell cultures followed by the addition of 1 ml of growth media. For the inoculated cell cultures in which CPE developed, cells adherent to the coverslips were washed once in PBS, then fixed in 100% acetone for 10 min at RT followed by air drying. Appropriate FITC-conjugated anti-viral antibodies were added to the cells and incubated at 37°C for at least 30 min. Unless otherwise indicated, the following anti-viral antibodies were obtained from VMRD (Pullman, WA) and diluted 1∶4 in VMRD 4× rinse buffer: bovine adenovirus types 1 and 5; bluetongue; bovine viral diarrhea virus; bovine coronavirus; bovine herpesvirus types 1, 2, and 5; bovine parainfluenza 3 virus; bovine parvovirus; bovine reovirus; bovine respiratory syncytial virus; and bovine rotavirus (1∶30 in PBS; USDA, NVSL, Ames, IA). After incubation, the cells were washed in PBS, then overlaid with mounting media (VMRD FA mounting fluid or PBS:glycerol [1∶1]) and coverslipped. Fluorescence indicating the presence of a virus was visually determined using a Nikon Optiphot Labophot Episcopic fluorescence microscope.

Negative staining electron microscopy

MDBK cells were grown to 90% confluency in 25 cm2 flasks (Midwest Scientific) containing 10 ml of growth media then inoculated with the 0.5 ml of the filtered virus isolate. After 90% CPE development, cells were harvested by freeze and thawing one time. The cell culture fluid was then clarified by centrifugation at 930× g for 20 min at 4°C. Four ml of the clarified supernatant containing 1.25% (final) neutral buffered formalin was subjected to ultracentrifugation (76 k× g for 1 hour at 4°C) to pellet the virus particles. The pellet subsequently was suspended in 100 µl of supernatant, stirred with a wooden stick and vortexed until thoroughly mixed. The sample was then placed as a drop on parafilm and a formvar plastic and carbon-coated copper grid was placed on top of the specimen droplet for 15 minutes. Excess sample was then removed with filter paper and the grid placed on 2% ammonium molybdate for 2 minutes. The grid was then dried by removing the excess fluid with filter paper, placed into a grid box and covered with drierite crystals for 10 minutes. The grid was then examined in the transmission electron microscope at 20 k and 100 k magnifications.

Degenerate Oligonucleotide Primer- PCR (DOP-PCR)

Viral capsid enrichment.

One mL aliquots of infected and non-infected Madin-Darby bovine kidney (MDBK) cell culture supernatant were digested with DNase I at a final concentration of 80 U/mL DNase I (2000 units lyophilized DNase I; Sigma). The lyophilized DNase I was resuspended in 200 µL DNase I dilution buffer (10 mM HEPES-KOH pH 7.9, 30 mM CaCl2, 30 mM MgCl2, 50% v/v glycerol). The digestion was carried out in 83 mM NaC2H3O2, 4.2 mM Mg2SO4, and 25 mM NaCl for 30 min at 37°C. The DNase digestion was stopped by adding EDTA to a final concentration of 41.5 mM. The digested samples were purified by ultracentrifugation through a 2 ml 1 M NaCl, 10 mM Tris/HCl pH 7.5 solution cushion that was overlaid with 1 ml of 10 mM Tris/HCl pH 7.5 and the digested sample. Centrifugation was carried out for 1.5 hrs, at 4°C and 64,000× g in a SW60 rotor (Beckman Coulter Inc., Fullerton, CA). The supernatant was carefully removed and the air-dried pellet was resuspended in 350 µL RLT+-buffer (Qiagen, Valencia, CA) containing 143 mM β-mercaptoethanol (MP Biomedicals Inc., Solon, OH). DNA and RNA were isolated in parallel from the same sample using the AllPrep DNA/RNA kit (Qiagen). DNA was eluted in 100 µL elution buffer; RNA was eluted in 50 µL H2O. RNA was transcribed in to cDNA with random hexamer primers according to the manufacturer's instructions (First Strand Kit, Invitrogen, Carlsbad, CA). Following the cDNA synthesis the RNA was digested away with 2 µL RNase H for 20 min at 37°C (First Strand Kit, Invitrogen).

DOP-PCR.

reactions contained 1.5 mM MgCl2, 10 mM KCl, 10 mM Tris pH 8.4, 200 µM dNTP, 2.4 µM DOP primer 5′-CCGACTCGAGINNNNNNTGTGG-3′ (Oligos Etc., Wilsonville, OR) and 2.5 U Low DNA Taq polymerase (Applied Biosystems, Foster City, CA). The reactions were carried out with 10 µL of template. Cycling conditions consisted of initial denaturation for 5 min at 95°C; followed by 5 cycles of 1 min at 94°C, 5 min at 25°C, slow ramping at 0.1°C/sec to 30°C, 4 min at 30°C, slow ramping at 0.1°C/sec to 37°C, 3 min at 37°C, slow ramping at 0.1°C/sec to 42°C, 2 min at 42°C, slow ramping at 0.1°C/sec to 55°C, 55°C for 1 min, 72°C for 2 min; 35 cycles as follows: 94°C for 20 sec, 55°C for 1 min, 72°C for 1 min with the addition of 1 second per cycle to the extension step; final extension at 72°C for 10 min.

The DOP-PCR products were analyzed and purified by agarose gel electrophoresis. Distinct bands were excised, purified (GenElute, Sigma, St. Louis, MO), and ligated into the pCR4-TOPO® vector (Invitrogen). The ligation products were used to transform competent One Shot® TOP10 bacteria (Invitrogen) according to the manufacturer's instructions. Colony PCR was performed; these PCR products were purified (QIAquick PCR Purification, Qiagen) and sequenced using M13 primers. Sequencing was performed using a 3130×l Genetic Analyzer (Applied Biosystems, Foster City, CA). Sequences from the obtained clones were compared with the non-redundant (nr) database in GenBank using TBLASTX (NCBI, Bethesda, MD).

Virus genome sequencing.

The sequence data obtained from DOP-PCR products as well as data for EV-E and EV-F genomes available in GenBank (NC_001859.1, AY508697.1, AY508696.1, D00214.1, AF123433.1, AF123432.1) were used to develop multiple primer pairs (Table 1) to amplify and sequence the full genome of the novel virus. PCR conditions were: initial denaturation 2 min 95°C, 40 cycles of 94°C for 30 sec, 54°C for 1 min, 72°C for 1 min, concluding with a final extension step of 72°C for 1 min. PCR products were cloned and sequenced as described above.

The 3′ end of the viral genome was sequenced using rapid amplification of cDNA ends (RACE) with a commercially available kit (Smart RACE cDNA amplification kit, Clontech, Mountain View, CA) according to the manufacturer's instructions. Viral RNA was extracted from infected cell cultures and cDNA was synthesized using MMLV reverse transcriptase and an oligo dT primer provided with the commercial kit. The cDNA was then PCR-amplified with a gene specific primer (GSP) and the universal primer provided in the kit. The GSP (SM-11) was designed based on the cloned alpaca virus cDNA sequence and the EV-E and EV-F sequences in the GenBank database and was approximately 650-bp upstream of the 3′ poly A tail. The resulting PCR products were gel purified cloned into a TA cloning vector. The 3′ viral ends were sequenced directly from the purified PCR products and from the cloned cDNA to verify the correct sequence.

Formalin-fixed tissues

Lung tissues from alpaca and horses involved in the initial outbreak were formalin-fixed and paraffin embedded (FFPE) for future analysis and histopathology. Following the identification of enteroviruses sequences from MDBK cell cultures by DOP-PCR, we attempted to test these FFPE tissues by specific PCR assays for enterovirus sequences. Five sections of 10 µm thickness of each tissue were removed from each block with a microtome and RNA was extracted with the RNeasy FFPE Kit according to the manufacturer's instructions (Qiagen). RNA was reverse transcribed into cDNA according to the manufacturer's instructions (First strand synthesis kit, Invitrogen).

Several specific PCRs for enterovirus sequences were performed using primers designed for this alpaca sourced enterovirus (SM 9–32, Table 1) ranging from 100–500-bp and published bovine enterovirus primers [24], [39] (Beld, BEV, N-BEV; Table 1). PCRs for two bovine housekeeping genes, β–actin and glyceraldehyde 3-phosphate dehydrogenase (GAPDH) [40], [41], were also performed to verify the quality of the RNA (Table 1). Additionally DOP-PCR was performed on the cDNA as described above. PCR products were cloned and sequenced to verify the identity of the nucleic acid amplified.

Phylogenetic analysis

The deduced amino acid sequences from the alpaca virus polyprotein containing the capsid, polymerase, and protease genes were aligned with other homologous enteroviruses available in the GenBank database. The sequences were aligned with Clustal W and neighbor-joining phylogenetic trees were constructed and viewed using Treeview [42] and Phylip software [43] with bootstrap confidence values determined by 1000 replications. GenBank accession numbers used in the analyses were: Enterovirus F, strain IL/Alpaca KC748420, BEV-A, now EV-E (serotypes 1–4) and BEV-B, now EV-F (serotypes 1–4) AF123433, D00214, AF123432, DQ092770, DQ092795, AY508696, AY508697, DQ092794, NC001859, AY462106, DQ092786, DQ092787, JQ690748, EU886967, HQ917060, NC07767, JQ690741, AY724745, and JX538037; HEV A FM955278, HEV B AF029859, HEV C U05876, and HEV D D00820; polio virus Mahoney strain N002058; and porcine enterovirus B NC_004441. Not all of the viral sequences used for analyses contained the complete polyprotein sequences, and therefore some do not appear in all of the panels of Figures 4 to 8. The capsid sequences for the capped langur (JX538037) and unclassified AY24745 sequences are partial. EV-F serotypes 5 and 6 [16] do not have nucleic acid sequences available and do not appear in these analyses.

Acknowledgments

We acknowledge the contributions of Debbie Cassout (University of Illinois), who performed the virus isolation and direct fluorescent antibody work, and Lou Ann Miller (University of Illinois), who performed the negative staining EM work.

Author Contributions

Conceived and designed the experiments: CU GS SM LB RF PK. Performed the experiments: CU GS SM LB RF. Analyzed the data: CU GS SM LB RF PK. Contributed reagents/materials/analysis tools: CU GS SM LB RF PK. Wrote the paper: CU GS SM LB RF PK.

References

  1. 1. Kapil S, Yeary T, Evermann JF (2009) Viral diseases of new world camelids. Vet Clin North Am Food Anim Pract 25(2): 323–37.
  2. 2. United States Department of Agriculture (2009) 2007 Census of Agriculture. Available: http://www.agcensus.usda.gov/Publications/2007/Full_Report/usv1.pdf. Accessed 2013 Jun 14.
  3. 3. Barrington GM, Parish SM, Tyler JW (1997) Chronic weight loss in an immunodeficient adult llama. J Am Vet Med Assoc 211(3): 294–5; discussion 6–8.
  4. 4. Wernery U (2001) Camelid immunoglobulins and their importance for the new-born–a review. J Vet Med B Infect Dis Vet Public Health 48(8): 561–8.
  5. 5. Wentz PA, Belknap EB, Brock KV, Collins JK, Pugh DG (2003) Evaluation of bovine viral diarrhea virus in New World camelids. J Am Vet Med Assoc 223(2): 223–8.
  6. 6. Belknap EB, Collins JK, Larsen RS, Conrad KP (2000) Bovine viral diarrhea virus in New World camelids. J Vet Diagn Invest 12(6): 568–70.
  7. 7. Goyal SM, Bouljihad M, Haugerud S, Ridpath JF (2002) Isolation of bovine viral diarrhea virus from an alpaca. J Vet Diagn Invest 14(6): 523–5.
  8. 8. Barnett J, Twomey DF, Millar MF, Bell S, Bradshaw J, et al. (2008) BVDV in British alpacas. Vet Rec 162(24): 795.
  9. 9. Kim SG, Anderson RR, Yu JZ, Zylich NC, Kinde H, et al. (2009) Genotyping and phylogenetic analysis of bovine viral diarrhea virus isolates from BVDV infected alpacas in North America. Vet Microbiol 136(3–4): 209–16.
  10. 10. Topliff CL, Smith DR, Clowser SL, Steffen DJ, Henningson JN, et al. (2009) Prevalence of bovine viral diarrhea virus infections in alpacas in the United States. J Am Vet Med Assoc 234(4): 519–29.
  11. 11. Mattson DE, Baker RJ, Catania JE, Imbur SR, Wellejus KM, et al. (2006) Persistent infection with bovine viral diarrhea virus in an alpaca. J Am Vet Med Assoc 228(11): 1762–5.
  12. 12. Mattson DE (1994) Update on llama medicine. Viral diseases. Vet Clin North Am Food Anim Pract 10(2): 345–51.
  13. 13. Byers SR, Snekvik KR, Righter DJ, Evermann JF, Bradway DS, et al. (2009) Disseminated Bovine viral diarrhea virus in a persistently infected alpaca (Vicugna pacos) cria. J Vet Diagn Invest 21(1): 145–8.
  14. 14. Knowles NJ, Hovi T, Hyypia T, King AMQ, Lindberg AM, et al.. (2012) Picornaviridae. In: Virus taxonomy: classification and nomenclature of viruses: Ninth Report of the International Committee on Taxonomy of Viruses. San Diego: Elsevier.
  15. 15. Zell R, Krumbholz A, Dauber M, Hoey E, Wutzler P (2006) Molecular-based reclassification of the bovine enteroviruses. J Gen Virol 87(Pt 2): 375–85.
  16. 16. Knowles NJ (2013) The Picornavirus Pages Available: www.picornaviridae.com. Accessed 2013 Jun 14.
  17. 17. Uhlenhaut C, Cohen JI, Pavletic S, Illei G, Gea-Banacloche JC, et al. (2012) Use of a novel virus detection assay to identify coronavirus HKU1 in the lungs of a hematopoietic stem cell transplant recipient with fatal pneumonia. Transpl Infect Dis 14(1): 79–85.
  18. 18. Nanda S, Jayan G, Voulgaropoulou F, Sierra-Honigmann AM, Uhlenhaut C, et al. (2008) Universal virus detection by degenerate-oligonucleotide primed polymerase chain reaction of purified viral nucleic acids. J Virol Methods 152(1–2): 18–24.
  19. 19. Uhlenhaut C, Cohen JI, Fedorko D, Nanda S, Krause PR (2009) Use of a universal virus detection assay to identify human metapneumovirus in a hematopoietic stem cell transplant recipient with pneumonia of unknown origin. J Clin Virol 44(4): 337–9.
  20. 20. Goens SD, Botero S, Zemla A, Zhou CE, Perdue ML (2004) Bovine enterovirus 2: complete genomic sequence and molecular modelling of a reference strain and a wild-type isolate from endemically infected US cattle. J Gen Virol 85(Pt 11): 3195–203.
  21. 21. Oberste MS, Feeroz MM, Maher K, Nix WA, Engel GA, et al. (2013) Characterizing the picornavirus landscape among synanthropic nonhuman primates in Bangladesh, 2007 to 2008. J Virol 87(1): 558–71.
  22. 22. Pallansch MA, Roos RP (2007) Enteroviruses: Polioviruses, Coxsackieviruses, Echoviruses, and Newer Enteroviruses. In: Knipe DM, Howley PM, editors. Fields Virology: Lippincott Williams & Wilkins, 723–75.
  23. 23. Wernery U, Kaaden O (2004) Foot-and-mouth disease in camelids: a review. Vet J 168(2): 134–42.
  24. 24. Jimenez-Clavero MA, Escribano-Romero E, Mansilla C, Gomez N, Cordoba L, et al. (2005) Survey of bovine enterovirus in biological and environmental samples by a highly sensitive real-time reverse transcription-PCR. Appl Environ Microbiol 71(7): 3536–43.
  25. 25. Phillip JI, Darbyshire JH (1971) Respiratory viruses of cattle. Adv Vet Sci Comp Med 15: 159–99.
  26. 26. Moll T, Ulrich MI (1963) Biologic characteristics of certain bovine enteric viruses. Am J Vet Res 24: 545–50.
  27. 27. Huck RA, Cartwright SF, Yamanouchi K, Bankowski RA, Howarth JA (1963) Polioencephalomyelitis of pigs- the identification of viruses related to the Teschen and T80 groups in the United States. Am J Vet Res 24: 1207–11.
  28. 28. Oberste MS, Feeroz MM, Maher K, Nix WA, Engel GA, et al. (2013) Naturally acquired picornavirus infections in primates at the Dhaka zoo. J Virol 87(1): 572–80.
  29. 29. Zheng T (2007) Characterisation of two enteroviruses isolated from Australian brushtail possums (Trichosurus vulpecula) in New Zealand. Arch Virol 152(1): 191–8.
  30. 30. Middleton DJ, Weingartl HM (2012) Henipaviruses in Their Natural Animal Hosts. Curr Top Microbiol Immunol 359: 105–121.
  31. 31. Halpin K, Young P, Field H, Mackenzie J (2000) Isolation of Hendra virus from pteropid bats: a natural reservoir of Hendra virus. J Gen Virol 81: 1927–32.
  32. 32. Shi Z, Hu Z (2008) A review of studies on animal reservoirs of the SARS coronavirus. Virus Res 133: 74–87.
  33. 33. Lednicky JA, Dubach J, Kinsel MJ, Meehan TP, Bocchetta M, et al. (2004) Genetically distant American Canine distemper virus lineages have recently caused epizootics with somewhat different characteristics in raccoons living around a large suburban zoo in the USA. Virol J 1: 2.
  34. 34. Kameo Y, Nagao Y, Nishio Y, Shimoda H, Nakano H, et al. (2012) Epizootic canine distemper virus infection among wild mammals. Vet Microbiol 154(3–4): 222–9.
  35. 35. Wood SL, Thomson GW, Haines DM (1995) Canine distemper virus-like infection in a captive African lioness. Can Vet J 36(1): 34–5.
  36. 36. Roelke-Parker ME, Munson L, Packer C, Kock R, Cleaveland S, et al. (1996) A canine distemper virus epidemic in Serengeti lions (Panthera leo). Nature 379(6564): 441–5.
  37. 37. Guiserix M, Bahi-Jaber N, Fouchet D, Sauvage F, Pontier D (2007) The canine distemper epidemic in Serengeti: are lions victims of a new highly virulent canine distemper virus strain, or is pathogen circulation stochasticity to blame? J R Soc Interface 4(17): 1127–34.
  38. 38. Munson L, Terio KA, Kock R, Mlengeya T, Roelke ME, et al. (2008) Climate extremes promote fatal co-infections during canine distemper epidemics in African lions. PLoS One 3(6): e2545.
  39. 39. Beld M, Minnaar R, Weel J, Sol C, Damen M, et al. (2004) Highly sensitive assay for detection of enterovirus in clinical specimens by reverse transcription-PCR with an armored RNA internal control. J Clin Microbiol 42(7): 3059–64.
  40. 40. Moniwa M, Clavijo A, Li M, Collignon B, Kitching R (2007) Performance of a foot-and-mouth disease virus reverse transcription-polymerase chain reaction with amplification controls between three real-time instruments. J Vet Diagn Invest 19(1): 9–20.
  41. 41. Walker C, Meier S, Mitchell M, Roche J, Littlejohn M (2009) Evaluation of real-time PCR endogenous control genes for analysis of gene expression in bovine endometrium. BMC Mol Biol 10: 100
  42. 42. Page RD (1996) TreeView: an application to display phylogenetic trees on personal computers. Comput Appl Biosci 12(4): 357–8.
  43. 43. Felsenstein J (1997) An alternating least squares approach to inferring phylogenies from pairwise distances. Syst Biol 46(1): 101–11.