Bluetongue virus spread in Europe is a consequence of climatic, landscape and vertebrate host factors as revealed by phylogeographic inference

Reference Material Production and Milk Protein Concentration as Elements to Improve Bluetongue Serological Diagnosis in Bulk Tank Milk

The serological surveillance of bluetongue in bulk tank milk is an efficient and cost-effective method for the early detection of bluetongue virus incursions in unvaccinated free areas of the disease. In addition, the availability of standardized and reliable reagents and refined diagnostic procedures with high sensitivity and specificity are essential for surveillance purposes. However, no available reference materials for bluetongue virus serological surveillance in bulk tank milk exist. This study shows the production and characterization of reference material for the implementation of a commercially available bluetongue milk ELISA test in official laboratories, as well as the evaluation of a procedure to increase the sensitivity in samples with low levels of antibodies. This procedure, based on milk protein concentration, allowed us to notably increase the ELISA test's analytical sensitivity, which is useful for milk samples from farms with low within-herd prevalence or pools of bulk tank milk samples. The standardized milk reference material produced here, together with the evaluated procedure to improve analytical sensitivity, could be applied as tools to ensure an accurate diagnosis by official laboratories in bluetongue unvaccinated free areas.

How fast are viruses spreading in the wild?

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

Effects of climate change and human activities on vector-borne diseases

Vector-borne diseases are transmitted by haematophagous arthropods (for example, mosquitoes, ticks and sandflies) to humans and wild and domestic animals, with the largest burden on global public health disproportionately affecting people in tropical and subtropical areas. Because vectors are ectothermic, climate and weather alterations (for example, temperature, rainfall and humidity) can affect their reproduction, survival, geographic distribution and, consequently, ability to transmit pathogens. However, the effects of climate change on vector-borne diseases can be multifaceted and complex, sometimes with ambiguous consequences. In this Review, we discuss the potential effects of climate change, weather and other anthropogenic factors, including land use, human mobility and behaviour, as possible contributors to the redistribution of vectors and spread of vector-borne diseases worldwide.

Estimation of the economic impact of a bluetongue serotype 4 outbreak in Tunisia

Introduction

Since 1999, Tunisia has experienced multiple occurrences of Bluetongue (BT) outbreaks, leading to numerous reported cases of infection and mortality in flocks. The re-emergence of the disease in 2020 caused substantial economic losses in cattle, attributed to the incursion of serotype BTV-4.

Methods

To evaluate the economic impact of the recent BT episode, we conducted a retrospective study on outbreaks that occurred in Tunisia between August and November 2020, focusing on the impact at the owner's level and its effects on both small ruminants and cattle. A total of 234 ruminant farms (sheep, cattle, and mixed) were randomly selected across Tunisian governorates and included in the study to estimate both the direct and indirect costs of these outbreaks.

Results

Total costs were calculated as the sum of losses and expenditures resulting from the BT outbreaks. At the animal level, total losses were estimated to range between 116.280 and 207.086 TND for one infected ewe (€33.721 and 60.055). For one lactating cow, costs varied between 2,590.724 and 3,171.107 TND (€751.310 and 919.621). In cattle, exposure to BTV led to a daily unit milk yield decrease of 12.50 to 14.66 L over an average period of 5 months. Diseased sheep experienced weight loss ranging between 4 and 10 kg during the BT outbreaks. The total mean cost of the 2020 BT outbreak in Tunisian investigated farms was estimated at 1,935 million TND (million €561.15) (range: 1,489 and 2,474 million TND; 431.81 and million €717.46). The most influential costs of the total BT outbreaks were the decrease in milk yield, mortality, and veterinary treatment.

Discussion

This study gives valuable insights on the economic impact of the incursion of a new serotype of BT in a naive population in Tunisia. Considering the substantial costs incurred, it is imperative that this disease receives increased attention from stakeholders, including animal owners, veterinary services, practitioners, and decision-makers.

Mathematical modeling at the livestock-wildlife interface: scoping review of drivers of disease transmission between species

Modeling of infectious diseases at the livestock-wildlife interface is a unique subset of mathematical modeling with many innate challenges. To ascertain the characteristics of the models used in these scenarios, a scoping review of the scientific literature was conducted. Fifty-six studies qualified for inclusion. Only 14 diseases at this interface have benefited from the utility of mathematical modeling, despite a far greater number of shared diseases. The most represented species combinations were cattle and badgers (for bovine tuberculosis, 14), and pigs and wild boar [for African (8) and classical (3) swine fever, and foot-and-mouth and disease (1)]. Assessing control strategies was the overwhelming primary research objective (27), with most studies examining control strategies applied to wildlife hosts and the effect on domestic hosts (10) or both wild and domestic hosts (5). In spatially-explicit models, while livestock species can often be represented through explicit and identifiable location data (such as farm, herd, or pasture locations), wildlife locations are often inferred using habitat suitability as a proxy. Though there are innate assumptions that may not be fully accurate when using habitat suitability to represent wildlife presence, especially for wildlife the parsimony principle plays a large role in modeling diseases at this interface, where parameters are difficult to document or require a high level of data for inference. Explaining observed transmission dynamics was another common model objective, though the relative contribution of involved species to epizootic propagation was only ascertained in a few models. More direct evidence of disease spill-over, as can be obtained through genomic approaches based on pathogen sequences, could be a useful complement to further inform such modeling. As computational and programmatic capabilities advance, the resolution of the models and data used in these models will likely be able to increase as well, with a potential goal being the linking of modern complex ecological models with the depth of dynamics responsible for pathogen transmission. Controlling diseases at this interface is a critical step toward improving both livestock and wildlife health, and mechanistic models are becoming increasingly used to explore the strategies needed to confront these diseases.

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.

Widespread Reassortment Contributes to Antigenic Shift in Bluetongue Viruses from South Africa

Bluetongue (BT), a viral disease of ruminants, is endemic throughout South Africa, where outbreaks of different serotypes occur. The predominant serotypes can differ annually due to herd immunity provided by annual vaccinations using a live attenuated vaccine (LAV). This has led to both wild-type and vaccine strains co-circulating in the field, potentially leading to novel viral strains due to reassortment and recombination. Little is known about the molecular evolution of the virus in the field in South Africa. The purpose of this study was to investigate the genetic diversity of field strains of BTV in South Africa and to provide an initial assessment of the evolutionary processes shaping BTV genetic diversity in the field. Complete genomes of 35 field viruses belonging to 11 serotypes, collected from different regions of the country between 2011 and 2017, were sequenced. The sequences were phylogenetically analysed in relation to all the BTV sequences available from GenBank, including the LAVs and reference strains, resulting in the analyses and reassortment detection of 305 BTVs. Phylogenomic analysis indicated a geographical selection of the genome segments, irrespective of the serotype. Based on the initial assessment of the current genomic clades that circulate in South Africa, the selection for specific clades is prevalent in directing genome segment reassortment, which seems to exclude the vaccine strains and in multiple cases involves Segment-2 resulting in antigenic shift.

VectorNet: collaborative mapping of arthropod disease vectors in Europe and surrounding areas since 2010

BackgroundArthropod vectors such as ticks, mosquitoes, sandflies and biting midges are of public and veterinary health significance because of the pathogens they can transmit. Understanding their distributions is a key means of assessing risk. VectorNet maps their distribution in the EU and surrounding areas.AimWe aim to describe the methodology underlying VectorNet maps, encourage standardisation and evaluate output.Methods: Vector distribution and surveillance activity data have been collected since 2010 from a combination of literature searches, field-survey data by entomologist volunteers via a network facilitated for each participating country and expert validation. Data were collated by VectorNet members and extensively validated during data entry and mapping processes.ResultsAs of 2021, the VectorNet archive consisted of ca 475,000 records relating to > 330 species. Maps for 42 species are routinely produced online at subnational administrative unit resolution. On VectorNet maps, there are relatively few areas where surveillance has been recorded but there are no distribution data. Comparison with other continental databases, namely the Global Biodiversity Information Facility and VectorBase show that VectorNet has 5-10 times as many records overall, although three species are better represented in the other databases. In addition, VectorNet maps show where species are absent. VectorNet's impact as assessed by citations (ca 60 per year) and web statistics (58,000 views) is substantial and its maps are widely used as reference material by professionals and the public.ConclusionVectorNet maps are the pre-eminent source of rigorously validated arthropod vector maps for Europe and its surrounding areas.

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.

An integrated eco-evolutionary framework to predict population-level responses of climate-sensitive pathogens

It is critical to gain insight into how climate change impacts evolutionary responses within climate-sensitive pathogen populations, such as increased resilience, opportunistic responses and the emergence of dominant variants from highly variable genomic backgrounds and subsequent global dispersal. This review proposes a framework to support such analysis, by combining genomic evolutionary analysis with climate time-series data in a novel spatiotemporal dataframe for use within machine learning applications, to understand past and future evolutionary pathogen responses to climate change. Recommendations are presented to increase the feasibility of interdisciplinary applications, including the importance of robust spatiotemporal metadata accompanying genome submission to databases. Such workflows will inform accessible public health tools and early-warning systems, to aid decision-making and mitigate future human health threats.

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.

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.

Risk-based serological survey of bluetongue and the first evidence of bluetongue virus serotype 26 circulation in Tunisia

BACKGROUND

Bluetongue (BT), a vector-borne disease of wild and domestic ruminants, is responsible for severe economic losses in flocks. To reduce this impact, a surveillance and control plan was implemented in Tunisia. However, the epidemiological situation of BT remains incompletely understood, especially for the circulating serotypes.

OBJECTIVE

The aim of this survey was to determine the seroprevalence, to identify the circulating serotypes and to identify the associated risk factors for bluetongue virus (BTV) circulation in Tunisia using risk-based sampling (RBS).

METHODS

A total of 3314 blood samples were randomly collected from 67 sectors using risk-based sampling and screened by competitive enzyme-linked immunosorbent assays (c-ELISAs). Out of the 1330 positive samples, 200 samples were analysed by serum neutralization test (SNT) to identify circulating BTV serotypes.

RESULTS

Of 3314 sera, 1330 were c-ELISA-positive (40.1%) for antibodies against the BTV structural protein VP7. The result of SNT showed the presence of BTV-1, BTV-2, BTV-3, BTV-4 and, for the first time in Tunisia, BTV-26. The logistic regression model revealed that older animals had nearly two times the odds of being infected with BTV compared to younger animals. Flocks with a history of BT were almost 1.5 times more likely to be at risk for contracting BTV infection. The flock size, housing indoors and intensive production system were significant protective factors.

CONCLUSIONS

High seroprevalence of BTV among sheep was highlighted in Tunisia. The neutralization test showed the presence of the following BTV serotypes: BTV-1, BTV-2, BTV-3, BTV-4 and, for the first time in Tunisia, BTV-26. Age, production system and flock size were important variables associated with BTV infection in sheep. This finding is crucial, as it will allow the adjustment of the BT control programme in Tunisia.

Analysis of bluetongue disease epizootics in sheep of Andhra Pradesh, India using spatial and temporal autocorrelation

Bluetongue (BT) disease poses a constant risk to the livestock population around the world. A better understanding of the risk factors will enable a more accurate prediction of the place and time of high-risk events. Mapping the disease epizootics over a period in a particular geographic area will identify the spatial distribution of disease occurrence. A Geographical Information System (GIS) based methodology to analyze the relationship between bluetongue epizootics and spatial-temporal patterns was used for the years 2000 to 2015 in sheep of Andhra Pradesh, India. Autocorrelation (ACF), partial autocorrelation (PACF), and cross-correlation (CCF) analyses were carried out to find the self-dependency between BT epizootics and their dependencies on environmental factors and livestock population. The association with climatic or remote sensing variables at different months lag, including wind speed, temperature, rainfall, relative humidity, normalized difference vegetation index (NDVI), normalized difference water index (NDWI), land surface temperature (LST), was also examined. The ACF & PACF of BT epizootics with its lag showed a significant positive autocorrelation with a month's lag (r = 0.41). Cross-correlations between the environmental variables and BT epizootics indicated the significant positive correlations at 0, 1, and 2 month's lag of rainfall, relative humidity, normalized difference water index (NDWI), and normalized difference vegetation index (NDVI). Spatial autocorrelation analysis estimated the univariate global Moran's I value of 0.21. Meanwhile, the local Moran's I value for the year 2000 (r = 0.32) showed a high degree of spatial autocorrelation. The spatial autocorrelation analysis revealed that the BT epizootics in sheep are having considerable spatial association among the outbreaks in nearby districts, and have to be taken care of while making any forecasting or disease prediction with other risk factors.

Red deer reveal spatial risks of Crimean-Congo haemorrhagic fever virus infection

Crimean-Congo haemorrhagic fever virus (CCHFV) continues to cause new human cases in Iberia while its spatial distribution and ecological determinants remain unknown. The virus remains active in a silent tick-animal cycle to which animals contribute maintaining the tick populations and the virus itself. Wild ungulates, in particular red deer, are essential hosts for Hyalomma ticks in Iberia, which are the principal competent vector of CCHFV. Red deer could be an excellent model to understand the ecological determinants of CCHFV as well as to predict infection risks for humans because it is large, gregarious, abundant and the principal host for Hyalomma lusitanicum. We designed a cross-sectional study, analysed the presence of CCHFV antibodies in 1444 deer from 82 populations, and statistically modelled exposure risk with host and environmental predictors. The best-fitted statistical model was projected for peninsular Spain to map infection risks. Fifty out of 82 deer populations were seropositive, with individual population prevalence as high as 88%. The highest prevalence of exposure to CCHFV occurred in the southwest of the Iberian Peninsula. Climate and ungulate abundance were the most influential predictors of the risk of exposure to the virus. The highest risk regions were those where H. lusitanicum is most abundant. Eight of the nine primary human cases occurred in or bordering these regions, demonstrating that the model predicts human infection risk accurately. A recent human case of CCHF occurred in northwestern Spain, a region that the model predicted as low risk, pointing out that it needs improvement to capture all determinants of the CCHFV infection risk. In this study, we have been able to identify the main ecological determinants of CCHFV, and we have also managed to create an accurate model to assess the risk of CCHFV infection.

Environment, vector, or host? Using machine learning to untangle the mechanisms driving arbovirus outbreaks

Climatic, landscape, and host features are critical components in shaping outbreaks of vector-borne diseases. However, the relationship between the outbreaks of vector-borne pathogens and their environmental drivers is typically complicated, nonlinear, and may vary by taxonomic units below the species level (e.g., strain or serotype). Here, we aim to untangle how these complex forces shape the risk of outbreaks of Bluetongue virus (BTV); a vector-borne pathogen that is continuously emerging and re-emerging across Europe, with severe economic implications. We tested if the ecological predictors of BTV outbreak risk were serotype-specific by examining the most prevalent serotypes recorded in Europe (1, 4, and 8). We used a robust machine learning (ML) pipeline and 23 relevant environmental features to fit predictive models to 24,245 outbreaks reported in 25 European countries between 2000 and 2019. Our ML models demonstrated high predictive performance for all BTV serotypes (accuracies > 0.87) and revealed strong nonlinear relationships between BTV outbreak risk and environmental and host features. Serotype-specific analysis suggests, however, that each of the major serotypes (1, 4, and 8) had a unique outbreak risk profile. For example, temperature and midge abundance were as the most important characteristics shaping serotype 1, whereas for serotype 4 goat density and temperature were more important. We were also able to identify strong interactive effects between environmental and host characteristics that were also serotype specific. Our ML pipeline was able to reveal more in-depth insights into the complex epidemiology of BTVs and can guide policymakers in intervention strategies to help reduce the economic implications and social cost of this important pathogen.

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.

Serological Cross-Reactions between Expressed VP2 Proteins from Different Bluetongue Virus Serotypes

Bluetongue (BT) is a severe and economically important disease of ruminants that is widely distributed around the world, caused by the bluetongue virus (BTV). More than 28 different BTV serotypes have been identified in serum neutralisation tests (SNT), which, along with geographic variants (topotypes) within each serotype, reflect differences in BTV outer-capsid protein VP2. VP2 is the primary target for neutralising antibodies, although the basis for cross-reactions and serological variations between and within BTV serotypes is poorly understood. Recombinant BTV VP2 proteins (rVP2) were expressed in Nicotiana benthamiana, based on sequence data for isolates of thirteen BTV serotypes (primarily from Europe), including three 'novel' serotypes (BTV-25, -26 and -27) and alternative topotypes of four serotypes. Cross-reactions within and between these viruses were explored using rabbit anti-rVP2 sera and post BTV-infection sheep reference-antisera, in I-ELISA (with rVP2 target antigens) and SNT (with reference strains of BTV-1 to -24, -26 and -27). Strong reactions were generally detected with homologous rVP2 proteins or virus strains/serotypes. The sheep antisera were largely serotype-specific in SNT, but more cross-reactive by ELISA. Rabbit antisera were more cross-reactive in SNT, and showed widespread, high titre cross-reactions against homologous and heterologous rVP2 proteins in ELISA. Results were analysed and visualised by antigenic cartography, showing closer relationships in some, but not all cases, between VP2 topotypes within the same serotype, and between serotypes belonging to the same 'VP2 nucleotype'.

Blood meal analysis: host-feeding patterns of biting midges (Diptera, Ceratopogonidae, Culicoides Latreille) in Slovakia

Biting midges of the genus Culicoides are vectors of important pathogens affecting domestic and wild animals and have played a major role in the re-emergence of new outbreaks of bluetongue (BTV) and Schmallenberg (SBV) viruses in Europe. To determine vector-host specificity, trophic preference from blood meal analysis is of major importance in the surveillance of arthropod-borne diseases. Of 28,752 specimens collected, we identified 17 Culicoides species and investigated a total of 48 host sequences from the blood meals. Culicoides obsoletus/C. scoticus, C. dewulfi, C. pulicaris, C. lupicaris, C. punctatus, C. newsteadi, C. riethi, and C. furcillatus were found to feed on mammals (cattle, horses, and humans), birds (domestic chickens), small rodents (Apodemus flavicollis), and hares (Lepus europaeus). To our knowledge, this is the first study investigating trophic preferences of Culicoides spp. in Slovakia. This study demonstrated that Culicoides species are able to feed on domesticated host vertebrates as well as birds, rodents, and humans.

Phylogeographic analysis of foot-and-mouth disease virus serotype O dispersal and associated drivers in East Africa

The continued endemicity of foot and mouth disease virus (FMDV) in East Africa has significant implications for livestock production and poverty reduction, yet its complex epidemiology in endemic settings remains poorly understood. Identifying FMDV dispersal routes and drivers of transmission is key to improved control strategies. Environmental heterogeneity and anthropogenic drivers (e.g., demand for animal products) can impact viral spread by influencing host movements. Here, we utilized FMDV serotype O VP1 genetic sequences and corresponding spatiotemporal data in order to (i) infer the recent dispersal history, and (II) investigate the impact of external factors (cattle density, human population density, proximity to livestock markets, and drought) on dispersal velocity, location, and direction of FMDV serotype O in East Africa. We identified statistical evidence of long-distance transmission events, and we found that FMDV serotype O tends to remain circulating in areas of high cattle density, high human population density, and in close proximity to livestock markets. The latter two findings highlight the influence of anthropogenic factors on FMDV serotype O spread in this region. These findings contribute to the understanding of FMDV epidemiology in East Africa and can help guide improved control measures.

Landscape Seroprevalence of Three Hemorrhagic Disease-Causing Viruses in a Wild Cervid

Disease plays a major role in shaping wildlife populations worldwide, and changes in landscape conditions can significantly influence risk of pathogen exposure, a threat to vulnerable wild species. Three viruses that cause hemorrhagic disease affect cervid populations in the USA (Odocoileus hemionus adenovirus, bluetongue virus, and epizootic hemorrhagic disease virus), but little is known of their distribution and prevalence in wild populations. We explored the distribution and co-occurrence of seroprevalence of these three pathogens in southern mule deer (Odocoileus hemionus fuliginatus), a subspecies of conservation concern and a harvested species native to southern California, to evaluate the distribution of exposure to these pathogens relative to landscape attributes. We found that habitat type, level of development, and proximity to livestock may affect hemorrhagic disease seroprevalence in southern mule deer. Continued monitoring of hemorrhagic disease-causing viruses in areas where deer are in proximity to cattle and human development is needed to better understand the implications of future outbreaks in wild populations and to identify opportunities to mitigate disease impacts in southern mule deer and other cervid species.

High dispersal capacity of Culicoides obsoletus (Diptera: Ceratopogonidae), vector of bluetongue and Schmallenberg viruses, revealed by landscape genetic analyses

BACKGROUND

In the last two decades, recurrent epizootics of bluetongue virus and Schmallenberg virus have been reported in the western Palearctic region. These viruses affect domestic cattle, sheep, goats and wild ruminants and are transmitted by native hematophagous midges of the genus Culicoides (Diptera: Ceratopogonidae). Culicoides dispersal is known to be stratified, i.e. due to a combination of dispersal processes occurring actively at short distances and passively or semi-actively at long distances, allowing individuals to jump hundreds of kilometers.

METHODS

Here, we aim to identify the environmental factors that promote or limit gene flow of Culicoides obsoletus, an abundant and widespread vector species in Europe, using an innovative framework integrating spatial, population genetics and statistical approaches. A total of 348 individuals were sampled in 46 sites in France and were genotyped using 13 newly designed microsatellite markers.

RESULTS

We found low genetic differentiation and a weak population structure for C. obsoletus across the country. Using three complementary inter-individual genetic distances, we did not detect any significant isolation by distance, but did detect significant anisotropic isolation by distance on a north-south axis. We employed a multiple regression on distance matrices approach to investigate the correlation between genetic and environmental distances. Among all the environmental factors that were tested, only cattle density seems to have an impact on C. obsoletus gene flow.

CONCLUSIONS

The high dispersal capacity of C. obsoletus over land found in the present study calls for a re-evaluation of the impact of Culicoides on virus dispersal, and highlights the urgent need to better integrate molecular, spatial and statistical information to guide vector-borne disease control.

Bluetongue virus in South America: current status based on phylogenetic analysis

Bluetongue (BT) is an insect-borne disease affecting domestic and wild ruminants. Bluetongue virus (BTV) is the causative agent of the BT disease. BT outbreaks have been widely recorded worldwide. However, in the South American subcontinent, accurate information about the disease and molecular epidemiology is still lacking because little effort has been made to cover the region. This study comprises an exhaustive phylogenetic analysis including all BTV sequences available in databases and reports new Argentinean sequences for Seg 8 and Seg 9. Maximum-likelihood phylogenetic analyses were conducted for Seg 2, Seg 3, Seg 6, Seg 7, Seg 8, Seg 9 and Seg 10. Throughout the study, wide circulation and genetic continuity along the American continent were detected. Also, reassortment events are reported, and the historical virus introduction path into and through South America is suggested.

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.

Global emergence and evolutionary dynamics of bluetongue virus

Bluetongue virus (BTV) epidemics are responsible for worldwide economic losses of up to US$ 3 billion. Understanding the global evolutionary epidemiology of BTV is critical in designing intervention programs. Here we employed phylodynamic models to quantify the evolutionary characteristics, spatiotemporal origins, and multi-host transmission dynamics of BTV across the globe. We inferred that goats are the ancestral hosts for BTV but are less likely to be important for cross-species transmission, sheep and cattle continue to be important for the transmission and maintenance of infection between other species. Our models pointed to China and India, countries with the highest population of goats, as the likely ancestral country for BTV emergence and dispersal worldwide over 1000 years ago. However, the increased diversification and dispersal of BTV coincided with the initiation of transcontinental livestock trade after the 1850s. Our analysis uncovered important epidemiological aspects of BTV that may guide future molecular surveillance of BTV.

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.

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.

Continuous Cell Lines from the European Biting Midge Culicoides nubeculosus (Meigen, 1830)

Culicoides biting midges (Diptera: Ceratopogonidae) transmit arboviruses of veterinary or medical importance, including bluetongue virus (BTV) and Schmallenberg virus, as well as causing severe irritation to livestock and humans. Arthropod cell lines are essential laboratory research tools for the isolation and propagation of vector-borne pathogens and the investigation of host-vector-pathogen interactions. Here we report the establishment of two continuous cell lines, CNE/LULS44 and CNE/LULS47, from embryos of Culicoides nubeculosus, a midge distributed throughout the Western Palearctic region. Species origin of the cultured cells was confirmed by polymerase chain reaction (PCR) amplification and sequencing of a fragment of the cytochrome oxidase 1 gene, and the absence of bacterial contamination was confirmed by bacterial 16S rRNA PCR. Both lines have been successfully cryopreserved and resuscitated. The majority of cells examined in both lines had the expected diploid chromosome number of 2n = 6. Transmission electron microscopy of CNE/LULS44 cells revealed the presence of large mitochondria within cells of a diverse population, while arrays of virus-like particles were not seen. CNE/LULS44 cells supported replication of a strain of BTV serotype 1, but not of a strain of serotype 26 which is not known to be insect-transmitted. These new cell lines will expand the scope of research on Culicoides-borne pathogens.

A model for the assessment of bluetongue virus serotype 1 persistence in Spain

Bluetongue virus (BTV) is an arbovirus of ruminants that has been circulating in Europe continuously for more than two decades and has become endemic in some countries such as Spain. Spain is ideal for BTV epidemiological studies since BTV outbreaks from different sources and serotypes have occurred continuously there since 2000; BTV-1 has been reported there from 2007 to 2017. Here we develop a model for BTV-1 endemic scenario to estimate the risk of an area becoming endemic, as well as to identify the most influential factors for BTV-1 persistence. We created abundance maps at 1-km2 spatial resolution for the main vectors in Spain, Culicoides imicola and Obsoletus and Pulicaris complexes, by combining environmental satellite data with occurrence models and a random forest machine learning algorithm. The endemic model included vector abundance and host-related variables (farm density). The three most relevant variables in the endemic model were the abundance of C. imicola and Obsoletus complex and density of goat farms (AUC 0.86); this model suggests that BTV-1 is more likely to become endemic in central and southwestern regions of Spain. It only requires host- and vector-related variables to identify areas at greater risk of becoming endemic for bluetongue. Our results highlight the importance of suitable Culicoides spp. prediction maps for bluetongue epidemiological studies and decision-making about control and eradication measures.

Identification and Characterization of Bluetongue Virus Serotype 14 in Russia

This paper reports a case of bluetongue virus (BTV) infection in the Smolensk and Kaluga regions of Russia in 2011-2012. The virus was initially detected in heifers transferred in Russia from Germany through Poland and Belarus in 2011. On day 27 of quarantine, RNA and infectious viruses of BTV were detected in four heifers, but five were serologically positive. However, on day 3 before shipment, all heifers were seronegative and PCR-negative for BTV. Thus, a few animals from this consignment were viremic without any evident subclinical infection. Based on Seg-2 (VP2 gene) and Seg-5 (NS1 gene) sequencing, the recovered virus had 99.86-100% nucleotide identity with BTV-14-like viruses such as the vaccine BTV-14 strain RSArrrr/BTV 14 and the BTV-14 isolates detected in Lithuania and Poland in 2012. Subsequently, BTV-14 was also reported in local animals in two regions of Russia. During the monitoring survey, 1623 local animals within a 300-km radius were tested, of which 471 tested positive by ELISA and 183 by PCR for BTV-14 RNA. No other serotypes were identified in either imported or aboriginal animals within that radius. The Culicoides midges trapped at the site of the outbreak in May 2012 tested positive for the BTV-14 genome, indicating that the possible mechanism of spread most likely occurs via vector bites. However, further investigation is required to confirm this hypothesis, which would provide an improved understanding of the circulation and overwintering of BTV in northern latitudes.

Agricultural and geographic factors shaped the North American 2015 highly pathogenic avian influenza H5N2 outbreak

The 2014-2015 highly pathogenic avian influenza (HPAI) H5NX outbreak represents the largest and most expensive HPAI outbreak in the United States to date. Despite extensive traditional and molecular epidemiological studies, factors associated with the spread of HPAI among midwestern poultry premises remain unclear. To better understand the dynamics of this outbreak, 182 full genome HPAI H5N2 sequences isolated from commercial layer chicken and turkey production premises were analyzed using evolutionary models able to accommodate epidemiological and geographic information. Epidemiological compartmental models embedded in a phylogenetic framework provided evidence that poultry type acted as a barrier to the transmission of virus among midwestern poultry farms. Furthermore, after initial introduction, the propagation of HPAI cases was self-sustainable within the commercial poultry industries. Discrete trait diffusion models indicated that within state viral transitions occurred more frequently than inter-state transitions. Distance and sample size were very strongly supported as associated with viral transition between county groups (Bayes Factor > 30.0). Together these findings indicate that the different types of midwestern poultry industries were not a single homogenous population, but rather, the outbreak was shaped by poultry industries and geographic factors.

Using phylogeographic approaches to analyse the dispersal history, velocity and direction of viral lineages - Application to rabies virus spread in Iran

Recent years have seen the extensive use of phylogeographic approaches to unveil the dispersal history of virus epidemics. Spatially explicit reconstructions of viral spread represent valuable sources of lineage movement data that can be exploited to investigate the impact of underlying environmental layers on the dispersal of pathogens. Here, we performed phylogeographic inference and applied different post hoc approaches to analyse a new and comprehensive data set of viral genomes to elucidate the dispersal history and dynamics of rabies virus (RABV) in Iran, which have remained largely unknown. We first analysed the association between environmental factors and variations in dispersal velocity among lineages. Second, we present, test and apply a new approach to study the link between environmental conditions and the dispersal direction of lineages. The statistical performance (power of detection, false-positive rate) of this new method was assessed using simulations. We performed phylogeographic analyses of RABV genomes, allowing us to describe the large diversity of RABV in Iran and to confirm the cocirculation of several clades in the country. Overall, we estimate a relatively high lineage dispersal velocity, similar to previous estimates for dog rabies virus spread in northern Africa. Finally, we highlight a tendency for RABV lineages to spread in accessible areas associated with high human population density. Our analytical workflow illustrates how phylogeographic approaches can be used to investigate the impact of environmental factors on several aspects of viral dispersal dynamics.

Comparing patterns and scales of plant virus phylogeography: Rice yellow mottle virus in Madagascar and in continental Africa

Rice yellow mottle virus (RYMV) in Madagascar Island provides an opportunity to study the spread of a plant virus disease after a relatively recent introduction in a large and isolated country with a heterogeneous host landscape ecology. Here, we take advantage of field survey data on the occurrence of RYMV disease throughout Madagascar dating back to the 1970s, and of virus genetic data from ninety-four isolates collected since 1989 in most regions of the country to reconstruct the epidemic history. We find that the Malagasy isolates belong to a unique recombinant strain that most likely entered Madagascar through a long-distance introduction from the most eastern part of mainland Africa. We infer the spread of RYMV as a continuous process using a Bayesian statistical framework. In order to calibrate the time scale in calendar time units in this analysis, we pool the information about the RYMV evolutionary rate from several geographical partitions. Whereas the field surveys and the phylogeographic reconstructions both point to a rapid southward invasion across hundreds of kilometers throughout Madagascar within three to four decades, they differ on the inferred origin location and time of the epidemic. The phylogeographic reconstructions suggest a lineage displacement and unveil a re-invasion of the northern regions that may have remained unnoticed otherwise. Despite ecological differences that could affect the transmission potential of RYMV in Madagascar and in mainland Africa, we estimate similar invasion and dispersal rates. We could not identify environmental factors that have a relevant impact on the lineage dispersal velocity of RYMV in Madagascar. This study highlights the value and complementarity of (historical) nongenetic and (more contemporaneous) genetic surveillance data for reconstructing the history of spread of plant viruses.

A cross-sectional serosurvey in a sheep population in central Italy following a bluetongue epidemic

Bluetongue (BT) is a viral disease that affects ruminants and is transmitted by midges of the genus Culicoides spp. The seroprevalence, the clinical form and the occurrence rates significantly differ in relation to several factors such as bluetongue virus (BTV) serotype, host species, breed susceptibility, specific previous exposure, vector ecology, husbandry and health status. Following the 2001-2006 BTV2 and BTV16 epidemics in central Italy, a new epidemic caused by BTV1 occurred in 2013-2015 causing 398 outbreaks in a susceptible population of about 1 million ruminants. The present study assessed the BTV1 seroprevalence in the sheep population of central Italy by conducting two cross-sectional surveys, in the proximity of and within BT outbreak farms. A total of 2,984 sheep from 437 farms were sampled. The animal-level prevalence was 19% (95% CI: 17-21%), the between-herd prevalence was 46% (95% CI: 41-51%) and the within-herd prevalence was 21% (95% CI: 16-26%). Risk factors were investigated by logistic regression models. Living on a farm where an outbreak occurred and the number of outbreaks in proximity of the farm were identified as risk factors, while herd size was identified as a protective factor. This study represents the first BT survey in southern Europe and reports valuable findings on BTV epidemiology. Despite intensive virus circulation, the estimated seroprevalences were low. The assessment of the population immunity level is crucial for defining an efficient vaccination strategy and for predicting the impact of future virus circulation. In view of the low seroprevalence detected albeit an extensive BTV1 circulation, the population immunity was likely to be inadequate in preventing new BTV1 epidemics. Moreover, considering the recurrent introduction of new serotypes from North Africa and the Balkans, the control of multi-serotype BTV infections will continue to present a challenge in the near future.

Transmission dynamics of re-emerging rabies in domestic dogs of rural China

Despite ongoing efforts to control transmission, rabies prevention remains a challenge in many developing countries, especially in rural areas of China where re-emerging rabies is under-reported due to a lack of sustained animal surveillance. By taking advantage of detailed genomic and epidemiological data for the re-emerging rabies outbreak in Yunnan Province, China, collected between 1999 and 2015, we reconstruct the demographic and dispersal history of domestic dog rabies virus (RABV) as well as the dynamics of dog-to-dog and dog-to-human transmission. Phylogeographic analyses reveal a lower diffusion coefficient than previously estimated for dog RABV dissemination in northern Africa. Furthermore, epidemiological analyses reveal transmission rates between dogs, as well as between dogs and humans, lower than estimates for Africa. Finally, we show that reconstructed epidemic history of RABV among dogs and the dynamics of rabid dogs are consistent with the recorded human rabies cases. This work illustrates the benefits of combining phylogeographic and epidemic modelling approaches for uncovering the spatiotemporal dynamics of zoonotic diseases, with both approaches providing estimates of key epidemiological parameters.

Pathogens in space: Advancing understanding of pathogen dynamics and disease ecology through landscape genetics

Landscape genetics has provided many insights into how heterogeneous landscape features drive processes influencing spatial genetic variation in free-living organisms. This rapidly developing field has focused heavily on vertebrates, and expansion of this scope to the study of infectious diseases holds great potential for landscape geneticists and disease ecologists alike. The potential application of landscape genetics to infectious agents has garnered attention at formative stages in the development of landscape genetics, but systematic examination is lacking. We comprehensively review how landscape genetics is being used to better understand pathogen dynamics. We characterize the field and evaluate the types of questions addressed, approaches used and systems studied. We also review the now established landscape genetic methods and their realized and potential applications to disease ecology. Lastly, we identify emerging frontiers in the landscape genetic study of infectious agents, including recent phylogeographic approaches and frameworks for studying complex multihost and host-vector systems. Our review emphasizes the expanding utility of landscape genetic methods available for elucidating key pathogen dynamics (particularly transmission and spread) and also how landscape genetic studies of pathogens can provide insight into host population dynamics. Through this review, we convey how increasing awareness of the complementarity of landscape genetics and disease ecology among practitioners of each field promises to drive important cross-disciplinary advances.

On the importance of negative controls in viral landscape phylogeography

Phylogeographic reconstructions are becoming an established procedure to evaluate the factors that could impact virus spread. While a discrete phylogeographic approach can be used to test predictors of transition rates among discrete locations, alternative continuous phylogeographic reconstructions can also be exploited to investigate the impact of underlying environmental layers on the dispersal velocity of a virus. The two approaches are complementary tools for studying pathogens' spread, but in both cases, care must be taken to avoid misinterpretations. Here, we analyse rice yellow mottle virus (RYMV) sequence data from West and East Africa to illustrate how both approaches can be used to study the impact of environmental factors on the virus’ dispersal frequency and velocity. While it was previously reported that host connectivity was a major determinant of RYMV spread, we show that this was a false positive result due to the lack of appropriate negative controls. We also discuss and compare the phylodynamic tools currently available for investigating the impact of environmental factors on virus spread.