Quasispecies dynamics and fixation of a synonymous mutation in hantavirus transmission

Puumala Virus Variants Circulating in Forests of Ardennes, France: Ten Years of Genetic Evolution

In Europe, Puumala virus (PUUV) transmitted by the bank vole (Myodes glareolus) is the causative agent of nephropathia epidemica (NE), a mild form of haemorrhagic fever with renal syndrome. In France, very little is known about the spatial and temporal variability of the virus circulating within bank vole populations. The present study involved monitoring of bank vole population dynamics and PUUV microdiversity over a ten-year period (2000–2009) in two forests of the Ardennes region: Elan and Croix-Scaille. Ardennes region is characterised by different environmental conditions associated with different NE epidemiology. Bank vole density and population parameters were estimated using the capture/marking/recapture method, and blood samples were collected to monitor the overall seroprevalence of PUUV in rodent populations. Phylogenetic analyses of fifty-five sequences were performed to illustrate the genetic diversity of PUUV variants between forests. The pattern of the two forests differed clearly. In the Elan forest, the rodent survival was higher, and this limited turn-over resulted in a lower seroprevalence and diversity of PUUV sequences than in the Croix-Scaille forest. Uncovering the links between host dynamics and virus microevolution is improving our understanding of PUUV distribution in rodents and the NE risk.

How Bank Vole-PUUV Interactions Influence the Eco-Evolutionary Processes Driving Nephropathia Epidemica Epidemiology-An Experimental and Genomic Approach

In Europe, Puumala virus (PUUV) is responsible for nephropathia epidemica (NE), a mild form of hemorrhagic fever with renal syndrome (HFRS). Despite the presence of its reservoir, the bank vole, on most of French territory, the geographic distribution of NE cases is heterogeneous and NE endemic and non-endemic areas have been reported. In this study we analyzed whether bank vole-PUUV interactions could partly shape these epidemiological differences. We performed crossed-experimental infections using wild bank voles from French endemic (Ardennes) and non-endemic (Loiret) areas and two French PUUV strains isolated from these areas. The serological response and dynamics of PUUV infection were compared between the four cross-infection combinations. Due to logistical constraints, this study was based on a small number of animals. Based on this experimental design, we saw a stronger serological response and presence of PUUV in excretory organs (bladder) in bank voles infected with the PUUV endemic strain. Moreover, the within-host viral diversity in excretory organs seemed to be higher than in other non-excretory organs for the NE endemic cross-infection but not for the NE non-endemic cross-infection. Despite the small number of rodents included, our results showed that genetically different PUUV strains and in a lesser extent their interaction with sympatric bank voles, could affect virus replication and diversity. This could impact PUUV excretion/transmission between rodents and to humans and in turn at least partly shape NE epidemiology in France.

Detection and Genetic Characterization of Puumala Orthohantavirus S-Segment in Areas of France Non-Endemic for Nephropathia Epidemica

Puumala virus (PUUV) in Europe causes nephropathia epidemica (NE), a mild form of hemorrhagic fever with renal syndrome (HFRS). The incidence of NE is highly heterogeneous spatially, whereas the geographic distribution of the wild reservoir of PUUV, the bank vole, is essentially homogeneous. Our understanding of the processes driving this heterogeneity remains incomplete due to gaps in knowledge. Little is known about the current distribution and genetic variation of PUUV in the areas outside the well-identified zones of NE endemicity. We trapped bank voles in four forests in French regions in which NE is considered non-endemic, but sporadic NE cases have been reported recently. We tested bank voles for anti-PUUV IgG and characterized the S segment sequences of PUUV from seropositive animals. Phylogenetic analyses revealed specific amino-acid signatures and genetic differences between PUUV circulating in non-endemic and nearby NE-endemic areas. We also showed, in temporal surveys, that the amino-acid sequences of PUUV had undergone fewer recent changes in areas non-endemic for NE than in endemic areas. The evolutionary history of the current French PUUV clusters was investigated by phylogeographic approaches, and the results were considered in the context of the history of French forests. Our findings highlight the need to monitor the circulation and genetics of PUUV in a larger array of bank vole populations, to improve our understanding of the risk of NE.

The Needs for Developing Experiments on Reservoirs in Hantavirus Research: Accomplishments, Challenges and Promises for the Future

Due to their large geographic distribution and potential high mortality rates in human infections, hantaviruses constitute a worldwide threat to public health. As such, they have been the subject of a large array of clinical, virological and eco-evolutionary studies. Many experiments have been conducted in vitro or on animal models to identify the mechanisms leading to pathogenesis in humans and to develop treatments of hantavirus diseases. Experimental research has also been dedicated to the understanding of the relationship between hantaviruses and their reservoirs. However, these studies remain too scarce considering the diversity of hantavirus/reservoir pairs identified, and the wide range of issues that need to be addressed. In this review, we present a synthesis of the experimental studies that have been conducted on hantaviruses and their reservoirs. We aim at summarizing the knowledge gathered from this research, and to emphasize the gaps that need to be filled. Despite the many difficulties encountered to carry hantavirus experiments, we advocate for the need of such studies in the future, at the interface of evolutionary ecology and virology. They are critical to address emerging areas of research, including hantavirus evolution and the epidemiological consequences of individual variation in infection outcomes.

An in vivo system for directed experimental evolution of rabbit haemorrhagic disease virus

The calicivirus Rabbit haemorrhagic disease virus (RHDV) is widely used in Australia as a biocontrol agent to manage wild European rabbit (Oryctolagus cuniculus) populations. However, widespread herd immunity limits the effectiveness of the currently used strain, CAPM V-351. To overcome this, we developed an experimental platform for the selection and characterisation of novel RHDV strains. As RHDV does not replicate in cell culture, variant viruses were selected by serially passaging a highly virulent RHDV field isolate in immunologically naïve laboratory rabbits that were passively immunised 18–24 hours post-challenge with a neutralising monoclonal antibody. After seven passages, two amino acid substitutions in the P2 domain of the capsid protein became fixed within the virus population. Furthermore, a synonymous substitution within the coding sequence of the viral polymerase appeared and was also maintained in all subsequent passages. These findings demonstrate proof-of-concept that RHDV evolution can be experimentally manipulated to select for virus variants with altered phenotypes, in this case partial immune escape.

Spatiotemporal dynamics of Puumala hantavirus associated with its rodent host, Myodes glareolus

Many viruses significantly impact human and animal health. Understanding the population dynamics of these viruses and their hosts can provide important insights for epidemiology and virus evolution. Puumala virus (PUUV) is a European hantavirus that may cause regional outbreaks of hemorrhagic fever with renal syndrome in humans. Here, we analyzed the spatiotemporal dynamics of PUUV circulating in local populations of its rodent reservoir host, the bank vole (Myodes glareolus) during eight years. Phylogenetic and population genetic analyses of all three genome segments of PUUV showed strong geographical structuring at a very local scale. There was a high temporal turnover of virus strains in the local bank vole populations, but several virus strains persisted through multiple years. Phylodynamic analyses showed no significant changes in the local effective population sizes of PUUV, although vole numbers and virus prevalence fluctuated widely. Microsatellite data demonstrated also a temporally persisting subdivision between local vole populations, but these groups did not correspond to the subdivision in the virus strains. We conclude that restricted transmission between vole populations and genetic drift play important roles in shaping the genetic structure and temporal dynamics of PUUV in its natural host which has several implications for zoonotic risks of the human population.

Silent mutations at codons 65 and 66 in reverse transcriptase alleviate indel formation and restore fitness in subtype B HIV-1 containing D67N and K70R drug resistance mutations

Resistance to combined antiretroviral therapy (cART) in HIV-1-infected individuals is typically due to nonsynonymous mutations that change the protein sequence; however, the selection of synonymous or ‘silent’ mutations in the HIV-1 genome with cART has been reported. These silent K65K and K66K mutations in the HIV-1 reverse transcriptase (RT) occur in over 35% of drug-experienced individuals and are highly associated with the thymidine analog mutations D67N and K70R, which confer decreased susceptibility to most nucleoside and nucleotide RT inhibitors. However, the basis for selection of these silent mutations under selective drug pressure is unknown. Using Illumina next-generation sequencing, we demonstrate that the D67N/K70R substitutions in HIV-1 RT increase indel frequency by 100-fold at RT codons 65–67, consequently impairing viral fitness. Introduction of either K65K or K66K into HIV-1 containing D67N/K70R reversed the error-prone DNA synthesis at codons 65–67 in RT and improved viral replication fitness, but did not impact RT inhibitor drug susceptibility. These data provide new mechanistic insights into the role of silent mutations selected during antiretroviral therapy and have broader implications for the relevance of silent mutations in the evolution and fitness of RNA viruses.

Changes in diversification patterns and signatures of selection during the evolution of murinae-associated hantaviruses

In the last 50 years, hantaviruses have significantly affected public health worldwide, but the exact extent of the distribution of hantavirus diseases, species and lineages and the risk of their emergence into new geographic areas are still poorly known. In particular, the determinants of molecular evolution of hantaviruses circulating in different geographical areas or different host species are poorly documented. Yet, this understanding is essential for the establishment of more accurate scenarios of hantavirus emergence under different climatic and environmental constraints. In this study, we focused on Murinae-associated hantaviruses (mainly Seoul Dobrava and Hantaan virus) using sequences available in GenBank and conducted several complementary phylogenetic inferences. We sought for signatures of selection and changes in patterns and rates of diversification in order to characterize hantaviruses’ molecular evolution at different geographical scales (global and local). We then investigated whether these events were localized in particular geographic areas. Our phylogenetic analyses supported the assumption that RNA virus molecular variations were under strong evolutionary constraints and revealed changes in patterns of diversification during the evolutionary history of hantaviruses. These analyses provide new knowledge on the molecular evolution of hantaviruses at different scales of time and space.

Comparative genomic and phylogenetic analysis of the first Usutu virus isolate from a human patient presenting with neurological symptoms

Usutu virus (USUV) is a mosquito-borne flavivirus, belonging to the Japanese encephalitis antigenic complex, that circulates among mosquitoes and birds. We describe and analyze the complete genome sequence of the first USUV strain isolated from an immunocompromised patient with neuroinvasive disease. This USUV isolate showed an overall nucleotide identity of 99% and 96%, respectively, with the genomes of isolates from Europe and Africa. Comparison of the human USUV complete polyprotein sequence with bird-derived strains, showed two unique amino acid substitutions. In particular, one substitution (S595G) was situated in the DIII domain of the viral Envelope protein that is recognized by flavivirus neutralizing antibodies. An additional amino acid substitution (D3425E) was identified in the RNA-dependent RNA polymerase (RdRp) domain of the NS5 protein. This substitution is remarkable since E3425 is highly conserved among the other USUV isolates that were not associated with human infection. However, a similar substitution was observed in Japanese encephalitis and in West Nile viruses isolated from humans. Phylogenetic analysis of the human USUV strain revealed a close relationship with an Italian strain isolated in 2009. Analysis of synonymous nucleotide substitutions (SNSs) among the different USUV genomes showed a specific evolutionary divergence among different countries. In addition, 15 SNSs were identified as unique in the human isolate. We also identified four specific nucleotide substitutions in the 5' and 3' untranslated regions (UTRs) in the human isolate that were not present in the other USUV sequences. Our analyses provide the basis for further experimental studies aimed at defining the effective role of these mutations in the USUV genome, their potential role in the development of viral variants pathogenic for humans and their evolution and dispersal out of Africa.

Genome-wide patterns of intrahuman dengue virus diversity reveal associations with viral phylogenetic clade and interhost diversity

Analogous to observations in RNA viruses such as human immunodeficiency virus, genetic variation associated with intrahost dengue virus (DENV) populations has been postulated to influence viral fitness and disease pathogenesis. Previous attempts to investigate intrahost genetic variation in DENV characterized only a few viral genes or a limited number of full-length genomes. We developed a whole-genome amplification approach coupled with deep sequencing to capture intrahost diversity across the entire coding region of DENV-2. Using this approach, we sequenced DENV-2 genomes from the serum of 22 Nicaraguan individuals with secondary DENV infection and captured ∼75% of the DENV genome in each sample (range, 40 to 98%). We identified and quantified variants using a highly sensitive and specific method and determined that the extent of diversity was considerably lower than previous estimates. Significant differences in intrahost diversity were detected between genes and also between antigenically distinct domains of the Envelope gene. Interestingly, a strong association was discerned between the extent of intrahost diversity in a few genes and viral clade identity. Additionally, the abundance of viral variants within a host, as well as the impact of viral mutations on amino acid encoding and predicted protein function, determined whether intrahost variants were observed at the interhost level in circulating Nicaraguan DENV-2 populations, strongly suggestive of purifying selection across transmission events. Our data illustrate the value of high-coverage genome-wide analysis of intrahost diversity for high-resolution mapping of the relationship between intrahost diversity and clinical, epidemiological, and virological parameters of viral infection.

A global perspective on hantavirus ecology, epidemiology, and disease

Hantaviruses are enzootic viruses that maintain persistent infections in their rodent hosts without apparent disease symptoms. The spillover of these viruses to humans can lead to one of two serious illnesses, hantavirus pulmonary syndrome and hemorrhagic fever with renal syndrome. In recent years, there has been an improved understanding of the epidemiology, pathogenesis, and natural history of these viruses following an increase in the number of outbreaks in the Americas. In this review, current concepts regarding the ecology of and disease associated with these serious human pathogens are presented. Priorities for future research suggest an integration of the ecology and evolution of these and other host-virus ecosystems through modeling and hypothesis-driven research with the risk of emergence, host switching/spillover, and disease transmission to humans.

Associations between MHC genes and Puumala virus infection in Myodes glareolus are detected in wild populations, but not from experimental infection data

We analysed the influence of MHC class II Dqa and Drb genes on Puumala virus (PUUV) infection in bank voles (Myodes glareolus). We considered voles sampled in five European localities or derived from a previous experiment that showed variable infection success of PUUV. The genetic variation observed in the Dqa and Drb genes was assessed by using single-strand conformation polymorphism and pyrosequencing methods, respectively. Patterns were compared with those obtained from 13 microsatellites. We revealed significant genetic differentiation between PUUV-seronegative and -seropositive bank voles sampled in wild populations, at the Drb gene only. The absence of genetic differentiation observed at neutral microsatellites confirmed the important role of selective pressures in shaping these Drb patterns. Also, we found no significant associations between infection success and MHC alleles among laboratory-colonized bank voles, which is explained by a loss of genetic variability that occurred during the captivity of these voles.

Rapid accumulation of virulent rift valley Fever virus in mice from an attenuated virus carrying a single nucleotide substitution in the m RNA

Background

Rift Valley fever virus (RVFV), a member of the genus Phlebovirus within the family Bunyaviridae, is a negative-stranded RNA virus with a tripartite genome. RVFV is transmitted by mosquitoes and causes fever and severe hemorrhagic illness among humans, while in livestock it causes fever and high abortion rates.

Methodology/Principal Findings

Sequence analysis showed that a wild-type RVFV ZH501 preparation consisted of two major viral subpopulations, with a single nucleotide heterogeneity at nucleotide 847 of M segment (M847); one had a G residue at M847 encoding glycine in a major viral envelope Gn protein, while the other carried A residue encoding glutamic acid at the corresponding site. Two ZH501-derived viruses, rZH501-M847-G and rZH501-M847-A, carried identical genomic sequences, except that the former and the latter had G and A, respectively, at M847 were recovered by using a reverse genetics system. Intraperitoneal inoculation of rZH501-M847-A into mice caused a rapid and efficient viral accumulation in the sera, livers, spleens, kidneys and brains, and killed most of the mice within 8 days, whereas rZH501-M847-G caused low viremia titers, did not replicate as efficiently as did rZH501-M847-A in these organs, and had attenuated virulence to mice. Remarkably, as early as 2 days postinfection with rZH501-M847-G, the viruses carrying A at M847 emerged and became the major virus population thereafter, while replicating viruses retained the input A residue at M847 in rZH501-M847-A-infected mice.

Conclusions/Significance

These data demonstrated that the single nucleotide substitution in the Gn protein substantially affected the RVFV mouse virulence and that a virus population carrying the virulent viral genotype quickly emerged and became the major viral population within a few days in mice that were inoculated with the attenuated virus.

Hantavirus evolution in relation to its rodent and insectivore hosts: no evidence for codivergence

Hantaviruses are considered one of the best examples of a long-term association between RNA viruses and their hosts. Based on the appearance of strong host specificity, it has been suggested that hantaviruses cospeciated with the rodents and insectivores they infect since these mammals last shared a common ancestor, approximately 100 million years ago. We tested this hypothesis of host-virus codivergence in two ways: 1) we used cophylogenetic reconciliation analysis to assess the fit of the virus tree onto that of the host and 2) we estimated the evolutionary rates and divergence times for the Hantavirus genus using a Bayesian Markov Chain Monte Carlo method and similarly compared these with those of their hosts. Our reconciliation analysis provided no evidence for a history of codivergence between hantaviruses and their hosts. Further, the divergence times for the Hantavirus genus were many orders of magnitude too recent to correspond with the timescale of their hosts' speciation. We therefore propose that apparent similarities between the phylogenies of hantaviruses and their mammalian hosts are the result of a more recent history of preferential host switching and local adaptation. Based on the presence of clade-defining amino acids in all genomic segments, we propose that the patterns of amino acid replacement in these viruses are also compatible with a history of host-specific adaptation.

Accumulation of point mutations and reassortment of genomic RNA segments are involved in the microevolution of Puumala hantavirus in a bank vole (Myodes glareolus) population

The genetic diversity of Puumala hantavirus (PUUV) was studied in a local population of its natural host, the bank vole (Myodes glareolus). The trapping area (2.5 x 2.5 km) at Konnevesi, Central Finland, included 14 trapping sites, at least 500 m apart; altogether, 147 voles were captured during May and October 2005. Partial sequences of the S, M and L viral genome segments were recovered from 40 animals. Seven, 12 and 17 variants were detected for the S, M and L sequences, respectively; these represent new wild-type PUUV strains that belong to the Finnish genetic lineage. The genetic diversity of PUUV strains from Konnevesi was 0.2-4.9 % for the S segment, 0.2-4.8 % for the M segment and 0.2-9.7 % for the L segment. Most nucleotide substitutions were synonymous and most deduced amino acid substitutions were conservative, probably due to strong stabilizing selection operating at the protein level. Based on both sequence markers and phylogenetic clustering, the S, M and L sequences could be assigned to two groups, 'A' and 'B'. Notably, not all bank voles carried S, M and L sequences belonging to the same group, i.e. S(A)M(A)L(A) or S(B)M(B)L(B). A substantial proportion (8/40, 20 %) of the newly characterized PUUV strains possessed reassortant genomes such as S(B)M(A)L(A), S(A)M(B)L(B) or S(B)M(A)L(B). These results suggest that at least some of the PUUV reassortants are viable and can survive in the presence of their parental strains.