Malpais spring virus is a new species in the genus vesiculovirus
© Vasilakis et al; licensee BioMed Central Ltd. 2013
Received: 9 January 2013
Accepted: 21 February 2013
Published: 4 March 2013
Malpais Spring virus (MSPV) is a mosquito-borne rhabdovirus that infects a variety of wild and feral ungulates in New Mexico, including horses and deer. Although, initial serologic tests and electron microscopy at the time of isolation nearly 25 years ago provided evidence that MSPV is a novel virus, possibly related to vesiculoviruses, the virus still has not been approved as a new species.
Use of the illumina platform allowed us to obtain the complete genome of MSPV. Analysis of the complete 11019 nt genome sequence of the prototype 85-488NM strain of MSPV indicates that it encodes the five common rhabdovirus structural proteins (N, P, M, G and L) with alternative ORFs (> 180 nt) in the N, M and G genes, including a 249 nt ORF in the G gene predicted to encode a 9.26 kDa highly basic transmembrane protein. Although antigenically very distant, phylogenetic analysis of the L gene indicates that MSPV is most closely related to Jurona virus, also isolated from mosquitoes in Brazil, as well as a number of other vesiculoviruses.
In sum, our analysis indicates MSPV should be classified as a member of the genus Vesiculovirus, family Rhabdoviridae. The complete genome sequence of MSPV will be helpful in the development of a reverse genetics system to study the unique aspects of this vesiculovirus in vivo and in vitro, and will assist development of specific diagnostic tests to study the epidemiology of MSPV infection.
KeywordsMalpais spring virus (MSPV) Family rhabdoviridae Genus vesiculovirus Genome sequence Phylogeny
Malpais Spring virus
Nearest neighbor interchange
RNA dependent RNA polymerase
Vesicular stomatitis Indiana virus
Malpais Spring virus (MSPV) is a negative-sense RNA virus currently assigned to the family Rhabdoviridae in the order Mononegavirales. MSPV was first isolated from pools of Aedes campestris and Psorophora signipennis mosquitoes collected during an investigation of unexplained deaths among feral horses living near a spring in the Malpais Lava Flow on the White Sands Missile Range in New Mexico, USA . Initial serologic tests and electron microscopy provided evidence that MSPV is a novel rhabdovirus, possibly related to vesiculoviruses . MSPV is lethal to newborn mice following intracerebral (i.c.) inoculation and causes a cytopathic effect (CPE) in Vero cells, forming sizeable (3–4 mm in diameter) plaques six days post-inoculation . Plaque reduction neutralization tests conducted using sera from wild mammals living in the area where the virus was recovered, indicated a relatively high prevalence (30-100%) of specific MSPV neutralizing antibodies among feral horses, coyotes, mule deer, gemsbok (Oryx gazella) and pronghorn, although the virus could not be clearly associated with the horse deaths .
The prototype strain of MSPV (85-488NM) was obtained from the World Reference Center for Emerging Viruses and Arboviruses at the University of Texas Medical Branch. The virus stock used had been passaged four times by i.c. inoculation of newborn mice and three times in Vero cells.
Next generation sequencing
The complete genome of MSPV (85-488NM) was determined by de novo sequencing using the Illumina platform.
Viral RNA (0.05-1.7 μg) was fragmented by incubation at 94o C for eight (8) minutes in 19.5 μl of fragmentation buffer (Illumina 15016648). First and second strand synthesis, adapter ligation and amplification of the library were performed using the Illumina TruSeq RNA Samplec Preparation kit under conditions prescribed by the manufacturer (Illumina). Samples were tracked using the “index tags” incorporated into the adapters as defined by the manufacturer.
Cluster formation of the library DNA templates was performed using the TruSeq PE Cluster Kit v3 (Illumina) and the Illumina cBot workstation using conditions recommended by the manufacturer. Paired end 50 base sequencing by synthesis was performed using TruSeq SBS kit v3 (Illumina) on an Illumina HiSeq 1000 using protocols defined by the manufacturer. Cluster density per lane was 645–980 k/mm2 and post filter reads ranged from 148–178 million per lane. Base call conversion to sequence reads was performed using CASAVA-1.8.2. Virus assembly was performed using SeqMan Lasergene software (DNASTAR). In certain cases, pre-filtering of reads to remove host sequence enhanced the assembly process. Assembly was carried out using a fasta file of hamster ribosomal RNA sequences to remove host DNA from the assembly, thus reducing the number of contigs present.
Nucleotide sequence accession number
The GenBank accession number of MSPV prototype 85-488NM strain is KC412247.
The L protein sequence of MSPV was compared with those of 44 other rhabdoviruses downloaded from GenBank (members of the genera Cytorhabdovirus, Novirhabdovirus and Nucleorhabdovirus were excluded because their excessive divergence reduced phylogenetic resolution). All protein sequences were aligned using MUSCLE  under default settings, and ambiguously aligned regions were then removed using the Gblocks program . This resulted in a final sequence alignment of 1136 amino acid residues. The phylogenetic relationships among these sequences were determined using the maximum likelihood (ML) method available in PhyML 3.0 , employing the WAG+Γ model of amino acid substitution and subtree pruning and regrafting (SPR) branch-swapping. The robustness of each node was determined using 1,000 bootstrap replicates utilizing nearest neighbor interchange (NNI) branch-swapping.
Antigens and immune reagents
Antigenic relationships of MSPV with confirmed and tentative members of the genus Vesiculovirus
VSV Indiana (VSIV)
VSV New Jersey (VSNJV)
Gray Lodge (GLOV)
La Joya (LJV)
Jug Bogdanovac (YBV)
Malpais Spring (MSPV)
Complement fixation tests were conducted according to a microtechnique described previously , using 2 full units of guinea-pig complement. Titers were recorded as the highest dilutions giving 3+ or 4+ fixation of complement on a scale of 0–4+.
Amino acid sequence alignments (MUSCLE) using available rhabdovirus N, G and L proteins indicated that MSPV shares highest sequence identity to vesiculoviruses. Phylogenetic analysis based on an alignment of 1136 amino acid residues of the L protein sequence of 45 rhabdoviruses indicated that MSPV is most closely related (although with weak bootstrap support) to JURV, which has been isolated from Haemagogus spp. mosquitoes in the state of Para, Brazil in 1962 [10, 11], as well as to Isfahan, CHPV and Perinet viruses (Figure 2). These five viruses form a well-supported (100% bootstrap support) phylogenetic cluster within the genus Vesiculovirus.
The genome organization and phylogenetic relationships indicate that Malpais Spring virus should be classified as a new species in the genus Vesiculovirus. This is also consistent with the previous report of a distant antigenic relationship with JURV, which also falls phylogenetically within the vesiculovirus cluster. Evidence of alternative ORFs in the MSPV N, M and G genes is unique amongst vesiculovirus genomes sequenced to date but overlapping ORFs appear to occur commonly in rhabdoviruses , including VSIV in which an alternative ORF in the P gene encodes the 55-aa and 65-aa C and C’ proteins which are known to be expressed in infected cells from alternative initiation codons [12, 13]. The complete genome sequence of MSPV will be helpful in the development of a reverse genetics system to study the unique aspects of this vesiculovirus in vivo and in vitro, and will assist development of specific diagnostic tests to study the epidemiology of MSPV infection.
This work is supported in part by the Department of Pathology start up funds and a grant from the Institute for Human Infections and Immunity, University of Texas Medical Branch (NV), NIH contract HHSN272201000040I/HHSN27200004/D04 (RBT).
- Dietzgen RG, Calisher CH, Kurath G, Kuzmin IV, Rodriguez LL, Stone DM, Tesh RB, Tordo N, Walker PJ, Wetzel T, Whitfield AE: Rhabdoviridae. In Virus Taxonomy, Ninth Report of the International Committe on Taxonomy of Viruses. Edited by: King AMQ, Adams MJ, Carstens EB, Lefkowitz EJ. London: Elsevier/Academic Press; 2012:686–713.
- Clark GG, Calisher CH, Crabbs CL, Canestorp KM, Tesh RB, Bowen RA, Taylor DE: Malpais spring virus: a new vesiculovirus from mosquitoes collected in New Mexico and evidence of infected indigenous and exotic ungulates. Am J Trop Med Hyg 1988, 39:586–592.PubMed
- Edgar RC: MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinforma 2004, 5:113.View Article
- Talavera G, Castresana J: Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. Syst Biol 2007, 56:564–577.PubMedView Article
- Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O: New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol 2010, 59:307–321.PubMedView Article
- Clarke DH, Casals J: Techniques for hemagglutination and hemagglutination-inhibition with arthropod-borne viruses. Am J Trop Med Hyg 1958, 7:561–573.PubMed
- Tesh RB, Travassos Da Rosa AP, Travassos Da Rosa JS: Antigenic relationship among rhabdoviruses infecting terrestrial vertebrates. J Gen Virol 1983, 64:169–176.PubMedView Article
- Walker PJ, Dietzgen RG, Joubert DA, Blasdell KR: Rhabdovirus accessory genes. Virus Res 2011, 162:110–125.PubMedView Article
- Walker PJ, Kongsuwan K: Deduced structural model for animal rhabdovirus glycoproteins. J Gen Virol 1999, 80:1211–1220.PubMed
- Karabatsos N: Jurona virus. 3rd edition. San Antonio: American Society of Tropical Medicine and Hygiene; 1985:527–528. [International catalogue of arboviruses including certain other viruses of vertebrates]
- Shope RE, de Andrade AH, Bensabath G, Causey OR, Humphrey PS: The epidemiology of EEE WEE, SLE and Turlock viruses, with special reference to birds, in a tropical rain forest near Belem, Brazil. Am J Epidemiol 1966, 84:467–477.PubMed
- Peluso RW, Richardson JC, Talon J, Lock M: Identification of a set of proteins (C' and C) encoded by the bicistronic P gene of the Indiana serotype of vesicular stomatitis virus and analysis of their effect on transcription by the viral RNA polymerase. Virology 1996, 218:335–342.PubMedView Article
- Spiropoulou CF, Nichol ST: A small highly basic protein is encoded in overlapping frame within the P gene of vesicular stomatitis virus. J Virol 1993, 67:3103–3110.PubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.