Explaining the geographic spread of emerging epidemics: a framework for comparing viral phylogenies and environmental landscape data

Genesis and Spread of Novel Highly Pathogenic Avian Influenza A(H5N1) Clade 2.3.4.4b Virus Genotype EA-2023-DG Reassortant, Western Europe

In Europe, highly pathogenic avian influenza (HPAI) virus circulates in avian wildlife, undergoing frequent reassortment, sporadic introductions in domestic birds, and spillover to mammals. An H5N1 clade 2.3.4.4b reassortant, EA-2023-DG, affecting wild and domestic birds was detected in western Europe in November 2023. Six of its RNA segments came from the EA-2021-AB genotype, but the polymerase basic 2 and polymerase acidic segments originated from low pathogenicity avian influenza viruses. Discrete phylogeographic analyses of concatenated genomes and single polymerase basic 2 and polymerase acidic segments suggested reassortment in summer 2023 near the southwestern Baltic Sea. Subsequent continuous phylogeographic analysis of all concatenated EA-2023-DG genomes highlighted circulation in northwestern Europe until June 2024 and long-distance dispersal toward France, Norway, England, Slovakia, Switzerland, and Austria. Those results illustrate the value of phylodynamic approaches to investigate emergence of novel avian influenza virus variants, trace their subsequent dispersal history, and provide vital clues for informing outbreak prevention and intervention policies.

Tracing more than two decades of Japanese encephalitis virus circulation in mainland China

Japanese encephalitis is a viral disease caused by the Japanese encephalitis virus (JEV), primarily affecting rural areas of Asia and western Pacific region. China remains one of the main epicenters, experiencing a significant burden of human and animal cases despite vaccination efforts. The ecology of this arbovirus is complex, involving Culex mosquitoes as primary vectors, wading birds as natural reservoirs, and pigs as amplifying hosts. Given the virus's epidemiological importance in China, combined with the country's expanding pig farming industry and diverse climates, investigating the virus spread and its environmental drivers is needed to address its persistent burden. In this study, we conducted phylogeographic analyses by combining publicly available JEV envelope gene sequences from China and other regions. Our reconstructions revealed multiple introduction events leading to various circulating JEV clades in China, with one predominant clade. Additionally, our analyses showed a diffusion capacity of JEV exceeding previous estimates for co-circulating arboviruses. These differences could be attributed to pig trade or bird migration, calling for further investigations into the drivers of JEV spread.

IMPORTANCE

Japanese encephalitis virus (JEV) is the cause of Japanese encephalitis, a significant health concern in China. Despite being one of the most studied mosquito-borne viruses, no previous studies have combined genomic and geographic data to investigate the spatial epidemiology and dispersal capacity of the virus. In this study, we analyzed genomic, geographic, and environmental data to trace the dispersal history of JEV in China and explore the environmental factors influencing its distribution. Our findings show that JEV circulates predominantly in areas with higher temperatures, dense human and pig populations, and favorable conditions for Culex mosquitoes. Notably, our analyses showed a higher diffusion capacity of JEV compared to co-circulating viruses, possibly driven by factors like pig trade and bird migration. Our analysis calls for improved genomic surveillance and establishes a baseline for future studies on the effects of climate change, agricultural practices, and bird migration on JEV circulation.

Influence of Host and Landscape-Associated Factors in the Infection and Transmission of Pathogens: The Case of Directly Transmitted Virus in Mammals

BACKGROUND

Among pathogens associated with mammals, numerous viruses with a direct transmission route impact human, domestic and wild species health. Host and landscape factors affect viral infection and transmission dynamics of these viruses, along with barriers to host dispersal and gene exchange. However, studies show biases toward certain locations, hosts and detected pathogens, with regional variations in similar host-virus associations.

METHODS

Using a systematic review, in two electronic repositories for articles published until December 2022, we analysed the available information on host- and landscape-associated factors influencing the infection and transmission of directly transmitted viruses in mammals.

RESULTS

In the analysis, about 50% of papers examined either host traits, landscape composition or configuration measures, while approximately 24% combined host and landscape-associated factors. Additionally, approximately 17% of the articles included climatic data and 30% integrated factors related to anthropogenic impact, as these variables have a role in host density, distribution and virus persistence. The most significant and frequent host traits used as predictor variables were sex, age, body weight, host density and species identity. Land cover was the most evaluated landscape attribute, while some explored configuration variables like edge density and fragmentation indexes. Finally, temperature, precipitation and features such as human population density and human footprint index were also typically measured and found impactful.

CONCLUSION

Given the many contributions host- and landscape-related factors have in pathogen dynamics, this systematic study contributes to a better knowledge of host-virus dynamics and the identification of variables and gaps that can be used for disease prevention.

Detection of divergent Orthohantavirus tulaense provides insight into wide host range and viral evolutionary patterns

Orthohantavirus tulaense (TULV) is a member of the orthohantavirus genus and distributed in Europe and Asia. To shed light on TULV epidemiology and evolution, we trapped wild rodents from eastern Turkiye and found 15 TULV positive rodents. Sequencing and phylogenetic analyses confirmed the presence of diverse TULV strains. Global phylogenetic characterization suggested 5 distinct TULV lineages. Global phylogeographic reconstruction estimated different rooting times for each three segments, a potential ancestor location in Eastern Black Sea region, and strongly supported phylogeographic structure with 11 clusters. Dispersal velocity of TULV was estimated to be much faster than some other orthohantaviruses. Eastern Black Sea seemed to have lineages evolving faster and genetically closer to proto-Tula virus. Host switching estimates suggested potential switching events from Microtus arvalis to M. obscurus to M. irani with host-dependent sub-clustering within geographic clusters and suggested substantial evidence for no clear virus jumps from M. arvalis to M. irani.

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 that can be informative of the dispersal dynamics and the capacity to spread among hosts. Heterogeneous sampling efforts of genomic sequences can however impact the accuracy of phylogeographic dispersal metrics. While the impact of spatial sampling bias on the outcomes of continuous phylogeographic inference has previously been explored, the impact of sampling intensity (i.e., sampling size) when aiming to characterise dispersal patterns through continuous phylogeographic reconstructions has not yet been thoroughly evaluated. In our study, we use simulations to evaluate the robustness of 3 dispersal metrics - a lineage dispersal velocity, a diffusion coefficient, and an isolation-by-distance (IBD) signal metric - to the sampling intensity. Our results reveal that both the diffusion coefficient and IBD signal metrics appear to be the most robust to the number of samples considered for the phylogeographic reconstruction. We then use these 2 dispersal metrics to compare the dispersal pattern and capacity of various viruses spreading in animal populations. Our comparative analysis reveals a broad range of IBD 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 use of lineage dispersal metrics to consider in future studies and illustrates their application to compare the spread of viruses in various settings.

Spatio-temporal spread of Lassa virus and a new rodent host in the Mano River Union area, West Africa

The spread of Lassa virus (LASV) in Guinea, Liberia and Sierra Leone, which together are named the Mano River Union (MRU) area, was examined phylogeographically. To provide a reliable evolutionary scenario, new rodent-derived, whole LASV sequences were included. These were generated by metatranscriptomic next-generation sequencing from rodents sampled between 2003 and 2020 in 21 localities of Guinea and Sierra Leone. An analysis was performed using BEAST to perform continuous phylogeographic inference and EvoLaps v36 to visualize spatio-temporal spread. LASV was identified as expected in its primary host reservoir, the Natal multimammate mouse (Mastomys natalensis), and also in two Guinean multimammate mice (Mastomys erythroleucus) in northern Sierra Leone and two rusty-bellied brush-furred mice (Lophuromys sikapusi) in southern Sierra Leone. This finding is consistent with the latter two species being secondary host reservoirs. The strains in these three species were very closely related in LASV lineage IV. Phylogenetic analysis indicated that the most recent common ancestor of lineage IV existed 316-374 years ago and revealed distinct, well-supported clades from Sierra Leone (Bo, Kabala and Kenema), Guinea (Faranah, Kissidougou-Guekedou and Macenta) and Liberia (Phebe-Ganta). The phylogeographic scenario suggests southern Guinea as the point of origin of LASV in the MRU area, with subsequent spread to towards Mali, Liberia and Sierra Leone at a mean speed of 1.6 to 1.1 km/year.

Ecological Changes Exacerbating the Spread of Invasive Ticks has Driven the Dispersal of Severe Fever with Thrombocytopenia Syndrome Virus Throughout Southeast Asia

Severe fever with thrombocytopenia syndrome virus (SFTSV) is a tick-borne virus recognized by the World Health Organization as an emerging infectious disease of growing concern. Utilizing phylodynamic and phylogeographic methods, we have reconstructed the origin and transmission patterns of SFTSV lineages and the roles demographic, ecological, and climatic factors have played in shaping its emergence and spread throughout Asia. Environmental changes and fluctuations in tick populations, exacerbated by the widespread use of pesticides, have contributed significantly to its geographic expansion. The increased adaptability of Lineage L2 strains to the Haemaphysalis longicornis vector has facilitated the dispersal of SFTSV through Southeast Asia. Increased surveillance and proactive measures are needed to prevent further spread to Australia, Indonesia, and North America.

Ecological dynamics of moa extinctions reveal convergent refugia that today harbour flightless birds

Human settlement of islands across the Pacific Ocean was followed by waves of faunal extinctions that occurred so rapidly that their dynamics are difficult to reconstruct in space and time. These extinctions included large, wingless birds called moa that were endemic to New Zealand. Here we reconstructed the range and extinction dynamics of six genetically distinct species of moa across New Zealand at a fine spatiotemporal resolution, using hundreds of thousands of process-explicit simulations of climate-human-moa interactions, which were validated against inferences of occurrence and range contraction from an extensive fossil record. These process-based simulations revealed important interspecific differences in the ecological and demographic attributes of moa and established how these differences influenced likely trajectories of geographic and demographic declines of moa following Polynesian colonization of New Zealand. We show that despite these interspecific differences in extinction dynamics, the spatial patterns of geographic range collapse of moa species were probably similar. It is most likely that the final populations of all moa species persisted in suboptimal habitats in cold, mountainous areas that were generally last and least impacted by people. We find that these refugia for the last populations of moa continue to serve as isolated sanctuaries for New Zealand's remaining flightless birds, providing fresh insights for conserving endemic species in the face of current and future threats.

Spatiotemporal and Species-Crossing Transmission Dynamics of Subclade 2.3.4.4b H5Nx HPAIVs

Subclade 2.3.4.4b H5Nx highly pathogenic avian influenza (HPAI) viruses, emerged in 2013 with multiple subtypes of H5N8, H5N1, and H5N6, had unprecedently caused a global epizootic by H5N1 since 2021, which had devasted multiple species of wild birds, poultry, and wild mammals (terrestrial and marine) with a high mortality, causing severe ecological damage. The infected wild mammals may become new "mixers" for influenza viruses, posing the potential transmission to human. Frequent outbreaks of subclade 2.3.4.4b H5Nx viruses among wild birds and poultry had exposed major gaps in our knowledge on their evolution, spatiotemporal diffusion, and species-crossing transmission. Here, we integrated the phylogenetic and epidemiological data of subclade 2.3.4.4b H5Nx viruses in public database and used Bayesian phylodynamic analysis to reveal the pattern of the global large-scale transmission. Phylogenic analysis demonstrated that the HA gene of these viruses diverged into two dominant clusters around 2015 and 2016. The Bayesian phylodynamic analysis illustrated that the viruses presented spatiotemporally complex transmission network with geographical and host relative expansion and recombination with different subtypes of NA segment. Spatially, the Russian Federation (Siberia) was identified as the primary hub for virus transmission, which was further facilitated by the establishment of strong epidemiological linkages between West Europe and broader regions, such as North America. As for hosts, wild Anseriformes were the primary species for the virus spillover, contributing to the spatial expansion and rapid diffusion globally of subclade 2.3.4.4b viruses. We investigated the phylogeny of subclade 2.3.4.4b H5Nx viruses and the spatiotemporal pattern of transmission with initial location and the primary host, which could provide comprehensive insights for subclade 2.3.4.4b H5Nx viruses. Due to the wild birds involved the widespread of subclade 2.3.4.4b H5Nx viruses, the epizootics in poultry are inevitable, so we highly recommend to apply the policy of culling plus with vaccination to protect the poultry industry and potentially protect the public health.

Investigating the cause of a 2021 winter wave of COVID-19 in a border region in eastern Germany: a mixed-methods study, August to November 2021

It is so far unclear how the COVID-19 winter waves started and what should be done to prevent possible future waves. In this study, we deciphered the dynamic course of a winter wave in 2021 in Saxony, a state in Eastern Germany neighbouring the Czech Republic and Poland. The study was carried out through the integration of multiple virus genomic epidemiology approaches to track transmission chains, identify emerging variants and investigate dynamic changes in transmission clusters. For identified local variants of interest, functional evaluations were performed. Multiple long-lasting community transmission clusters have been identified acting as driving force for the winter wave 2021. Analysis of the dynamic courses of two representative clusters indicated a similar transmission pattern. However, the transmission cluster caused by a locally occurring new Delta variant AY.36.1 showed a distinct transmission pattern, and functional analyses revealed a replication advantage of it. This study indicated that long-lasting community transmission clusters starting since early autumn caused by imported or locally occurring variants all contributed to the development of the 2021 winter wave. The information we achieved might help future pandemic prevention.

Complete genome reconstruction of the global and European regional dispersal history of the lumpy skin disease virus

IMPORTANCE

Lumpy skin disease virus (LSDV) has a complex epidemiology involving multiple strains, recombination, and vaccination. Its DNA genome provides limited genetic variation to trace outbreaks in space and time. Sequencing of LSDV whole genomes has also been patchy at global and regional scales. Here, we provide the first fine-grained whole genome sequence sampling of a constrained LSDV outbreak (southeastern Europe, 2015-2017), which we analyze along with global publicly available genomes. We formally evaluate the past occurrence of recombination events as well as the temporal signal that is required for calibrating molecular clock models and subsequently conduct a time-calibrated spatially explicit phylogeographic reconstruction. Our study further illustrates the importance of accounting for recombination events before reconstructing global and regional dynamics of DNA viruses. More LSDV whole genomes from endemic areas are needed to obtain a comprehensive understanding of global LSDV dispersal dynamics.

Uncovering the endemic circulation of rabies in Cambodia

In epidemiology, endemicity characterizes sustained pathogen circulation in a geographical area, which involves a circulation that is not being maintained by external introductions. Because it could potentially shape the design of public health interventions, there is an interest in fully uncovering the endemic pattern of a disease. Here, we use a phylogeographic approach to investigate the endemic signature of rabies virus (RABV) circulation in Cambodia. Cambodia is located in one of the most affected regions by rabies in the world, but RABV circulation between and within Southeast Asian countries remains understudied. Our analyses are based on a new comprehensive data set of 199 RABV genomes collected between 2014 and 2017 as well as previously published Southeast Asian RABV sequences. We show that most Cambodian sequences belong to a distinct clade that has been circulating almost exclusively in Cambodia. Our results thus point towards rabies circulation in Cambodia that does not rely on external introductions. We further characterize within-Cambodia RABV circulation by estimating lineage dispersal metrics that appear to be similar to other settings, and by performing landscape phylogeographic analyses to investigate environmental factors impacting the dispersal dynamic of viral lineages. The latter analyses do not lead to the identification of environmental variables that would be associated with the heterogeneity of viral lineage dispersal velocities, which calls for a better understanding of local dog ecology and further investigations of the potential drivers of RABV spread in the region. Overall, our study illustrates how phylogeographic investigations can be performed to assess and characterize viral endemicity in a context of relatively limited data.

Identifying potential significant factors impacting zero-inflated proportion data

Classical supervised methods like linear regression and decision trees are not completely adapted for identifying impacting factors on a response variable corresponding to zero-inflated proportion data (ZIPD) that are dependent, continuous and bounded. In this article we propose a within-block permutation-based methodology to identify factors (discrete or continuous) that are significantly correlated with ZIPD, we propose a performance indicator quantifying the percentage of correlation explained by the subset of significant factors, and we show how to predict the ranks of the response variables conditionally on the observation of these factors. The methodology is illustrated on simulated data and on two real data sets dealing with epidemiology. In the first data set, ZIPD correspond to probabilities of transmission of Influenza between horses. In the second data set, ZIPD correspond to probabilities that geographic entities (eg, states and countries) have the same COVID-19 mortality dynamics.

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.

When Everything Becomes Bigger: Big Data for Big Poultry Production

In future decades, the demand for poultry meat and eggs is predicted to considerably increase in pace with human population growth. Although this expansion clearly represents a remarkable opportunity for the sector, it conceals a multitude of challenges. Pollution and land erosion, competition for limited resources between animal and human nutrition, animal welfare concerns, limitations on the use of growth promoters and antimicrobial agents, and increasing risks and effects of animal infectious diseases and zoonoses are several topics that have received attention from authorities and the public. The increase in poultry production must be achieved mainly through optimization and increased efficiency. The increasing ability to generate large amounts of data ("big data") is pervasive in both modern society and the farming industry. Information accessibility-coupled with the availability of tools and computational power to store, share, integrate, and analyze data with automatic and flexible algorithms-offers an unprecedented opportunity to develop tools to maximize farm profitability, reduce socio-environmental impacts, and increase animal and human health and welfare. A detailed description of all topics and applications of big data analysis in poultry farming would be infeasible. Therefore, the present work briefly reviews the application of sensor technologies, such as optical, acoustic, and wearable sensors, as well as infrared thermal imaging and optical flow, to poultry farming. The principles and benefits of advanced statistical techniques, such as machine learning and deep learning, and their use in developing effective and reliable classification and prediction models to benefit the farming system, are also discussed. Finally, recent progress in pathogen genome sequencing and analysis is discussed, highlighting practical applications in epidemiological tracking, and reconstruction of microorganisms' population dynamics, evolution, and spread. The benefits of the objective evaluation of the effectiveness of applied control strategies are also considered. Although human-artificial intelligence collaborations in the livestock sector can be frightening because they require farmers and employees in the sector to adapt to new roles, challenges, and competencies-and because several unknowns, limitations, and open-ended questions are inevitable-their overall benefits appear to be far greater than their drawbacks. As more farms and companies connect to technology, artificial intelligence (AI) and sensing technologies will begin to play a greater role in identifying patterns and solutions to pressing problems in modern animal farming, thus providing remarkable production-based and commercial advantages. Moreover, the combination of diverse sources and types of data will also become fundamental for the development of predictive models able to anticipate, rather than merely detect, disease occurrence. The increasing availability of sensors, infrastructures, and tools for big data collection, storage, sharing, and analysis-together with the use of open standards and integration with pathogen molecular epidemiology-have the potential to address the major challenge of producing higher-quality, more healthful food on a larger scale in a more sustainable manner, thereby protecting ecosystems, preserving natural resources, and improving animal and human welfare and health.

Unraveling the epidemiology of Mycobacterium bovis using whole-genome sequencing combined with environmental and demographic data

When studying the dynamics of a pathogen in a host population, one crucial question is whether it transitioned from an epidemic (i.e., the pathogen population and the number of infected hosts are increasing) to an endemic stable state (i.e., the pathogen population reached an equilibrium). For slow-growing and slow-evolving clonal pathogens such as Mycobacterium bovis, the causative agent of bovine (or animal) and zoonotic tuberculosis, it can be challenging to discriminate between these two states. This is a result of the combination of suboptimal detection tests so that the actual extent of the pathogen prevalence is often unknown, as well as of the low genetic diversity, which can hide the temporal signal provided by the accumulation of mutations in the bacterial DNA. In recent years, the increased availability, efficiency, and reliability of genomic reading techniques, such as whole-genome sequencing (WGS), have significantly increased the amount of information we can use to study infectious diseases, and therefore, it has improved the precision of epidemiological inferences for pathogens such as M. bovis. In this study, we use WGS to gain insights into the epidemiology of M. bovis in Cameroon, a developing country where the pathogen has been reported for decades. A total of 91 high-quality sequences were obtained from tissue samples collected in four abattoirs, 64 of which were with complete metadata. We combined these with environmental, demographic, ecological, and cattle movement data to generate inferences using phylodynamic models. Our findings suggest M. bovis in Cameroon is slowly expanding its epidemiological range over time; therefore, endemic stability is unlikely. This suggests that animal movement plays an important role in transmission. The simultaneous prevalence of M. bovis in co-located cattle and humans highlights the risk of such transmission being zoonotic. Therefore, using genomic tools as part of surveillance would vastly improve our understanding of disease ecology and control strategies.

Phylogeographic reconstruction of the emergence and spread of Powassan virus in the northeastern United States

Powassan virus is an emerging tick-borne virus of concern for public health, but very little is known about its transmission patterns and ecology. Here, we expanded the genomic dataset by sequencing 279 Powassan viruses isolated from Ixodes scapularis ticks from the northeastern United States. Our phylogeographic reconstructions revealed that Powassan virus lineage II was likely introduced or emerged from a relict population in the Northeast between 1940 and 1975. Sequences strongly clustered by sampling location, suggesting a highly focal geographical distribution. Our analyses further indicated that Powassan virus lineage II emerged in the northeastern United States mostly following a south-to-north pattern, with a weighted lineage dispersal velocity of ~3 km/y. Since the emergence in the Northeast, we found an overall increase in the effective population size of Powassan virus lineage II, but with growth stagnating during recent years. The cascading effect of population expansion of white-tailed deer and I. scapularis populations likely facilitated the emergence of Powassan virus in the northeastern United States.

Effects of culling vampire bats on the spatial spread and spillover of rabies virus

Controlling pathogen circulation in wildlife reservoirs is notoriously challenging. In Latin America, vampire bats have been culled for decades in hopes of mitigating lethal rabies infections in humans and livestock. Whether culls reduce or exacerbate rabies transmission remains controversial. Using Bayesian state-space models, we show that a 2-year, spatially extensive bat cull in an area of exceptional rabies incidence in Peru failed to reduce spillover to livestock, despite reducing bat population density. Viral whole genome sequencing and phylogeographic analyses further demonstrated that culling before virus arrival slowed viral spatial spread, but reactive culling accelerated spread, suggesting that culling-induced changes in bat dispersal promoted viral invasions. Our findings question the core assumptions of density-dependent transmission and localized viral maintenance that underlie culling bats as a rabies prevention strategy and provide an epidemiological and evolutionary framework to understand the outcomes of interventions in complex wildlife disease systems.

Early Genomic Surveillance and Phylogeographic Analysis of Getah Virus, a Reemerging Arbovirus, in Livestock in China

Getah virus (GETV) mainly causes disease in livestock and may pose an epidemic risk due to its expanding host range and the potential of long-distance dispersal through animal trade. Here, we used metagenomic next-generation sequencing (mNGS) to identify GETV as the pathogen responsible for reemerging swine disease in China and subsequently estimated key epidemiological parameters using phylodynamic and spatially-explicit phylogeographic approaches. The GETV isolates were able to replicate in a variety of cell lines, including human cells, and showed high pathogenicity in a mouse model, suggesting the potential for more mammal hosts. We obtained 16 complete genomes and 79 E2 gene sequences from viral strains collected in China from 2016 to 2021 through large-scale surveillance among livestock, pets, and mosquitoes. Our phylogenetic analysis revealed that three major GETV lineages are responsible for the current epidemic in livestock in China. We identified three potential positively selected sites and mutations of interest in E2, which may impact the transmissibility and pathogenicity of the virus. Phylodynamic inference of the GETV demographic dynamics identified an association between livestock meat consumption and the evolution of viral genetic diversity. Finally, phylogeographic reconstruction of GETV dispersal indicated that the sampled lineages have preferentially circulated within areas associated with relatively higher mean annual temperature and pig population density. Our results highlight the importance of continuous surveillance of GETV among livestock in southern Chinese regions associated with relatively high temperatures. IMPORTANCE Although livestock is known to be the primary reservoir of Getah virus (GETV) in Asian countries, where identification is largely based on serology, the evolutionary history and spatial epidemiology of GETV in these regions remain largely unknown. Through our sequencing efforts, we provided robust support for lineage delineation of GETV and identified three major lineages that are responsible for the current epidemic in livestock in China. We further analyzed genomic and epidemiological data to reconstruct the recent demographic and dispersal history of GETV in domestic animals in China and to explore the impact of environmental factors on its genetic diversity and its diffusion. Notably, except for livestock meat consumption, other pig-related factors such as the evolution of live pig transport and pork production do not show a significant association with the evolution of viral genetic diversity, pointing out that further studies should investigate the potential contribution of other host species to the GETV outbreak. Our analysis of GETV demonstrates the need for wider animal species surveillance and provides a baseline for future studies of the molecular epidemiology and early warning of emerging arboviruses in China.

Phylogeographic analysis of Tula hantavirus highlights a single introduction to central Europe

Orthohantaviruses are zoonotic pathogens of humans, unique among the bunyaviruses in not being transmitted by an arthropod vector. Tula orthohantavirus (TULV) is an old-world hantavirus, of yet unclear human pathogenicity, with few reported cases of clinically relevant human infection. So far, phylogeographic studies exploring the global pathways of hantaviral migration are scarce and generally do not focus on a specific hantavirus species. The aim of the present study was to reconstruct the dispersal history of TULV lineages across Eurasia based on S segment sequences sampled from different geographic areas. Maximum-likelihood and Bayesian inference methods were used to perform the phylogenetic analysis and phylogeographic reconstructions. Sampling time and trapping localities were obtained for a total of 735 TULV S segment sequences available in public databases at the time of the study. The estimated substitution rate of the analyzed partial S segment alignment was 2.26 × 10-3 substitutions/site/year (95 per cent highest posterior density interval: 1.79 × 10-3 to 2.75 × 10-3). Continuous phylogeography of TULV S segment sequences placed the potential root and origin of TULV spread in the Black Sea region. In our study, we detect a single-lineage introduction of TULV to Europe, followed by local viral circulation further on.

Spatial and temporal dynamic analysis of rabies: A review of current methodologies

Rabies continues to be one of the deadliest, high risk diseases worldwide, posing a severe threat to public health. The lack of human-to-human transmission means that the spread of rabies is not significantly affected by the distribution of humans or migra- tion. Thus, the spatiotemporal dynamic of cases in both wild and domestic animals is an important issue that can result in human cases. This paper gives an overview of the methodologies for the spatial and temporal dynamic analysis of this disease. It introduces the most representative research progress of spatial aggregation, dynamic transmission, spatiotemporal distribution, epidemiologi- cal analysis and application of modelling in the study of rabies transmission in recent years. This overview should be useful for investigating the spatial and temporal dynamics of rabies, as it could help understanding the spread of cases as well as contribute to the development of better prevention and control strategies in ecology and epidemiology.

Predicting the evolution of the Lassa virus endemic area and population at risk over the next decades

Lassa fever is a severe viral hemorrhagic fever caused by a zoonotic virus that repeatedly spills over to humans from its rodent reservoirs. It is currently not known how climate and land use changes could affect the endemic area of this virus, currently limited to parts of West Africa. By exploring the environmental data associated with virus occurrence using ecological niche modelling, we show how temperature, precipitation and the presence of pastures determine ecological suitability for virus circulation. Based on projections of climate, land use, and population changes, we find that regions in Central and East Africa will likely become suitable for Lassa virus over the next decades and estimate that the total population living in ecological conditions that are suitable for Lassa virus circulation may drastically increase by 2070. By analysing geotagged viral genomes using spatially-explicit phylogeography and simulating virus dispersal, we find that in the event of Lassa virus being introduced into a new suitable region, its spread might remain spatially limited over the first decades.

Clustering and mapping the first COVID-19 outbreak in France

BACKGROUND

With more than 160 000 confirmed COVID-19 cases and about 30 000 deceased people at the end of June 2020, France was one of the countries most affected by the coronavirus crisis worldwide. We aim to assess the efficiency of global lockdown policy in limiting spatial contamination through an in-depth reanalysis of spatial statistics in France during the first lockdown and immediate post-lockdown phases.

METHODS

To reach that goal, we use an integrated approach at the crossroads of geography, spatial epidemiology, and public health science. To eliminate any ambiguity relevant to the scope of the study, attention focused at first on data quality assessment. The data used originate from official databases (Santé Publique France) and the analysis is performed at a departmental level. We then developed spatial autocorrelation analysis, thematic mapping, hot spot analysis, and multivariate clustering.

RESULTS

We observe the extreme heterogeneity of local situations and demonstrate that clustering and intensity are decorrelated indicators. Thematic mapping allows us to identify five "ghost" clusters, whereas hot spot analysis detects two positive and two negative clusters. Our re-evaluation also highlights that spatial dissemination follows a twofold logic, zonal contiguity and linear development, thus determining a "metastatic" propagation pattern.

CONCLUSIONS

One of the most problematic issues about COVID-19 management by the authorities is the limited capacity to identify hot spots. Clustering of epidemic events is often biased because of inappropriate data quality assessment and algorithms eliminating statistical-spatial outliers. Enhanced detection techniques allow for a better identification of hot and cold spots, which may lead to more effective political decisions during epidemic outbreaks.

Genetic diversity and connectivity of the Ostreid herpesvirus 1 populations in France: A first attempt to phylogeographic inference for a marine mollusc disease

The genetic diversity of viral populations is a key driver of the spatial and temporal diffusion of viruses; yet, studying the diversity of whole genomes from natural populations still remains a challenge. Phylodynamic approaches are commonly used for RNA viruses harboring small genomes but have only rarely been applied to DNA viruses with larger genomes. Here, we used the Pacific oyster mortality syndrome (a disease that affects oyster farms around the world) as a model to study the genetic diversity of its causative agent, the Ostreid herpesvirus 1 (OsHV-1) in the three main French oyster-farming areas. Using ultra-deep sequencing on individual moribund oysters and an innovative combination of bioinformatics tools, we de novo assembled twenty-one OsHV-1 new genomes. Combining quantification of major and minor genetic variations, phylogenetic analysis, and ancestral state reconstruction of discrete traits approaches, we assessed the connectivity of OsHV-1 viral populations between the three oyster-farming areas. Our results suggest that the Marennes-Oléron Bay represents the main source of OsHV-1 diversity, from where the virus has dispersed to other farming areas, a scenario consistent with current practices of oyster transfers in France. We demonstrate that phylodynamic approaches can be applied to aquatic DNA viruses to determine how epidemiological, immunological, and evolutionary processes act and potentially interact to shape their diversity patterns.

OUP accepted manuscript

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) evolution plays a significant role in shaping the dynamics of the coronavirus disease 2019 pandemic. To monitor the evolution of SARS-CoV-2 variants, through international collaborations, we performed genomic epidemiology analyses on a weekly basis with SARS-CoV-2 samples collected from a border region between Germany, Poland, and the Czech Republic in a global background. For identified virus mutant variants, active viruses were isolated and functional evaluations were performed to test their replication fitness and neutralization sensitivity against vaccine-elicited serum neutralizing antibodies. Thereby we identified a new B.1.1.7 sub-lineage carrying additional mutations of nucleoprotein G204P and open-reading-frame-8 K68stop. Of note, this B.1.1.7 sub-lineage is the predominant B.1.1.7 variant in several European countries such as Czech Republic, Austria, and Slovakia. The earliest samples belonging to this sub-lineage were detected in November 2020 in a few countries in the European continent, but not in the UK. We have also detected its further evolution with extra spike mutations D138Y and A701V, which are signature mutations shared with the Gamma and Beta variants, respectively. Antibody neutralization assay of virus variant isolations has revealed that the variant with extra spike mutations is 3.2-fold less sensitive to vaccine-elicited antibodies as compared to the other B.1.1.7 variants tested, indicating potential for immune evasion, but it also exhibited reduced replication fitness, suggesting lower transmissibility. The wide spread of this B.1.1.7 sub-lineage was related to the pandemic waves in early 2021 in various European countries. These findings about the emergence, spread, evolution, infection, and transmission abilities of this B.1.1.7 sub-lineage add to our understanding about the pandemic development in Europe and highlight the importance of international collaboration on virus mutant surveillance.

Porcine Reproductive and Respiratory Syndrome (PRRS) Epidemiology in an Integrated Pig Company of Northern Italy: A Multilevel Threat Requiring Multilevel Interventions

Porcine reproductive and respiratory syndrome (PRRS) is probably the most relevant viral disease affecting pig farming. Despite the remarkable efforts paid in terms of vaccination administration and biosecurity, eradication and long-term control have often been frustrated. Unfortunately, few studies are currently available that objectively link, using a formal statistical approach, viral molecular epidemiology to the risk factors determining the observed scenario. The purpose of the present study is to contribute to filling this knowledge gap taking advantage of the advancements in the field of phylodynamics. Approximately one-thousand ORF7 sequences were obtained from strains collected between 2004 and 2021 from the largest Italian pig company, which implements strict compartmentalization among independent three-sites (i.e., sow herds, nurseries and finishing units) pig flows. The history and dynamics of the viral population and its evolution over time were reconstructed and linked to managerial choices. The viral fluxes within and among independent pig flows were evaluated, and the contribution of other integrated pig companies and rurally risen pigs in mediating such spreading was investigated. Moreover, viral circulation in Northern Italy was reconstructed using a continuous phylogeographic approach, and the impact of several environmental features on PRRSV strain persistence and spreading velocity was assessed. The results demonstrate that PRRSV epidemiology is shaped by a multitude of factors, including pig herd management (e.g., immunization strategy), implementation of strict-independent pig flows, and environmental features (e.g., climate, altitude, pig density, road density, etc.) among the others. Small farms and rurally raised animals also emerged as a potential threat for larger, integrated companies. These pieces of evidence suggest that none of the implemented measures can be considered effective alone, and a multidimensional approach, ranging from individual herd management to collaboration and information sharing among different companies, is mandatory for effective infection control.

Exploiting genomic surveillance to map the spatio-temporal dispersal of SARS-CoV-2 spike mutations in Belgium across 2020

At the end of 2020, several new variants of SARS-CoV-2—designated variants of concern—were detected and quickly suspected to be associated with a higher transmissibility and possible escape of vaccine-induced immunity. In Belgium, this discovery has motivated the initiation of a more ambitious genomic surveillance program, which is drastically increasing the number of SARS-CoV-2 genomes to analyse for monitoring the circulation of viral lineages and variants of concern. In order to efficiently analyse the massive collection of genomic data that are the result of such increased sequencing efforts, streamlined analytical strategies are crucial. In this study, we illustrate how to efficiently map the spatio-temporal dispersal of target mutations at a regional level. As a proof of concept, we focus on the Belgian province of Liège that has been consistently sampled throughout 2020, but was also one of the main epicenters of the second European epidemic wave. Specifically, we employ a recently developed phylogeographic workflow to infer the regional dispersal history of viral lineages associated with three specific mutations on the spike protein (S98F, A222V and S477N) and to quantify their relative importance through time. Our analytical pipeline enables analysing large data sets and has the potential to be quickly applied and updated to track target mutations in space and time throughout the course of an epidemic.

A Review of Environmental Risk Factors for African Swine Fever in European Wild Boar

A detailed understanding of environmental risk factors for African swine fever (ASF) in wild boar will be not only essential for risk assessments but also for timely and spatially informed allocation of resources in order to manage wild boar-targeted ASF control measures efficiently. Here, we review currently known environmental risk factors that can influence the occurrence of ASF virus infection in wild boar when compared to disease occurrence in wild boar of a non-exposed reference scenario. Accordingly, the exposure of wild boar to environmental risk factors related to (1) climate, (2) land cover, (3) human activity, (4) wild boar and (5) ASF were evaluated. As key environmental risk factors in this review, increased ASF occurrence in wild boar was associated with seasonal patterns, forest coverage, presence of water, human presence, farming activities, wild boar density and ASF nearness. The review highlights inconsistencies in some of these risk factor associations with disease detection in space and time and may provide valuable insights for the investigation of ASF transmission dynamics. The examined risk information was applied to consider potential improvements of the ASF control strategy in wild boar regarding disease surveillance, hunting, wild boar carcass searches and ASF barrier implementation.

On the Use of Phylogeographic Inference to Infer the Dispersal History of Rabies Virus: A Review Study

Rabies is a neglected zoonotic disease which is caused by negative strand RNA-viruses belonging to the genus Lyssavirus. Within this genus, rabies viruses circulate in a diverse set of mammalian reservoir hosts, is present worldwide, and is almost always fatal in non-vaccinated humans. Approximately 59,000 people are still estimated to die from rabies each year, leading to a global initiative to work towards the goal of zero human deaths from dog-mediated rabies by 2030, requiring scientific efforts from different research fields. The past decade has seen a much increased use of phylogeographic and phylodynamic analyses to study the evolution and spread of rabies virus. We here review published studies in these research areas, making a distinction between the geographic resolution associated with the available sequence data. We pay special attention to environmental factors that these studies found to be relevant to the spread of rabies virus. Importantly, we highlight a knowledge gap in terms of applying these methods when all required data were available but not fully exploited. We conclude with an overview of recent methodological developments that have yet to be applied in phylogeographic and phylodynamic analyses of rabies virus.

Rivers and landscape ecology of a plant virus, Rice yellow mottle virus along the Niger Valley

To investigate the spread of Rice yellow mottle virus (RYMV) along the Niger River, regular sampling of virus isolates was conducted along 500 km of the Niger Valley in the Republic of Niger and was complemented by additional sampling in neighbouring countries in West Africa and Central Africa. The spread of RYMV into and within the Republic of Niger was inferred as a continuous process using a Bayesian statistical framework applied previously to reconstruct its dispersal history in West Africa, East Africa, and Madagascar. The spatial resolution along this section of the Niger River was the highest implemented for RYMV and possibly for any plant virus. We benefited from the results of early field surveys of the disease for the validation of the phylogeographic reconstruction and from the well-documented history of rice cultivation changes along the Niger River for their interpretation. As a prerequisite, the temporal signal of the RYMV data sets was revisited in the light of recent methodological advances. The role of the hydrographic network of the Niger Basin in RYMV spread was examined, and the link between virus population dynamics and the extent of irrigated rice was assessed. RYMV was introduced along the Niger River in the Republic of Niger in the early 1980s from areas to the southwest of the country where rice was increasingly grown. Viral spread was triggered by a major irrigation scheme made of a set of rice perimeters along the river valley. The subsequent spatial and temporal host continuity and the inoculum build-up allowed for a rapid spread of RYMV along the Niger River, upstream and downstream, over hundreds of kilometres, and led to the development of severe epidemics. There was no evidence of long-distance dissemination of the virus through natural water. Floating rice in the main meanders of the Middle Niger did not contribute to virus dispersal from West Africa to Central Africa. RYMV along the Niger River is an insightful example of how agricultural intensification favours pathogen emergence and spread.

Relax, Keep Walking - A Practical Guide to Continuous Phylogeographic Inference with BEAST

Spatially explicit phylogeographic analyses can be performed with an inference framework that employs relaxed random walks to reconstruct phylogenetic dispersal histories in continuous space. This core model was first implemented 10 years ago and has opened up new opportunities in the field of phylodynamics, allowing researchers to map and analyze the spatial dissemination of rapidly evolving pathogens. We here provide a detailed and step-by-step guide on how to set up, run, and interpret continuous phylogeographic analyses using the programs BEAUti, BEAST, Tracer, and TreeAnnotator.

Towards a more healthy conservation paradigm: integrating disease and molecular ecology to aid biological conservation†

Parasites, and the diseases they cause, are important from an ecological and evolutionary perspective because they can negatively affect host fitness and can regulate host populations. Consequently, conservation biology has long recognized the vital role that parasites can play in the process of species endangerment and recovery. However, we are only beginning to understand how deeply parasites are embedded in ecological systems, and there is a growing recognition of the important ways in which parasites affect ecosystem structure and function. Thus, there is an urgent need to revisit how parasites are viewed from a conservation perspective and broaden the role that disease ecology plays in conservation-related research and outcomes. This review broadly focusses on the role that disease ecology can play in biological conservation. Our review specifically emphasizes on how the integration of tools and analytical approaches associated with both disease and molecular ecology can be leveraged to aid conservation biology. Our review first concentrates on disease-mediated extinctions and wildlife epidemics. We then focus on elucidating how host–parasite interactions has improved our understanding of the eco-evolutionary dynamics affecting hosts at the individual, population, community and ecosystem scales. We believe that the role of parasites as drivers and indicators of ecosystem health is especially an exciting area of research that has the potential to fundamentally alter our view of parasites and their role in biological conservation. The review concludes with a broad overview of the current and potential applications of modern genomic tools in disease ecology to aid biological conservation.

Landscape epidemiology of Batrachochytrium salamandrivorans: reconciling data limitations and conservation urgency

Starting in 2010, rapid fire salamander (Salamandra salamandra) population declines in northwestern Europe heralded the emergence of Batrachochytrium salamandrivorans (Bsal), a salamander-pathogenic chytrid fungus. Bsal poses an imminent threat to global salamander diversity owing to its wide host range, high pathogenicity, and long-term persistence in ecosystems. While there is a pressing need to develop further research and conservation actions, data limitations inherent to recent pathogen emergence obscure necessary insights into Bsal disease ecology. Here, we use a hierarchical modeling framework to describe Bsal landscape epidemiology of outbreak sites in light of these methodological challenges. Using model selection and machine learning, we find that Bsal presence is associated with humid and relatively cool, stable climates. Outbreaks are generally located in areas characterized by low landscape heterogeneity and low steepness of slope. We further find an association between Bsal presence and high trail density, suggesting that human-mediated spread may increase risk for spillover between populations. We then use distribution modeling to show that favorable conditions occur in lowlands influenced by the North Sea, where increased survey effort is needed to determine how Bsal impacts local newt populations, but also in hill- and mountain ranges in northeastern France and the lower half of Germany. Finally, connectivity analyses suggest that these hill- and mountain ranges may act as stepping stones for further spread southward. Our results provide initial insight into regional environmental conditions underlying Bsal epizootics, present updated invasibility predictions for northwestern Europe, and lead us to discuss a wide variety of potential survey and research actions needed to advance future conservation and mitigation efforts.

The roles of environmental variation and parasite survival in virulence-transmission relationships

Disease outbreaks are a consequence of interactions among the three components of a host-parasite system: the infectious agent, the host and the environment. While virulence and transmission are widely investigated, most studies of parasite life-history trade-offs are conducted with theoretical models or tractable experimental systems where transmission is standardized and the environment controlled. Yet, biotic and abiotic environmental factors can strongly affect disease dynamics, and ultimately, host-parasite coevolution. Here, we review research on how environmental context alters virulence-transmission relationships, focusing on the off-host portion of the parasite life cycle, and how variation in parasite survival affects the evolution of virulence and transmission. We review three inter-related 'approaches' that have dominated the study of the evolution of virulence and transmission for different host-parasite systems: (i) evolutionary trade-off theory, (ii) parasite local adaptation and (iii) parasite phylodynamics. These approaches consider the role of the environment in virulence and transmission evolution from different angles, which entail different advantages and potential biases. We suggest improvements to how to investigate virulence-transmission relationships, through conceptual and methodological developments and taking environmental context into consideration. By combining developments in life-history evolution, phylogenetics, adaptive dynamics and comparative genomics, we can improve our understanding of virulence-transmission relationships across a diversity of host-parasite systems that have eluded experimental study of parasite life history.

Development of a Dissemination Platform for Spatiotemporal and Phylogenetic Analysis of Avian Infectious Bronchitis Virus

Infecting large portions of the global poultry populations, the avian infectious bronchitis virus (IBV) remains a major economic burden in North America. With more than 30 serotypes globally distributed, Arkansas, Connecticut, Delaware, Georgia, and Massachusetts are among the most predominant serotypes in the United States. Even though vaccination is widely used, the high mutation rate exhibited by IBV is continuously triggering the emergence of new viral strains and hindering control and prevention measures. For that reason, targeted strategies based on constantly updated information on the IBV circulation are necessary. Here, we sampled IBV-infected farms from one US state and collected and analyzed 65 genetic sequences coming from three different lineages along with the immunization information of each sampled farm. Phylodynamic analyses showed that IBV dispersal velocity was 12.3 km/year. The majority of IBV infections appeared to have derived from the introduction of the Arkansas DPI serotype, and the Arkansas DPI and Georgia 13 were the predominant serotypes. When analyzed against IBV sequences collected across the United States and deposited in the GenBank database, the most likely viral origin of our sequences was from the states of Alabama, Georgia, and Delaware. Information about vaccination showed that the MILDVAC-MASS+ARK vaccine was applied on 26% of the farms. Using a publicly accessible open-source tool for real-time interactive tracking of pathogen spread and evolution, we analyzed the spatiotemporal spread of IBV and developed an online reporting dashboard. Overall, our work demonstrates how the combination of genetic and spatial information could be used to track the spread and evolution of poultry diseases, providing timely information to the industry. Our results could allow producers and veterinarians to monitor in near-real time the current IBV strain circulating, making it more informative, for example, in vaccination-related decisions.

Mathematical modelling and phylodynamics for the study of dog rabies dynamics and control: A scoping review

BACKGROUND

Rabies is a fatal yet vaccine-preventable disease. In the last two decades, domestic dog populations have been shown to constitute the predominant reservoir of rabies in developing countries, causing 99% of human rabies cases. Despite substantial control efforts, dog rabies is still widely endemic and is spreading across previously rabies-free areas. Developing a detailed understanding of dog rabies dynamics and the impact of vaccination is essential to optimize existing control strategies and developing new ones. In this scoping review, we aimed at disentangling the respective contributions of mathematical models and phylodynamic approaches to advancing the understanding of rabies dynamics and control in domestic dog populations. We also addressed the methodological limitations of both approaches and the remaining issues related to studying rabies spread and how this could be applied to rabies control.

METHODOLOGY/PRINCIPAL FINDINGS

We reviewed how mathematical modelling of disease dynamics and phylodynamics have been developed and used to characterize dog rabies dynamics and control. Through a detailed search of the PubMed, Web of Science, and Scopus databases, we identified a total of n = 59 relevant studies using mathematical models (n = 30), phylodynamic inference (n = 22) and interdisciplinary approaches (n = 7). We found that despite often relying on scarce rabies epidemiological data, mathematical models investigated multiple aspects of rabies dynamics and control. These models confirmed the overwhelming efficacy of massive dog vaccination campaigns in all settings and unraveled the role of dog population structure and frequent introductions in dog rabies maintenance. Phylodynamic approaches successfully disentangled the evolutionary and environmental determinants of rabies dispersal and consistently reported support for the role of reintroduction events and human-mediated transportation over long distances in the maintenance of rabies in endemic areas. Potential biases in data collection still need to be properly accounted for in most of these analyses. Finally, interdisciplinary studies were determined to provide the most comprehensive assessments through hypothesis generation and testing. They also represent new avenues, especially concerning the reconstruction of local transmission chains or clusters through data integration.

CONCLUSIONS/SIGNIFICANCE

Despite advances in rabies knowledge, substantial uncertainty remains regarding the mechanisms of local spread, the role of wildlife in dog rabies maintenance, and the impact of community behavior on the efficacy of control strategies including vaccination of dogs. Future integrative approaches that use phylodynamic analyses and mechanistic models within a single framework could take full advantage of not only viral sequences but also additional epidemiological information as well as dog ecology data to refine our understanding of rabies spread and control. This would represent a significant improvement on past studies and a promising opportunity for canine rabies research in the frame of the One Health concept that aims to achieve better public health outcomes through cross-sector collaboration.

Metapopulation Structure Predicts Population Dynamics in the Cakile maritima-Alternaria brassicicola Host-Pathogen Interaction

AbstractIn symbiotic interactions, spatiotemporal variation in the distribution or population dynamics of one species represents spatial and temporal heterogeneity of the landscape for the other. Such interdependent demographic dynamics result in situations where the relative importance of biotic and abiotic factors in determining ecological processes is complicated to decipher. Using a detailed survey of three metapopulations of the succulent plant Cakile maritima and the necrotrophic fungus Alternaria brassicicola located along the southeastern Australian coast, we developed a series of statistical analyses-namely, synchrony analysis, patch occupancy dynamics, and a spatially explicit metapopulation model-to understand how habitat quality, weather conditions, dispersal, and spatial structure determine metapopulation dynamics. Climatic conditions are important drivers, likely explaining the high synchrony among populations. Host availability, landscape features facilitating dispersal, and habitat conditions also impact the occurrence and spread of disease. Overall, we show that the collection of extensive data on host and pathogen population dynamics, in combination with spatially explicit epidemiological modeling, makes it possible to accurately predict disease dynamics-even when there is extreme variability in host population dynamics. Finally, we discuss the importance of genetic information for predicting demographic dynamics in this pathosystem.

Host relatedness and landscape connectivity shape pathogen spread in the puma, a large secretive carnivore

Urban expansion can fundamentally alter wildlife movement and gene flow, but how urbanization alters pathogen spread is poorly understood. Here, we combine high resolution host and viral genomic data with landscape variables to examine the context of viral spread in puma (Puma concolor) from two contrasting regions: one bounded by the wildland urban interface (WUI) and one unbounded with minimal anthropogenic development (UB). We found landscape variables and host gene flow explained significant amounts of variation of feline immunodeficiency virus (FIV) spread in the WUI, but not in the unbounded region. The most important predictors of viral spread also differed; host spatial proximity, host relatedness, and mountain ranges played a role in FIV spread in the WUI, whereas roads might have facilitated viral spread in the unbounded region. Our research demonstrates how anthropogenic landscapes can alter pathogen spread, providing a more nuanced understanding of host-pathogen relationships to inform disease ecology in free-ranging species.

Spatiotemporal reconstruction and transmission dynamics during the 2016-17 H5N8 highly pathogenic avian influenza epidemic in Italy

Effective control of avian diseases in domestic populations requires understanding of the transmission dynamics facilitating viral emergence and spread. In 2016-17, Italy experienced a significant avian influenza epidemic caused by a highly pathogenic A(H5N8) virus, which affected domestic premises housing around 2.7 million birds, primarily in the north-eastern regions with the highest density of poultry farms (Lombardy, Emilia-Romagna and Veneto). We perform integrated analyses of genetic, spatiotemporal and host data within a Bayesian phylogenetic framework. Using continuous and discrete phylogeography, we estimate the locations of movements responsible for the spread and persistence of the epidemic. The information derived from these analyses on rates of transmission between regions through time can be used to assess the success of control measures. Using an approach based on phylogenetic-temporal distances between domestic cases, we infer the presence of cryptic wild bird-mediated transmission, information that can be used to complement existing epidemiological methods for distinguishing transmission within the domestic population from incursions across the wildlife-domestic interface, a common challenge in veterinary epidemiology. Spatiotemporal reconstruction of the epidemic reveals a highly skewed distribution of virus movements with a high proportion of shorter distance local movements interspersed with occasional long-distance dispersal events associated with wild birds. We also show how such inference be used to identify possible instances of human-mediated movements where distances between phylogenetically linked domestic cases are unusually high.

Comparative Circulation Dynamics of the Five Main HIV Types in China

The HIV epidemic in China accounts for 3% of the global HIV incidence. We compared the patterns and determinants of interprovincial spread of the five most prevalent circulating types. HIV pol sequences sampled across China were used to identify relevant transmission networks of the five most relevant HIV-1 types (B and circulating recombinant forms [CRFs] CRF01_AE, CRF07_BC, CRF08_BC, and CRF55_01B) in China. From these, the dispersal history across provinces was inferred. A generalized linear model (GLM) was used to test the association between migration rates among provinces and several measures of human mobility. A total of 10,707 sequences were collected between 2004 and 2017 across 26 provinces, among which 1,962 are newly reported here. A mean of 18 (minimum and maximum, 1 and 54) independent transmission networks involving up to 17 provinces were identified. Discrete phylogeographic analysis largely recapitulates the documented spread of the HIV types, which in turn, mirrors within-China population migration flows to a large extent. In line with the different spatiotemporal spread dynamics, the identified drivers thereof were also heterogeneous but are consistent with a central role of human mobility. The comparative analysis of the dispersal dynamics of the five main HIV types circulating in China suggests a key role of large population centers and developed transportation infrastructures as hubs of HIV dispersal. This advocates for coordinated public health efforts in addition to local targeted interventions.IMPORTANCE While traditional epidemiological studies are of great interest in describing the dynamics of epidemics, they struggle to fully capture the geospatial dynamics and factors driving the dispersal of pathogens like HIV as they have difficulties capturing linkages between infections. To overcome this, we used a discrete phylogeographic approach coupled to a generalized linear model extension to characterize the dynamics and drivers of the across-province spread of the five main HIV types circulating in China. Our results indicate that large urbanized areas with dense populations and developed transportation infrastructures are facilitators of HIV dispersal throughout China and highlight the need to consider harmonized country-wide public policies to control local HIV epidemics.

Epidemiological hypothesis testing using a phylogeographic and phylodynamic framework

Computational analyses of pathogen genomes are increasingly used to unravel the dispersal history and transmission dynamics of epidemics. Here, we show how to go beyond historical reconstructions and use spatially-explicit phylogeographic and phylodynamic approaches to formally test epidemiological hypotheses. We illustrate our approach by focusing on the West Nile virus (WNV) spread in North America that has substantially impacted public, veterinary, and wildlife health. We apply an analytical workflow to a comprehensive WNV genome collection to test the impact of environmental factors on the dispersal of viral lineages and on viral population genetic diversity through time. We find that WNV lineages tend to disperse faster in areas with higher temperatures and we identify temporal variation in temperature as a main predictor of viral genetic diversity through time. By contrasting inference with simulation, we find no evidence for viral lineages to preferentially circulate within the same migratory bird flyway, suggesting a substantial role for non-migratory birds or mosquito dispersal along the longitudinal gradient.

Incorporating heterogeneous sampling probabilities in continuous phylogeographic inference - Application to H5N1 spread in the Mekong region

MOTIVATION

The potentially low precision associated with the geographic origin of sampled sequences represents an important limitation for spatially explicit (i.e. continuous) phylogeographic inference of fast-evolving pathogens such as RNA viruses. A substantial proportion of publicly available sequences is geo-referenced at broad spatial scale such as the administrative unit of origin, rather than more precise locations (e.g. geographic coordinates). Most frequently, such sequences are either discarded prior to continuous phylogeographic inference or arbitrarily assigned to the geographic coordinates of the centroid of their administrative area of origin for lack of a better alternative.

RESULTS

We here implement and describe a new approach that allows to incorporate heterogeneous prior sampling probabilities over a geographic area. External data, such as outbreak locations, are used to specify these prior sampling probabilities over a collection of sub-polygons. We apply this new method to the analysis of highly pathogenic avian influenza H5N1 clade data in the Mekong region. Our method allows to properly include, in continuous phylogeographic analyses, H5N1 sequences that are only associated with large administrative areas of origin and assign them with more accurate locations. Finally, we use continuous phylogeographic reconstructions to analyse the dispersal dynamics of different H5N1 clades and investigate the impact of environmental factors on lineage dispersal velocities.

AVAILABILITY AND IMPLEMENTATION

Our new method allowing heterogeneous sampling priors for continuous phylogeographic inference is implemented in the open-source multi-platform software package BEAST 1.10.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genomic Epidemiology, Evolution, and Transmission Dynamics of Porcine Deltacoronavirus

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has shown once again that coronavirus (CoV) in animals are potential sources for epidemics in humans. Porcine deltacoronavirus (PDCoV) is an emerging enteropathogen of swine with a worldwide distribution. Here, we implemented and described an approach to analyze the epidemiology of PDCoV following its emergence in the pig population. We performed an integrated analysis of full genome sequence data from 21 newly sequenced viruses, along with comprehensive epidemiological surveillance data collected globally over the last 15 years. We found four distinct phylogenetic lineages of PDCoV, which differ in their geographic circulation patterns. Interestingly, we identified more frequent intra- and interlineage recombination and higher virus genetic diversity in the Chinese lineages compared with the USA lineage where pigs are raised in different farming systems and ecological environments. Most recombination breakpoints are located in the ORF1ab gene rather than in genes encoding structural proteins. We also identified five amino acids under positive selection in the spike protein suggesting a role for adaptive evolution. According to structural mapping, three positively selected sites are located in the N-terminal domain of the S1 subunit, which is the most likely involved in binding to a carbohydrate receptor, whereas the other two are located in or near the fusion peptide of the S2 subunit and thus might affect membrane fusion. Finally, our phylogeographic investigations highlighted notable South-North transmission as well as frequent long-distance dispersal events in China that could implicate human-mediated transmission. Our findings provide new insights into the evolution and dispersal of PDCoV that contribute to our understanding of the critical factors involved in CoVs emergence.

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.

Porcine reproductive and respiratory syndrome virus dissemination across pig production systems in the United States

Porcine reproductive and respiratory syndrome virus (PRRSV) remains widespread in the North American pig population. Despite improvements in virus characterization, it is unclear whether PRRSV infections are a product of viral circulation within production systems (local) or across production systems (external). Here, we examined the local and external dissemination dynamics of PRRSV and the processes facilitating its spread in three production systems. Overall, PRRSV genetic diversity has declined since 2018, while phylodynamic results support frequent external transmission. We found that PRRSV dissemination predominantly occurred mostly through transmission between farms of different production companies for several months, especially from November until May, a timeframe already established as PRRSV season. Although local PRRSV dissemination occurred mainly through regular pig flow (from sow to nursery and then to finisher farms), an important flux of PRRSV dissemination also occurred in the opposite direction, from finisher to sow and nursery farms, highlighting the importance of downstream farms as sources of the virus. Our results also showed that farms with pig densities of 500 to 1,000 pig/km² and farms located at a range within 0.5 km and 0.7 km from major roads were more likely to be infected by PRRSV, whereas farms at an elevation of 41 to 61 meters and surrounded by denser vegetation were less likely to be infected, indicating their role as dissemination barriers. In conclusion, our results demonstrate that external dissemination was intense, and reinforce the importance of farm proximity on PRRSV spread. Thus, consideration of farm location, geographic characteristics and animal densities across production systems may help to forecast PRRSV collateral dissemination.

The macroecology of the COVID-19 pandemic in the Anthropocene

Severe acute respiratory syndrome coronavirus 2, the virus that causes coronavirus disease 2019 (COVID-19), has expanded rapidly throughout the world. Thus, it is important to understand how global factors linked with the functioning of the Anthropocene are responsible for the COVID-19 outbreak. We tested hypotheses that the number of COVID-19 cases, number of deaths and growth rate of recorded infections: (1) are positively associated with population density as well as (2) proportion of the human population living in urban areas as a proxies of interpersonal contact rate, (3) age of the population in a given country as an indication of that population's susceptibility to COVID-19; (4) net migration rate and (5) number of tourists as proxies of infection pressure, and negatively associated with (5) gross domestic product which is a proxy of health care quality. Data at the country level were compiled from publicly available databases and analysed with gradient boosting regression trees after controlling for confounding factors (e.g. geographic location). We found a positive association between the number of COVID-19 cases in a given country and gross domestic product, number of tourists, and geographic longitude. The number of deaths was positively associated with gross domestic product, number of tourists in a country, and geographic longitude. The effects of gross domestic product and number of tourists were non-linear, with clear thresholds above which the number of COVID-19 cases and deaths increased rapidly. The growth rate of COVID-19 cases was positively linked to the number of tourists and gross domestic product. The growth rate of COVID-19 cases was negatively associated with the mean age of the population and geographic longitude. Growth was slower in less urbanised countries. This study demonstrates that the characteristics of the human population and high mobility, but not population density, may help explain the global spread of the virus. In addition, geography, possibly via climate, may play a role in the pandemic. The unexpected positive and strong association between gross domestic product and number of cases, deaths, and growth rate suggests that COVID-19 may be a new civilisation disease affecting rich economies.

Phylodynamic analysis and evaluation of the balance between anthropic and environmental factors affecting IBV spreading among Italian poultry farms

Infectious bronchitis virus (IBV) control is mainly based on wide vaccine administration. Although effective, its efficacy is not absolute, the viral circulation is not prevented and some side effects cannot be denied. Despite this, the determinants of IBV epidemiology and the factors affecting its circulation are still largely unknown and poorly investigated. In the present study, 361 IBV QX (the most relevant field genotype in Italy) sequences were obtained between 2012 and 2016 from the two main Italian integrated poultry companies. Several biostatistical and bioinformatics approaches were used to reconstruct the history of the QX genotype in Italy and to assess the effect of different environmental, climatic and social factors on its spreading patterns. Moreover, two structured coalescent models were considered in order to investigate if an actual compartmentalization occurs between the two integrated poultry companies and the role of a third "ghost" deme, representative of minor industrial poultry companies and the rural sector. The obtained results suggest that the integration of the poultry companies is an effective barrier against IBV spreading, since the strains sampled from the two companies formed two essentially-independent clades. Remarkably, the only exceptions were represented by farms located in the high densely populated poultry area of Northern Italy. The inclusion of a third deme in the model revealed the likely role of other poultry companies and rural farms (particularly concentrated in Northern Italy) as sources of strain introduction into one of the major poultry companies, whose farms are mainly located in the high densely populated poultry area of Northern Italy. Accordingly, when the effect of different environmental and urban parameters on IBV geographic spreading was investigated, no factor seems to contribute to IBV dispersal velocity, being poultry population density the only exception. Finally, the different viral population pattern observed in the two companies over the same time period supports the pivotal role of management and control strategies on IBV epidemiology. Overall, the present study results stress the crucial relevance of human action rather than environmental factors, highlighting the direct benefits that could derive from improved management and organization of the poultry sector on a larger scale.