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

Molecular evolution of Hokkaido virus, a genotype of Orthohantavirus puumalaense, among Myodes rodents

Viruses in the genus Orthohantavirus within the family Hantaviridae cause human hantavirus infections and represent a threat to public health. Hokkaido virus (HOKV), a genotype of Orthohantavirus puumalaense (Puumala virus; PUUV), was first identified in Tobetsu, Hokkaido, Japan. Although it is genetically related to the prototype of PUUV, the evolutionary pathway of HOKV is unclear. We conducted a field survey in a forest in Tobetsu in 2022 and captured 44 rodents. Complete coding genome sequences of HOKVs were obtained from five viral-RNA-positive rodents (four Myodes rufocanus bedfordiae and one Apodemus speciosus). Phylogenetic analysis revealed a close relationship between the phylogenies and geographical origins of M. rufocanus-related orthohantaviruses. Comparison of the phylogenetic trees of the S segments of orthohantaviruses and the cytochrome b genes of Myodes species suggested that Myodes-related orthohantaviruses evolved in Myodes rodent species as a result of genetic isolation and host switching.

Cassava brown streak virus evolves with a nucleotide-substitution rate that is typical for the family Potyviridae

The ipomoviruses (family Potyviridae) that cause cassava brown streak disease (cassava brown streak virus [CBSV] and Uganda cassava brown streak virus [UCBSV]) are damaging plant pathogens that affect the sustainability of cassava production in East and Central Africa. However, little is known about the rate at which the viruses evolve and when they emerged in Africa - which inform how easily these viruses can host shift and resist RNAi approaches for control. We present here the rates of evolution determined from the coat protein gene (CP) of CBSV (Temporal signal in a UCBSV dataset was not sufficient for comparable analysis). Our BEAST analysis estimated the CBSV CP evolves at a mean rate of 1.43 × 10-3 nucleotide substitutions per site per year, with the most recent common ancestor of sampled CBSV isolates existing in 1944 (95% HPD, between years 1922 - 1963). We compared the published measured and estimated rates of evolution of CPs from ten families of plant viruses and showed that CBSV is an average-evolving potyvirid, but that members of Potyviridae evolve more quickly than members of Virgaviridae and the single representatives of Betaflexiviridae, Bunyaviridae, Caulimoviridae and Closteroviridae.

RNA virus diversity in rodents

Many zoonotic disease emergencies are associated with RNA viruses in rodents that substantially impact public health. With the widespread application of meta-genomics and meta-transcriptomics for virus discovery over the last decade, viral sequences deposited in public databases have expanded rapidly, and the number of novel viruses discovered in rodents has increased. As important reservoirs of zoonotic viruses, rodents have attracted increasing attention for the risk of potential spillover of rodent-borne viruses. However, knowledge of rodent viral diversity and the major factors contributing to the risk of zoonotic epidemic outbreaks remains limited. Therefore, this study analyzes the diversity and composition of rodent RNA viruses using virus records from the Database of Rodent-associated Viruses (DRodVir/ZOVER), which covers the published literatures and records in GenBank database, reviews the main rodent RNA virus-induced human infectious diseases, and discusses potential challenges in this field.

HantaNet: A New MicrobeTrace Application for Hantavirus Classification, Genomic Surveillance, Epidemiology and Outbreak Investigations

Hantaviruses zoonotically infect humans worldwide with pathogenic consequences and are mainly spread by rodents that shed aerosolized virus particles in urine and feces. Bioinformatics methods for hantavirus diagnostics, genomic surveillance and epidemiology are currently lacking a comprehensive approach for data sharing, integration, visualization, analytics and reporting. With the possibility of hantavirus cases going undetected and spreading over international borders, a significant reporting delay can miss linked transmission events and impedes timely, targeted public health interventions. To overcome these challenges, we built HantaNet, a standalone visualization engine for hantavirus genomes that facilitates viral surveillance and classification for early outbreak detection and response. HantaNet is powered by MicrobeTrace, a browser-based multitool originally developed at the Centers for Disease Control and Prevention (CDC) to visualize HIV clusters and transmission networks. HantaNet integrates coding gene sequences and standardized metadata from hantavirus reference genomes into three separate gene modules for dashboard visualization of phylogenetic trees, viral strain clusters for classification, epidemiological networks and spatiotemporal analysis. We used 85 hantavirus reference datasets from GenBank to validate HantaNet as a classification and enhanced visualization tool, and as a public repository to download standardized sequence data and metadata for building analytic datasets. HantaNet is a model on how to deploy MicrobeTrace-specific tools to advance pathogen surveillance, epidemiology and public health globally.

Host phylogeny shapes viral transmission networks in an island ecosystem

Virus transmission between host species underpins disease emergence. Both host phylogenetic relatedness and aspects of their ecology, such as species interactions and predator-prey relationships, may govern rates and patterns of cross-species virus transmission and hence zoonotic risk. To address the impact of host phylogeny and ecology on virus diversity and evolution, we characterized the virome structure of a relatively isolated island ecological community in Fiordland, New Zealand, that are linked through a food web. We show that phylogenetic barriers that inhibited cross-species virus transmission occurred at the level of host phyla (between the Chordata, Arthropoda and Streptophyta) as well as at lower taxonomic levels. By contrast, host ecology, manifest as predator-prey interactions and diet, had a smaller influence on virome composition, especially at higher taxonomic levels. The virus-host community comprised a 'small world' network, in which hosts with a high diversity of viruses were more likely to acquire new viruses, and generalist viruses that infect multiple hosts were more likely to infect additional species compared to host specialist viruses. Such a highly connected ecological community increases the likelihood of cross-species virus transmission, particularly among closely related species, and suggests that host generalist viruses present the greatest risk of disease emergence.

A framework for understanding and predicting orthohantavirus functional traits

Orthohantaviruses present a global public health threat; there are 58 distinct viruses currently recognized and case fatality of pathogenic orthohantaviruses ranges from <0.1% to 50%. An Old World versus New World dichotomy is frequently applied to distinguish human diseases caused by orthohantaviruses. However, this geographic grouping masks the importance of phylogeny and virus-host ecology in shaping orthohantavirus traits, especially since related arvicoline rodents and their orthohantaviruses are found in both regions. We argue that orthohantaviruses can be separated into three phylogenetically based rodent host groups with differences in key functional traits, including human disease, transmission route, and virus-host fidelity. This framework can help understand and predict traits of under-studied and newly discovered orthohantaviruses and guide public health and biosafety policy.

Deep Sequencing to Reveal Phylo-Geographic Relationships of Juquitiba Virus in Paraguay

Several hantaviruses result in zoonotic infections of significant public health concern, causing hemorrhagic fever with renal syndrome (HFRS) or hantavirus cardiopulmonary syndrome (HCPS) in the Old and New World, respectively. Given a 35% case fatality rate, disease-causing New World hantaviruses require a greater understanding of their biology, genetic diversity, and geographical distribution. Juquitiba hantaviruses have been identified in Oligoryzomys nigripes in Brazil, Paraguay, and Uruguay. Brazil has reported the most HCPS cases associated with this virus. We used a multiplexed, amplicon-based PCR strategy to screen and deep-sequence the virus harbored within lung tissues collected from Oligoryzomys species during rodent field collections in southern (Itapúa) and western (Boquerón) Paraguay. No Juquitiba-like hantaviruses were identified in Boquerón. Herein, we report the full-length S and M segments of the Juquitiba hantaviruses identified in Paraguay from O. nigripes. We also report the phylogenetic relationships of the Juquitiba hantaviruses in rodents collected from Itapúa with those previously collected in Canindeyú. We showed, using the TN93 nucleotide substitution model, the coalescent (constant-size) population tree model, and Bayesian inference implemented in the Bayesian evolutionary analysis by sampling trees (BEAST) framework, that the Juquitiba virus lineage in Itapúa is distinct from that in Canindeyú. Our spatiotemporal analysis showed significantly different time to the most recent ancestor (TMRA) estimates between the M and S segments, but a common geographic origin. Our estimates suggest the additional geographic diversity of the Juquitiba virus within the Interior Atlantic Forest and highlight the need for more extensive sampling across this biome.

A Brief History of Bunyaviral Family Hantaviridae

The discovery of Hantaan virus as an etiologic agent of hemorrhagic fever with renal syndrome in South Korea in 1978 led to identification of related pathogenic and nonpathogenic rodent-borne viruses in Asia and Europe. Their global distribution was recognized in 1993 after connecting newly discovered relatives of these viruses to hantavirus pulmonary syndrome in the Americas. The 1971 description of the shrew-infecting Hantaan-virus-like Thottapalayam virus was long considered an anomaly. Today, this virus and many others that infect eulipotyphlans, bats, fish, rodents, and reptiles are classified among several genera in the continuously expanding family Hantaviridae.

Lions and brown bears colonized North America in multiple synchronous waves of dispersal across the Bering Land Bridge

The Bering Land Bridge connecting North America and Eurasia was periodically exposed and inundated by oscillating sea levels during the Pleistocene glacial cycles. This land connection allowed the intermittent dispersal of animals, including humans, between Western Beringia (far northeast Asia) and Eastern Beringia (northwest North America), changing the faunal community composition of both continents. The Pleistocene glacial cycles also had profound impacts on temperature, precipitation and vegetation, impacting faunal community structure and demography. While these palaeoenvironmental impacts have been studied in many large herbivores from Beringia (e.g., bison, mammoths, horses), the Pleistocene population dynamics of the diverse guild of carnivorans present in the region are less well understood, due to their lower abundances. In this study, we analyse mitochondrial genome data from ancient brown bears (Ursus arctos; n = 103) and lions (Panthera spp.; n = 39), two megafaunal carnivorans that dispersed into North America during the Pleistocene. Our results reveal striking synchronicity in the population dynamics of Beringian lions and brown bears, with multiple waves of dispersal across the Bering Land Bridge coinciding with glacial periods of low sea levels, as well as synchronous local extinctions in Eastern Beringia during Marine Isotope Stage 3. The evolutionary histories of these two taxa underline the crucial biogeographical role of the Bering Land Bridge in the distribution, turnover and maintenance of megafaunal populations in North America.

The dinucleotide composition of sugarcane mosaic virus is shaped more by protein coding regions than by host species

Sugarcane mosaic virus (SCMV), which belongs to the Potyvirus genus of the family Potyviridae, causes mosaic diseases in canna, sugarcane and maize worldwide. Previously, the genetic variations, timescale, codon usage patterns and host adaptions of SCMV were determined. However, the dinucleotide composition and the dinucleotide bias from hosts or the protein coding regions of the virus have yet to be investigated. In this study, comprehensive analyses of the dinucleotide composition and dinucleotide bias from hosts, lineages and protein coding regions of SCMV were performed using 131 complete genomic sequences. We found that UpG and CpA were largely overrepresented while UpA, CpC, and CpG were largely underrepresented in the polyprotein and 11 protein coding region data sets. SCMV dinucleotide composition bias is more strongly dependent on the protein coding regions than on hosts. A weak association between the dinucleotide composition and SCMV lineages was also observed. Our analysis provides a novel perspective on the molecular evolutionary mechanisms of SCMV and may provide a better understanding of future research on the origin and evolutionary patterns of SCMV.

Identifying co-phylogenetic hotspots for zoonotic disease

The incidence of zoonotic diseases is increasing worldwide, which makes identifying parasites likely to become zoonotic and hosts likely to harbour zoonotic parasites a critical concern. Prior work indicates that there is a higher risk of zoonotic spillover accruing from closely related hosts and from hosts that are infected with a high phylogenetic diversity of parasites. This suggests that host and parasite evolutionary history may be important drivers of spillover, but identifying whether host-parasite associations are more strongly structured by the host, parasite or both requires co-phylogenetic analyses that combine host-parasite association data with host and parasite phylogenies. Here, we use host-parasite datasets containing associations between helminth taxa and free-range mammals in combination with phylogenetic models to explore whether host, parasite, or both host and parasite evolutionary history influences host-parasite associations. We find that host phylogenetic history is most important for driving patterns of helminth-mammal association, indicating that zoonoses are most likely to come from a host's close relatives. More broadly, our results suggest that co-phylogenetic analyses across broad taxonomic scales can provide a novel perspective for surveying potential emerging infectious diseases. This article is part of the theme issue 'Infectious disease macroecology: parasite diversity and dynamics across the globe'.

Discovery and Genetic Characterization of Novel Paramyxoviruses Related to the Genus Henipavirus in Crocidura Species in the Republic of Korea

Paramyxoviruses, negative-sense single-stranded RNA viruses, pose a critical threat to human public health. Currently, 78 species, 17 genera, and 4 subfamilies of paramyxoviruses are harbored by multiple natural reservoirs, including rodents, bats, birds, reptiles, and fish. Henipaviruses are critical zoonotic pathogens that cause severe acute respiratory distress and neurological diseases in humans. Using reverse transcription-polymerase chain reaction, 115 Crocidura species individuals were examined for the prevalence of paramyxovirus infections. Paramyxovirus RNA was observed in 26 (22.6%) shrews collected at five trapping sites, Republic of Korea. Herein, we report two genetically distinct novel paramyxoviruses (genus: Henipavirus): Gamak virus (GAKV) and Daeryong virus (DARV) isolated from C. lasiura and C. shantungensis, respectively. Two GAKVs and one DARV were nearly completely sequenced using next-generation sequencing. GAKV and DARV contain six genes (3′-N-P-M-F-G-L-5′) with genome sizes of 18,460 nucleotides and 19,471 nucleotides, respectively. The phylogenetic inference demonstrated that GAKV and DARV form independent genetic lineages of Henipavirus in Crocidura species. GAKV-infected human lung epithelial cells elicited the induction of type I/III interferons, interferon-stimulated genes, and proinflammatory cytokines. In conclusion, this study contributes further understandings of the molecular prevalence, genetic characteristics and diversity, and zoonotic potential of novel paramyxoviruses in shrews.

Phylodynamic and phylogeographic analysis of the complete genome of the West Nile virus lineage 2 (WNV-2) in the Mediterranean basin

Background

The West Nile virus is a highly contagious agent for a wide range of hosts. Its spread in the Mediterranean region raises several questions about its origin and the risk factors underlying the virus’s dispersal.

Materials and methods

The present study aims to reconstruct the temporal and spatial phylodynamics of West Nile virus lineage 2 in the Mediterranean region using 75 complete genome sequences from different host species retrieved from international databases.

Results

This data set suggests that current strains of WNV-2 began spreading in South Africa or nearby regions in the early twentieth century, and it migrated northwards via at least one route crossing the Mediterranean to reach Hungary in the early 2000s, before spreading throughout Europe. Another introduction event, according to the data set collected and analyses performed, is inferred to have occurred in around 1978. Migratory birds constitute, among others, additional risk factors that enhance the geographical transmission of the infection.

Conclusion

Our data underline the importance of the spatial–temporal tracking of migratory birds and phylodynamic reconstruction in setting up an efficient surveillance system for emerging and reemerging zoonoses in the Mediterranean region.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01902-w.

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.

Rodent-Borne Orthohantaviruses in Vietnam, Madagascar and Japan

Hantaviruses are harbored by multiple small mammal species in Asia, Europe, Africa, and the Americas. To ascertain the geographic distribution and virus-host relationships of rodent-borne hantaviruses in Japan, Vietnam, Myanmar, and Madagascar, RNAlater™-preserved lung tissues of 981 rodents representing 40 species, collected in 2011-2017, were analyzed for hantavirus RNA by RT-PCR. Our data showed Hantaan orthohantavirus Da Bie Shan strain in the Chinese white-bellied rat (Niviventer confucianus) in Vietnam, Thailand; orthohantavirus Anjo strain in the black rat (Rattus rattus) in Madagascar; and Puumala orthohantavirus Hokkaido strain in the grey-sided vole (Myodes rufocanus) in Japan. The Hokkaido strain of Puumala virus was also detected in the large Japanese field mouse (Apodemus speciosus) and small Japanese field mouse (Apodemus argenteus), with evidence of host-switching as determined by co-phylogeny mapping.

Edging on Mutational Bias, Induced Natural Selection From Host and Natural Reservoirs Predominates Codon Usage Evolution in Hantaan Virus

The molecular evolutionary dynamics that shape hantaviruses’ evolution are poorly understood even now, besides the contribution of virus-host interaction to their evolution remains an open question. Our study aimed to investigate these two aspects in Hantaan virus (HTNV)—the prototype of hantaviruses and an emerging zoonotic pathogen that infects humans, causing hemorrhagic fever with renal syndrome (HFRS): endemic in Far East Russia, China, and South Korea—via a comprehensive, phylogenetic-dependent codon usage analysis. We found that host- and natural reservoir-induced natural selection is the primary determinant of its biased codon choices, exceeding the mutational bias effect. The phylogenetic analysis of HTNV strains resulted in three distinct clades: South Korean, Russian, and Chinese. An effective number of codon (ENC) analysis showed a slightly biased codon usage in HTNV genomes. Nucleotide composition and RSCU analyses revealed a significant bias toward A/U nucleotides and A/U-ended codons, indicating the potential influence of mutational bias on the codon usage patterns of HTNV. Via ENC-plot, Parity Rule 2 (PR2), and neutrality plot analyses, we would conclude the presence of both mutation pressure and natural selection effect in shaping the codon usage patterns of HTNV; however, natural selection is the dominant factor influencing its codon usage bias. Codon adaptation index (CAI), Relative codon deoptimization index (RCDI), and Similarity Index (SiD) analyses uncovered the intense selection pressure from the host (Human) and natural reservoirs (Striped field mouse and Chinese white-bellied rat) in shaping HTNV biased codon choices. Our study clearly revealed the evolutionary processes in HTNV and the role of virus-host interaction in its evolution. Moreover, it opens the door for a more comprehensive codon usage analysis for all hantaviruses species to determine their molecular evolutionary dynamics and adaptability to several hosts and environments. We believe that our research will help in a better and deep understanding of HTNV evolution that will serve its future basic research and aid live attenuated vaccines design.

Identification of a novel hantavirus strain in the root vole (Microtus oeconomus) in Lithuania, Eastern Europe

Hantaviruses are zoonotic pathogens that can cause subclinical to lethal infections in humans. In Europe, five orthohantaviruses are present in rodents: Myodes-associated Puumala orthohantavirus (PUUV), Microtus-associated Tula orthohantavirus, Traemmersee hantavirus (TRAV)/ Tatenale hantavirus (TATV)/ Kielder hantavirus, rat-borne Seoul orthohantavirus, and Apodemus-associated Dobrava-Belgrade orthohantavirus (DOBV). Human PUUV and DOBV infections were detected previously in Lithuania, but the presence of Microtus-associated hantaviruses is not known. For this study we screened 234 Microtus voles, including root voles (Microtus oeconomus), field voles (Microtus agrestis) and common voles (Microtus arvalis) from Lithuania for hantavirus infections. This initial screening was based on reverse transcription-polymerase chain reaction (RT-PCR) targeting the S segment and serological analysis. A novel hantavirus was detected in eight of 79 root voles tentatively named "Rusne virus" according to the capture location and complete genome sequences were determined. In the coding regions of all three genome segments, Rusne virus showed high sequence similarity to TRAV and TATV and clustered with Kielder hantavirus in phylogenetic analyses of partial S and L segment sequences. Pairwise evolutionary distance analysis confirmed Rusne virus as a strain of the species TRAV/TATV. Moreover, we synthesized the entire nucleocapsid (N) protein of Rusne virus in Saccharomyces cerevisiae. We observed cross-reactivity of antibodies raised against other hantaviruses, including PUUV, with this new N protein. ELISA investigation of all 234 voles detected Rusne virus-reactive antibodies exclusively in four of 79 root voles, all being also RNA positive, but not in any other vole species. In conclusion, the detection of Rusne virus RNA in multiple root voles at the same trapping site during three years and its absence in sympatric field voles suggests root voles as the reservoir host of this novel virus. Future investigations should evaluate host association of TRAV, TATV, Kielder virus and the novel Rusne virus and their evolutionary relationships.

Bayesian Phylogeographic Analysis Incorporating Predictors and Individual Travel Histories in BEAST

Advances in sequencing technologies have tremendously reduced the time and costs associated with sequence generation, making genomic data an important asset for routine public health practices. Within this context, phylogenetic and phylogeographic inference has become a popular method to study disease transmission. In a Bayesian context, these approaches have the benefit of accommodating phylogenetic uncertainty, and popular implementations provide the possibility to parameterize the transition rates between locations as a function of epidemiological and ecological data to reconstruct spatial spread while simultaneously identifying the main factors impacting the spatial spread dynamics. Recent developments enable researchers to make use of travel history data of infected individuals in the reconstruction of pathogen spread, offering increased inference accuracy and mitigating sampling bias. Here, we describe a detailed workflow to reconstruct the spatial spread of a pathogen through Bayesian phylogeographic analysis in discrete space using these novel approaches, implemented in BEAST. The individual protocols focus on how to incorporate molecular data, covariates of spread, and individual travel history data into the analysis. © 2021 Wiley Periodicals LLC. Basic Protocol 1: Creating a SARS-CoV-2 MSA using sequences from GISAID Basic Protocol 2: Setting up a discrete trait phylogeographic reconstruction in BEAUti Basic Protocol 3: Phylogeographic reconstruction incorporating travel history information Basic Protocol 4: Visualizing ancestral spatial trajectories for specific taxa.

Bayesian Evaluation of Temporal Signal in Measurably Evolving Populations

Phylogenetic methods can use the sampling times of molecular sequence data to calibrate the molecular clock, enabling the estimation of evolutionary rates and timescales for rapidly evolving pathogens and data sets containing ancient DNA samples. A key aspect of such calibrations is whether a sufficient amount of molecular evolution has occurred over the sampling time window, that is, whether the data can be treated as having come from a measurably evolving population. Here, we investigate the performance of a fully Bayesian evaluation of temporal signal (BETS) in sequence data. The method involves comparing the fit to the data of two models: a model in which the data are accompanied by the actual (heterochronous) sampling times, and a model in which the samples are constrained to be contemporaneous (isochronous). We conducted simulations under a wide range of conditions to demonstrate that BETS accurately classifies data sets according to whether they contain temporal signal or not, even when there is substantial among-lineage rate variation. We explore the behavior of this classification in analyses of five empirical data sets: modern samples of A/H1N1 influenza virus, the bacterium Bordetella pertussis, coronaviruses from mammalian hosts, ancient DNA from Hepatitis B virus, and mitochondrial genomes of dog species. Our results indicate that BETS is an effective alternative to other tests of temporal signal. In particular, this method has the key advantage of allowing a coherent assessment of the entire model, including the molecular clock and tree prior which are essential aspects of Bayesian phylodynamic analyses.

Immunological Responses to Seoul Orthohantavirus in Experimentally and Naturally Infected Brown Rats (Rattus norvegicus)

To clarify the mechanism of Seoul orthohantavirus (SEOV) persistence, we compared the humoral and cell-mediated immune responses to SEOV in experimentally and naturally infected brown rats. Rats that were experimentally infected by the intraperitoneal route showed transient immunoglobulin M (IgM) production, followed by an increased anti-SEOV immunoglobulin G (IgG) antibody response and maturation of IgG avidity. The level of SEOV-specific cytotoxic T lymphocytes (CTLs) peaked at 6 days after inoculation and the viral genome disappeared from serum. In contrast, naturally infected brown rats simultaneously had a high rate of SEOV-specific IgM and IgG antibodies (28/43). Most of the IgM-positive rats (24/27) had the SEOV genome in their lungs, suggesting that chronic SEOV infection was established in those rats. In female rats with IgG avidity maturation, the viral load in the lungs was decreased. On the other hand, there was no relationship between IgG avidity and viral load in the lungs in male rats. A CTL response was not detected in naturally infected rats. The difference between immune responses in the experimentally and naturally infected rats is associated with the establishment of chronic infection in natural hosts.

Phylodynamics and Codon Usage Pattern Analysis of Broad Bean Wilt Virus 2

Broad bean wilt virus 2 (BBWV-2), which belongs to the genus Fabavirus of the family Secoviridae, is an important pathogen that causes damage to broad bean, pepper, yam, spinach and other economically important ornamental and horticultural crops worldwide. Previously, only limited reports have shown the genetic variation of BBWV2. Meanwhile, the detailed evolutionary changes, synonymous codon usage bias and host adaptation of this virus are largely unclear. Here, we performed comprehensive analyses of the phylodynamics, reassortment, composition bias and codon usage pattern of BBWV2 using forty-two complete genome sequences of BBWV-2 isolates together with two other full-length RNA1 sequences and six full-length RNA2 sequences. Both recombination and reassortment had a significant influence on the genomic evolution of BBWV2. Through phylogenetic analysis we detected three and four lineages based on the ORF1 and ORF2 nonrecombinant sequences, respectively. The evolutionary rates of the two BBWV2 ORF coding sequences were 8.895 × 10⁻⁴ and 4.560 × 10⁻⁴ subs/site/year, respectively. We found a relatively conserved and stable genomic composition with a lower codon usage choice in the two BBWV2 protein coding sequences. ENC-plot and neutrality plot analyses showed that natural selection is the key factor shaping the codon usage pattern of BBWV2. Strong correlations between BBWV2 and broad bean and pepper were observed from similarity index (SiD), codon adaptation index (CAI) and relative codon deoptimization index (RCDI) analyses. Our study is the first to evaluate the phylodynamics, codon usage patterns and adaptive evolution of a fabavirus, and our results may be useful for the understanding of the origin of this virus.

Orthohantaviruses, Emerging Zoonotic Pathogens

Orthohantaviruses give rise to the emerging infections such as of hemorrhagic fever with renal syndrome (HFRS) and hantavirus pulmonary syndrome (HPS) in Eurasia and the Americas, respectively. In this review we will provide a comprehensive analysis of orthohantaviruses distribution and circulation in Eurasia and address the genetic diversity and evolution of Puumala orthohantavirus (PUUV), which causes HFRS in this region. Current data indicate that the geographical location and migration of the natural hosts can lead to the orthohantaviruses genetic diversity as the rodents adapt to the new environmental conditions. The data shows that a high level of diversity characterizes the genome of orthohantaviruses, and the PUUV genome is the most divergent. The reasons for the high genome diversity are mainly caused by point mutations and reassortment, which occur in the genome segments. However, it still remains unclear whether this diversity is linked to the disease's severity. We anticipate that the information provided in this review will be useful for optimizing and developing preventive strategies of HFRS, an emerging zoonosis with potentially very high mortality rates.

Genetic Diversity and Phylogeography of Thottapalayam thottimvirus (Hantaviridae) in Asian House Shrew (Suncus murinus) in Eurasia

Murid and cricetid rodents were previously believed to be the principal reservoir hosts of hantaviruses. Recently, however, multiple newfound hantaviruses have been discovered in shrews, moles, and bats, suggesting a complex evolutionary history. Little is known about the genetic diversity and geographic distribution of the prototype shrew-borne hantavirus, Thottapalayam thottimvirus (TPMV), carried by the Asian house shrew (Suncus murinus), which is widespread in Asia, Africa, and the Middle East. Comparison of TPMV genomic sequences from two Asian house shrews captured in Myanmar and Pakistan with TPMV strains in GenBank revealed that the Myanmar TPMV strain (H2763) was closely related to the prototype TPMV strain (VRC66412) from India. In the L-segment tree, on the other hand, the Pakistan TPMV strain (PK3629) appeared to be the most divergent, followed by TPMV strains from Nepal, then the Indian-Myanmar strains, and finally TPMV strains from China. The Myanmar strain of TPMV showed sequence similarity of 79.3-96.1% at the nucleotide level, but the deduced amino acid sequences showed a high degree of conservation of more than 94% with TPMV strains from Nepal, India, Pakistan, and China. Cophylogenetic analysis of host cytochrome b and TPMV strains suggested that the Pakistan TPMV strain was mismatched. Phylogenetic trees, based on host cytochrome b and cytochrome c oxidase subunit I genes of mitochondrial DNA, and on host recombination activating gene 1 of nuclear DNA, suggested that the Asian house shrew and Asian highland shrew (Suncus montanus) comprised a species complex. Overall, the geographic-specific clustering of TPMV strains in Asian countries suggested local host-specific adaptation. Additional in-depth studies are warranted to ascertain if TPMV originated in Asian house shrews on the Indian subcontinent.

Genetic diversity and evolution of Hantaan virus in China and its neighbors

BACKGROUND

Hantaan virus (HTNV; family Hantaviridae, order Bunyavirales) causes hemorrhagic fever with renal syndrome (HFRS), which has raised serious concerns in Eurasia, especially in China, Russia, and South Korea. Previous studies reported genetic diversity and phylogenetic features of HTNV in different parts of China, but the analyses from the holistic perspective are rare.

METHODOLOGY AND PRINCIPAL FINDINGS

To better understand HTNV genetic diversity and gene evolution, we analyzed all available complete sequences derived from the small (S) and medium (M) segments with bioinformatic tools. Eleven phylogenetic groups were defined and showed geographic clustering; 42 significant amino acid variant sites were found, and 19 of them were located in immune epitopes; nine recombinant events and eight reassortments with highly divergent sequences were found and analyzed. We found that sequences from Guizhou showed high genetic divergence, contributing to multiple lineages of the phylogenetic tree and also to the recombination and reassortment events. Bayesian stochastic search variable selection analysis revealed that Heilongjiang, Shaanxi, and Guizhou played important roles in HTNV evolution and migration; the virus may originate from Zhejiang Province in the eastern part of China; and the virus population size expanded from the 1980s to 1990s.

CONCLUSIONS/SIGNIFICANCE

These findings revealed the original and evolutionary features of HTNV, which will help to illustrate hantavirus epidemic trends, thus aiding in disease control and prevention.

A seventeenth-century Mycobacterium tuberculosis genome supports a Neolithic emergence of the Mycobacterium tuberculosis complex

BACKGROUND

Although tuberculosis accounts for the highest mortality from a bacterial infection on a global scale, questions persist regarding its origin. One hypothesis based on modern Mycobacterium tuberculosis complex (MTBC) genomes suggests their most recent common ancestor followed human migrations out of Africa approximately 70,000 years before present. However, studies using ancient genomes as calibration points have yielded much younger dates of less than 6000 years. Here, we aim to address this discrepancy through the analysis of the highest-coverage and highest-quality ancient MTBC genome available to date, reconstructed from a calcified lung nodule of Bishop Peder Winstrup of Lund (b. 1605-d. 1679).

RESULTS

A metagenomic approach for taxonomic classification of whole DNA content permitted the identification of abundant DNA belonging to the human host and the MTBC, with few non-TB bacterial taxa comprising the background. Genomic enrichment enabled the reconstruction of a 141-fold coverage M. tuberculosis genome. In utilizing this high-quality, high-coverage seventeenth-century genome as a calibration point for dating the MTBC, we employed multiple Bayesian tree models, including birth-death models, which allowed us to model pathogen population dynamics and data sampling strategies more realistically than those based on the coalescent.

CONCLUSIONS

The results of our metagenomic analysis demonstrate the unique preservation environment calcified nodules provide for DNA. Importantly, we estimate a most recent common ancestor date for the MTBC of between 2190 and 4501 before present and for Lineage 4 of between 929 and 2084 before present using multiple models, confirming a Neolithic emergence for the MTBC.

Tilapia lake virus (TiLV): Genomic epidemiology and its early origin

Tilapia lake virus (TiLV) is an emerging virus that is rapidly spreading across the world. Over the past 6 years (2014-2020), TiLV outbreaks had been reported in at least 16 countries, spanning three continents, including Asia, Africa, and America. Despite its enormous economic impact, its origin, evolution and epidemiology are still largely poorly characterized. Here, we report eight TiLV whole-genome sequences from Thailand sampled between 2014 and 2019. Together with publicly available sequences from various regions of the world, we estimated the origin of TiLV to be between 2003 and 2009, 5-10 years before the first report of the virus in Israel in 2014. Our analyses consistently showed that TiLV started to spread in 2000s, and reached its peak in 2014-2016, matching well with the timing of its first report. From 2016 onwards, the global TiLV population declined steadily. This could be a result of herd immunity building up in the fish population, and/or a reflection of a better awareness of the virus coupled with a better and more cautious protocol of Tilapia importation. Despite the fact that we included all publicly available sequences, our analyses revealed long unsampled histories of TiLVs in many countries, especially towards its basal diversification. This result highlights the lack and the need for systematic surveillance of TiLV in fish.

Genetic changes and host adaptability in sugarcane mosaic virus based on complete genome sequences

Sugarcane mosaic virus (SCMV), a member of the genus Potyvirus in the family Potyviridae, is an important pathogen that causes mosaic diseases in maize, sugarcane, canna and other graminaceous species worldwide. Previously, several reports have showed the genetic variation and population structure of SCMV. However, the evolutionary dynamics, synonymous codon usage pattern and adaptive evolution of the virus is unclear. In this study, we performed comprehensive analyses of phylodynamics, composition bias and codon usage of SCMV using 108 complete genomic sequences. Our phylogenetic analysis found six host- and geographically confined phylogenetic lineages within the SCMV non-recombinant isolates. We found a relatively stable and conserved genomic composition with a lower codon usage choice in the SCMV protein coding sequences. Mutation pressure and natural selection have shaped the codon usage patterns of the SCMV protein coding sequences with natural selection being the dominant factor. The codon adaptation index (CAI), relative codon deoptimization index (RCDI) and similarity index (SiD) analyses revealed a stronger correlation between SCMV and maize than between SCMV and sugarcane or canna. Our study is the first to evaluate the codon usage pattern of SCMV based on complete sequences and may provide a better understanding of the origin of SCMV and its evolutionary patterns for future research.

Neutralizing Monoclonal Antibodies against the Gn and the Gc of the Andes Virus Glycoprotein Spike Complex Protect from Virus Challenge in a Preclinical Hamster Model

Hantaviruses are the etiological agent of hemorrhagic fever with renal syndrome (HFRS) and hantavirus cardiopulmonary syndrome (HCPS). The latter is associated with case fatality rates ranging from 30% to 50%. HCPS cases are rare, with approximately 300 recorded annually in the Americas. Recently, an HCPS outbreak of unprecedented size has been occurring in and around Epuyén, in the southwestern Argentinian state of Chubut. Since November of 2018, at least 29 cases have been laboratory confirmed, and human-to-human transmission is suspected. Despite posing a significant threat to public health, no treatment or vaccine is available for hantaviral disease. Here, we describe an effort to identify, characterize, and develop neutralizing and protective antibodies against the glycoprotein complex (Gn and Gc) of Andes virus (ANDV), the causative agent of the Epuyén outbreak. Using murine hybridoma technology, we generated 19 distinct monoclonal antibodies (MAbs) against ANDV GnGc. When tested for neutralization against a recombinant vesicular stomatitis virus expressing the Andes glycoprotein (GP) (VSV-ANDV), 12 MAbs showed potent neutralization and 8 showed activity in an antibody-dependent cellular cytotoxicity reporter assay. Escape mutant analysis revealed that neutralizing MAbs targeted both the Gn and the Gc. Four MAbs that bound different epitopes were selected for preclinical studies and were found to be 100% protective against lethality in a Syrian hamster model of ANDV infection. These data suggest the existence of a wide array of neutralizing antibody epitopes on hantavirus GnGc with unique properties and mechanisms of action.IMPORTANCE Infections with New World hantaviruses are associated with high case fatality rates, and no specific vaccine or treatment options exist. Furthermore, the biology of the hantaviral GnGc complex, its antigenicity, and its fusion machinery are poorly understood. Protective monoclonal antibodies against GnGc have the potential to be developed into therapeutics against hantaviral disease and are also great tools to elucidate the biology of the glycoprotein complex.

Evolution and Global Transmission of a Multidrug-Resistant, Community-Associated Methicillin-Resistant Staphylococcus aureus Lineage from the Indian Subcontinent

The evolution and global transmission of antimicrobial resistance have been well documented for Gram-negative bacteria and health care-associated epidemic pathogens, often emerging from regions with heavy antimicrobial use. However, the degree to which similar processes occur with Gram-positive bacteria in the community setting is less well understood. In this study, we traced the recent origins and global spread of a multidrug-resistant, community-associated Staphylococcus aureus lineage from the Indian subcontinent, the Bengal Bay clone (ST772). We generated whole-genome sequence data of 340 isolates from 14 countries, including the first isolates from Bangladesh and India, to reconstruct the evolutionary history and genomic epidemiology of the lineage. Our data show that the clone emerged on the Indian subcontinent in the early 1960s and disseminated rapidly in the 1990s. Short-term outbreaks in community and health care settings occurred following intercontinental transmission, typically associated with travel and family contacts on the subcontinent, but ongoing endemic transmission was uncommon. Acquisition of a multidrug resistance integrated plasmid was instrumental in the emergence of a single dominant and globally disseminated clade in the early 1990s. Phenotypic data on biofilm, growth, and toxicity point to antimicrobial resistance as the driving force in the evolution of ST772. The Bengal Bay clone therefore combines the multidrug resistance of traditional health care-associated clones with the epidemiological transmission of community-associated methicillin-resistant S. aureus (MRSA). Our study demonstrates the importance of whole-genome sequencing for tracking the evolution of emerging and resistant pathogens. It provides a critical framework for ongoing surveillance of the clone on the Indian subcontinent and elsewhere.IMPORTANCE The Bengal Bay clone (ST772) is a community-associated and multidrug-resistant Staphylococcus aureus lineage first isolated from Bangladesh and India in 2004. In this study, we showed that the Bengal Bay clone emerged from a virulent progenitor circulating on the Indian subcontinent. Its subsequent global transmission was associated with travel or family contact in the region. ST772 progressively acquired specific resistance elements at limited cost to its fitness and continues to be exported globally, resulting in small-scale community and health care outbreaks. The Bengal Bay clone therefore combines the virulence potential and epidemiology of community-associated clones with the multidrug resistance of health care-associated S. aureus lineages. This study demonstrates the importance of whole-genome sequencing for the surveillance of highly antibiotic-resistant pathogens, which may emerge in the community setting of regions with poor antibiotic stewardship and rapidly spread into hospitals and communities across the world.

Evolutionary Rates and Phylogeographical Analysis of Odontoglossum Ringspot Virus Based on the 166 Coat Protein Gene Sequences

Odontoglossum ringspot virus (ORSV) is a member of the genus Tobamovirus. It is one of the most prevalent viruses infecting orchids worldwide. Earlier studies reported the genetic variability of ORSV isolates from Korea and China. However, the evolutionary rate, timescale, and phylogeographical analyses of ORSV were unclear. Twenty-one coat protein (CP) gene sequences of ORSV were determined in this study, and used them together with 145 CP sequences obtained from GenBank to infer the genetic diversities, evolutionary rate, timescale and migration of ORSV populations. Evolutionary rate of ORSV populations was 1.25 × 10^-3 nucleotides/site/y. The most recent common ancestors came from 30 year ago (95% confidence intervals, 26-40). Based on CP gene, ORSV migrated from mainland China and South Korea to Taiwan island, Germany, Australia, Singapore, and Indonesia, and it also circulated within east Asia. Our study is the first attempt to evaluate the evolutionary rates, timescales and migration dynamics of ORSV.

Proteomics Computational Analyses Suggest that the Antennavirus Glycoprotein Complex Includes a Class I Viral Fusion Protein (α-Penetrene) with an Internal Zinc-Binding Domain and a Stable Signal Peptide

A metatranscriptomic study of RNA viruses in cold-blooded vertebrates identified two related viruses from frogfish (Antennarius striatus) that represent a new genus Antennavirus in the family Arenaviridae (Order: Bunyavirales). Computational analyses were used to identify features common to class I viral fusion proteins (VFPs) in antennavirus glycoproteins, including an N-terminal fusion peptide, two extended alpha-helices, an intrahelical loop, and a carboxyl terminal transmembrane domain. Like mammarenavirus and hartmanivirus glycoproteins, the antennavirus glycoproteins have an intracellular zinc-binding domain and a long virion-associated stable signal peptide (SSP). The glycoproteins of reptarenaviruses are also class I VFPs, but do not contain zinc-binding domains nor do they encode SSPs. Divergent evolution from a common progenitor potentially explains similarities of antennavirus, mammarenavirus, and hartmanivirus glycoproteins, with an ancient recombination event resulting in a divergent reptarenavirus glycoprotein.

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.

Đakrông virus, a novel mobatvirus (Hantaviridae) harbored by the Stoliczka's Asian trident bat (Aselliscus stoliczkanus) in Vietnam

The recent discovery of genetically distinct shrew- and mole-borne viruses belonging to the newly defined family Hantaviridae (order Bunyavirales) has spurred an extended search for hantaviruses in RNAlater®-preserved lung tissues from 215 bats (order Chiroptera) representing five families (Hipposideridae, Megadermatidae, Pteropodidae, Rhinolophidae and Vespertilionidae), collected in Vietnam during 2012 to 2014. A newly identified hantavirus, designated Đakrông virus (DKGV), was detected in one of two Stoliczka’s Asian trident bats (Aselliscus stoliczkanus), from Đakrông Nature Reserve in Quảng Trị Province. Using maximum-likelihood and Bayesian methods, phylogenetic trees based on the full-length S, M and L segments showed that DKGV occupied a basal position with other mobatviruses, suggesting that primordial hantaviruses may have been hosted by ancestral bats.

Paleomicrobiology: Diagnosis and Evolution of Ancient Pathogens

The last century has witnessed progress in the study of ancient infectious disease from purely medical descriptions of past ailments to dynamic interpretations of past population health that draw upon multiple perspectives. The recent adoption of high-throughput DNA sequencing has led to an expanded understanding of pathogen presence, evolution, and ecology across the globe. This genomic revolution has led to the identification of disease-causing microbes in both expected and unexpected contexts, while also providing for the genomic characterization of ancient pathogens previously believed to be unattainable by available methods. In this review we explore the development of DNA-based ancient pathogen research, the specialized methods and tools that have emerged to authenticate and explore infectious disease of the past, and the unique challenges that persist in molecular paleopathology. We offer guidelines to mitigate the impact of these challenges, which will allow for more reliable interpretations of data in this rapidly evolving field of investigation.

Potato virus Y; the Andean connection

Potato virus Y (PVY) causes disease in potatoes and other solanaceous crops. The appearance of its necrogenic strains in the 1980s made it the most economically important virus of potatoes. We report the isolation and genomic sequences of 32 Peruvian isolates of PVY which, together with 428 published PVY genomic sequences, gave an alignment of 460 sequences. Of these 190 (41%) were non-recombinant, and 162 of these provided a dated phylogeny, that corresponds well with the likely history of PVY, and show that PVY originated in South America which is where potatoes were first domesticated. The most basal divergences of the PVY population produced the N and C: O phylogroups; the origin of the N phylogroup is clearly Andean, but that of the O and C phylogroups is unknown, although they may have been first to establish in European crops. The current PVY population originated around 156 CE. PVY was probably first taken from South America to Europe in the 16th century in tubers. Most of the present PVY diversity emerged in the second half of the 19th century, after the Phytophthora infestans epidemics of the mid-19th century destroyed the European crop and stimulated potato breeding. Imported breeding lines were shared, and there was no quarantine. The early O population was joined later by N phylogroup isolates and their recombinants generated the R1 and R2 populations of damaging necrogenic strains. Our dating study has confirmed that human activity has dominated the phylodynamics of PVY for the last two millennia.

Bayesian inference of ancestral dates on bacterial phylogenetic trees

The sequencing and comparative analysis of a collection of bacterial genomes from a single species or lineage of interest can lead to key insights into its evolution, ecology or epidemiology. The tool of choice for such a study is often to build a phylogenetic tree, and more specifically when possible a dated phylogeny, in which the dates of all common ancestors are estimated. Here, we propose a new Bayesian methodology to construct dated phylogenies which is specifically designed for bacterial genomics. Unlike previous Bayesian methods aimed at building dated phylogenies, we consider that the phylogenetic relationships between the genomes have been previously evaluated using a standard phylogenetic method, which makes our methodology much faster and scalable. This two-step approach also allows us to directly exploit existing phylogenetic methods that detect bacterial recombination, and therefore to account for the effect of recombination in the construction of a dated phylogeny. We analysed many simulated datasets in order to benchmark the performance of our approach in a wide range of situations. Furthermore, we present applications to three different real datasets from recent bacterial genomic studies. Our methodology is implemented in a R package called BactDating which is freely available for download at https://github.com/xavierdidelot/BactDating.

Disease management at the wildlife-livestock interface: Using whole-genome sequencing to study the role of elk in Mycobacterium bovis transmission in Michigan, USA

The role of wildlife in the persistence and spread of livestock diseases is difficult to quantify and control. These difficulties are exacerbated when several wildlife species are potentially involved. Bovine tuberculosis (bTB), caused by Mycobacterium bovis, has experienced an ecological shift in Michigan, with spillover from cattle leading to an endemically infected white-tailed deer (deer) population. It has potentially substantial implications for the health and well-being of both wildlife and livestock and incurs a significant economic cost to industry and government. Deer are known to act as a reservoir of infection, with evidence of M. bovis transmission to sympatric elk and cattle populations. However, the role of elk in the circulation of M. bovis is uncertain; they are few in number, but range further than deer, so may enable long distance spread. Combining Whole Genome Sequences (WGS) for M. bovis isolates from exceptionally well-observed populations of elk, deer and cattle with spatiotemporal locations, we use spatial and Bayesian phylogenetic analyses to show strong spatiotemporal admixture of M. bovis isolates. Clustering of bTB in elk and cattle suggests either intraspecies transmission within the two populations, or exposure to a common source. However, there is no support for significant pathogen transfer amongst elk and cattle, and our data are in accordance with existing evidence that interspecies transmission in Michigan is likely only maintained by deer. This study demonstrates the value of whole genome population studies of M. bovis transmission at the wildlife-livestock interface, providing insights into bTB management in an endemic system.

Temporal analysis and adaptive evolution of the global population of potato virus M

Potato virus M (PVM), which is a member of the genus Carlavirus in the family Betaflexviridae, causes critical economic losses of nightshade crops. PVM is transmitted by aphids in a non-persistent manner, by sap inoculation and also transmitted in tubers. Previously, several reports described the genetic structure of PVM. However, the evolutionary rate, timescale, spread and adaptation evolution of the virus have not been examined. In this study, we investigated the phylodynamics of PVM using 145 nucleotide sequences of the coat protein gene and 117 sequences of the cysteine-rich nucleic acid-binding protein (NABP) gene, which were sampled between 1985 and 2013. We found that at least three lineages with isolates that were defined geographically but not by the original host were clustered. The evolutionary rate of the NABP (1.06 × 10^-2) was faster than that of the CP (4.12 × 10^-3). The time to the most recent common ancestors (TMRCAs) is similar between CP (CIs 31-110) and NABP (CIs 28-33) genes. Based on CP and NABP genes, PVM migrated from China to Canada, Iran, India and European countries, and it circulated within China. Our study is the first attempt to evaluate the evolutionary rates, timescales and migration dynamics of PVM.

The mystery of the phylogeographic structural pattern in rodent-borne hantaviruses

Hantaviruses (order Bunyavirales, family Hantaviridae) are important zoonotic pathogens. Because of the great diversity of their reservoir hosts, hantaviruses are excellent models to evaluate the dynamics of virus-host co-evolution. To understand the mechanisms behind the evolutionary history of hantaviruses through virus-reservoir interactions, it is important to know how the radiation and diversity of hantaviruses occurred. In this paper, we evaluate the pattern of hantavirus diversification based on a complete S segment representing major groups of hantaviruses found in the Americas. Phylogenetic analyses revealed a high degree of phylogeographic structure and a surprising pattern of geographical distribution of New World hantaviruses. The available data suggest that hantaviruses related to the Arvicolinae rodent subfamily in North America probably emerged and initially adapted from a shared common ancestor of the Tula virus. The first clade of hantaviruses associated with Neotominae occupied a stem lineage, especially those that emerged in Central America or Mexico. Hantaviruses from Central America and Mexico found in Neotominae rodents spread northward and probably gave rise to the first phylogroup of hantaviruses associated with Sigmodontinae in North America. Two preferential host-switching transmissions in hantaviruses apparently gave rise to two different paraphyletic group in Neotominae and Sigmodontinae. Our study supports a probable epicenter of diversification in Central America and/or Mexico for hantaviruses related to both the Neotominae and Sigmodontinae subfamilies.

Orthohantavirus genotype Lechiguanas in Oligoryzomys nigripes (Rodentia: Cricetidae): New evidence of host-switching

To identify and predict situations of increased risk of orthohantavirus infection in humans, it is necessary to study the relationships between the virus and its rodent hosts. The present study investigated orthohantavirus infection in an assemblage of wild Sigmodontinae rodents of the Paraná Delta, Argentina, and providing new evidence of host-switching events. Rodents belonging to the species Oxymycterus rufus (n = 187), Akodon azarae (n = 82), Oligoryzomys flavescens (n = 80), Oligoryzomys nigripes (n = 47), Scapteromys aquaticus (n = 38), Deltamys kempi (n = 7) and Holochilus brasiliensis (n = 2) were captured at 4 sampling sites during 20 trapping sessions. Blood samples were analyzed by IgG ELISA and livers by a nested reverse transcription PCR for the diagnosis of orthohantavirus infection. The amplified products of the S and M orthohantavirus genomes were sequenced and analyzed to determine similarities with species of the Orthohantavirus genus. The species of the Oligoryzomys positive to the virus were confirmed by amplifying and sequencing the complete cyt b gene. Of the 443 serum samples analyzed by IgG ELISA, A. azarae presented the highest host-specific prevalence value (10/82, 12.2%) followed by Ol. nigripes (4/47, 8.5%) and Ox. rufus (1/187, 0.5%). All the sero-positive Ol. nigripes (n = 4) were positive to the amplification of the S and M segments of the Lechiguanas genotype (98% nucleotide identity for both segments). This is surprising given that Ol. nigripes has been previously associated with Juquitiba genotype, not Lechiguanas. The latter is generally associated with Ol. flavescens, which in our study were all sero-negative. In addition, the association Ox. rufus - Pergamino genotype found here is, to our knowledge, novel and another potential evidence of host-switching considering that Pergamino has been originally associated with A. azarae. These findings contribute to the building evidence that contradicts the one-genotype-one-reservoir species premise in the association between rodent reservoirs and orthohantaviruses, and supports the hypothesis that the community structure of sympatric host species may contribute to orthohantavirus dynamics.

A Novel Hantavirus of the European Mole, Bruges Virus, Is Involved in Frequent Nova Virus Coinfections

Hantaviruses are zoonotic viruses with a complex evolutionary history of virus–host coevolution and cross-species transmission. Although hantaviruses have a broad reservoir host range, virus–host relationships were previously thought to be strict, with a single virus species infecting a single host species. Here, we describe Bruges virus, a novel hantavirus harbored by the European mole (Talpa europaea), which is the well-known host of Nova virus. Phylogenetic analyses of all three genomic segments showed tree topology inconsistencies, suggesting that Bruges virus has emerged from cross-species transmission and ancient reassortment events. A high number of coinfections with Bruges and Nova viruses was detected, but no evidence was found for reassortment between these two hantaviruses. These findings highlight the complexity of hantavirus evolution and the importance of further investigation of hantavirus–reservoir relationships.

Ancient Evolution of Mammarenaviruses: Adaptation via Changes in the L Protein and No Evidence for Host-Virus Codivergence

The Mammarenavirus genus includes several pathogenic species of rodent-borne viruses. Old World (OW) mammarenaviruses infect rodents in the Murinae subfamily and are mainly transmitted in Africa and Asia; New World (NW) mammarenaviruses are found in rodents of the Cricetidae subfamily in the Americas. We applied a selection-informed method to estimate that OW and NW mammarenaviruses diverged less than ∼45,000 years ago (ya). By incorporating phylogeographic inference, we show that NW mammarenaviruses emerged in the Latin America-Caribbean region ∼41,400–3,300 ya, whereas OW mammarenaviruses originated ∼23,100–1,880 ya, most likely in Southern Africa. Cophylogenetic analysis indicated that cospeciation did not contribute significantly to mammarenavirus–host associations. Finally, we show that extremely strong selective pressure on the viral polymerase accompanied the speciation of NW viruses. These data suggest that the evolutionary history of mammarenaviruses was not driven by codivergence with their hosts. The viral polymerase should be regarded as a major determinant of mammarenavirus adaptation.

Ancient human parvovirus B19 in Eurasia reveals its long-term association with humans

Human parvovirus B19 (B19V) is a ubiquitous human pathogen associated with a number of conditions, such as fifth disease in children and arthritis and arthralgias in adults. B19V is thought to evolve exceptionally rapidly among DNA viruses, with substitution rates previously estimated to be closer to those typical of RNA viruses. On the basis of genetic sequences up to ∼70 years of age, the most recent common ancestor of all B19V has been dated to the early 1800s, and it has been suggested that genotype 1, the most common B19V genotype, only started circulating in the 1960s. Here we present 10 genomes (63.9-99.7% genome coverage) of B19V from dental and skeletal remains of individuals who lived in Eurasia and Greenland from ∼0.5 to ∼6.9 thousand years ago (kya). In a phylogenetic analysis, five of the ancient B19V sequences fall within or basal to the modern genotype 1, and five fall basal to genotype 2, showing a long-term association of B19V with humans. The most recent common ancestor of all B19V is placed ∼12.6 kya, and we find a substitution rate that is an order of magnitude lower than inferred previously. Further, we are able to date the recombination event between genotypes 1 and 3 that formed genotype 2 to ∼5.0-6.8 kya. This study emphasizes the importance of ancient viral sequences for our understanding of virus evolution and phylogenetics.

Inferring demographic parameters in bacterial genomic data using Bayesian and hybrid phylogenetic methods

Background

Recent developments in sequencing technologies make it possible to obtain genome sequences from a large number of isolates in a very short time. Bayesian phylogenetic approaches can take advantage of these data by simultaneously inferring the phylogenetic tree, evolutionary timescale, and demographic parameters (such as population growth rates), while naturally integrating uncertainty in all parameters. Despite their desirable properties, Bayesian approaches can be computationally intensive, hindering their use for outbreak investigations involving genome data for a large numbers of pathogen isolates. An alternative to using full Bayesian inference is to use a hybrid approach, where the phylogenetic tree and evolutionary timescale are estimated first using maximum likelihood. Under this hybrid approach, demographic parameters are inferred from estimated trees instead of the sequence data, using maximum likelihood, Bayesian inference, or approximate Bayesian computation. This can vastly reduce the computational burden, but has the disadvantage of ignoring the uncertainty in the phylogenetic tree and evolutionary timescale.

Results

We compared the performance of a fully Bayesian and a hybrid method by analysing six whole-genome SNP data sets from a range of bacteria and simulations. The estimates from the two methods were very similar, suggesting that the hybrid method is a valid alternative for very large datasets. However, we also found that congruence between these methods is contingent on the presence of strong temporal structure in the data (i.e. clocklike behaviour), which is typically verified using a date-randomisation test in a Bayesian framework. To reduce the computational burden of this Bayesian test we implemented a date-randomisation test using a rapid maximum likelihood method, which has similar performance to its Bayesian counterpart.

Conclusions

Hybrid approaches can produce reliable inferences of evolutionary timescales and phylodynamic parameters in a fraction of the time required for fully Bayesian analyses. As such, they are a valuable alternative in outbreak studies involving a large number of isolates.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1210-5) contains supplementary material, which is available to authorized users.

Exploring the Robustness of the Parsimonious Reconciliation Method in Host-Symbiont Cophylogeny

The aim of this paper is to explore the robustness of the parsimonious host-symbiont tree reconciliation method under editing or small perturbations of the input. The editing involves making different choices of unique symbiont mapping to a host in the case where multiple associations exist. This is made necessary by the fact that the tree reconciliation model is currently unable to handle such associations. The analysis performed could however also address the problem of errors. The perturbations are re-rootings of the symbiont tree to deal with a possibly wrong placement of the root specially in the case of fast-evolving species. In order to do this robustness analysis, we introduce a simulation scheme specifically designed for the host-symbiont cophylogeny context, as well as a measure to compare sets of tree reconciliations, both of which are of interest by themselves.

A comparison of methods for estimating substitution rates from ancient DNA sequence data

Background

Phylogenetic analysis of DNA from modern and ancient samples allows the reconstruction of important demographic and evolutionary processes. A critical component of these analyses is the estimation of evolutionary rates, which can be calibrated using information about the ages of the samples. However, the reliability of these rate estimates can be negatively affected by among-lineage rate variation and non-random sampling. Using a simulation study, we compared the performance of three phylogenetic methods for inferring evolutionary rates from time-structured data sets: regression of root-to-tip distances, least-squares dating, and Bayesian inference. We also applied these three methods to time-structured mitogenomic data sets from six vertebrate species.

Results

Our results from 12 simulation scenarios show that the three methods produce reliable estimates when the substitution rate is high, rate variation is low, and samples of similar ages are not all grouped together in the tree (i.e., low phylo-temporal clustering). The interaction of these factors is particularly important for least-squares dating and Bayesian estimation of evolutionary rates. The three estimation methods produced consistent estimates of rates across most of the six mitogenomic data sets, with sequence data from horses being an exception.

Conclusions

We recommend that phylogenetic studies of ancient DNA sequences should use multiple methods of inference and test for the presence of temporal signal, among-lineage rate variation, and phylo-temporal clustering in the data.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1192-3) contains supplementary material, which is available to authorized users.

A geography-aware reconciliation method to investigate diversification patterns in host/parasite interactions

Cospeciation studies aim at investigating whether hosts and symbionts speciate simultaneously or whether the associations diversify through host shifts. This problem is often tackled through reconciliation analyses that map the symbiont phylogeny onto the host phylogeny by mixing different types of diversification events. These reconciliations can be difficult to interpret and are not always biologically realistic. Researchers have underlined that the biogeographic histories of both hosts and symbionts influence the probability of cospeciation and host switches, but up to now no reconciliation software integrates geographic data. We present a new functionality in the Mowgli software that bridges this gap. The user can provide geographic information on both the host and symbiont extant and ancestral taxa. Constraints in the reconciliation algorithm have been implemented to generate biologically realistic codiversification scenarios. We apply our method to the fig/fig wasp association and infer diversification scenarios that differ from reconciliations ignoring geographic information. In addition, we updated the reconciliation viewer SylvX to visualize ancestral character states on the phylogenetic trees and highlight parts of reconciliations that are geographically inconsistent when not accounting for geographic constraints. We suggest that the comparison of reconciliations obtained with and without such constraints can help solving ambiguities in the biogeographic histories of the partners. With the development of robust methods in historical biogeography, and the advent of next-generation sequencing that leads to better-resolved trees, a geography-aware reconciliation method represents a substantial advance that is likely to be useful to researchers studying the evolution of biotic interactions and biogeography.

Host shifts result in parallel genetic changes when viruses evolve in closely related species

Host shifts, where a pathogen invades and establishes in a new host species, are a major source of emerging infectious diseases. They frequently occur between related host species and often rely on the pathogen evolving adaptations that increase their fitness in the novel host species. To investigate genetic changes in novel hosts, we experimentally evolved replicate lineages of an RNA virus (Drosophila C Virus) in 19 different species of Drosophilidae and deep sequenced the viral genomes. We found a strong pattern of parallel evolution, where viral lineages from the same host were genetically more similar to each other than to lineages from other host species. When we compared viruses that had evolved in different host species, we found that parallel genetic changes were more likely to occur if the two host species were closely related. This suggests that when a virus adapts to one host it might also become better adapted to closely related host species. This may explain in part why host shifts tend to occur between related species, and may mean that when a new pathogen appears in a given species, closely related species may become vulnerable to the new disease.

Host shifts, where a pathogen jumps from one host species to another, are a major source of infectious disease. Hosts shifts are more likely to occur between related host species and often rely on the pathogen evolving adaptations that increase their fitness in the novel host. Here we have investigated how viruses evolve in different host species, by experimentally evolving replicate lineages of an RNA virus in 19 different host species that shared a common ancestor 40 million years ago. We then deep sequenced the genomes of these viruses to examine the genetic changes that have occurred in different host species that vary in their relatedness. We found that parallel mutations–that are indicative of selection–were significantly more likely to occur within viral lineages from the same host, and between viruses evolved in closely related species. This suggests that a mutation that may adapt a virus to a given host, may also adapt it to closely related host species.

Global Diversity and Distribution of Hantaviruses and Their Hosts

Rodents represent 42% of the world's mammalian biodiversity encompassing 2,277 species populating every continent (except Antarctica) and are reservoir hosts for a wide diversity of disease agents. Thus, knowing the identity, diversity, host-pathogen relationships, and geographic distribution of rodent-borne zoonotic pathogens, is essential for predicting and mitigating zoonotic disease outbreaks. Hantaviruses are hosted by numerous rodent reservoirs. However, the diversity of rodents harboring hantaviruses is likely unknown because research is biased toward specific reservoir hosts and viruses. An up-to-date, systematic review covering all known rodent hosts is lacking. Herein, we document gaps in our knowledge of the diversity and distribution of rodent species that host hantaviruses. Of the currently recognized 681 cricetid, 730 murid, 61 nesomyid, and 278 sciurid species, we determined that 11.3, 2.1, 1.6, and 1.1%, respectively, have known associations with hantaviruses. The diversity of hantaviruses hosted by rodents and their distribution among host species supports a reassessment of the paradigm that each virus is associated with a single-host species. We examine these host-virus associations on a global taxonomic and geographical scale with emphasis on the rodent host diversity and distribution. Previous reviews have been centered on the viruses and not the mammalian hosts. Thus, we provide a perspective not previously addressed.