Spatio-temporal modelling of Culicoides Latreille (Diptera: Ceratopogonidae) populations on Reunion Island (Indian Ocean).

BACKGROUND: Reunion Island regularly faces outbreaks of bluetongue and epizootic hemorrhagic diseases, two insect-borne orbiviral diseases of ruminants. Hematophagous midges of the genus Culicoides (Diptera: Ceratopogonidae) are the vectors of bluetongue (BTV) and epizootic hemorrhagic disease (EHDV) viruses. In a previous study, statistical models based on environmental and meteorological data were developed for the five Culicoides species present in the island to provide a better understanding of their ecology and predict their presence and abundance. The purpose of this study was to couple these statistical models with a Geographic Information System (GIS) to produce dynamic maps of the distribution of Culicoides throughout the island. METHODS: Based on meteorological data from ground weather stations and satellite-derived environmental data, the abundance of each of the five Culicoides species was estimated for the 2214 husbandry locations on the island for the period ranging from February 2016 to June 2018. A large-scale Culicoides sampling campaign including 100 farms was carried out in March 2018 to validate the model. RESULTS: According to the model predictions, no husbandry location was free of Culicoides throughout the study period. The five Culicoides species were present on average in 57.0% of the husbandry locations for C. bolitinos Meiswinkel, 40.7% for C. enderleini Cornet & Brunhes, 26.5% for C. grahamii Austen, 87.1% for C. imicola Kieffer and 91.8% for C. kibatiensis Goetghebuer. The models also showed high seasonal variations in their distribution. During the validation process, predictions were acceptable for C. bolitinos, C. enderleini and C. kibatiensis, with normalized root mean square errors (NRMSE) of 15.4%, 13.6% and 16.5%, respectively. The NRMSE was 27.4% for C. grahamii. For C. imicola, the NRMSE was acceptable (11.9%) considering all husbandry locations except in two specific areas, the Cirque de Salazie-an inner mountainous part of the island-and the sea edge, where the model overestimated its abundance. CONCLUSIONS: Our model provides, for the first time to our knowledge, an operational tool to better understand and predict the distribution of Culicoides in Reunion Island. As it predicts a wide spatial distribution of the five Culicoides species throughout the year and taking into consideration their vector competence, our results suggest that BTV and EHDV can circulate continuously on the island. As further actions, our model could be coupled with an epidemiological model of BTV and EHDV transmission to improve risk assessment of Culicoides-borne diseases on the island.

Modelling temporal dynamics of Culicoides Latreille (Diptera: Ceratopogonidae) populations on Reunion Island (Indian Ocean), vectors of viruses of veterinary importance.

BACKGROUND: Reunion Island regularly faces outbreaks of epizootic haemorrhagic disease (EHD) and bluetongue (BT), two viral diseases transmitted by haematophagous midges of the genus Culicoides (Diptera: Ceratopogonidae) to ruminants. To date, five species of Culicoides are recorded in Reunion Island in which the first two are proven vector species: Culicoides bolitinos, C. imicola, C. enderleini, C. grahamii and C. kibatiensis. Meteorological and environmental factors can severely constrain Culicoides populations and activities and thereby affect dispersion and intensity of transmission of Culicoides-borne viruses. The aim of this study was to describe and predict the temporal dynamics of all Culicoides species present in Reunion Island. METHODS: Between 2016 and 2018, 55 biweekly Culicoides catches using Onderstepoort Veterinary Institute traps were set up in 11 sites. A hurdle model (i.e. a presence/absence model combined with an abundance model) was developed for each species in order to determine meteorological and environmental drivers of presence and abundance of Culicoides. RESULTS: Abundance displayed very strong heterogeneity between sites. Average Culicoides catch per site per night ranged from 4 to 45,875 individuals. Culicoides imicola was dominant at low altitude and C. kibatiensis at high altitude. A marked seasonality was observed for the three other species with annual variations. Twelve groups of variables were tested. It was found that presence and/or abundance of all five Culicoides species were driven by common parameters: rain, temperature, vegetation index, forested environment and host density. Other parameters such as wind speed and farm building opening size governed abundance level of some species. In addition, Culicoides populations were also affected by meteorological parameters and/or vegetation index with different lags of time, suggesting an impact on immature stages. Taking into account all the parameters for the final hurdle model, the error rate by Normalized Root mean Square Error ranged from 4.4 to 8.5%. CONCLUSIONS: To our knowledge, this is the first study to model Culicoides population dynamics in Reunion Island. In the absence of vaccination and vector control strategies, determining periods of high abundance of Culicoides is a crucial first step towards identifying periods at high risk of transmission for the two economically important viruses they transmit.