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

Quantitative risk assessment for the introduction of bluetongue virus into mainland Europe by long-distance wind dispersal of Culicoides spp.: A case study from Sardinia

Europe faces regular introductions and reintroductions of bluetongue virus (BTV) serotypes, most recently exemplified by the incursion of serotype 3 in the Netherlands. Although the long-distance wind dispersal of the disease vector, Culicoides spp., is recognized as a virus introduction pathway, it remains understudied in risk assessments. A Quantitative Risk Assessment framework was developed to estimate the risk of BTV-3 incursion into mainland Europe from Sardinia, where the virus has been present since 2018. We used an atmospheric transport model (HYbrid Single-Particle Lagrangian Integrated Trajectory) to infer the probability of airborne dispersion of the insect vector. Epidemiological disease parameters quantified the virus prevalence in vector population in Sardinia and its potential first transmission after introduction in a new area. When assuming a 24h maximal flight duration, the risk of BTV introduction from Sardinia is limited to the Mediterranean Basin, mainly affecting the southwestern area of the Italian Peninsula, Sicily, Malta, and Corsica. The risk extends to the northern and central parts of Italy, Balearic archipelago, and mainland France and Spain, mostly when maximal flight duration is longer than 24h. Additional knowledge on vector flight conditions and Obsoletus complex-specific parameters could improve the robustness of the model. Providing both spatial and temporal insights into BTV introduction risks, our framework is a key tool to guide global surveillance and preparedness against epizootics.

A study of the composition of the Obsoletus complex and genetic diversity of Culicoides obsoletus populations in Spain

Background

The Culicoides obsoletus species complex (henceforth ‘Obsoletus complex’) is implicated in the transmission of several arboviruses that can cause severe disease in livestock, such as bluetongue, African horse sickness, epizootic hemorrhagic disease and Schmallenberg disease. Thus, this study aimed to increase our knowledge of the composition and genetic diversity of the Obsoletus complex by partial sequencing of the cytochrome c oxidase I (cox1) gene in poorly studied areas of Spain.

Methods

A study of C. obsoletus populations was carried out using a single-tube multiplex polymerase chain reaction (PCR) assay that was designed to differentiate the Obsoletus complex sibling species Culicoides obsoletus and Culicoides scoticus, based on the partial amplification of the cox1 gene, as well as cox1 georeferenced sequences from Spain available at GenBank. We sampled 117 insects of the Obsoletus complex from six locations and used a total of 238 sequences of C. obsoletus (ss) individuals (sampled here, and from GenBank) from 14 sites in mainland Spain, the Balearic Islands and the Canary Islands for genetic diversity and phylogenetic analyses.

Results

We identified 90 C. obsoletus (ss), 19 Culicoides scoticus and five Culicoides montanus midges from the six collection sites sampled, and found that the genetic diversity of C. obsoletus (ss) were higher in mainland Spain than in the Canary Islands. The multiplex PCR had limitations in terms of specificity, and no cryptic species within the Obsoletus complex were identified.

Conclusions

Within the Obsoletus complex, C. obsoletus (ss) was the predominant species in the analyzed sites of mainland Spain. Information about the species composition of the Obsoletus complex could be of relevance for future epidemiological studies when specific aspects of the vector competence and capacity of each species have been identified. Our results indicate that the intraspecific divergence is higher in C. obsoletus (ss) northern populations, and demonstrate the isolation of C. obsoletus (ss) populations of the Canary Islands.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-04841-z.