A small focus of hemorrhagic fever (HF) cases occurred near Cochabamba, Bolivia, in December 2003 and January 2004. Specimens were available from only one fatal case, which had a clinical course that included fever, headache, arthralgia, myalgia, and vomiting with subsequent deterioration and multiple hemorrhagic signs. A non-cytopathic virus was isolated from two of the patient serum samples, and identified as an arenavirus by IFA staining with a rabbit polyvalent antiserum raised against South American arenaviruses known to be associated with HF (Guanarito, Machupo, and Sabiá). RT-PCR analysis and subsequent analysis of the complete virus S and L RNA segment sequences identified the virus as a member of the New World Clade B arenaviruses, which includes all the pathogenic South American arenaviruses. The virus was shown to be most closely related to Sabiá virus, but with 26% and 30% nucleotide difference in the S and L segments, and 26%, 28%, 15% and 22% amino acid differences for the L, Z, N, and GP proteins, respectively, indicating the virus represents a newly discovered arenavirus, for which we propose the name Chapare virus. In conclusion, two different arenaviruses, Machupo and Chapare, can be associated with severe HF cases in Bolivia.
Four rodent-borne arenaviruses are known to cause hemorrhagic fever (HF) in the New World. These include Junín, Machupo, Guanarito, and Sabiá viruses, which are found in rural areas of Argentina, Bolivia, Venezuela, and Brazil, respectively. In December 2003 and January 2004, a small number of HF cases were reported in rural Bolivia in an area outside the known Machupo HF endemic zone, and sera from one fatal case was available for laboratory testing. The man had symptoms similar to those seen with other arenaviral HF cases—acute febrile illness beginning with headache, joint and muscle pain, and vomiting—and rapidly progressed to shock, bleeding, and death at 14 days post onset of illness. Virus was isolated from two of the patient's serum samples and identified as an arenavirus by reaction of virus infected cells with arenavirus-specific antibodies and by genetic detection techniques (PCR). Subsequent complete genome analysis of the virus showed the virus to be a distinct newly discovered member of the arenavirus family, and the name Chapare virus was proposed. The virus is phylogenetically related to other arenaviruses that naturally cause hemorrhagic fever in South America, particularly Sabiá virus. Physicians should consider Chapare virus as a potential etiologic agent when encountering HF cases in the region.
Citation: Delgado S, Erickson BR, Agudo R, Blair PJ, Vallejo E, Albariño CG, et al. (2008) Chapare Virus, a Newly Discovered Arenavirus Isolated from a Fatal Hemorrhagic Fever Case in Bolivia. PLoS Pathog 4(4): e1000047. https://doi.org/10.1371/journal.ppat.1000047
Editor: Michael J. Buchmeier, University of California Irvine, United States of America
Received: January 18, 2008; Accepted: March 20, 2008; Published: April 18, 2008
This is an open-access article distributed under the terms of the Creative Commons Public Domain declaration which stipulates that, once placed in the public domain, this work may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose.
Funding: These investigations were funded by the Centers for Disease Control and Prevention and the United States Navy. The findings and conclusions in this report are those of the authors and do not necessarily represent the views of the funding agencies.
Competing interests: The authors have declared that no competing interests exist.
The family Arenaviridae is composed of largely rodent-borne viruses which are divided into Old World and New World complexes ,. Lassa and lymphocytic choriomeningitis (LCM) viruses are considered the most important of Old World arenaviruses due to their association with severe disease. The New World complex is divided into 3 major Clades (A, B and C), with Clade B containing all the hemorrhagic fever (HF) associated viruses ,,,. These are Junín, Machupo, Guanarito and Sabiá viruses, the cause of Argentine, Bolivian, Venezuelan, and Brazilian HF, respectively . Three of these viruses, Junín, Machupo, and Guanarito, can be associated with large HF outbreaks and untreated case fatalities can be in excess of 30%. The clinical picture is similar for each of these diseases. Onset of symptoms follows an incubation period of 1–2 weeks. Initial symptoms often include fever, malaise, myalgia and anorexia, followed approx. 3–4 days later by headache, back pain, dizziness, nausea, vomiting, and severe prostration. Hemorrhagic and neurologic symptoms, including petechiae and bleeding gums, tremors, and lethargy are common. About a third of untreated cases go on to develop more severe neurologic and/or hemorrhagic symptoms, with diffuse echymoses, and bleeding from mucous membranes or puncture sites, and/or delirium, coma and convulsions. Machupo virus, vectored by Calomys callosus rodents , is the only known pathogenic arenavirus found in Bolivia, although another arenavirus, Latino virus, has also been isolated from Calomys callosus in Bolivia . Despite broad distribution of this rodent host, which is thought to include the lowlands of Bolivia, east-central Brazil, Paraguay and northern Argentina , Machupo virus-associated HF cases have originated only in the Beni department in northeastern Bolivia (Figure 1). We report here the investigation of a fatal HF case which occurred near Cochabamba, Cochabamba Department, Bolivia in December, 2003, and identify the associated arenavirus as a unique newly discovered virus, Chapare virus.
In late 2003 reports were received of a small cluster of HF cases in a rural area near the Chapare River, close to Cochabamba, Bolivia in the eastern foothills of the Andes (Figure 1). Exact details of the number of cases and verification of symptoms were difficult to obtain. However, a clinical description and blood specimens were available for one fatal case. This patient, a 22 year old male, had lived for the last 4 years in Samuzabeti, a small town located 35 km northeast of Villa Tunari. He was a tailor and also a farmer. Coca is the main crop in this area. He had no history of travel and no contact with any case with compatible illness for at least 4 weeks prior to his disease onset on January 3rd, 2004. In addition, no members of case household or other close contacts were affected. His clinical course included fever, headache, arthralgia, myalgia and vomiting with subsequent deterioration and multiple hemorrhagic signs and death on January 17th, 2004 (14 days post onset). Based on these symptoms, the patient was initially suspected of having yellow fever or dengue HF. However, initial tests for these agents were negative. Initial IgM, IgG, antigen capture and RT-PCR testing for Machupo virus or related arenaviruses were also negative.
Patient specimens were sent to the biosafety level 4 (BSL4) containment laboratory at the Special Pathogens Branch in Atlanta where virus isolation attempts could be performed. These specimens consisted of 4 acute phase sera, collected on days 4, 7, 9 and 14 post onset of disease. Both the day 7 and day 9 sera yielded a non-cytopathic virus by growth in Vero E6 cells. These were identified by immunofluorescent antibody (IFA) staining with rabbit polyvalent hyperimmune serum raised against South American arenaviruses previously known to be associated with HF (Guanarito, Machupo, and Sabiá). RT-PCR analysis of the virus isolate RNAs amplified a 481 bp fragment which yielded nucleotide sequence related to known New World Clade B arenaviruses (which includes all the South American HF associated arenaviruses). Full length virus genome sequences were successfully determined for the virus isolated from the day 9 post onset bleed (designated strain 810419) by RT-PCR and sequence analysis followed by primer walking utilizing newly derived sequence information. The full length S segment was amplified by using the 19C primer designed based on the conserved RNA termini of New World Arenaviruses , whereas the L segment was amplified in multiple sections using a variety of primers (sequences available on request). Sequence analysis of the complete S and L segments confirmed that this virus, proposed name Chapare, was a unique member of the Clade B New World arenaviruses ,,,,. The virus was found to be most closely related to Sabiá virus, but with 26 and 30% nucleotide difference in the complete S and L segments, and 26, 28, 15 and 22% amino acid differences for the L, Z, N and GP proteins, respectively (Tables 1 and S1). The genetic differences between Chapare virus and other Clade B viruses range from 36–40% for the complete S segment and 39–40% for the complete L segment (data not shown). These nucleotide and amino acid sequence divergence levels are in excess of those seen among strains of the same species of New World arenavirus (Tables 2 and S1) ,,. For instance, the greatest difference seen between complete S segments of virus strains is 14% (within Allpahuayo virus strains) and 10% for the complete L segment (among Machupo virus strains) .
Chapare virus was found to be monophyletic with Sabiá virus on phylogenetic analysis of the nucleotide or encoded amino acid sequences of the complete S or L segment (Figure 2), or NP, GP, L or Z ORFs (data not shown). No evidence of reassortment or recombination between Chapare virus and other arenaviruses was found. There is no overall change in the structure of the trees except for the previously described , switch of the Clade A/Rec viruses from Clade A for the NP gene to Clade B for the GP gene (data not shown).
Complete S and L segments for New World arenaviruses were analyzed by Bayesian inference of phylogeny (MrBayes3.1.2) using the sequence of Pichindé virus as the outgroup. Multiple strains are grouped with small brackets and large brackets group the arenavirus Clades: A, A/Rec, B, and C. The Genbank accession numbers for the S segment analysis include: Allpahuayo (AY012686), Amapari (AF485256), Bear Canyon (AY924392), Catarina (DQ865245), Chapare (EU260463), Cupixi (AF512832), Flexal (AF485257), Guanarito (NC_005077), Junín (AY619641, NC_005081, AY746353), Latino (AF512830), Machupo (AY924208, AY619645, AY924202, NC_005078), Oliveros (U34248), Paraná (AF485261), Pichindé (NC_006447), Pirital (NC_005894), Sabiá (NC_006317), Tacaribe (NC_004293), Tamiami (AF485263), and Whitewater Arroyo (AF485264). The Genbank accession numbers used for the L segment analysis include: Allpahuayo (NC_010249), Amapari (AY924389), Bear Canyon (AY924390), Chapare (EU260464), Cupixi (NC_010252), Guanarito (NC_005082), Junín (NC_005080, AY819707, AY619640), Machupo (AY624354, NC_005079, AY619644), Oliveros (NC_010250), Pichindé (NC_006439), Pirital (NC_005897), Sabiá (NC_006313), Tacaribe (NC_004292), Tamiami (AY924393), and Whitewater Arroyo (AY924395).
The pathogenicity of the New World Clade B viruses correlates with the efficient interaction of their GP1 surface glycoproteins with the human cellular receptor, transferring receptor 1 (TfR1) ,. Presumably, Chapare virus will be found to have similar TfR1 binding properties, but this remains to be confirmed. Even assuming this to be true, the diversity of the GP1 amino acid sequences of Junín, Machupo, Guanarito, Sabiá and Chapare viruses is such that one cannot easily discern the GP1 domain involved in high efficiency binding to TfR1 solely on the basis of amino acid sequence alignments.
The relationship of Chapare virus from Bolivia to Sabiá virus from Brazil is intriguing. Both these virus clearly cause HF similar to that seen with Junín, Machupo and Guanarito viruses. The single HF case associated with a naturally acquired Sabiá virus infection was reported in the community of Sabiá, in Sao Paulo, Brazil in 1990 . The exact site of exposure was unclear, as was the rodent reservoir. Yellow fever was the initial suspicion in the Sabiá case and that associated with the Chapare virus infection as both had associated extensive liver necrosis. More extensive liver involvement may be a feature shared between these viruses, as it is not commonly observed with HF associated with the other New World arenaviruses (although it is occasionally seen).
Due to the difficulties of working in this resource poor rural region, initial follow up efforts in the Chapare area, did not yield a more precise description of the reported cluster of cases with similar illness, and a limited ecological study did not reveal the rodent reservoir of this virus. It is hoped that more extensive studies in the area will reveal the extent to which Chapare virus poses a public health problem in this area, and shed light on the source of human infection. In summary, three arenaviruses are now known to be present in Bolivia, namely Machupo and Latino viruses (both hosted by Calomys callosus) and Chapare viruses (reservoir unknown). Furthermore, both Machupo and Chapare viruses are agents of fatal hemorrhagic fever in Bolivia.
Materials and Methods
Diagnostic Amplification & Identification
Initial virus genetic detection and analysis was conducted on total RNA extracted from infected Vero E6 cells, using Tripure Isolation Reagent (Roche Applied Science, Indianapolis, IN) in a ratio of 1∶5 and incubated at room temperature for a minimum of 10 min. Total RNA was isolated by using the RNaid Kit following the manufacturer's recommendations (Qbiogene Inc., Carlsbad, CA), and the extracted RNA was reconstituted in 50 µL H2O. Broadly reactive Arenavirus primers used for initial identification were designed for the L polymerase gene on the L segment (L4160F, GCA GAR TTY AAA TCI AGA TT; L4393R, CCR TYI ASC CAR TCT ITI ACA TC; L4292F, GAT CAT TCI RTY GCI AAT GG; L4841R, CAI AII CCT ATA AAI CCW GAT G)  and the glycoprotein gene on the S segment (GP878+, GAC RTG CCW GGI GGI TAY TG; GP1126-, TAC CAA AAT TTG TGT ART TRC ART AIG G; GP1153+, CCT TAY TGY AAY TAC ACI AAA TTT TGG T; GP1396-, ATG TGY CTR TGI GTI GGI AW).
Reverse Transcription (RT) was done using 2.5 µL of RNA in a 25 µL total reaction volume and AMV RT (Promega Biosciences, San Luis Obispo, CA) at 42°C for 1 hr. Subsequent PCR amplification using FastStart Taq DNA Polymerase with GC-rich solution (Roche) was performed using 5 µL of the completed RT reaction in a 25 µL reaction volume with the following cycling conditions: 2 min at 95°C, (36 cycles of 1 min at 95°C, 1 min at 45°C, 2 min at 72°C), and a final elongation of 10 min at 72°C. Resulting DNA products were visualized and purified using a 1% agarose gel, and the Qiagen Gel Extraction Kit (Qiagen, Valencia, CA). PCR products were sequenced directly (without cloning) using the corresponding primers in a BigDye Terminator v3.1 reaction on a 3100 Genetic Analyzer (Applied Biosystems, Foster City, CA). Sequence was further analyzed using Sequencher (Gene Codes Corporation, Ann Arbor, MI).
Complete Genome Amplification & Analysis
To obtain full length sequence for each segment, alignments of all New World arenavirus complete genomes were used to design primers for the conserved regions (available upon request). The full length S segment was generated following the Thermoscript RT-PCR system's directions (Invitrogen, Carlsbad, CA) and using the 19C primer . Reverse transcription was conducted at 55°C, while the PCR profile was the same as stated above with an increased extension time of 4 minutes.
Different fractions of the full-length L RNA were amplified using 2-step or 1-step RT-PCR protocols and following the manufacturer recommendations. Briefly, cDNA was synthesized in the first approach using 10 µl of purified RNA, specific primers, dNTPs and Superscript III (Invitrogen) in 20 µl reactions. Amplification reactions were done using 5 µl of cDNA, specific primers, dNTPs and Platinum Taq DNA polymerase High Fidelity (Invitrogen) in 50 ul reactions. Alternatively, 1-step RT-PCR were performed using 5 µl of RNA, dNTPs and the enzyme blend provided by the SuperScript III One-Step RT-PCR System with Platinum Taq High Fidelity (Invitrogen) in a 50 µl reactions. Amplification reactions were analyzed in TBE/agarose gels and DNA bands purified using QIAquick Gel Extraction Kit (Qiagen). Sequencing reactions were done as described above.
All full length S and L segment sequences available in Genbank were used to compute pairwise uncorrected genetic distances using PAUP 4.0b10 (Sinauer Associates) for the following viruses: Allpahuayo, Amapari, Chapare, Flexal, Guanarito, Junín, Machupo, Paraná, Pichindé, Pirital, Sabiá, Tacaribe, Tamiami, and Whitewater Arroyo.
A representative sub-set of full length sequences (omitting multiple near identical variants of the same virus) were included in a Bayesian phylogenetic analysis. Sequence alignments were done with ClustalX  with manual adjustments and phylogenetic analysis was done with MrBayes3.1.2  using the GTR+I+G model in 2 runs of 500,000 generations using the sequence of Pichindé virus as the outgroup.
The authors would like to thank Debbie L. Cannon and Jennifer B. Oliver for technical assistance, and Dr. Serena Reeder for her excellent assistance with the phylogenetic analysis presented in this work.
Conceived and designed the experiments: CA JO SN. Performed the experiments: BE CA JC. Analyzed the data: BE RA PB EV CA JV JC PR TK JO SN. Contributed reagents/materials/analysis tools: SD RA PB EV JV JC PR JO. Wrote the paper: BE SN. Initial clinician that reported the case, described the clinical course and provided specimens: SD. Clinicians involved in the case investigation: RA and EV. Participated in field investigation: PB. Participated in case investigation: JV. Edited the paper: TK. Led NMRCD investigation: JO.
- 1. Buchmeier MJ (2007) Arenaviridae. In: Knipe DMaH PM, editor. Fields Virology. 5th ed. Philadelphia: Lippincott, Williams and Wilkins. pp. 1791–1828.
- 2. Gonzalez JP, Emonet S, de Lamballerie X, Charrel R (2007) Arenaviruses. Curr Top Microbiol Immunol 315: 253–288.
- 3. Bowen MD, Peters CJ, Nichol ST (1997) Phylogenetic analysis of the Arenaviridae: patterns of virus evolution and evidence for cospeciation between arenaviruses and their rodent hosts. Molecular Phylogenetics & Evolution 8: 301–316.
- 4. Archer AM, Rico-Hesse R (2002) High genetic divergence and recombination in Arenaviruses from the Americas. Virology 304: 274–281.
- 5. Charrel RN, Lemasson JJ, Garbutt M, Khelifa R, De Micco P, et al. (2003) New insights into the evolutionary relationships between arenaviruses provided by comparative analysis of small and large segment sequences. Virology 317: 191–196.
- 6. Salvato MS, Clegg JCS, Buchmeier MJ, Charrel RN, Gonzalez JP, Lukashevich IS, Peters CJ, Rico-Hesse R, Romanowski V (2005) Family Arenaviridae. In: Fauquet CM, Mayo MA, Maniloff J, Desselberger U, Ball LA, editors. Virus Taxonomy: Eighth Report of the International Committee on Taxonomy of Viruses. Elsevier Academic Press. pp. 725–733.
- 7. Johnson KM, Kuns ML, Mackenzie RB, Webb PA, Yunker CE (1966) Isolation of Machupo virus from wild rodent Calomys callosus. Am J Trop Med Hyg 15: 103–106.
- 8. Murphy FA, Webb PA, Johnson KM, Whitfield SG, Chappell WA (1970) Arenoviruses in Vero cells: ultrastructural studies. J Virol 6: 507–518.
- 9. Salazar-Bravo J, Dragoo JW, Bowen MD, Peters CJ, Ksiazek TG, et al. (2002) Natural nidality in Bolivian hemorrhagic fever and the systematics of the reservoir species. Infection, Genetics & Evolution 1: 191–199.
- 10. Bowen MD, Peters CJ, Nichol ST (1996) The phylogeny of New World (Tacaribe complex) arenaviruses. Virology 219: 285–290.
- 11. Weaver SC, Salas RA, de Manzione N, Fulhorst CF, Duno G, et al. (2000) Guanarito virus (Arenaviridae) isolates from endemic and outlying localities in Venezuela: sequence comparisons among and within strains isolated from Venezuelan hemorrhagic fever patients and rodents. Virology 266: 189–195.
- 12. Cajimat MN, Fulhorst CF (2004) Phylogeny of the Venezuelan arenaviruses. Virus Research 102: 199–206.
- 13. Weaver SC, Salas RA, de Manzione N, Fulhorst CF, Travasos da Rosa AP, et al. (2001) Extreme genetic diversity among Pirital virus (Arenaviridae) isolates from western Venezuela. Virology 285: 110–118.
- 14. Cajimat MN, Milazzo ML, Hess BD, Rood MP, Fulhorst CF (2007) Principal host relationships and evolutionary history of the North American arenaviruses. Virology 367: 235–243.
- 15. Moncayo AC, Hice CL, Watts DM, Travassos de Rosa AP, Guzman H, et al. (2001) Allpahuayo virus: a newly recognized arenavirus (arenaviridae) from arboreal rice rats (oecomys bicolor and oecomys paricola) in northeastern peru. Virology 284: 277–286.
- 16. Flanagan ML, Oldenburg J, Reignier T, Holt N, Hamilton GA, et al. (2008) New world Clade B arenaviruses can use transferrin receptor 1 (TfR1)-dependent and -independent entry pathways, and glycoproteins from human pathogenic strains are associated with the use of TfR1. J Virol 82: 938–948.
- 17. Radoshitzky SR, Abraham J, Spiropoulou CF, Kuhn JH, Nguyen D, et al. (2007) Transferrin receptor 1 is a cellular receptor for New World haemorrhagic fever arenaviruses. Nature 446: 92–96.
- 18. Lisieux T, Coimbra M, Nassar ES, Burattini MN, de Souza LT, et al. (1994) New arenavirus isolated in Brazil. Lancet 343: 391–392.
- 19. Fischer SA, Graham MB, Kuehnert MJ, Kotton CN, Srinivasan A, et al. (2006) Transmission of lymphocytic choriomeningitis virus by organ transplantation. N Engl J Med 354: 2235–2249.
- 20. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG (1997) The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res 25: 4876–4882.
- 21. Ronquist F, Huelsenbeck JP (2003) MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19: 1572–1574.