Identifying the patterns and drivers of Puumala hantavirus enzootic dynamics using reservoir sampling

Fine-scale spatial and social patterns of SARS-CoV-2 transmission from identical pathogen sequences

Pathogen genomics can provide insights into disease transmission patterns, but new methods are needed to handle modern large-scale pathogen genome datasets. Genetically proximal viruses indicate epidemiological linkage and are informative about transmission events. Here, we leverage pairs of identical sequences using 114,298 SARS-CoV-2 genomes collected via sentinel surveillance from March 2021 to December 2022 in Washington State, USA, with linked age and residence information to characterize fine-scale transmission. The location of pairs of identical sequences is highly consistent with expectations from mobility and social contact data. Outliers in the relationship between genetic and mobility data can be explained by SARS-CoV-2 transmission between postal codes with male prisons, consistent with transmission between prison facilities. Transmission patterns between age groups vary across spatial scales. Finally, we use the timing of sequence collection to understand the age groups driving transmission. This work improves our ability to characterize transmission from large pathogen genome datasets.

How fast are viruses spreading in the wild?

Genomic data collected from viral outbreaks can be exploited to reconstruct the dispersal history of viral lineages in a two-dimensional space using continuous phylogeographic inference. These spatially explicit reconstructions can subsequently be used to estimate dispersal metrics allowing to unveil the dispersal dynamics and evaluate the capacity to spread among hosts. Heterogeneous sampling intensity of genomic sequences can however impact the accuracy of dispersal insights gained through phylogeographic inference. In our study, we implement a simulation framework to evaluate the robustness of three dispersal metrics - a lineage dispersal velocity, a diffusion coefficient, and an isolation-by-distance signal metric - to the sampling effort. Our results reveal that both the diffusion coefficient and isolation-by-distance signal metrics appear to be robust to the number of samples considered for the phylogeographic reconstruction. We then use these two dispersal metrics to compare the dispersal pattern and capacity of various viruses spreading in animal populations. Our comparative analysis reveals a broad range of isolation-by-distance patterns and diffusion coefficients mostly reflecting the dispersal capacity of the main infected host species but also, in some cases, the likely signature of rapid and/or long-distance dispersal events driven by human-mediated movements through animal trade. Overall, our study provides key recommendations for the lineage dispersal metrics to consider in future studies and illustrates their application to compare the spread of viruses in various settings.

Phylostems: a new graphical tool to investigate temporal signal of heterochronous sequences datasets

Motivation

Molecular tip-dating of phylogenetic trees is a growing discipline that uses DNA sequences sampled at different points in time to co-estimate the timing of evolutionary events with rates of molecular evolution. Importantly, such inferences should only be performed on datasets displaying sufficient temporal signal, a feature important to test prior to any tip-dating inference. For this purpose, the most popular method considered to-date has been the 'root-to-tip regression' which consist in fitting a linear regression of the number of substitutions accumulated from the root to the tips of a phylogenetic tree as a function of sampling times. The main limitation of the regression method, in its current implementation, relies in the fact that the temporal signal can only be tested at the whole-tree scale (i.e. its root).

Results

To overcome this limitation we introduce Phylostems, a new graphical user-friendly tool developed to investigate temporal signal within every clade of a phylogenetic tree. We provide a 'how to' guide by running Phylostems on an empirical dataset and supply guidance for results interpretation.

Availability and implementation

Phylostems is freely available at https://pvbmt-apps.cirad.fr/apps/phylostems.

Serum Markers Associated with Disease Severity in a Bosnian Hemorrhagic Fever with Renal Syndrome Cohort

Puumala orthohantavirus (PUUV) is endemic in Europe and can cause hemorrhagic fever with renal syndrome (nephropathia epidemica). Disease features include fever, thrombocytopenia, and acute kidney injury (AKI). This retrospective cohort study of forty PUUV patients aims to characterize associations of serum immunological, hemostatic or kidney injury markers to disease severity. While interleukin-18 (IL-18) was significantly increased in severely thrombocytopenic patients (<100 × 109 platelets/L) compared to patients with higher platelet counts, RANTES was significantly decreased in these patients. These data suggest that patients with significant thrombocytopenia might have experienced pronounced Th1 immune responses. When kidney dysfunction was used as the primary disease outcome, recently identified AKI biomarkers (Cystatin C, insulin-like growth factor-binding protein 7, Nephrin, and trefoil factor 3) were significantly upregulated in patients with severe PUUV infection, defined as the estimated glomerular filtration rate (eGFR) below 30 m/min/1.73 m2. The increased expression of these markers specifically indicates pathology in glomeruli and proximal tubuli. Furthermore, E-selectin was significantly higher while interferon gamma-induced protein 10 (IP-10) was significantly lower in PUUV patients with more severe kidney dysfunction compared to patients with higher eGFR-values. Increased E-selectin illustrates the central role of endothelial cell activation, whereas decreased IP-10 could indicate a less important role of this cytokine in the pathogenesis of kidney dysfunction.

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.

Impact of Landscape on Host-Parasite Genetic Diversity and Distribution Using the Puumala orthohantavirus-Bank Vole System

In nature, host specificity has a strong impact on the parasite's distribution, prevalence, and genetic diversity. The host's population dynamics is expected to shape the distribution of host-specific parasites. In turn, the parasite's genetic structure is predicted to mirror that of the host. Here, we study the tandem Puumala orthohantavirus (PUUV)-bank vole system. The genetic diversity of 310 bank voles and 33 PUUV isolates from 10 characterized localities of Northeast France was assessed. Our findings show that the genetic diversity of both PUUV and voles, was positively correlated with forest coverage and contiguity of habitats. While the genetic diversity of voles was weakly structured in space, that of PUUV was found to be strongly structured, suggesting that the dispersion of voles was not sufficient to ensure a broad PUUV dissemination. Genetic diversity of PUUV was mainly shaped by purifying selection. Genetic drift and extinction events were better reflected than local adaptation of PUUV. These contrasting patterns of microevolution have important consequences for the understanding of PUUV distribution and epidemiology.

Genetic Diversity of Puumala orthohantavirus in Rodents and Human Patients in Austria, 2012-2019

Puumala orthohantavirus (PUUV) has a wide distribution throughout Europe. Distinctive temporal patterns of spillover into the human population are related to population dynamics of the reservoir host, the bank vole (Clethrionomys glareolus). As the rodent host is tied to specific habitats with small individual ranges, PUUV genetic diversity is also highly correlated with geographic distance. Using sequenced portions of viral S and M segments, we determined whether geographic clusters were supported. Human cases of PUUV infections are concentrated in southeastern Austria. We detected four distinct genotypes: two genotypes of the Alpe-Adria (ALAD) lineage typically associated with southeast Europe, and two sublineages of the Central Europe (CE) lineage. One cluster of CE genotypes represents a phylogenetically distinct sublineage compared to previously reported CE clades, and extends the boundary of the CE lineage further south than previously reported.

Circuitscape in Julia: Empowering Dynamic Approaches to Connectivity Assessment

The conservation field is experiencing a rapid increase in the amount, variety, and quality of spatial data that can help us understand species movement and landscape connectivity patterns. As interest grows in more dynamic representations of movement potential, modelers are often limited by the capacity of their analytic tools to handle these datasets. Technology developments in software and high-performance computing are rapidly emerging in many fields, but uptake within conservation may lag, as our tools or our choice of computing language can constrain our ability to keep pace. We recently updated Circuitscape, a widely used connectivity analysis tool developed by Brad McRae and Viral Shah, by implementing it in Julia, a high-performance computing language. In this initial re-code (Circuitscape 5.0) and later updates, we improved computational efficiency and parallelism, achieving major speed improvements, and enabling assessments across larger extents or with higher resolution data. Here, we reflect on the benefits to conservation of strengthening collaborations with computer scientists, and extract examples from a collection of 572 Circuitscape applications to illustrate how through a decade of repeated investment in the software, applications have been many, varied, and increasingly dynamic. Beyond empowering continued innovations in dynamic connectivity, we expect that faster run times will play an important role in facilitating co-production of connectivity assessments with stakeholders, increasing the likelihood that connectivity science will be incorporated in land use decisions.

Genomic and phylodynamic analysis of sapoviruses isolated in Henan Province, China

In this study, we determined the near-complete and partial genome sequences of ten SaV isolates. Phylogenetic analysis based on full-length VP1 and RdRp nucleotide sequences indicated that nine isolates were of GI.1 and one was GII.3. Evolutionary dynamics analysis indicated that GI.1 and GII.3 SaVs evolved at different rates, the latter evolving more rapidly. Cluster analysis indicated that distantly related GI.1 SaVs were more similar in their amino acid compositions than were GII.3 SaVs. The data provided in this study may facilitate studies on SaV genomic diversity and epidemiological patterns in China and worldwide.

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.

Genomic epidemiology and evolutionary dynamics of respiratory syncytial virus group B in Kilifi, Kenya, 2015-17

Respiratory syncytial virus (RSV) circulates worldwide, occurring seasonally in communities, and is a leading cause of acute respiratory illness in young children. There is paucity of genomic data from purposively sampled populations by which to investigate evolutionary dynamics and transmission patterns of RSV. Here we present an analysis of 295 RSV group B (RSVB) genomes from Kilifi, coastal Kenya, sampled from individuals seeking outpatient care in nine health facilities across a defined geographical area (∼890 km²), over two RSV epidemics between 2015 and 2017. RSVB diversity was characterized by multiple virus introductions into the area and co-circulation of distinct genetic clusters, which transmitted and diversified locally with varying frequency. Increase in relative genetic diversity paralleled seasonal virus incidence. Importantly, we identified a cluster of viruses that emerged in the 2016/17 epidemic, carrying distinct amino-acid signatures including a novel nonsynonymous change (K68Q) in antigenic site ∅ in the Fusion protein. RSVB diversity was additionally marked by signature nonsynonymous substitutions that were unique to particular genomic clusters, some under diversifying selection. Our findings provide insights into recent evolutionary and epidemiological behaviors of RSVB, and highlight possible emergence of a novel antigenic variant, which has implications on current prophylactic strategies in development.

Assessing Genome-Wide Diversity in European Hantaviruses through Sequence Capture from Natural Host Samples

Research on the ecology and evolution of viruses is often hampered by the limitation of sequence information to short parts of the genomes or single genomes derived from cultures. In this study, we use hybrid sequence capture enrichment in combination with high-throughput sequencing to provide efficient access to full genomes of European hantaviruses from rodent samples obtained in the field. We applied this methodology to Tula (TULV) and Puumala (PUUV) orthohantaviruses for which analyses from natural host samples are typically restricted to partial sequences of their tri-segmented RNA genome. We assembled a total of ten novel hantavirus genomes de novo with very high coverage (on average >99%) and sequencing depth (average >247×). A comparison with partial Sanger sequences indicated an accuracy of >99.9% for the assemblies. An analysis of two common vole (Microtus arvalis) samples infected with two TULV strains each allowed for the de novo assembly of all four TULV genomes. Combining the novel sequences with all available TULV and PUUV genomes revealed very similar patterns of sequence diversity along the genomes, except for remarkably higher diversity in the non-coding region of the S-segment in PUUV. The genomic distribution of polymorphisms in the coding sequence was similar between the species, but differed between the segments with the highest sequence divergence of 0.274 for the M-segment, 0.265 for the S-segment, and 0.248 for the L-segment (overall 0.258). Phylogenetic analyses showed the clustering of genome sequences consistent with their geographic distribution within each species. Genome-wide data yielded extremely high node support values, despite the impact of strong mutational saturation that is expected for hantavirus sequences obtained over large spatial distances. We conclude that genome sequencing based on capture enrichment protocols provides an efficient means for ecological and evolutionary investigations of hantaviruses at an unprecedented completeness and depth.

Spatial and Temporal Evolutionary Patterns in Puumala Orthohantavirus (PUUV) S Segment

The S segment of bank vole (Clethrionomys glareolus)-associated Puumala orthohantavirus (PUUV) contains two overlapping open reading frames coding for the nucleocapsid (N) and a non-structural (NSs) protein. To identify the influence of bank vole population dynamics on PUUV S segment sequence evolution and test for spillover infections in sympatric rodent species, during 2010–2014, 883 bank voles, 357 yellow-necked mice (Apodemus flavicollis), 62 wood mice (A. sylvaticus), 149 common voles (Microtus arvalis) and 8 field voles (M. agrestis) were collected in Baden-Wuerttemberg and North Rhine-Westphalia, Germany. In total, 27.9% and 22.3% of bank voles were positive for PUUV-reactive antibodies and PUUV-specific RNA, respectively. One of eight field voles was PUUV RNA-positive, indicating a spillover infection, but none of the other species showed evidence of PUUV infection. Phylogenetic and isolation-by-distance analyses demonstrated a spatial clustering of PUUV S segment sequences. In the hantavirus outbreak years 2010 and 2012, PUUV RNA prevalence was higher in our study regions compared to non-outbreak years 2011, 2013 and 2014. NSs amino acid and nucleotide sequence types showed temporal and/or local variation, whereas the N protein was highly conserved in the NSs overlapping region and, to a lower rate, in the N alone coding part.

Meeting report: Eleventh International Conference on Hantaviruses

The 2019 11th International Conference on Hantaviruses (ICH 2019) was organized by the International Society for Hantaviruses (ISH), and held on September 1-4, 2019, at the Irish College, in Leuven, Belgium. These ICHs have been held every three years since 1989. ICH 2019 was attended by 158 participants from 33 countries. The current report summarizes research presented on all aspects of hantavirology: ecology; pathogenesis and immune responses; virus phylogeny, replication and morphogenesis; epidemiology; vaccines, therapeutics and prevention; and clinical aspects and diagnosis.

Synanthropic rodents as virus reservoirs and transmitters

This review focuses on reports of hepatitis E virus, hantavirus, rotavirus, coronavirus, and arenavirus in synanthropic rodents (Rattus rattus, Rattus norvegicus, and Mus musculus) within urban environments. Despite their potential impact on human health, relatively few studies have addressed the monitoring of these viruses in rodents. Comprehensive control and preventive activities should include actions such as the elimination or reduction of rat and mouse populations, sanitary education, reduction of shelters for the animals, and restriction of the access of rodents to residences, water, and food supplies.

Phylogeography of Puumala orthohantavirus in Europe

Puumala virus is an RNA virus hosted by the bank vole (Myodes glareolus) and is today present in most European countries. Whilst it is generally accepted that hantaviruses have been tightly co-evolving with their hosts, Puumala virus (PUUV) evolutionary history is still controversial and so far has not been studied at the whole European level. This study attempts to reconstruct the phylogeographical spread of modern PUUV throughout Europe during the last postglacial period in the light of an upgraded dataset of complete PUUV small (S) segment sequences and by using most recent computational approaches. Taking advantage of the knowledge on the past migrations of its host, we identified at least three potential independent dispersal routes of PUUV during postglacial recolonization of Europe by the bank vole. From the Alpe-Adrian region (Balkan, Austria, and Hungary) to Western European countries (Germany, France, Belgium, and Netherland), and South Scandinavia. From the vicinity of Carpathian Mountains to the Baltic countries and to Poland, Russia, and Finland. The dissemination towards Denmark and North Scandinavia is more hypothetical and probably involved several independent streams from south and north Fennoscandia.