Updating the global occurrence of Culicoides imicola, a vector for emerging viral diseases

Epidemiological Analyses of the First Incursion of the Epizootic Hemorrhagic Disease Virus Serotype 8 in Tunisia, 2021-2022

Epizootic hemorrhagic disease (EHD) is a non-contagious arthropod-transmitted viral disease and a World Organization for Animal Health (WOAH)-listed disease of domestic and wild ruminants since 2008. EHDV is transmitted among susceptible animals by a few species of midges of genus Culicoides. During the fall of 2021, a large outbreak caused by the epizootic hemorrhagic disease virus (EHDV), identified as serotype 8, was reported in Tunisian dairy and beef farms with Bluetongue virus (BTV)-like clinical signs. The disease was detected later in the south of Italy, in Spain, in Portugal and, more recently, in France, where it caused severe infections in cattle. This was the first evidence of EHDV-8 circulation outside Australia since 1982. In this study, we analyzed the epidemiological situation of the 2021-2022 EHDV outbreaks reported in Tunisia, providing a detailed description of the spatiotemporal evolution of the disease. We attempted to identify the eco-climatic factors associated with infected areas using generalized linear models (GLMs). Our results demonstrated that environmental factors mostly associated with the presence of C. imicola, such as digital elevation model (DEM), slope, normalized difference vegetation index (NDVI), and night-time land surface temperature (NLST)) were by far the most explanatory variables for EHD repartition cases in Tunisia that may have consequences in neighboring countries, both in Africa and Europe through the spread of infected vectors. The risk maps elaborated could be useful for disease control and prevention strategies.

Viral diversity and blood-feeding patterns of Afrotropical Culicoides biting midges (Diptera: Ceratopogonidae)

Introduction

Culicoides biting midges (Diptera: Ceratopogonidae) are vectors of arboviral pathogens that primarily affect livestock represented by Schmallenberg virus (SBV), epizootic hemorrhagic disease virus (EHDV) and bluetongue virus (BTV). In Kenya, studies examining the bionomic features of Culicoides including species diversity, blood-feeding habits, and association with viruses are limited.

Methods

Adult Culicoides were surveyed using CDC light traps in two semi-arid ecologies, Baringo and Kajiado counties, in Kenya. Blood-fed specimens were analysed through polymerase chain reaction (PCR) and sequencing of cytochrome oxidase subunit 1 (cox1) barcoding region. Culicoides pools were screened for virus infection by generic RT-PCR and next-generation sequencing (NGS).

Results

Analysis of blood-fed specimens confirmed that midges had fed on cattle, goats, sheep, zebra, and birds. Cox1 barcoding of the sampled specimens revealed the presence of known vectors of BTV and epizootic hemorrhagic disease virus (EHDV) including species in the Imicola group (Culicoides imicola) and Schultzei group (C. enderleni, C. kingi, and C. chultzei). Culicoides leucostictus and a cryptic species distantly related to the Imicola group were also identified. Screening of generated pools (11,006 individuals assigned to 333 pools) by generic RT-PCR revealed presence of seven phylogenetically distinct viruses grouping in the genera Goukovirus, Pacuvirus and Orthobunyavirus. The viruses showed an overall minimum infection rate (MIR) of 7.0% (66/333, 95% confidence interval (CI) 5.5-8.9). In addition, full coding sequences of two new iflaviruses, tentatively named Oloisinyai_1 and Oloisinyai_2, were generated by next-generation sequencing (NGS) from individual homogenate of Culicoides pool.

Conclusion

The results indicate a high genetic diversity of viruses in Kenyan biting midges. Further insights into host-vector-virus interactions as well as investigations on the potential clinical significance of the detected viruses are warranted.

Biting Midges (Diptera: Ceratopogonidae) as Vectors of Viruses

Biting midges of the genus Culicoides occur almost globally and can regionally and seasonally reach high abundances. Most species are hematophagous, feeding on all groups of vertebrates, including humans. In addition to being nuisance pests, they are able to transmit disease agents, with some viruses causing high morbidity and/or mortality in ruminants, horses and humans. Despite their impact on animal husbandry, public health and tourism, knowledge on the biology and ecology of culicoid biting midges and their interactions with ingested pathogens or symbiotic microorganisms is limited. Research is challenging due to unknown larval habitats, the insects' tiny size, the inability to establish and breed most species in the laboratory and the laborious maintenance of colonies of the few species that can be reared in the laboratory. Consequently, the natural transmission of pathogens has experimentally been demonstrated for few species while, for others, only indirect evidence of vector potential exists. Most experimental data are available for Culicoides sonorensis and C. nubeculosus, the only species kept in western-world insectaries. This contribution gives an overview on important biting midge vectors, transmitted viruses, culicoid-borne viral diseases and their epidemiologies and summarizes the little knowledge on interactions between biting midges, their microflora and culicoid-borne arboviruses.

Modeling the spatial distribution of Culicoides species (Diptera: Ceratopogonidae) as vectors of animal diseases in Ethiopia

Culicoides biting midges (Diptera: Ceratopogonidae) are the major vectors of bluetongue, Schmallenberg, and African horse sickness viruses. This study was conducted to survey Culicoides species in different parts of Ethiopia and to develop habitat suitability for the major Culicoides species in Ethiopia. Culicoides traps were set in different parts of the country from December 2018 to April 2021 using UV light Onderstepoort traps and the collected Culicoides were sorted to species level. To develop the species distribution model for the two predominant Culicoides species, namely Culicoides imicola and C. kingi, an ensemble modeling technique was used with the Biomod2 package of R software. KAPPA True skill statistics (TSS) and ROC curve were used to evaluate the accuracy of species distribution models. In the ensemble modeling, models which score TSS values greater than 0.8 were considered. Negative binomialregression models were used to evaluate the relationship between C. imicola and C. kingi catch and various environmental and climatic factors. During the study period, a total of 9148 Culicoides were collected from 66 trapping sites. Of the total 9148, 8576 of them belongs to seven species and the remaining 572 Culicoides were unidentified. The predominant species was C. imicola (52.8%), followed by C. kingi (23.6%). The abundance of these two species was highly influenced by the agro-ecological zone of the capture sites and the proximity of the capture sites to livestock farms. Climatic variables such as mean annual minimum and maximum temperature and mean annual rainfall were found to influence the catch of C. imicola at the different study sites. The ensemble model performed very well for both species with KAPPA (0.9), TSS (0.98), and ROC (0.999) for C. imicola and KAPPA (0.889), TSS (0.999), and ROC (0.999) for C. kingi. Culicoides imicola has a larger suitability range compared to C. kingi. The Great Rift Valley in Ethiopia, the southern and eastern parts of the country, and the areas along the Blue Nile and Lake Tana basins in northern Ethiopia were particularly suitable for C. imicola. High suitability for C. kingi was found in central Ethiopia and the Southern Nations, Nationalities and Peoples Region (SNNPR). The habitat suitability model developed here could help researchers better understand where the above vector-borne diseases are likely to occur and target surveillance to high-risk areas.

Dissecting the Species-Specific Virome in Culicoides of Thrace

Biting midges (Culicoides) are vectors of arboviruses of both veterinary and medical importance. The surge of emerging and reemerging vector-borne diseases and their expansion in geographical areas affected by climate change has increased the importance of understanding their capacity to contribute to novel and emerging infectious diseases. The study of Culicoides virome is the first step in the assessment of this potential. In this study, we analyzed the RNA virome of 10 Culicoides species within the geographical area of Thrace in the southeastern part of Europe, a crossing point between Asia and Europe and important path of various arboviruses, utilizing the Ion Torrent next-generation sequencing (NGS) platform and a custom bioinformatics pipeline based on TRINITY assembler and alignment algorithms. The analysis of the RNA virome of 10 Culicoides species resulted in the identification of the genomic signatures of 14 novel RNA viruses, including three fully assembled viruses and four segmented viruses with at least one segment fully assembled, most of which were significantly divergent from previously identified related viruses from the Solemoviridae, Phasmaviridae, Phenuiviridae, Reoviridae, Chuviridae, Partitiviridae, Orthomyxoviridae, Rhabdoviridae, and Flaviviridae families. Each Culicoides species carried a species-specific set of viruses, some of which are related to viruses from other insect vectors in the same area, contributing to the idea of a virus-carrier web within the ecosystem. The identified viruses not only expand our current knowledge on the virome of Culicoides but also set the basis of the genetic diversity of such viruses in the area of southeastern Europe. Furthermore, our study highlights that such metagenomic approaches should include as many species as possible of the local virus-carrier web that interact and share the virome of a geographical area.

Modeling the current distribution suitability and future dynamics of Culicoides imicola under climate change scenarios

BACKGROUND

African horse sickness, a transboundary and non-contagious arboviral infectious disease of equids, has spread without any warning from sub-Saharan Africa towards the Southeast Asian countries in 2020. It is imperative to predict the global distribution of Culicoides imicola (C. imicola), which was the main vector of African horse sickness virus.

METHODS

The occurrence records of C. imicola were mainly obtained from the published literature and the Global Biodiversity Information Facility database. The maximum entropy algorithm was used to model the current distribution suitability and future dynamics of C. imicola under climate change scenarios.

RESULTS

The modeling results showed that the currently suitable habitats for C. imicola were distributed in most of the southern part areas of America, southwestern Europe, most of Africa, the coastal areas of the Middle East, almost all regions of South Asia, southern China, a few countries in Southeast Asia, and the whole Australia. Our model also revealed the important environmental variables on the distribution of C. imicola were temperature seasonality, precipitation of coldest quarter, and mean temperature of wettest quarter. Representative Concentration Pathways (RCPs) is an assumption of possible greenhouse gases emissions in the future. Under future climate change scenarios, the area of habitat suitability increased and decreased with time, and RCP 8.5 in the 2070s gave the worst prediction. Moreover, the habitat suitability of C. imicola will likely expand to higher latitudes. The prediction of this study is of strategic significance for vector surveillance and the prevention of vector-borne diseases.

Identification of a novel bluetongue virus 1 specific B cell epitope using monoclonal antibodies against the VP2 protein

Bluetongue (BT) is a non-contact infectious disease caused by Bluetongue virus (BTV), which can be transmitted by vector insects such as Culicoides and Aedes mosquitoes. The BTV VP2 protein encoded by the L2 gene is located at the outermost layer of the virus particle, plays a key role on mediating the adsorption and entry of virus, and it is also a main antigenic protein widely used for vaccine development. In this study, the BTV1 VP2 gene was cloned into pFastBac™Dual vector, and expressed in insect Sf21 cells. Immunized mice with purified recombinant VP2 protein can induce higher levels of antibodies. Three anti BTV1 VP2 monoclonal antibodies (mAbs) were generated (17E9C6, 17E9C8, 17E9H12), and showed high specific reactivity with recombinant VP2 protein and inactivated BTV1 virus. Finally, a novel linear B-cell epitope ^296-KEPAD^-300 on recombinant VP2 protein was identified by using three mAbs react with a series of continue-truncated peptides. The results of this study may provide new information on the structure and function of BTV1 VP2 protein and lay a foundation for the development of BTV1 diagnostic and prophylactic methods.

Spatial and temporal distribution of Culicoides species in the New England region of New South Wales, Australia between 1990 and 2018

Culicoides are one of the smallest hematophagous flies measuring 1–5 mm in size with only females seeking blood for egg development. The present study investigated spatio-temporal distribution of Culicoides species trapped between 1990 and 2018 at 13 sites in the New England region of NSW, Australia using automated light traps. Trapping locations were divided into three subregions (tablelands, slopes and plains). Nineteen Culicoides species were identified. Culicoides marksi and C. austropalpalis were the most abundant and widespread species. Culicoides brevitarsis, the principal vector of livestock diseases in New South Wales comprised 2.9% of the total catch and was detected in 12 of the 13 locations in the study. Abundance as determined by Log₁₀Culicoides count per trapping event for the eight most abundant species did not vary significantly with season but trended towards higher counts in summer for C. marksi (P = 0.09) and C. austropalpalis (P = 0.05). Significant geographic variation in abundance was observed for C. marksi, C. austropalpalis and C. dycei with counts decreasing with increasing altitude from the plains to the slopes and tablelands. Culicoides victoriae exhibited the reverse trend in abundance (P = 0.08). Greater abundance during the warmer seasons and at lower altitudes for C. marksi and C. austropalpalis was indicative of temperature and rainfall dependence in this region with moderate summer dominance in rainfall. The Shannon-Wiener diversity index of species was higher on the tablelands (H = 1.59) than the slopes (H = 1.33) and plains (H = 1.08) with evenness indices of 0.62, 0.46 and 0.39 respectively. Culicoides species on the tablelands were more diverse than on the slopes and plains where C. marksi and C. austropalpalis dominated. The temporal and spatial variation in abundance, diversity and evenness of species reported in this diverse region of Australia provides additional insight into Culicoides as pests and disease vectors and may contribute to future modelling studies.

Scientific Opinion on the assessment of the control measures of the category A diseases of Animal Health Law: African Horse Sickness

EFSA received a mandate from the European Commission to assess the effectiveness of some of the control measures against diseases included in the Category A list according to Regulation (EU) 2016/429 on transmissible animal diseases (‘Animal Health Law’). This opinion belongs to a series of opinions where these control measures will be assessed, with this opinion covering the assessment of control measures for African Horse Sickness (AHS). In this opinion, EFSA and the AHAW Panel of experts review the effectiveness of: (i) clinical and laboratory sampling procedures, (ii) monitoring period and (iii) the minimum radius of the protection and surveillance zone, and the minimum duration of measures in these zones. The general methodology used for this series of opinions has been published elsewhere; nonetheless, specific details of the transmission kernels used for the assessment of the minimum radius of the protection and surveillance zones are shown. Several scenarios for which these control measures were assessed were designed and agreed prior to the start of the assessment. In summary, sampling procedures described in the diagnostic manual for AHS were considered efficient for all Equidae considering the high case fatality rate expected. The monitoring period (14 days) was assessed as effective in every scenario, except for those relating to the epidemiological enquiry where the risk manager should consider increasing the monitoring period, based on the awareness of keepers, environmental conditions and the vector abundance in the region. The current protection zone (100 km) comprises more than 95% of the infections from an affected establishment. Both the radius and duration of the zones could be reduced, based on local environmental conditions and the time of year of the first index case. Recommendations provided for each of the scenarios assessed aim to support the European Commission in the drafting of further pieces of legislation relating to AHS.

Inhibition of Orbivirus Replication by Aurintricarboxylic Acid

Bluetongue virus (BTV) and African horse sickness virus (AHSV) are vector-borne viruses belonging to the Orbivirus genus, which are transmitted between hosts primarily by biting midges of the genus Culicoides. With recent BTV and AHSV outbreaks causing epidemics and important economy losses, there is a pressing need for efficacious drugs to treat and control the spread of these infections. The polyanionic aromatic compound aurintricarboxylic acid (ATA) has been shown to have a broad-spectrum antiviral activity. Here, we evaluated ATA as a potential antiviral compound against Orbivirus infections in both mammalian and insect cells. Notably, ATA was able to prevent the replication of BTV and AHSV in both cell types in a time- and concentration-dependent manner. In addition, we evaluated the effect of ATA in vivo using a mouse model of infection. ATA did not protect mice against a lethal challenge with BTV or AHSV, most probably due to the in vivo effect of ATA on immune system regulation. Overall, these results demonstrate that ATA has inhibitory activity against Orbivirus replication in vitro, but further in vivo analysis will be required before considering it as a potential therapy for future clinical evaluation.

Heterologous Combination of ChAdOx1 and MVA Vectors Expressing Protein NS1 as Vaccination Strategy to Induce Durable and Cross-Protective CD8+ T Cell Immunity to Bluetongue Virus

The sequence of non-structural protein NS1 of bluetongue virus (BTV), which contains immunodominant CD8+ T cell epitopes, is highly conserved among BTV serotypes, and has therefore become a major tool in the development of a universal BTV vaccine. In this work, we have engineered multiserotype BTV vaccine candidates based on recombinant chimpanzee adenovirus (ChAdOx1) and modified vaccinia virus Ankara (MVA) vectors expressing the NS1 protein of BTV-4 or its truncated form NS1-Nt. A single dose of ChAdOx1-NS1 or ChAdOx1-NS1-Nt induced a moderate CD8+ T cell response and protected IFNAR(-/-) mice against a lethal dose of BTV-4/MOR09, a reassortant strain between BTV-1 and BTV-4, although the animals showed low viremia after infection. Furthermore, IFNAR(-/-) mice immunized with a single dose of ChAdOx1-NS1 were protected after challenge with a lethal dose of BTV-8 in absence of viremia nor clinical signs. Additionally, the heterologous prime-boost ChAdOx1/MVA expressing NS1 or NS1-Nt elicited a robust NS1 specific CD8+ T cell response and protected the animals against BTV-4/MOR09 even 16 weeks after immunization, with undetectable levels of viremia at any time after challenge. Subsequently, the best immunization strategy based on ChAdOx1/MVA-NS1 was assayed in sheep. Non-immunized animals presented fever and viremia levels up to 10⁴ PFU/mL after infection. In contrast, although viremia was detected in immunized sheep, the level of virus in blood was 100 times lower than in non-immunized animals in absence of clinical signs.