Arenavirus Quasispecies and Their Biological Implications.

The family Arenaviridae currently comprises over 20 viral species, each of them associated with a main rodent species as the natural reservoir and in one case possibly phyllostomid bats. Moreover, recent findings have documented a divergent group of arenaviruses in captive alethinophidian snakes. Human infections occur through mucosal exposure to aerosols or by direct contact of abraded skin with infectious materials. Arenaviruses merit interest both as highly tractable experimental model systems to study acute and persistent infections and as clinically important human pathogens including Lassa (LASV) and Junin (JUNV) viruses, the causative agents of Lassa and Argentine hemorrhagic fevers (AHFs), respectively, for which there are no FDA-licensed vaccines, and current therapy is limited to an off-label use of ribavirin (Rib) that has significant limitations. Arenaviruses are enveloped viruses with a bi-segmented negative strand (NS) RNA genome. Each genome segment, L (ca 7.3 kb) and S (ca 3.5 kb), uses an ambisense coding strategy to direct the synthesis of two polypeptides in opposite orientation, separated by a noncoding intergenic region (IGR). The S genomic RNA encodes the virus nucleoprotein (NP) and the precursor (GPC) of the virus surface glycoprotein that mediates virus receptor recognition and cell entry via endocytosis. The L genome RNA encodes the viral RNA-dependent RNA polymerase (RdRp, or L polymerase) and the small (ca 11 kDa) RING finger protein Z that has functions of a bona fide matrix protein including directing virus budding. Arenaviruses were thought to be relatively stable genetically with intra- and interspecies amino acid sequence identities of 90-95 % and 44-63 %, respectively. However, recent evidence has documented extensive arenavirus genetic variability in the field. Moreover, dramatic phenotypic differences have been documented among closely related LCMV isolates. These data provide strong evidence of viral quasispecies involvement in arenavirus adaptability and pathogenesis. Here, we will review several aspects of the molecular biology of arenaviruses, phylogeny and evolution, and quasispecies dynamics of arenavirus populations for a better understanding of arenavirus pathogenesis, as well as for the development of novel antiviral strategies to combat arenavirus infections.

Arenaviruses as Both Important Human Pathogens and Highly Tractable Experimental Systems for the Investigation of Virus–Host Interactions

Arenaviruses and Their Impact in Human Health

Arenaviruses cause chronic infections of rodents across the world, and human infections occur through mucosal exposure to aerosols or by direct contact of abraded skin with infectious materials (Buchmeier et al. 2007). Both viral and host factors contribute to a variable outcome of arenavirus infection, ranging from virus control and clearance by the host defenses to chronic infection in the absence of clinical symptoms to severe disease. Several arenaviruses cause hemorrhagic fever (HF) disease in humans and pose a serious public health problem in their endemic regions. Thus, the Old World (OW) arenavirus, Lassa virus (LASV), infects several hundred thousand individuals yearly in West Africa resulting in a high number of Lassa fever (LF)cases associated with high morbidity and mortality (Bray 2005; Geisbert and Jahrling 2004). The current LF endemic areas cover approximately 80 % of Sierra Leone and Liberia; 50 % of Guinea; 40 % of Nigeria; 30 % of each of Côte d’Ivoire, Togo, and Benin; and 10 % of Ghana (Fichet-Calvet and Rogers 2009), and recent studies have found evidence for expansion of LASV endemic regions into Mali, Senegal, and the Democratic Republic of Congo (Sogoba et al. 2012). Thus, it is quite likely that population at risk of LF includes most of West Africa and can be as high as 200 million people (Richmond and Baglole 2003); hence, with dengue fever exception, the estimated global burden of LF is the highest among viral HF (Falzarano and Feldmann 2013). Notably, increased traveling to and from endemic regions has led to the importation of LF cases into non-endemic metropolitan areas around the globe (Freedman and Woodall 1999; Isaacson 2001). Likewise, the New World (NW) arenavirus, Junin virus (JUNV), causes Argentine HF (AHF), a disease endemic to the Argentinean Pampas with hemorrhagic and neurological manifestations and a case fatality of 15–30 % (Peters 2002). In addition, evidence indicates that the worldwide-distributed prototypic arenavirus, lymphocytic choriomeningitis virus (LCMV), is a neglected human pathogen of clinical significance (Jahrling and Peters 1992; Mets et al. 2000; Palacios et al. 2008; Peters 2006).

Public health concerns about arenavirus infections are exacerbated because of the lack of FDA-approved vaccines against HF arenaviral diseases and current therapy to combat arenavirus infections being limited to an off-label use of the nucleoside analogue ribavirin (Rib) (Huggins 1989; Kilgore et al. 1995, 1997; McCormick et al. 1986; McKee et al. 1988). However, Rib is only partially effective, and its use has several limitations including the need of intravenous (iv) and early administration for optimal efficacy, and significant side effects (Damonte and Coto 2002). Although the live-attenuated Candid1 strain of JUNV has been shown to be an effective vaccine against AHF (Enria and Barrera Oro 2002), outside Argentina Candid1 has only investigational new drug (IND) status and it does not protect against LF. Despite efforts dedicated to the development and testing of LASV vaccines, not a single LASV vaccine candidate has entered a clinical trial although several vaccine platforms, including the MOPV/LASV reassortant ML29 and recombinant VSV and vaccinia virus expressing viral antigens, have been successfully tested in non-human primates (Falzarano and Feldmann 2013; Lukashevich 2012).

LCMV Infection of the Mouse as the Rosetta Stone of Virus–Host Interaction

Studies using LCMV infection of its natural host, the mouse, have led to major advances in virology and immunology that apply universally to other microbial infections and viral infections of humans including virus-induced immunopathological disease, MHC restriction, T cell-mediated killing, and effectiveness of adoptive immune T cell therapy in clearing viral infection (Zinkernagel 2002; Oldstone 2002). Likewise, the initial description of the contribution to viral persistence of negative immune regulators such as PDL-1 was made first with LCMV in mice (Barber et al. 2006). The outcome of LCMV infection of its natural host, the mouse, varies dramatically depending on the age, immunocompetence and genetic background of the host, as well as the route of infection, and the strain and dose of infecting virus (Oldstone 2002; Zinkernagel 2002; Buchmeier et al. 2007). Thus, iv inoculation of adult immunocompetent mice with LCMV Armstrong (ARM) strain results in an acute infection that induces a protective immune response that mediates virus clearance in 10–14 days, a process predominantly mediated by virus-specific CD8+ cytotoxic T lymphocytes (CTL). In contrast, iv inoculation with a high (≥2 × 106 PFU) dose of the immunosuppressive LCMV clone 13 (Cl-13) strain, which differs from ARM only at three amino acid positions, causes persistent infection associated with generalized immune suppression (Ahmed and Oldstone 1988; Ahmed et al. 1984, 1988). These dramatic phenotypic differences between genetically very closely related viruses provide investigators with a superb experimental system to examine both host immune mechanisms of virus control and viral counteracting strategies.

Molecular and Cell Biology of Arenaviruses

Arenavirus Genome Organization

Arenaviruses are enveloped viruses with a bi-segmented negative-stranded (NS) RNA genome and a life cycle restricted to the cell cytoplasm (Buchmeier et al. 2007). Virions appear pleomorphic when examined by cryoelectron microscopy, ranging in size from 40 to more than 200 nm in diameter. Each genomic viral RNA segment, L (ca 7.3 kb) and S (ca 3.5 kb), uses an ambisense coding strategy to direct the synthesis of two polypeptides in opposite orientation, separated by a noncoding IGR that has been shown to act as a bona fide transcription termination signal (Fig. 1). There are significant differences in sequence and predicted folded structure between the S and L IGR, but among strains of the same arenavirus species, the S and L IGR sequences are highly conserved. The S RNA encodes the viral NP and the glycoprotein precursor (GPC) that is posttranslationally cleaved by the cellular site 1 protease (S1P) to yield the two mature virion glycoproteins GP1 and GP2 that form the spikes that decorate the virus surface (Beyer et al. 2003; Lenz et al. 2001; Pinschewer et al. 2003). GP1 mediates virus receptor recognition and cell entry via endocytosis, whereas GP2 mediates the pH-dependent fusion event required to release the virus ribonucleoprotein (RNP) core into the cytoplasm of infected cells (Buchmeier et al. 2007). The L RNA encodes the viral RNA-dependent RNA polymerase (RdRp, or L polymerase) and the small RING finger protein Z that has functions of a bona fide matrix protein including directing virus budding (Perez et al. 2003; Strecker et al. 2003).

Fig. 1

Arenavirus genome organization and virion structure. Arenaviruses are enveloped viruses with a bi-segmented NS RNA genome. Each genome segment uses an ambisense coding strategy to direct the synthesis of two viral polypeptides. The small (S, ca 3.5 kb) segment encodes for the glycoprotein precursor (GPC) and nucleoprotein (NP). Posttranslational processing of GPC by the cellular protease S1P results in the production of GP1 and GP2. The large (L, ca 7.3 kb) segment encodes for the virus RNA-dependent RNA polymerase (L) and a small RING finger protein (Z) that has the properties of the bona fide matrix proteins (M) found in many enveloped NS RNA viruses

Arenaviruses exhibit high degree of sequence conservation at the genome 3′-termini (17 out of 19 nt are identical), and, as with other NS RNA viruses, arenavirus genome termini exhibit terminal complementarity with the 5′- and 3′-ends of both L and S genome segments predicted to form panhandle structures (Buchmeier et al. 2007). This prediction is supported by electron microscopy (EM) data showing the presence within the virion particles of circular RNP complexes (Young and Howard 1983). It has been proposed, based on EM studies, that host ribosomes can be also incorporated into virions, but the biological significance of this remains to be determined (Muller et al. 1983).

Arenavirus Life Cycle

Cell entry

Arenavirus cell entry occurs via receptor-mediated endocytosis, which involves a pH-dependent fusion event between viral and cellular membranes within the acidic environment of the endosome (Rojek and Kunz 2008; Rojek et al. 2008a, b). The widely expressed cell surface receptor for extracellular matrix proteins (ECM), alpha-dystroglycan (αDG), has been identified as a main receptor for LCMV, LASV, and several other arenaviruses (Kunz et al. 2002, 2004). To function as an arenavirus receptor, αDG needs to undergo several posttranslational modifications including a critical like-acetylglucosaminyltransferase (LARGE)-dependent glycosylation (Kunz et al. 2005). Intriguingly, recent genome-wide studies have uncovered evidence of positive selection of specific LARGE alleles within the Yoruba people living in regions of West Africa where LASV is endemic (Sabeti et al. 2007). Endocytic vesicles of the OW arenavirus LCMV were found to be non-coated, and LCMV cell entry process was reported to be cholesterol-dependent but clathrin-, dynamin-, caveolin-, ARF6-, flotillin-, and actin-independent (Rojek et al. 2008b). OW arenavirus cell entry follows an endocytic pathway that bypasses classical Rab5/EEA-1 positive early endosomes, but merges into the classical endocytotic pathway at the level of the multivesicular body (MVB) and is then delivered to the late endosome in an endosomal sorting complex required for transport (ESCRT)-dependent process (Pasqual et al. 2011). The unusual endocytotic pathway followed by LCMV and LASV may reflect the natural cellular trafficking of their receptor αDG. Nevertheless, it cannot be ruled out that binding of GP to αDG could trigger a novel route of entry for αDG that is tailored specifically to benefit virus cell entry as previously documented for other pathogens. Differently to OW arenaviruses, many NW arenaviruses, including the HFA JUNV, use the human transferrin receptor to initiate cell entry via clathrin-mediated endocytosis (Rojek et al. 2008b).

Expression and replication of the viral genome

The fusion between viral and cellular membranes within the late endosome compartment releases the viral RNP into the cytoplasm where both RNA replication and gene transcription of the genome take place (Buchmeier et al. 2007). The NP and L coding regions are transcribed into a genomic complementary mRNA, whereas the GPC and Z coding regions are not translated directly from genomic RNA, but rather from genomic sense mRNAs that are transcribed using as templates the corresponding antigenome RNA species, which also function as replicative intermediates. Arenavirus mRNAs are non-polyadenylate and contain at their 5′-ends a short (4- to 5-nucleotide) non-templated sequence and a cap structure that are likely obtained from cellular mRNAs via a cap-snatching mechanism whose details remain to be determined, but the endonuclease activity associated with the arenavirus L polymerase likely contributes to this process (Morin et al. 2010). Transcription termination of viral mRNAs was mapped to multiple sites within the distal side of the IGR (Meyer et al. 2002; Meyer and Southern 1994; Tortorici et al. 2001), suggesting that the IGR acts as a bona fide transcription termination signal for the virus polymerase, which has been confirmed using cell-based minireplicon systems (Pinschewer et al. 2005).

NP and L are the only viral trans-acting factors required for efficient arenavirus RNA synthesis, whereas Z has been shown to exert a dose-dependent inhibitory effect on the RNA biosynthetic activities of the arenavirus polymerase complex (Cornu and de la Torre 2001, 2002). The arenavirus L protein has the characteristic sequence motifs conserved among the RdRp (L proteins) of NS RNA viruses (Poch et al. 1989). Bioinformatic analysis and functional assays have provided evidence that LASV L protein is organized into three distinct structural domains and that the N- and C-parts of L are able to functionally trans-complement each other (Brunotte et al. 2011). Notably, the recent EM characterization of a functional MACV L protein has revealed a core ring domain decorated by appendages, which may reflect a modular organization of the arenavirus polymerase (Kranzusch and Whelan 2011). The NP is the main structural element of the viral RNP that directs viral RNA synthesis. NP exhibits also a type I interferon (IFN-I) counteracting activity that was mapped to its C-terminus (Martinez-Sobrido et al. 2006, 2007, 2009). Two recently determined crystal structures of LASV NP at 1.80 Å (Qi et al. 2010) and 1.5 Å (Hastie et al. 2011) resolution identified distinct N- and C-terminal domains connected by a flexible linker. The N-terminal domain was proposed to contain a potential cap-binding site that could provide the host-derived primers to initiate transcription by the virus polymerase (Qi et al. 2010). However, analysis of the structure of the N-terminal domain of LASV NP in complex with ssRNA suggests that this originally proposed cap-binding site likely corresponds to a binding site for the viral genome (Hastie et al. 2011). The C-terminal of NP has a folding that mimics that of the DEDDH family of 3′–5′ exoribonucleases (Hastie et al. 2011; Qi et al. 2010). Functional studies confirmed the 3′–5′ exoribonuclease activity of LASV NP, which was proposed to be critical for the anti-IFN activity of NP but dispensable for the role of NP on replication and transcription of the viral genome (Hastie et al. 2011; Qi et al. 2010). However, LCMV with a mutant NP lacking the 3′–5′ exoribonuclease had a large decrease in fitness during its replication in IFN-deficient Vero cells (Martinez-Sobrido et al. 2009).

Assembly and buddingBudding

Assembly and cell release of infectious arenavirus progeny require both Z and correct processing of GPC into GP1 and GP2 by the cellular site 1 protease (S1P) (Kunz et al. 2003b; Perez et al. 2003; Pinschewer et al. 2003; Strecker et al. 2003; Urata et al. 2006). GP1 and GP2 together with a stable 58 amino acid long signal peptide (SSP) form the functional GP complex. SSP has been implicated in different aspects of the trafficking and function of the viral envelope glycoproteins (Saunders et al. 2007; York and Nunberg 2006, 2007; York et al. 2004). GP1 located at the top of the spike and held in place by ionic interactions with the N-terminus of the transmembrane GP-2 mediates virus interaction with host cell surface receptors, whereas GP2 exhibits similarities with fusion-active domains of other viral envelope proteins. The arenavirus RING finger protein Z is a structural component of the virion that has no homologue among other known NS RNA viruses. In LCMV-infected cells, Z has been shown to interact with several cellular proteins including the promyelocytic leukemia (PML) protein and the eukaryote translation initiation factor 4E (eIF4E), which have been proposed to contribute to the non-cytolytic nature of LCMV infection and repression of cap-dependent translation, respectively (Borden et al. 1998; Campbell Dwyer et al. 2000; Djavani et al. 2001; Volpon et al. 2010). For most enveloped viruses, a matrix (M) protein is involved in organizing the virion components prior to assembly. Interestingly, arenaviruses do not have an obvious counterpart of M. However, structural and functional studies indicated that Z was the arenavirus counterpart of the M protein found in many other enveloped NS RNA viruses that play a critical role in assembly and cell release of mature infectious virions (Perez et al. 2003; Strecker et al. 2003; Urata et al. 2006). Accordingly, Z has been shown to be the main driving force of arenavirus budding, a process mediated by the Z late (L) domain motifs: PTAP and PPPY similar to those known to control budding of several other viruses including HIV and Ebola virus, via interaction with specific host cell proteins (Perez et al. 2003; Strecker et al. 2003; Urata et al. 2006). Consistent with this observation, Z exhibited features characteristic of bona fide budding proteins: (1) ability to bud from cells by itself and (2) substituted efficiently for other L domain. Targeting of Z to the plasma membrane, the location of arenavirus budding, strictly required its myristoylation (Perez et al. 2004; Strecker et al. 2006). Results derived from cryoelectron microscopy of arenavirus particles (Neuman et al. 2005) were also consistent with the role of Z as a functional M protein.

Arenavirus Reverse Genetics

In contrast to positive-stranded RNA viruses, deproteinized genomic and antigenomic RNAs of NS RNA viruses, such as LCMV, cannot function as mRNAs and are not infectious. This reflects the fact that the template of the polymerases of NS RNA viruses is exclusively a nucleocapsid consisting of the genomic RNA tightly encapsidated by the NP, which associated with the virus polymerase proteins forms a RNP complex. Thus, generation of biologically active synthetic NS viruses from cDNA requires trans-complementation by all viral proteins involved in virus replication and transcription. These considerations hindered the application of recombinant DNA technology to the genetic analysis of these viruses. However, significant progress has been made in this area, and for members of all families of NS RNA viruses, it has become possible to develop cell-based minireplicon systems and to rescue infectious virus from cloned cDNAs (Neumann and Kawaoka 2004). The use of reverse genetics approaches has revolutionized the analysis of cis-acting sequences and trans-acting proteins required for virus replication, transcription, maturation, and budding. In addition, the possibility to generate predetermined specific mutations within the virus genome and analyze their phenotypic expression in vivo in the context of the virus natural infection is contributing very significantly to the elucidation of the molecular mechanisms underlying virus–host interactions at the cellular and molecular levels, which has provided investigators with novel and powerful approaches for the investigation of viral pathogenesis. In addition, these developments have also paved the way for engineering these viruses for vaccine and gene therapy purposes (Subbarao and Katz 2004; von Messling and Cattaneo 2004).

Arenavirus minireplicon (MR) systems

Cell-based minireplicon (MR) systems have been developed for a variety of arenaviruses (Emonet et al. 2011a), and their use provided conclusive evidence that NP and L are the minimal viral trans-acting factors required for efficient RNA synthesis mediated by the arenavirus polymerase complex (Hass et al. 2004; Lee et al. 2000; Lopez et al. 2001). Moreover, LCMV MR assays uncovered that oligomerization of L was required for the activity of the LCMV polymerase (Sanchez and de la Torre 2005). Consistent with this finding, LASV L protein was shown to contain both N- and C-termini domains that mediate L–L interaction (Brunotte et al. 2011). MR-based assays facilitated mutation-function studies that identified residues within the N-terminus of L that played a critical role in synthesis of viral mRNA but not in RNA replication (Lelke et al. 2010). Moreover, MR-based assays together with structural data uncovered within the N-terminus of LASV L protein, an endonuclease functional domain that was critically required for arenavirus transcription (Morin et al. 2010). Mutation-function analysis of the genome 5′- and 3′-termini using LCMV and LASV MR-based assays demonstrated that the activity of the arenavirus genomic promoter requires both sequence specificity within the highly conserved 3′-terminal 19 nt of arenavirus genomes and the integrity of the predicted panhandle structure formed between the 5′- and 3′-termini of viral genome RNAs (Hass et al. 2006; Perez and de la Torre 2002). Likewise, MR-based assays provided direct experimental confirmation that the IGR is a bona fide transcription termination signal (Pinschewer et al. 2005). These studies also uncovered that Z exhibited a dose-dependent inhibitory effect on both RNA replication and gene transcription by the arenavirus polymerase complex (Cornu and de la Torre 2001, 2002; Cornu et al. 2004; Hass et al. 2004; Lopez et al. 2001).

Generation of infectious recombinant arenaviruses

Generation of infectious virus from cloned cDNAs has been reported for LCMV (Flatz et al. 2006; Sanchez and de la Torre 2006), PICV (Lan et al. 2009), JUNV (Albarino et al. 2009; Emonet et al. 2011b), and LASV (Albarino et al. 2011). Both RNA polymerase I and T7 RNA polymerase (T7RP)-based systems haven been successfully used, and with similar efficiencies, to direct intracellular synthesis of L and S genome RNA species required for the rescue of infectious recombinant arenaviruses (Emonet et al. 2011a). Intriguingly, the rescue of PICV and JUNV using the T7RP-based system did not require plasmid-supplied viral NP and L proteins, indicating that T7RP-mediated RNA synthesis produced both viral antigenome RNA species that were substrate for encapsidation and replication, and mRNAs that serve as template to produce levels of NP and L sufficiently high to facilitate virus rescue. This phenomenon has been reported for several other negative-sense RNA viruses (Emonet et al. 2011a). Reverse genetics approaches similar to those used to rescue infectious LCMV, PICV, JUNV, and LASV from cloned cDNAs should be now applicable to other arenaviruses. The ability to generate recombinant arenaviruses with predetermined specific mutations and analyze their phenotypic expression in the context of the natural course of infection has opened new opportunities to investigate arenavirus–host interactions that influence a variable infection outcome, ranging from virus control and clearance by the host defenses to long-term chronic infection associated with subclinical disease, and severe acute disease including HF. Thus, the use of rLCMV/VSVG uncovered the arenavirus GP as a viral Achilles’ heal and provided the foundations for a novel strategy to develop safe and effective live-attenuated arenavirus vaccines (Pinschewer et al. 2003). Likewise, rLCMV/VSVG was very instrumental in facilitating studies aimed at understanding the regulation of CD8+ T cell function within the infected brain (Pinschewer et al. 2010), as well as how viruses can induce organ-specific immune disease in the absence of molecular mimicry and without disruption of self-tolerance (Merkler et al. 2006). Other engineered rLCM viruses have been used to study LASV cell entry pathway (Rojek et al. 2008b) and the role of NP in the inhibition of induction of IFN-I by LCMV (Martinez-Sobrido et al. 2009). Similarly, studies aimed at examining the critical role played by the S1P-mediated processing of arenavirus GPC in the virus life cycle were aided by the use of recombinant viruses where the S1P cleavage site within GPC was replaced by a canonical furin cleavage site (Albarino et al. 2009; Rojek et al. 2010; Urata et al. 2011). The rescue of attenuated and virulent forms of PICV (Lan et al. 2009; Liang et al. 2009) or the Docile and Aggressive strains of LCMV (Chen et al. 2008) has allowed for the identification of specific genetic determinants of virus virulence. Moreover, it has been possible to generate single-cycle infectious, reporter-expressing, recombinant LCMV in which the GPC open reading frame (ORF) was replaced by GFP (rLCMVDGP/GFP) (Rodrigo et al. 2011). Genetic complementation with plasmids or stable cell lines expressing arenavirus GPCs of interest produces the corresponding GPC-pseudotyped rLCMV∆GP/GFP that can be used to evaluate antibody responses to HF arenaviruses using a Biosafety Level 2 (BSL2) platform.

Despite recombinant arenaviruses could be generated very efficiently, the ability to rescue arenaviruses expressing additional genes of interest posed unexpected difficulties. This problem was overcome by the generation of tri-segmented rLCMV (r3LCMV) containing 1L and 2S segments. Each of the two S segments was altered to replace one of the two viral ORF, GPC, or NP, by a gene of interest (GOI) (Emonet et al. 2009b). The rationale behind this approach was that the physical separation of GP and NP into two different S segments would represent a strong selective pressure to select and maintain a virus capable of packaging 1L and 2S segments. A variety of r3LCMV has been rescued that express one or two additional GOI. Depending on the GOI expressed (protein function, size of the gene), these r3LCMV showed little or no attenuation in cultured cells and they exhibited long-term genetic stability as reflected by unaltered expression levels during serial virus passages of the GOI incorporated into the S segment. The use of r3LCMV expressing appropriate reporter genes should facilitate the development of chemical screens to identify antiviral drugs, as well as siRNA-based genetic screens to identify host cell genes contributing to the different steps of the arenavirus life cycle. However, the use of r3LCMV to study virus–host interactions in mice has encountered some limitations due to the observation that r3LCMV with WT growth properties in cultured cells usually exhibits attenuation in mice due to reasons that remain to be determined.

Genetic Variability of Arenaviruses

Organization and Phylogenetic Relationships of the Arenaviridae Family

Until recently, the family Arenaviridae has consisted of a single genus, Arenavirus, that includes 25 recognized species. Based on their antigenic properties, arenaviruses are divided into two distinct groups: The OW arenaviruses (Lassa–lymphocytic choriomeningitis serocomplex) include viruses indigenous to Africa and the ubiquitous LCMV, and the NW arenaviruses (Tacaribe serocomplex) include viruses indigenous to the Americas. Subsequent genetic analyses based on sequencing data supported very well the previously defined antigenic grouping (Bowen et al. 1997). The monogeneric phylogenetic view of arenaviruses has recently shifted dramatically with the discovery of a divergent group of arenaviruses in captive alethinophidian snakes (Bodewes et al. 2013; Hetzel et al. 2013; Stenglein et al. 2012). Preliminary phylogenetic studies indicate that these reptilian arenaviruses constitute a sister clade to mammalian arenaviruses, and have been proposed to represent a new genus within the family Arenaviridae.

OW arenaviruses form one monophyletic group that is deeply rooted to the three identified NW arenavirus phylogenetic groups, clades A, B, and C. Clade B contains the human pathogenic viruses Chapare (CHPV), Guanarito (GTOV), JUNV, Machupo (MACV), and Sabiá (SABV), as well as several nonpathogenic including Tacaribe virus (TCRV) (Fig. 2). Among OW arenaviruses, LASV, Mobala virus (MOBV), and Mopeia virus (MOPV) are monophyletic, while Ippy virus (IPPYV) and LCMV are more distantly related. The recently discovered Lujo virus (LUJV), endemic in Zambia, is most closely related to OW arenaviruses, but its GPC contains elements of NW arenavirus sequences. Both OW and NW arenaviruses have been evolved from a common ancestor, developing a wide range of antigenic specificity, host species range, pathogenicity, and virulence (Zapata and Salvato 2013).

Fig. 2

Rooted tree of arenavirus S segment. Phylogeny reconstruction by neighbor joining with 10,000 bootstrap replications under maximum composite likelihood substitution model was done with Mega 6 (Tamura et al. 2013). Bootstrap values below 80 % are hidden. A similar tree was obtained for the L segment. GenBank accession numbers of S segment (all complete sequences) used to generate the tree, as well as species names of rodent reservoirs were obtained from previously published data (King et al. 2012), with the following exceptions: Skinner tank virus (SKTV, EU123328.1) from Neotoma mexicana (USA); Tonto creek virus (TTCV, EF619034.1) from Neotoma albigula; Big brushy tank virus (BBTV, EF619035.1) from Neotoma albigula, Arizona (USA); Catarina virus (CTNV, DQ865244) from Neotoma micropus (Cajimat et al. 2013); and Luna virus (LUNV-LSK-1, AB693148.1) from Mastomys natalensis, Zambia (Ishii et al. 2012). Colored dots next virus isolates indicate reservoir species: green, neotominae subfamily; blue, sigmodontinae subfamily; red, murinae subfamily; cian, the bat Artibeus spp, proposed reservoir of Tacaribe virus (TCRV-p2b2); and magenta, the Boa constrictor snake, reservoir of the recently reported putative arenaviruses University of Helsinki virus (UHV-1) and inclusion body disease virus (IBDV-NL3) (Bodewes et al. 2013; Hetzel et al. 2013). Absence of colored dot indicates unknown reservoir

As with other RNA viruses, arenavirus replication is mediated by an error-prone polymerase complex, which together with reassortment events of its segmented genome generates the genetic variability that provides arenavirus with the potential for rapid evolution and high capability of adaptation to new environments. The contribution of recombination to arenavirus genetic variability has not been explored at length, but similarly to other NS RNA viruses, recombination events among arenavirus genomes appear to be extremely rare. Nevertheless, both reassortment and recombination events need to be taken into account to understand the evolution and population dynamics of arenaviruses (Albarino et al. 1998; Archer and Rico-Hesse 2002; Charrel and de Lamballerie 2003; Charrel et al. 2001, 2003). It is worth noting that arenaviruses exhibit a superinfection inhibition feature that diminishes the possibilities for recombination and reassortment events since both are dependent on two different arenaviruses infecting the same cell (Ellenberg et al. 2004, 2007; Moreno et al. 2012b). The mechanisms underlying virus superinfection inhibition vary among different viruses (Geleziunas et al. 1994; Horga et al. 2000; Huang et al. 2008; Morrison and McGinnes 1989). In the case of arenaviruses, high levels of Z protein inhibit de novo formation of a functional polymerase complex thus leading to inhibition of viral RNA synthesis (Cornu and de la Torre 2001, 2002; Cornu et al. 2004; Kranzusch and Whelan 2011; Loureiro et al. 2011). Hence, infected cells may express levels of Z that do not permit synthesis of RNA by the polymerase complex of a second infecting arenavirus (Ellenberg et al. 2004, 2007; Moreno et al. 2012b).

Mechanisms of Arenavirus Genetic Diversity

Arenavirus Mutation Frequencies

RNA viruses have high mutation rates, in the range of 10−3 to 10−5 (Drake and Holland 1999) due to the error-prone nature of their viral RdRp (Domingo 2007). Accordingly, early studies already suggested high genetic variability of arenavirus populations based on the observation that CTL escape mutants could be readily selected both in vivo (Pircher et al. 1990) and in cell culture (Aebischer et al. 1991; Lewicki et al. 1995). Likewise, emergence of neutralization-resistant LCMV mutants in persistently infected mice was found to be a rather frequent event associated with single point mutations (Ciurea et al. 2000, 2001). More recent studies have examined the mutation frequencies (the number of nucleotide substitutions and indels per nucleotide sequenced) in LCMV populations recovered after one or several passages in BHK-21 cells. These studies involved sequencing of multiple clones of several regions of the GP, NP, and L genes and found mutation frequencies in the range of 1.0 × 10−4 to 2.7 × 10−4 and 1.0 × 10−4 to 5.7 × 10−4 for populations after one and nine, respectively, cell culture passages (Grande-Perez et al. 2002). The highest mutation frequency detected corresponded to the NP region and the lowest to the L region. These values correlated with the corresponding values of Shannon entropy, parameter that measures the genetic heterogeneity of the viral population (Domingo 2006). Other studies have documented mutation frequencies for arenaviruses in the range of 2.6 × 10−4 to 5.5 × 10−4 mutations per nucleotide during serial passages in cell culture (Ruiz-Jarabo et al. 2003). In addition, mutation frequencies in the range of 1.2 × 10−4 to 3.5 × 10−4 substitutions per nucleotide were estimated for both the intracellular and released viral populations during LCMV persistence in BHK-21 cells (Grande-Perez et al. 2005b). These figures are consistent with a large body of literature documenting the genetic variability of both negative- and positive-stranded RNA viruses (Domingo 2007) and support the view that arenavirus populations are ensembles of genetically closely related mutant genomes, called quasispecies, subjected to natural selection and genetic drift (Sevilla and de la Torre 2006).

RNA Recombination in Arenaviruses

Discrepancies between NP- and GPC-based phylogenetic trees observed between clades A and B arenaviruses and North American (lineage A/Rec) arenaviruses including White Water Arroyo virus (WWAV) and the two newly discovered North America arenaviruses, Catarina and Skinner Tank viruses, led to hypothesize that a recombination event between members of clades A and B led to the emergence of lineage A/Rec (Archer and Rico-Hesse 2002; Charrel et al. 2001, 2002, 2003; Fulhorst et al. 2001). Additional recombination events have been suggested to explain the evolutionary relationships between Oliveros (OLIV) and Pichinde (PICV) viruses with OW arenaviruses (Albarino et al. 1998). However, the generation of recombinant arenaviruses requires cells to be simultaneously infected by two different arenaviruses, and although coinfection with two different arenaviruses has been reported in cultured cells (Ellenberg et al. 2004, 2007; Lukashevich et al. 2005; Moreno et al. 2012b; Riviere et al. 1986), in nature it may be less likely and further restricted by the superinfection exclusion feature associated with arenavirus infection (Ellenberg et al. 2004, 2007; Lukashevich et al. 2005; Moreno et al. 2012b; Riviere et al. 1986). Some arenaviruses share a common geographic localization, such as Guanarito and Pirital (PIRV) viruses in Venezuela (Fulhorst et al. 1999; Weaver et al. 2000); Junin and Oliveros viruses in Argentina (Mills et al. 1996); and Machupo, Latino, and Chapare viruses in Bolivia (Delgado et al. 2008). Likewise, different arenaviruses can some times infect the same rodent species, such as Guanarito and Pirital viruses, which were isolated both in cotton rats Sigmodon alstoni and in cane mice Zygodontomys brevicauda (Weaver et al. 2000). However, despite an extensive genetic characterization, the existence of a GTOV-PIRV recombinant virus has never been demonstrated (Cajimat and Fulhorst 2004; Fulhorst et al. 2008), which might be explained by superinfection exclusion among arenaviruses (Damonte et al. 1983; Ellenberg et al. 2004, 2007). Since no A/Rec arenavirus has been isolated in South America and neither clade A nor clade B arenavirus has been isolated in North America, it is difficult to explain how the recombination occurred. It is plausible that the two parental arenaviruses were originally present in North America but disappeared, either by natural extinction, migration or by being outcompeted by the recombinant arenavirus. The migration to North America of a recombinant arenavirus generated in South America seems improbable because South America arenaviruses most likely descend from an ancestral arenavirus located in North America (Cajimat et al. 2007). Moreover, North America arenaviruses have a long-term coevolutionary relationship with rodents from Neotominae subfamily, which are almost exclusively found in North America (Cajimat et al. 2007). It should be noted that clade C arenaviruses also exhibit a shift between GPC- and NP-based trees (Charrel et al. 2008), and it seems unlikely that such a rare event would be the origin of two different lineages. It seems more plausible that the phylogenetic information derived from arenavirus GPC is influenced by specific selection pressures (Cajimat et al. 2007). The reasons why GPC would be under a selection pressure different from the other arenavirus genes need to be determined. Likewise, the statistical significance of the recombination event found in North America and clade C arenaviruses using recombination detection softwares should be validated using an up-to-date arenavirus genetic database.

The recently discovered arenaviruses associated with cases of inclusion body disease (IBD) in snakes (Bodewes et al. 2013; Hetzel et al. 2013; Stenglein et al. 2012) might represent a very special case of heterologous recombination between an arenavirus and a filovirus within a coinfected host (snake). Phylogenetic analysis showed that snake-related arenaviruses constitute a monophyletic lineage separated from NW and OW arenaviruses (Bodewes et al. 2013; Stenglein et al. 2012). Intriguingly, the GPC of snake-related arenaviruses was more closely related to the GPC of filoviruses than to arenaviruses, but the natural history of events favoring this potential recombination event between filo- and arenaviruses remains to be elucidated.

Arenavirus Genomic Reassortments

Reassortment events require, like recombination, a coinfection event that can be readily recreated in cultured cells (Ellenberg et al. 2004, 2007; Lukashevich 1992; Lukashevich et al. 2005; Moreno et al. 2012b; Moshkoff et al. 2007; Riviere et al. 1986; Riviere and Oldstone 1986). However, in nature, there is only one documented case of an arenavirus reassortant event (Palacios et al. 2008). Extensive phylogenetic analysis of NW and OW arenaviruses revealed no proof of detectable genomic reassortment events during their evolution (Charrel et al. 2003), with the exception of an arenavirus isolated from three fatal cases of LCMV infection in transplantation patients who received organs from a single donor (Palacios et al. 2008). Phylogenetic analysis revealed a new arenavirus, in which L and S segments were related to LE and M1-M2 LCMV strains, respectively, suggesting that likely a reassortment event was involved in the generation of this arenavirus. The lower reassortant frequency in arenaviruses may be due to the superinfection exclusion exhibited by the members of this family (Ellenberg et al. 2004, 2007; Moreno et al. 2012b).

Since reassortment of arenaviruses in nature seems to be restricted to some specific combinations of L and S segments (Lukashevich 1992; Riviere et al. 1986), and recombination in arenaviruses is extremely infrequent (Archer and Rico-Hesse 2002; Charrel et al. 2002), mutations during genome replication by an error-prone polymerase arenaviruses are thought to be the main contributing factor to arenavirus genetic variability.

Arenavirus Hypermutation

Interferon-inducible adenosine deaminase that acts on double-stranded RNA (ADAR1) has been suggested to play a role in the innate immune response against LCMV (Zahn et al. 2007). ADARs are a family of dsRNA-binding enzymes present in animals that deaminate adenosine to inosine, and are involved in the regulation of a variety of posttranscriptional processes as well as in an antiviral response (Nishikura 2010). Since inosine and adenosine have different base-pairing properties, editing can alter protein-coding transcripts and noncoding targets including microRNA, small interfering RNA viral transcripts, and repeat elements (Hundley and Bass 2010; Mallela and Nishikura 2012). Upon LCMV infection of murine L929 cells, both expression and activity of ADAR1-L were upregulated and its editing activity was not inhibited by viral replication (Zahn et al. 2007). An adenosine-to-guanosine (A-to-G) transition bias was found in cDNA clones of the GPC gene derived from LCMV RNA isolated from L929 cells at seven days post-infection. Likewise, an A-to-G/U-to-C transition bias was observed in clones of LCMV genomic RNA species derived from RNA isolated from spleen of LCMV-infected C57BL/6 mice at four days post-infection, whereas no G-to-A/C-to-U transition was observed. Moreover, two of 54 analyzed clones showed an A-to-G or U-to-C hypermutation pattern in a very short stretch of the GP where 38 % of all adenosines or 14 % of uracils were mutated. In LCMV clones isolated from mice, no hypermutated LCMV clones were found. Approximately 50 % of the clones with an amino acid mutation were non-functional, suggesting that ADAR1-L-induced mutations in LCMV RNA reduce viral infectivity and protein function, thus supporting a contribution of ADAR1-L to the innate antiviral immune response. Nevertheless, multiplication of LCMV was not increased in mouse embryonic fibroblasts derived from ADAR1−/− mice, which raises some questions about a direct role of ADAR1 in restricting LCMV multiplication (Ward et al. 2011).

Hypermutation has been related to posttranscriptional editing by host adenosine deaminases of viral genomes from persistent infections of other RNA viruses such as measles virus (Cattaneo and Billeter 1992), human parainfluenza virus (Murphy et al. 1991; Rima et al. 2014), and respiratory syncytial viruses (Rueda et al. 1994). In these studies, the estimated mutation frequencies of U → C or A → G transitions were in the range of 3.3 × 10−2 to 1.0 × 10−1 mutations per nucleotide. The existence of hypermutated genomes had been hypothesized as proof of the genetic melting that occurs when the error catastrophe threshold is crossed as a result of chemical mutagenesis (Grande-Perez et al. 2002). However, evidence of hypermutated viral genomes has not been found in preextinction populations. In particular, hypermutated sequences of LCMV from 5-fluorouracil (FU) mutagenized populations were not amplified despite using a battery of oligonucleotide primers with different degrees of degeneracy targeting LCMV genomic RNA (Grande-Perez et al. 2005a). In this study, only one sequence was found that harbored nine mutations, six U → C transitions and the other three A → G, C → A, and G → A in an amplicon of 558 nucleotides (excluding primers), which renders a mutation frequency (relative to the consensus sequence of the sample from which it was obtained) of 1.6 × 10−2 mutations per nucleotide sequenced. No further analyses were conducted to find out whether this hypermutated molecule was the result of chemical mutagenesis or host adenosine deaminase. However, in a subsequent study (Martin et al. 2010), the participation of ADAR-like enzymes in the 5-FU-mediated increase of A → G transitions of LCMV genome was ruled out because the context in the immediate neighborhood of the observed A → G transitions was not preferred by the dsRNA adenosine deaminase (Polson and Bass 1994).

Contribution of Genetic Variability to Arenavirus–Host Interactions and Associated Disease

Coevolution of Arenaviruses and Their Natural Reservoirs: Intra- versus Inter-host Genetic Variation of Arenaviruses

Arenaviruses establish long-life persistent infections in their natural reservoirs in the absence of noticeable clinical symptoms (Childs and Peters 1993). Arenavirus natural reservoirs include rodent species from the family Muridae (Salazar-Bravo et al. 2002), but a notable exception is Tacaribe virus whose natural host is thought to be the fruit-eating bat Artibeus sp. (Downs et al. 1963). Coevolution of the arenaviruses with their natural hosts has likely been a major driving evolutionary force that has generated the currently existing arenavirus–host associations (Bowen et al. 1997), with groups of genetically closely related arenaviruses infecting closely related rodent species. Thus, OW arenaviruses are found in rodents from the subfamily Murinae, whereas NW arenavirus reservoirs belong to rodents from subfamilies Sigmodontinae and Neotominae (Fig. 2). Furthermore, inside the OW arenaviruses, the LASV complex (LASV, MOBV, MOPV, and IPPYV) is associated with Praomys rodents, while LCMV and Kodoko virus are found in rodents from the genus Mus (Lecompte et al. 2007). Likewise, within NW arenaviruses, the reservoirs of most North American arenaviruses are mainly rodents from the subfamily Neotominae, whereas the reservoirs of the South American arenaviruses belong to the subfamily Sigmodontinae (Emonet et al. 2009a). However, this coevolution hypothesis, which assumes also occasional host switching, could not be tested empirically due to the lack of detailed knowledge about rodent phylogeny. Studies trying to reconciliate the viral and host phylogenies for both OW (Hugot et al. 2001; Coulibaly-N’Golo et al. 2011) and NW (Jackson and Charleston 2004) arenaviruses did not find evidence of congruence between both phylogenies. Recent evidence has suggested that for OW arenaviruses, coevolutionary associations may have been concealed by the occurrence of multiple host switching and extinction events (Coulibaly-N’Golo et al. 2011). In addition, studies based on the recent availability of molecular data from both NW arenaviruses and their rodent reservoirs found evidence of arenavirus evolution via host switching, rather than codivergence between these arenaviruses and their hosts (Irwin et al. 2012).

Studies on the interaction between the surface glycoprotein GP1 of several NW arenaviruses and their cell entry receptor, the transferring receptor 1 (TfR1) (Abraham et al. 2009; Choe et al. 2011; Demogines et al. 2013; Flanagan et al. 2008; Radoshitzky et al. 2007, 2008), have suggested that arenaviruses might have exerted some degree of positive selection in their rodent reservoirs. These studies suggest that the dynamic of GP1–TfR1 interaction has had a direct impact in sequences of these viral and host genes. Interestingly, identified rapidly evolving positions within TfR1 were located within the TfR1 region that mediates interaction with the virus GP1. Mutations at these TfR1 residues altered the interaction with the viral GP1 without affecting TfR1 expression or function. As a consequence of the introduction of genetic divergence, this rodent receptor locus might have the potential to act as a species barrier to viral transmission limiting both the cross-species and zoonotic transmission of these viruses. Nevertheless, there is evidence that JUNV can enter cells very efficiently in a TfR1-independent manner, raising some questions about the true impact of the TfR1 gene as a key determinant of JUNV host range. NW arenaviruses, Junin virus, Machupo virus, and Guanarito virus, have acquired the ability to bind human TfR1 and are currently emerging in South America as zoonosis (Choe et al. 2011; Flanagan et al. 2008) associated with severe morbidity (Peters 2002). Intriguingly, a mutation within TfR1 predicted to affect infection by NW arenaviruses that use this receptor for cell entry was identified in Asian population (Demogines et al. 2013).

Genomic analysis of full-length genome sequences of 29 strains of LCMV collected from a variety of geographic regions at different times showed these viruses to be genetically highly diverse (Albarino et al. 2010). Several distinct lineages of LCMV could be identified, but there was little correlation with time or place of isolation. Bayesian analysis enabled estimation of the rate of evolution of the S and L genome segments of LCMV and Kodoko virus sequences and found them to be similar to other negative-stranded RNA viruses. The molecular evolutionary rate was estimated to be 3.3 × 10−4 and 3.7 × 10−4 substitutions per nucleotide and year for the S and L segments, respectively. The Bayesian analysis also indicated that the most recent common ancestor for LCMV S- and L genome RNA segments was estimated to be 1235 and 5142 years ago, respectively. These data indicate that LCMV is quite ancient, and the extensive diversity of the virus has accumulated over the past 1000–5000 years. Also, despite similar rates of evolution, the evolutionary history of the L segment appears to be more complex and can be traced back substantially longer than that for the S segment.

Arenaviruses coevolve with their natural hosts, but animal species other than the corresponding reservoir can be infected (Dalldorf 1939; Greenwood and Sanchez 2002; Salazar-Bravo et al. 2002). Information about possible arenavirus pathogenicity in their natural reservoir is rather limited, but data from LCMV infection of the mouse have shown that the same virus can behave differently depending on the host. For instance, C3H/st, BALB/WEHI, and SWR/J mice infected at birth with LCMV Armstrong, E-350, or Pasteur strains develop persistent infection, but only C3H/st mice develop growth hormone deficiency syndrome. In contrast, LCMV strains Traub and WE did not cause disease in those animals (Oldstone et al. 1984, 1985). Differential host-specific selection may be the reason why quasispecies diversity of Guanarito virus isolates from rodents in Venezuela had higher sequence variation than human isolates. One rodent isolate included a mixture of two phylogenetically distinct genotypes, suggesting a dual infection (Weaver et al. 2000). Likewise, the characterization by pyro-sequencing of the genetic stability and complexity of the Lassa vaccine candidate ML29 in vaccinated rhesus monkeys, marmosets, and mice (Zapata and Salvato 2013) revealed accumulation over time of some host-specific mutations. Protein prediction analysis showed that those mutations, at the amino acid level, induced potential structural changes in GPC and NP.

Contribution of Viral Quasispecies to Arenavirus Pathogenesis

RNA viruses within hosts exist as complex distribution of genetically closely related mutants, termed quasispecies, subjected to a process of genetic variation, competition, and selection (Domingo 1992; Domingo et al. 2001). Foremost, this swarm of mutants is possible because of the low copy-fidelity of viral RdRp. Studies in cultured cells found mutation frequencies of 1 × 10−4 to 5.7 × 10−4 substitutions per nucleotide for the Z, NP, GP, and L regions analyzed of LCMV (Grande-Perez et al. 2002, 2005a; Martin et al. 2010; Moreno et al. 2011), which are consistent with previous findings on genetic heterogeneity of LCMV (Sevilla et al. 2002) and other arenaviruses (Bowen et al. 1996, 1997; Bowen et al. 2000; Charrel and de Lamballerie 2003; Fulhorst et al. 2001; Garcia et al. 2000; Weaver et al. 2000, 2001) and also within the range commonly observed for riboviruses (10−3 to 10−5 substitutions per nucleotide) (Domingo et al. 1988; Drake and Holland 1999) (see Sect. 3.2.1).

Error-prone replication might contribute to spatial and temporal heterogeneities in RNA genome populations, favoring their rapid evolution and adaptation to different environments. The entire virus cloud contributes to the characteristics of the virus strain and will be the target of selection instead of individual variants (Perales et al. 2005). Cooperation and complementation events between variants of the quasispecies that permit the mutants carrying lethal mutations can be maintained within the population and even compete and interfere with fitter variants (de la Torre and Holland 1990; Gonzalez-Lopez et al. 2004; Martin et al. 2010; Mas et al. 2010; Perales et al. 2007; Zapata and Salvato 2013).

In addition to an error-prone polymerase, RNA recombination (Charrel et al. 2001; Fulhorst et al. 1996) and genome reassortment (Riviere 1987; Riviere et al. 1986; Riviere and Oldstone 1986) can also contribute to arenavirus genetic diversity. However, both recombination and reassortment events appear to be extremely rare in nature and only between phylogenetically closely related strains (Archer and Rico-Hesse 2002; Charrel et al. 2002; Emonet et al. 2006; Jay et al. 2005). The ability of the arenavirus to infect a wide variety of hosts could also contribute to their genetic variability (Zapata and Salvato 2013). The genetic diversity within and between arenavirus groups suggests that the spatial heterogeneity may be reflected in host range and pathogenicity (Blasdell et al. 2008). Biological functions of RNA virus distributions remain poorly understood, but evidence suggests that different arenavirus strains may exhibit different biological properties and pathogenesis attributable to differences in their quasispecies (Sevilla et al. 2002). Moreover, the genetic complexity of a quasispecies structure provides RNA viruses, including arenaviruses, with a great capacity to adapt to changing host environments (Novella et al. 1995). The quasispecies plasticity may facilitate that arenaviruses can overcome barriers to its spread, both within the host and between species (Oldstone and Campbell 2010; Sevilla et al. 2002). Viral variants, initially present at low frequency in the quasispecies, that become dominant under a selective pressure have been documented for different RNA viruses including rabies (Morimoto et al. 1998), coxsackie viruses (Beaucourt et al. 2011; Domingo et al. 2008), poliovirus (PV) (Pfeiffer and Kirkegaard 2005; Vignuzzi et al. 2006), foot-and-mouth disease virus (FMDV) (Haydon et al. 2001; Martin and Domingo 2008), and human parainfluenza virus (Prince et al. 2001). Arenaviruses are no exception, and selection of organ-specific LCMV variants with different pathogenic properties has been documented (Ahmed and Oldstone 1988; Ahmed et al. 1984). Likewise, as with other RNA viruses (Hotchin 1973; Hotchin et al. 1971; Hotchin and Sikora 1973), serial passages in cultured cells can yield LCMV variants that differ in pathogenicity with respect to the original isolate (Hotchin and Sikora 1973; Lukashevich 1992; Pulkkinen and Pfau 1970), further emphasizing the contribution of the quasispecies’ plasticity to arenavirus pathogenesis. These studies in cell culture provide strong support for the view that specific environmental conditions can influence virus–host interactions and promote selection of specific RNA species with distinct phenotypic features within the genetically complex quasispecies. The application of New Generation Sequencing to the field of virology is allowing scientists to start getting a detailed view of the genetic complexity and dynamic of RNA viral quasispecies of how they impact virus–host interactions and associated disease.

Selection of Immunosuppressive Variants During LCMV Persistence

Chronic infections provide investigators with an excellent experimental system for studying the contribution of virus quasispecies dynamics to virus–host interactions and associated disease as sustained virus multiplication in the persistently infected host facilitates virus evolution and adaptation resulting in selection of novel viral phenotypes that could influence the outcome of infection. Persistent infections can favor selection of viral variants with tissue-specific growth advantages that could be associated with different phenotypes and disease manifestations (Ahmed et al. 1984). Studies examining viral populations replicating in different tissues and organs of the same infected host have uncovered intra-host viral variability and selection of tissue-specific virus populations due to differences in selective pressures exerted by the different environments associated with these compartmentalized areas of the infected host (Ahmed and Oldstone 1988; Brown et al. 2011; Cabot et al. 2000; Cheng-Mayer et al. 1989; Deforges et al. 2004; Hall et al. 2001; Jelcic et al. 2004; Jridi et al. 2006; Sanjuan et al. 2004; Trivedi et al. 1994; Wright et al. 2011). LCMV can cause acute or persistent infections in mice, its natural host, and has been used as a model to study viral population dynamics in vivo (Ahmed et al. 1988; Evans et al. 1994). The age, genetics, and immune status of the host, as well as the composition of the viral quasispecies, route of infection, and dose of infecting virus, contribute to the outcome of the infection. Mice infected with LCMV within 24 h of birth do not develop an effective antiviral CTL response and become persistently infected (Matloubian et al. 1990; Salvato et al. 1991). On the other hand, iv inoculation of adult immunocompetent mice with a high (≥2 × 106 PFU) dose of Armstrong (ARM) strain of LCMV causes an acute infection that induces a protective immune response, which culminates in CD8+ CTL-mediated virus clearance. In contrast, a similar infection with the immunosuppressive clone 13 (Cl-13) strain of LCMV causes a persistent infection associated with generalized immune suppression (Ahmed et al. 1984, 1991; Ahmed and Oldstone 1988; Baranowski et al. 2000; Evans et al. 1994; Tishon et al. 1993; Wu-Hsieh et al. 1988). At early times, post-infection ARM exhibits tropism for macrophages of the red pulp of the spleen, while Cl-13 targets preferentially dendritic cells (DCs) of the spleen’s white pulp (Borrow et al. 1995; Sevilla et al. 2000). Notably, ARM and Cl-13 differ only at five nucleotide positions that involve three amino acid changes, two within GP1 [N(ARM)-176-D(Cl-13) and F(Arm)-260-L(Cl-13)] and one in the L polymerase [K(Arm)-1079-Q(C13)] (Matloubian et al. 1990, 1993; Salvato and Shimomaye 1989; Sevilla et al. 2000). Mutation F260L was found to be responsible for the strong binding affinity of Cl-13 to its receptor αDG that is expressed in DCs and facilitates over time high levels of DC infection by Cl-13 (Cao et al. 1998; Kunz et al. 2002), which is associated with Cl-13 persistence. However, at early stages of infection, ARM and Cl-13 infect similar numbers of DCs and macrophages indicating that differences between Cl-13 and ARM regarding early and yet-unknown interactions with host cell factors may determine whether LCMV persists (Cl-13) or is cleared ARM. Interestingly, mutation K1076Q in the L polymerase was found to play a critical role in the persisting phenotype of Cl-13 (Bergthaler et al. 2010; Sullivan et al. 2011).

Contribution of Viral Variants to LCMV-Induced Growth Hormone Deficiency Syndrome

Neonatal infection of C3He mice with certain strains of LCMV results in the development of a growth hormone disease syndrome (GHDS) characterized by growth retardation and severe hypoglycemia that result in death of the infected mice within 15–20 days of infection (Oldstone 2002; Oldstone et al. 1984, 1985). LCMV-induced GHSD was associated with decreasing levels of growth hormone (GH) mRNA that results in decreased GH protein synthesis in the anterior pituitary (Klavinskis and Oldstone 1989; Valsamakis et al. 1987). Virus-induced GHDS strictly correlated with the virus’s ability to replicate at high levels in the GH-producing cells in the anterior pituitary. Studies using reassortant viruses between strains of LCMV, which do (ARM) or do not (WE) cause GHSD, mapped the ability to cause this disorder to the S RNA (Riviere 1987). Notably, clonal analysis of the parental GHDS-nil WE population showed, as expected, that the majority (58/61) of the clones examined behaved as the parental WE clonal population and did not cause GHDS, but three clones isolated from the WE parental population did cause the characteristic GHDS (Buesa-Gomez et al. 1996). These results provided evidence that variants with the ability to cause GHDS (GHDS+) may be present within a GHDS-nil replicating WE parental population. Clonal isolates from the same parental virus population are genetically closely related. Accordingly, characterization of WEc54 (GHDS+) and WEc2.5 (GHDS−) clonal populations derived from the WE parental population identified a single amino acid at position 153 in GP1 as a critical viral genetic determinant required for induction of GHDS, illustrating how minor changes in the virus genome may have great impact in the biology of the infected host (Buesa-Gomez et al. 1996; Teng et al. 1996).

Implications of Arenavirus Genetic Variability for the Development of Antiviral Drugs and Vaccines to Combat Human Pathogenic Arenaviruses

Lethal Mutagenesis as a Novel Antiviral Strategy to Combat Arenavirus Infections

RNA viruses owe their potential for rapid evolution and adaptation to their genetic structure of quasispecies that is determined mainly by their low-fidelity replication machineries. Quasispecies is defined as the spectrum of closely related but nonidentical genomes (see Chaps. 1, 3, and 5). Within replicating viral quasispecies, many individual species may have the potential to exhibit distinct phenotypes including drug-resistant mutants that pose a main obstacle to the success of antiviral therapy. The high polymerase mutation rates responsible for the genetic heterogeneity displayed by RNA viruses are close to a limit, know as error threshold, beyond which the quasispecies theory predicts that biologically meaningful genetic information cannot be maintained. In the case of RNA viruses, violation of this limit is predicted to lead viral extinction, a process termed entry into error catastrophe (Eigen 2002). This model would predict that treatment of ribovirus-infected cells with mutagenic agents that increase the virus polymerase mutation rate could exert an antiviral effect by promoting viral entry into error catastrophe. Initial experimental evidence in support of this was obtained by Holland and colleagues when investigating the effect of mutagens on VSV and PV multiplication (Holland et al. 1990; Lee et al. 1997), and subsequently, the term “ lethal mutagenesis” was coined by Loeb and colleagues to describe the antiviral effect of the mutagen 5-hydrodeoxycytidine in HIV-1 infections (Loeb et al. 1999).

Because of the quasispecies genetic structure of RNA viruses, the interaction, either positive or negative, among all members of the quasispecies strongly influences the behavior of the viral population. Thereby, it is the quasispecies, rather than single independent genome species, that is subjected to selective pressures and evolution. Notably, low dose of a mutagenic agent can increase the number of mutations per genome (m) without reaching the error threshold (mc), which could result in the generation of defective-interfering genomes (DIG) that can interfere virus multiplication and production of infectious progeny, a process that has been termed “ lethal defection.” This model is supported by experimental results and in silico predictions (Ahmed et al. 1984, 1991; Ahmed and Oldstone 1988; Baranowski et al. 2000; Evans et al. 1994; Gonzalez-Lopez et al. 2004; Grande-Perez et al. 2005b; Manrubia et al. 2010; Moreno et al. 2012a; Perales et al. 2007; Tishon et al. 1993; Wu-Hsieh et al. 1988). Lethal defection can explain why viral populations subjected to lethal mutagenesis can be extinguished in the absence of the loss of the master sequence, the sequence present at the highest frequency within the quasispecies, and randomization of the virus genome sequence both predicted by the quasispecies theory (Grande-Perez et al. 2005a). Disappearance of the master sequence and melting of the genetic information of the population would be observed only in situations where m increases beyond m c (Fig. 3), a situation that is highly unlikely to be observed in viral populations as the inhibitory effects exerted by DIG together with the loss of essential viral functions due to single mutations within the genome can result in viral extinction without changes in the consensus sequence of the viral population.

Fig. 3

Lethal mutagenesis. The lethal mutagenesis concept predicts that increasing the dose of a mutagenic agent will result in an increasing number of mutations per genome (m) in the virus population. When m crosses an error threshold (m > m c), the virus population will enter into error catastrophe that is associated with the melting of the genetic information and viral extinction. The lethal defection model considers that under a moderate mutagenic dose, compatible with maintenance of the biological information, defective-interfering genomes (DIG) are generated and they interfere with the multiplication of the virus population, which could result in virus extinction in the absence of the sequence signatures of error catastrophe

Experimental evidences of lethal mutagenesis in cell culture have been achieved in several viral systems such as HIV-1 (Dapp et al. 2009, 2012; Harris et al. 2005), PV (Graci et al. 2007), or HCV (Ortega-Prieto et al. 2013). In coronavirus, the Nsp14 contains an 3′–5′ exonuclease domain that is responsible of the first described proofreading activity in riboviruses (Denison et al. 2011). Accordingly, viral mutants lacking this activity exhibit increased sensitivity to lethal mutagenesis (Smith et al. 2013). In arenaviruses, loss of infectivity and virus extinction through lethal mutagenesis have been achieved under laboratory conditions (Grande-Perez et al. 2005b; Martin et al. 2010; Moreno et al. 2011, 2012b; Ruiz-Jarabo et al. 2003) supporting the feasibility of entry into error catastrophe as a new antiviral strategy to combat arenavirus infections, a view that has been further reinforced by results from in vivo experiments (Ruiz-Jarabo et al. 2003; Sanz-Ramos et al. 2012).

Mutagenesis of LCMV with 100 µg/ml 5-FU achieved dose-dependent decreases in infectivity and systematic extinction after 48 h in one or two passages when BHK-21 cell cultures were inoculated with a moi of 0.01. However, nucleoside analogue 5-azacytidine (AZC) at 5 or 10 µg/ml was not able to eliminate LCMV after 13 passages in BHK-21 cells, not even when viral load was decreased tenfold (Grande-Perez et al. 2002). Increments in mutation frequencies in the L region reached 16.8-fold, whereas in the GP region, it was just 2.6-fold and 4.1-fold in the NP region in populations treated with a low dose (20 µg/ml) of 5-FU. However, in a preextinction population treated with 100 µg/ml of 5-FU, the mutation frequency of the polymerase increased 3.4-fold compared with the 1.6-fold found for the GP and NP regions (Grande-Perez et al. 2002). Notably, the largest increments in mutation frequencies did not correlate linearly with the lowest infectivity values (Grande-Perez et al. 2002). The type of mutations found differed between the two analogues being predominant the transitions U → C and A → G, followed by C → U in 5-FU-treated populations, whereas in AZC-treated populations, transversions C → G were the most frequent followed by C → A, G → C, and G → U. Despite the high diversity of LCMV mutagenized preextinction populations, their consensus genomic nucleotide sequences remained unaltered (Grande-Perez et al. 2005a). Furthermore, alternate virus passages in the presence and absence of 5-FU led to profound differences in fitness (measured as the capacity of LCMV to produce progeny), without modification of the consensus genomic sequence. LCMV viral populations in their way to extinction were enriched with mutants harboring a limited number of mutations called defectors (Grande-Perez et al. 2005a, b) (Fig. 3). Detailed analysis of the molecular scenario of LCMV persistent infections treated with 100 µg/ml of 5-FU showed an infectivity decrease that was paralleled by a significant increase of intracellular viral genomic L RNA. Preextinction populations recovered at different times post-infection showed specific infectivities (measured as plaque-forming units per genomic RNA molecule) lower than control populations. These experimental results together with silico simulations led to the lethal defection model of lethal mutagenesis that suggests the participation of defective-interfering genomes to explain virus extinction in cell culture when low-level mutagenesis is exerted by a mutagenic nucleoside analogue (Grande-Perez et al. 2005b). Interference caused by defectors is different to the LCMV natural defective-interfering particles, and it is produced as a consequence of 5-FU treatment (Martin et al. 2010). Supernatants of 5-FU-treated monolayers displayed interference capacity specific for LCMV, sensitive to UV irradiation and dose- and time-dependent (Martin et al. 2010). Interference did not involve significant increases of mutant spectrum complexity, as predicted by the lethal defection model (see also Chaps. 13 and 14).

Moreover, LCMV populations subjected to lethal mutagenesis exhibited a negative correlation between viral load and extinction, and infection at low moi in the presence of 5-FU was accompanied by a higher increase of mutant spectrum complexity and a higher antiviral effect than infection at high moi (Moreno et al. 2012b). This behavior was shared by VSV, another negative-stranded RNA virus, whereas the opposite was observed for the positive-stranded RNA viruses FMDV and EMCV. Increased mutagenesis caused by 5-FU treatment promoted the generation of DIG that inhibited viral RNA synthesis, which resulted in increased protection of the virus population against the mutagenic effect of 5-FU. This finding illustrates the biological consequences associated with the generation and activity of DIG. Nucleoside analogues are actively being pursued as candidate antiviral drugs that could exhibit a broad-spectrum mutagenic activity against a variety of riboviruses (Bonnac et al. 2013). The discovery of new nucleotide analogues that increase the viral mutation frequency while minimizing cell toxicity will represent a major advance for our understanding of lethal mutagenesis and its application as a novel antiviral strategy (Graci et al. 2008; Harki et al. 2002; Harris et al. 2005; Hicks et al. 2013).

Advantages of Combination Versus Single Drug Therapy to Inhibit Arenavirus Multiplication

Current antiviral drug therapy against arenavirus infections is limited to an off-label use of the nucleoside analogue ribavirin that has a limited prophylactic efficacy. In vitro and in vivo studies have documented the prophylactic and therapeutic value of RIb against several arenaviruses (Damonte and Coto 2002). Importantly, Rib reduced both morbidity and mortality in humans associated with LASV infection and experimentally in Machupo and Junin infections, if given early in the course of clinical disease (Damonte and Coto 2002). Despite its validated clinical benefits, the use of Rib has the limitations of its frequent and significant side effects, including anemia and congenital disorders, and the need of an early and intravenous administration for optimal efficacy. Several inhibitors of IMPDH (Andrei and De Clercq 1993), the S-adenosylhomocysteine (SAH) hydrolase and a variety of sulfated polysaccharides (Andrei and De Clercq 1990), phenotiazines compounds (Candurra et al. 1996), brassinosteroids (Wachsman et al. 2000), myristic acid (Cordo et al. 1999), and valproic acid (Vazquez-Calvo et al. 2013) have been reported to have antiarenaviral activity. However, most of these compounds displayed only modest and rather nonspecific effects associated with significant toxicity. Therefore, there is a pressing need for novel effective antiarenaviral drugs. In this regard, a recent high-throughput screening (HTS) using a virus-induced cytopathic effect (CPE)-based assay identified a potent small molecule inhibitor of Tacaribe virus and several other NW arenaviruses (Bolken et al. 2005). Likewise, cell-based HTS identified several small molecule inhibitors of virus cell entry mediated by LASV GP (Lee et al. 2008). These findings illustrate how complex chemical libraries used in the context of appropriate screening assays can be harnessed into a powerful tool to identify candidate antiviral drugs with highly specific activities.

Rib has been documented to be a broad-spectrum antiviral drug active against a variety of RNA viruses (Sidwell et al. 1974; Streeter et al. 1973). Several mechanisms of action have been proposed for the antiviral activity of Rib including direct incorporation into the RNA during replication, alteration of dNTP pools through IMPDH inhibition, direct inhibition of RNA synthesis, RNA capping inhibition, and immunomodulation (Graci and Cameron 2006). In LCMV-infected cells, treatment with high (100–400 µm) concentration of Rib did not increase mutation frequency values but inhibited dramatically both viral RNA synthesis and production of infectious progeny (Ruiz-Jarabo et al. 2003). However, treatment with low (20 µm) concentration of Rib increased 2.3-fold the mutation frequency of LCMV (Moreno et al. 2011), a finding that revealed a dose-dependent dual inhibitory and mutagenic activity of ribavirin for LCMV. The antiviral activity of Rib in LCMV-infected cells could be either totally or partially prevented by adding guanosine to the culture medium (Moreno et al. 2011). Inhibition of IMPDH by mycophenolic acid, which decreases intracellular GTP levels, did not increase but rather decreased LCMV mutation frequencies, thus raising some questions about whether targeting of IMPDH by Rib contributes to Rib-mediated inhibition of arenavirus multiplication. An alternative explanation for the protective effect of guanosine would be the dilution of Rib into the intracellular nucleotide pools, suggesting that Rib may be a competitive inhibitor (Moreno et al. 2011). The use of 5-FU as mutagenic agent and Rib either as inhibitor or as mutagen revealed that a treatment protocol based on a sequential inhibitor–mutagen administration exhibited the most potent antiviral activity compared to other combinations (Moreno et al. 2012a).

Rib-resistant escape mutants have been reported for several RNA viruses including PV, HCV, and FMDV (Pfeiffer and Kirkegaard 2003, 2005; Scheidel et al. 1987; Young et al. 2003) but are in general difficult to emerge. This likely reflects that Rib exerts its antiviral effect through different, mutually non-exclusive, mechanisms and thereby, a mutation conferring a virus with resistance to one mechanism may not allow viral escape from the overall antiviral activity of Rib. The observation that antiviral drugs exert their actions as a result of several mechanisms of action is not uncommon. Thus, 5-FU exerts its antiviral activity against FMDV through both its mutagenic activity and as inhibitor of initiation of viral RNA replication (Agudo et al. 2008). In coronavirus, the lack of the 3′–5′ exonuclease activity responsible of Nsp14 protein proofreading activity increases virus sensitivity to Rib and 5-FU mutagenesis, but evidence indicates that additional mechanisms of action associated with Rib treatment contribute to the overall observed reduction in virus infectivity (Smith et al. 2013). In the case of arenaviruses, Rib has been shown to exhibit a dose-dependent dual inhibitory and mutagenic activity (Moreno et al. 2011). At low concentrations (20 µm), Rib has been reported to increase the mutation frequency and mutant spectrum complexity of the population, whereas at higher concentrations (100 µm), Rib inhibited arenavirus replication without affecting the virus quasispecies complexity. The lack of a mutagenic effect of Rib when used at high concentration is not unexpected, because the corresponding strong inhibition of virus replication would limit the incorporation of Rib into the template RNA, which is required to exert its mutagenic effect.

The use of monotherapy-based treatments usually facilitates the selection of drug escape-resistant mutant within the population and subsequent virus rebound and therapy failure (Domingo 1989; Nijhuis et al. 2009; Young et al. 2003). Thus, inhibitor escape mutants have been isolated from virtually any virus for which a specific inhibitor has been developed, which poses a general problem in antimicrobial therapy (Domingo 2006). Combination therapy has proven effective to minimize and delay the selection of escape mutants (Pariente et al. 2003; Tapia et al. 2005). However, the emergence of variants resistant to more than one antiviral drug is not uncommon. A novel concept of combination therapy that is based on combining antiviral drugs targeting specific steps of the virus life cycle and lethal mutagenesis could contribute to overcome this problem. Importantly, recent evidence indicates that combination therapy of inhibitors of virus multiplication together with lethal mutagenesis can be dramatically more effective, while exhibiting much less toxicity, in controlling and clearing RNA virus infections than classic combination therapy (Perales et al. 2012). Moreover, evidence indicates that mutagen-mediated increase in virus mutation rate promotes the accumulation of non-functional viral gene products that often exhibit a strong dominant negative phenotype that further contributes to inhibition of virus multiplication, which has been referred to as lethal defection model of virus extinction (Grande-Perez et al. 2005b). Therefore, modest mutagenic intensities may be sufficient to help driving a virus to extinction, which further supports a combination antiviral therapy that incorporates lethal mutagenesis.

The interactions (neutral, positive, or negative) between drugs used in combination therapy are determined by the mechanism of action and pharmacokinetics of each drug, and thereby, the implementation of combination therapy requires a detailed understanding of the mechanisms of action of each drug. It is important noting that the mutagenic effect exerted by a drug could be masked or reduced by a replication inhibitor when used in combination, whereas the sequential administration of both drugs can result in a stronger antiviral effect than the sum of both drug effects (Moreno et al. 2012a; Perales et al. 2009). Therefore, combination therapy involving the use of a mutagenic agent should employ protocols involving the sequential administration of the drugs. In this regard, theoretical predictions (Eigen 1971, 2002) support the idea that reduced viral load promotes mutagen-induced viral extinction. Studies with FMDV provided the first experimental evidence for this prediction (Perales et al. 2009; Sierra et al. 2000). Likewise, studies with LCMV showed that administration of Rib at concentrations that favored its replication inhibition activity followed (sequential treatment) by administration of the mutagen 5-FU resulted in a stronger antiviral effect than when drugs were administered simultaneously or individually (Moreno et al. 2012a). This antiviral effect correlated with a low viral progeny production, low viral RNA synthesis, and a noticeable increase in quasispecies mutant spectrum complexity. The inhibition of viral replication by Rib potentiated the mutagenic effect of the subsequently added 5-FU, promoting the generation of DIG that contributed to the loss of infectivity. None of these effects were observed when mutagenic doses of Rib were used either in combination or sequentially with 5-FU. The lesser antiviral effect observed when inhibitory concentrations of Rib were administered together with the mutagen 5-FU could be explained due to a negative interaction between mutagen and replication inhibitor (Moreno et al. 2012a). The implications of these findings for the design of protocols to treat arenavirus infections in vivo remain to be investigated.

Arenavirus Genetic Variation and Development of Broad-Spectrum Vaccines

There are not FDA-approved vaccines against HF arenaviral diseases. The JUNV live-attenuated Candid1 strain has been shown to be an effective vaccine against Argentine hemorrhagic fever (AHF) disease (Enria and Barrera Oro 2002; Enria et al. 2008). However, outside Argentina, Candid1 has only IND status and studies addressing long-term immunity and safety have not been conducted. Moreover, Candid1 does not protect against LASV. Despite significant efforts dedicated to the development of LASV vaccines, not a single LASV vaccine candidate has entered a clinical trial although the MOPV/LASV reassortant ML29, as well as recombinant VSV and vaccinia virus expressing specific LASV antigens, has shown promising results (Falzarano and Feldmann 2013). Specifically, ML29 exhibited good safety and efficacy profiles in animal models of LASV infection. However, the high prevalence of HIV within LASV endemic regions raises safety concerns about the use of VSV- or vaccinia-based platforms.

Control of LASV infection seems to be mediated mainly by cellular immune responses, and significant titers of LASV neutralizing antibodies (NAbs) are usually observed only in patients who have clinically recovered (Jahrling and Peters 1984). However, passive antibody transfer has been shown to induce protection in animal models of LF (Jahrling 1983) and in limited human studies (Monath and Casals 1975) suggesting that an ideal vaccine should induce the right combination of cellular and humoral responses.

Sequence comparison of all proteins of seven pathogenic arenaviruses revealed a high degree of genetic variation among different geographical and temporal isolates of the same virus species (Bowen et al. 2000; Fulhorst et al. 2001; Garcia et al. 2000; Weaver et al. 2000). Overall, strain variation among 54 LASV strains was found to be as high as 27 and 15 % at the nucleotide and amino acid levels, respectively (Bowen et al. 2000). The high degree of LASV genetic diversity likely contributes to underestimating its prevalence and poses also a great challenge for vaccine development. In some endemic areas, up to 32 % of anti-LASV-positive samples would not have been detected, if the IFA had only been run with LASV-Josiah (Jos) strain, which is historically used as antigen in IFA or ELISA (Emmerich et al. 2008). This underscores the need for a multivalent or “universal” LASV vaccine able to protect against infections by the related isolates of the four lineages currently identified for LASV. The rML29/Jos expresses NP and GP antigens of Josiah strain, the prototype of LASV lineage IV, the largest sublineage, containing strains isolated in Guinea, Liberia, and Sierra Leone (Bowen et al. 2000). A single dose of ML29 vaccine effectively protected guinea pigs against challenge with homologous (lineage IV) or heterologous (lineage II) strains (Carrion et al. 2012). Lineage III and IV had a sister relationship, and the protection elicited by ML29 against LASV strains from lineage II strongly suggests that ML29 will also protect against members of lineage III. In contrast, sequence differences between LP (I) and Jos (IV) are 23.6 and 20.8 % for NP and GPC genes, respectively, and only one of the three recently predicted HL2-A2-restricted CTL epitopes was conserved in GP sequences of both viruses (Botten et al. 2006). This potential problem could be overcome by the use of arenavirus reverse genetics to generate rML29 expressing NP and GP antigens from lineages I and IV to generate a blended vaccine able to provide cross-protection against members of all LASV lineages. Despite the very promising features of the rML29 platform for the development of pan-LASV vaccine, it should be noted that the mechanisms of ML29 attenuation remain poorly understood and additional mutations, including reversions, in ML29 or reassortants between ML29 and circulating virulent LASV strains, could result in viruses with enhanced virulence.