Bluetongue in Spain: From the First Outbreak to 2012

The use of path analysis to determine effects of environmental factors on the adult seasonality of Culicoides (Diptera: Ceratopogonidae) vector species in Spain

Culicoides biting midges (Diptera: Ceratopogonidae) are the main vectors of livestock diseases such as bluetongue (BT) which mainly affect sheep and cattle. In Spain, bluetongue virus (BTV) is transmitted by several Culicoides taxa, including Culicoides imicola, Obsoletus complex, Culicoides newsteadi and Culicoides pulicaris that vary in seasonality and distribution, affecting the distribution and dynamics of BT outbreaks. Path analysis is useful for separating direct and indirect, biotic and abiotic determinants of species' population performance and is ideal for understanding the sensitivity of adult Culicoides dynamics to multiple environmental drivers. Start, end of season and length of overwintering of adult Culicoides were analysed across 329 sites in Spain sampled from 2005 to 2010 during the National Entomosurveillance Program for BTV with path analysis, to determine the direct and indirect effects of land use, climate and host factor variables. Culicoides taxa had species-specific responses to environmental variables. While the seasonality of adult C. imicola was strongly affected by topography, temperature, cover of agro-forestry and sclerophyllous vegetation, rainfall, livestock density, photoperiod in autumn and the abundance of Culicoides females, Obsoletus complex species seasonality was affected by land-use variables such as cover of natural grassland and broad-leaved forest. Culicoides female abundance was the most explanatory variable for the seasonality of C. newsteadi, while C. pulicaris showed that temperature during winter and the photoperiod in November had a strong effect on the start of the season and the length of overwinter period of this species. These results indicate that the seasonal vector-free period (SVFP) in Spain will vary between competent vector taxa and geographic locations, dependent on the different responses of each taxa to environmental conditions.

Development of microsatellite markers for population genetics of biting midges and a potential tool for species identification of Culicoides sonorensis Wirth & Jones

Background

Proper vector surveillance relies on the ability to identify species of interest accurately and efficiently, though this can be difficult in groups containing cryptic species. Culicoides Latreille is a genus of small biting flies responsible for the transmission of numerous pathogens to a multitude of vertebrates. Regarding pathogen transmission, the C. variipennis species complex is of particular interest in North America. Of the six species within this group, only C. sonorensis Wirth & Jones is a proven vector of bluetongue virus and epizootic hemorrhagic disease virus. Unfortunately, subtle morphological differences, cryptic species, and mitonuclear discordance make species identification in the C. variipennis complex challenging. Recently, single-nucleotide polymorphism (SNP) analysis enabled discrimination between the species of this group; however, this demanding approach is not practical for vector surveillance.

Methods

The aim of the current study was to develop a reliable and affordable way of distinguishing between the species within the C. variipennis complex, especially C. sonorensis. Twenty-five putative microsatellite markers were identified using the C. sonorensis genome and tested for amplification within five species of the C. variipennis complex. Machine learning was then used to determine which markers best explain the genetic differentiation between species. This led to the development of a subset of four and seven markers, which were also tested for species differentiation.

Results

A total of 21 microsatellite markers were successfully amplified in the species tested. Clustering analyses of all of these markers recovered the same species-level identification as the previous SNP data. Additionally, the subset of seven markers was equally capable of accurately distinguishing between the members of the C. variipennis complex as the 21 microsatellite markers. Finally, one microsatellite marker (C508) was found to be species-specific, only amplifying in the vector species C. sonorensis among the samples tested.

Conclusions

These microsatellites provide an affordable way to distinguish between the sibling species of the C. variipennis complex and could lead to a better understanding of the species dynamics within this group. Additionally, after further testing, marker C508 may allow for the identification of C. sonorensis with a single-tube assay, potentially providing a powerful new tool for vector surveillance in North America.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05189-8.

Safety and efficacy of a Bluetongue inactivated vaccine (serotypes 1 and 4) in sheep

A new inactivated vaccine against Bluetongue virus (BTV) serotypes 1 and 4, was developed from field isolates. Safety and efficacy of the vaccine were evaluated in sheep by serological monitoring and virus nucleic acid detection after experimental infection of vaccinated animals. Seroconversion was observed in vaccinated animals at day 14 post vaccination (pv) with neutralizing antibody titer of 1.9 and 1.8 for serotypes 1 and 4, respectively. The titer increase significantly after the booster reaching 2.7 and persist one year >1.5 for both serotypes. After challenge with virulent isolates, vireamia was recorded in control animals, as evident by q-PCR with threshold cycles (Ct) ranging from 24 to 31 and peaked at day 10 post challenge, while no vireamia was detected in vaccinated animals. Vaccinated sheep were fully protected against the disease and infection.

Environmental Drivers of Adult Seasonality and Abundance of Biting Midges Culicoides (Diptera: Ceratopogonidae), Bluetongue Vector Species in Spain

Bluetongue is a viral disease affecting wild and domestic ruminants transmitted by several species of biting midges Culicoides Latreille. The phenology of these insects were analyzed in relation to potential environmental drivers. Data from 329 sites in Spain were analyzed using Bayesian Generalized Linear Mixed Model (GLMM) approaches. The effects of environmental factors on adult female seasonality were contrasted. Obsoletus complex species (Diptera: Ceratopogonidae) were the most prevalent across sites, followed by Culicoides newsteadi Austen (Diptera: Ceratopogonidae). Activity of female Obsoletus complex species was longest in sites at low elevation, with warmer spring average temperatures and precipitation, as well as in sites with high abundance of cattle. The length of the Culicoides imicola Kieffer (Diptera: Ceratopogonidae) female adult season was also longest in sites at low elevation with higher coverage of broad-leaved vegetation. Long adult seasons of C. newsteadi were found in sites with warmer autumns and higher precipitation, high abundance of sheep. Culicoides pulicaris (Linnaeus) (Diptera: Ceratopogonidae) had longer adult periods in sites with a greater number of accumulated degree days over 10°C during winter. These results demonstrate the eco-climatic and seasonal differences among these four taxa in Spain, which may contribute to determining sites with suitable environmental circumstances for each particular species to inform assessments of the risk of Bluetongue virus outbreaks in this region.

Identifying Spanish Areas at More Risk of Monthly BTV Transmission with a Basic Reproduction Number Approach

Bluetongue virus (BTV) causes a disease that is endemic in Spain and its two major biological vector species, C. imicola and the Obsoletus complex species, differ greatly in their ecology and distribution. Understanding the seasonality of BTV transmission in risk areas is key to improving surveillance and control programs, as well as to better understand the pathogen transmission networks between wildlife and livestock. Here, monthly risk transmission maps were generated using risk categories based on well-known BTV R0 equations and predicted abundances of the two most relevant vectors in Spain. Previously, Culicoides spp. predicted abundances in mainland Spain and the Balearic Islands were obtained using remote sensing data and random forest machine learning algorithm. Risk transmission maps were externally assessed with the estimated date of infection of BTV-1 and BTV-4 historical outbreaks. Our results highlight the differences in risk transmission during April-October, June-August being the period with higher R0 values. Likewise, a natural barrier has been identified between northern and central-southern areas at risk that may hamper BTV spread between them. Our results can be relevant to implement risk-based interventions for the prevention, control and surveillance of BTV and other diseases shared between livestock and wildlife host populations.

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.

Cognition and behavior in sheep repetitively inoculated with aluminum adjuvant-containing vaccines or aluminum adjuvant only

Sheep health management strategies often include the use of aluminum (Al)-containing vaccines. These products were associated with the appearance of the ovine autoimmune/inflammatory syndrome induced by adjuvants (ASIA syndrome), which included an array of ethological changes in the affected animals. The aim of this pilot study was to investigate cognitive and behavioral changes in sheep subjected to a protocol of repetitive inoculation with Al-containing products. Twenty-one lambs were assigned to three groups (n = 7 each): Control, Adjuvant-only, and Vaccine. Vaccine group was inoculated with commercial Al- hydroxide containing vaccines; Adjuvant-only group received the equivalent dose of Al only (Alhydrogel®), and Control group received Phosphate-buffered saline. Sixteen inoculations were administered within a 349-day period. Ethological changes were studied in late summer (7 inoculations) and mid-winter (16 inoculations). Animals in Vaccine and Adjuvant-only groups exhibited individual and social behavioral changes. Affiliative interactions were significantly reduced, and aggressive interactions and stereotypies increased significantly. They also exhibited a significant increase in excitatory behavior and compulsive eating. There were increased levels of stress biomarkers in these two groups. In general, changes were more pronounced in the Vaccine group than they were in the Adjuvant-only group. Some changes were already significant in summer, after seven inoculations only. This study is the first to describe behavioral changes in sheep after having received repetitive injections of Al-containing products, and may explain some of the clinical signs observed in ovine ASIA syndrome.

Granulomas Following Subcutaneous Injection With Aluminum Adjuvant-Containing Products in Sheep

The use of vaccines including aluminum (Al)–based adjuvants is widespread among small ruminants and other animals. They are associated with the appearance of transient injection site nodules corresponding to granulomas. This study aims to characterize the morphology of these granulomas, to understand the role of the Al adjuvant in their genesis, and to establish the presence of the metal in regional lymph nodes. A total of 84 male neutered lambs were selected and divided into 3 treatment groups of 28 animals each: (1) vaccine (containing Al-based adjuvant), (2) adjuvant-only, and (3) control. A total of 19 subcutaneous injections were performed in a time frame of 15 months. Granulomas and regional lymph nodes were evaluated by clinicopathological means. All of the vaccine and 92.3% of the adjuvant-only lambs presented injection-site granulomas; the granulomas were more numerous in the group administered the vaccine. Bacterial culture in granulomas was always negative. Histologically, granulomas in the vaccine group presented a higher degree of severity. Al was specifically identified by lumogallion staining in granulomas and lymph nodes. Al median content was significantly higher ( P < .001) in the lymph nodes of the vaccine group (82.65 μg/g) compared with both adjuvant-only (2.53 μg/g) and control groups (0.96 μg/g). Scanning transmission electron microscopy demonstrated aggregates of Al within macrophages in vaccine and adjuvant-only groups. In these two groups, Al-based adjuvants induce persistent, sterile, subcutaneous granulomas with macrophage-driven translocation of Al to regional lymph nodes. Local translocation of Al may induce further accumulation in distant tissues and be related to the appearance of systemic signs.

An advection-deposition-survival model to assess the risk of introduction of vector-borne diseases through the wind: Application to bluetongue outbreaks in Spain

This work develops a methodology for estimating risk of wind-borne introduction of flying insects into a country, identifying areas and periods of high risk of vector-borne diseases incursion. This risk can be characterized by the role of suitable temperatures and wind currents in small insects' survival and movements, respectively. The model predicts the number density of introduced insects over space and time based on three processes: the advection due to wind currents, the deposition on the ground and the survival due to climatic conditions. Spanish livestock has suffered many bluetongue outbreaks since 2004 and numerous experts point to Culicoides transported by wind from affected areas in North Africa as a possible cause. This work implements numerical experiments simulating the introduction of Culicoides in 2004. The model identified southern and eastern Spain, particularly between June and November, as being at greatest risk of wind-borne Culicoides introduction, which matches field data on bluetongue outbreaks in Spain this year. This validation suggests that this model may be useful for predicting introduction of airborne pathogens of significance to animal productivity.

First molecular identification of the vertebrate hosts of Culicoides imicola in Europe and a review of its blood-feeding patterns worldwide: implications for the transmission of bluetongue disease and African horse sickness

Culicoides (Diptera: Ceratopogonidae) are vectors of pathogens that affect wildlife, livestock and, occasionally, humans. Culicoides imicola (Kieffer, 1913) is considered to be the main vector of the pathogens that cause bluetongue disease (BT) and African horse sickness (AHS) in southern Europe. The study of blood-feeding patterns in Culicoides is an essential step towards understanding the epidemiology of these pathogens. Molecular tools that increase the accuracy and sensitivity of traditional methods have been developed to identify the hosts of potential insect vectors. However, to the present group's knowledge, molecular studies that identify the hosts of C. imicola in Europe are lacking. The present study genetically characterizes the barcoding region of C. imicola trapped on farms in southern Spain and identifies its vertebrate hosts in the area. The report also reviews available information on the blood-feeding patterns of C. imicola worldwide. Culicoides imicola from Spain feed on blood of six mammals that include species known to be hosts of the BT and AHS viruses. This study provides evidence of the importance of livestock as sources of bloodmeals for C. imicola and the relevance of this species in the transmission of BT and AHS viruses in Europe.

Assessment of listing and categorisation of animal diseases within the framework of the Animal Health Law (Regulation (EU) No 2016/429): bluetongue

A specific concept of strain was developed in order to classify the BTV serotypes ever reported in Europe based on their properties of animal health impact: the genotype, morbidity, mortality, speed of spread, period and geographical area of occurrence were considered as classification parameters. According to this methodology the strain groups identified were (i) the BTV strains belonging to serotypes BTV‐1–24, (ii) some strains of serotypes BTV‐16 and (iii) small ruminant‐adapted strains belonging to serotypes BTV‐25, ‐27, ‐30. Those strain groups were assessed according to the criteria of the Animal Health Law (AHL), in particular criteria of Article 7, Article 5 on the eligibility of bluetongue to be listed, Article 9 for the categorisation according to disease prevention and control rules as in Annex IV and Article 8 on the list of animal species related to bluetongue. The assessment has been performed following a methodology composed of information collection, expert judgement at individual and collective level. The output is composed of the categorical answer, and for the questions where no consensus was reached, the different supporting views are reported. The strain group BTV (1–24) can be considered eligible to be listed for Union intervention as laid down in Article 5(3) of the AHL, while the strain group BTV‐25–30 and BTV‐16 cannot. The strain group BTV‐1–24 meets the criteria as in Sections 2 and 5 of Annex IV of the AHL, for the application of the disease prevention and control rules referred to in points (b) and (e) of Article 9(1) of the AHL. The animal species that can be considered to be listed for BTV‐1–24 according to Article 8(3) are several species of Bovidae, Cervidae and Camelidae as susceptible species; domestic cattle, sheep and red deer as reservoir hosts, midges insect of genus Culicoides spp. as vector species.

Understanding Spatio-Temporal Variability in the Reproduction Ratio of the Bluetongue (BTV-1) Epidemic in Southern Spain (Andalusia) in 2007 Using Epidemic Trees

Andalusia (Southern Spain) is considered one of the main routes of introduction of bluetongue virus (BTV) into Europe, evidenced by a devastating epidemic caused by BTV-1 in 2007. Understanding the pattern and the drivers of BTV-1 spread in Andalusia is critical for effective detection and control of future epidemics. A long-standing metric for quantifying the behaviour of infectious diseases is the case-reproduction ratio (R_t), defined as the average number of secondary cases arising from a single infected case at time t (for t>0). Here we apply a method using epidemic trees to estimate the between-herd case reproduction ratio directly from epidemic data allowing the spatial and temporal variability in transmission to be described. We then relate this variability to predictors describing the hosts, vectors and the environment to better understand why the epidemic spread more quickly in some regions or periods. The R_t value for the BTV-1 epidemic in Andalusia peaked in July at 4.6, at the start of the epidemic, then decreased to 2.2 by August, dropped below 1 by September (0.8), and by October it had decreased to 0.02. BTV spread was the consequence of both local transmission within established disease foci and BTV expansion to distant new areas (i.e. new foci), which resulted in a high variability in BTV transmission, not only among different areas, but particularly through time, which suggests that general control measures applied at broad spatial scales are unlikely to be effective. This high variability through time was probably due to the impact of temperature on BTV transmission, as evidenced by a reduction in the value of R_t by 0.0041 for every unit increase (day) in the extrinsic incubation period (EIP), which is itself directly dependent on temperature. Moreover, within the range of values at which BTV-1 transmission occurred in Andalusia (20.6°C to 29.5°C) there was a positive correlation between temperature and R_t values, although the relationship was not linear, probably as a result of the complex relationship between temperature and the different parameters affecting BTV transmission. R_t values for BTV-1 in Andalusia fell below the threshold of 1 when temperatures dropped below 21°C, a much higher threshold than that reported in other BTV outbreaks, such as the BTV-8 epidemic in Northern Europe. This divergence may be explained by differences in the adaptation to temperature of the main vectors of the BTV-1 epidemic in Andalusia (Culicoides imicola) compared those of the BTV-8 epidemic in Northern Europe (Culicoides obsoletus). Importantly, we found that BTV transmission (R_t value) increased significantly in areas with higher densities of sheep. Our analysis also established that control of BTV-1 in Andalusia was complicated by the simultaneous establishment of several distant foci at the start of the epidemic, which may have been caused by several independent introductions of infected vectors from the North of Africa. We discuss the implications of these findings for BTV surveillance and control in this region of Europe.

Spatial-temporal Trends and Factors Associated with the Bluetongue Virus Seropositivity in Large Game Hunting Areas from Southern Spain

An epidemiological study was carried out to determine the spatial-temporal trends and risk factors potentially involved in the seropositivity to bluetongue virus (BTV) in hunting areas with presence of red deer (Cervus elaphus). A total of 60 of 98 (61.2%; CI95% : 51.6-70.9) hunting areas sampled presented at least one seropositive red deer. Antibodies against BTV were detected in juvenile animals during the hunting seasons 2007/2008 to 2013/2014 in 15 of 98 (15.3%) hunting areas, which indicates an uninterrupted circulation of BTV in this period. A multivariate logistic regression model showed that the red deer density at hunting area level (>22 individuals/km(2) ), the annual abundance of Culicoides imicola (>1.4 mosquitoes/sampling) and the goat density at municipality level (>24.1 individuals/km(2) ) were factors significantly associated with BTV seropositivity in hunting areas. Control measures against BTV in the studied area include vaccination programmes in wild and domestic ruminants, movement control in areas with high densities and abundance of red deer and C. imicola, respectively. Considering the potential risk of BTV re-emergence, red deer should be included in the BT surveillance programmes in regions where these species share habitats with livestock.

Identification of Suitable Areas for African Horse Sickness Virus Infections in Spanish Equine Populations

African horse sickness (AHS) is one of the most important vector-borne viral infectious diseases of equines, transmitted mainly by Culicoides spp. The re-emergence of Culicoides-borne diseases in Europe, such as the recent bluetongue (BT) or Schmallenberg outbreaks, has raised concern about the potential re-introduction and further spread of AHS virus (AHSV) in Europe. Spain has one of the largest European equine populations. In addition, its geographical, environmental and entomological conditions favour AHSV infections, as shown by the historical outbreaks in the 1990s. The establishment of risk-based surveillance strategies would allow the early detection and rapid control of any potential AHSV outbreak. This study aimed to identify the areas and time periods that are suitable or at high risk for AHS occurrence in Spain using a GIS-based multicriteria decision framework. Specifically risk maps for AHS occurrence were produced using a weighted linear combination of the main risk factors of disease, namely extrinsic incubation period, equine density and distribution of competent Culicoides populations. Model results revealed that the south-western and north-central areas of Spain and the Balearic Islands are the areas at the highest risk for AHSV infections, particularly in late summer months. Conversely, Galicia, Castile and Leon and La Rioja can be considered as low-risk regions. This result was validated with historical AHS and BT outbreaks in Spain, and with the Culicoides vector distribution area. The model results, together with current Spanish equine production features, should provide the foundations to design risk-based and more cost-effective surveillance strategies for the early detection and rapid control potential of AHS outbreaks in Spain.

The role of wildlife in bluetongue virus maintenance in Europe: lessons learned after the natural infection in Spain

Bluetongue (BT) is a re-emergent vector-borne viral disease of domestic and wild ruminants caused by bluetongue virus (BTV), a member of the genus Orbivirus. A complex multi-host, multi-vector and multi-pathogen (26 serotypes) transmission and maintenance network has recently emerged in Europe, and wild ruminants are regarded as an important node in this network. This review analyses the reservoir role of wild ruminants in Europe, identifying gaps in knowledge and proposing actions. Wild ruminant species are indicators of BTV circulation. Excepting the mouflon (Ovis aries musimon), European wild ungulates do not develop clinical disease. Diagnostic techniques used in wildlife do not differ from those used in domestic ruminants provided they are validated. Demographic, behavioural and physiological traits of wild hosts modulate their relationship with BTV vectors and with the virus itself. While BTV has been eradicated from central and northern Europe, it is still circulating in the Mediterranean Basin. We propose that currently two BTV cycles coexist in certain regions of the Mediterranean Basin, a wild one largely driven by deer of the subfamily Cervinae and a domestic one. These are probably linked through shared Culicoides vectors of several species. We suggest that wildlife might be contributing to this situation through vector maintenance and virus maintenance. Additionally, differences in temperature and other environmental factors add complexity to the Mediterranean habitats as compared to central and northern European ones. Intervention options in wildlife populations are limited. There is a need to know the role of wildlife in maintaining Culicoides populations, and to know which Culicoides species mediate the wildlife-livestock-BTV transmission events. There is also a clear need to study more in depth the links between Cervinae deer densities, environmental factors and BTV maintenance. Regarding disease control, we suggest that research efforts should be focused on wildlife population and wildlife disease monitoring.