Carsten Kirkeby

Culicoids as Vectors of Schmallenberg Virus

To the Editor: In autumn 2011, an unidentified disease of livestock was reported on both sides of the Dutch–Germany border. By using metagenomics, the etiologic agent of this disease was identified as a novel orthobunyavirus and named Schmallenberg virus (SBV) (1). Other members of the genus Orthobunyavirus (e.g., Akabane virus) are widespread in Africa and Asia; biting midges (Culicoides spp.) and mosquitoes are responsible for transmitting these viruses. Hence, we reasonably assumed that European culicoids might be responsible for transmitting SBV within Europe. We present evidence that culicoids captured October 2011 in Denmark contained SBV RNA and most likely are vectors for this agent. In autumn 2011, culicoids were collected from several sites within Denmark. One site, a chicken farm in Hokkerup (Figure A1), was selected for study because of its location close (6 km) to the German border and proximity (<10 km) to an SBV-infected sheep farm in Germany, as reported on March 9, 2012, by the Friedrich Loeffler Institute surveillance website (www.fli.bund.de). The culicoids were collected during October 14–16 by using a Mosquito Magnet Independence trap (Mosquito Magnet, Lititz, PA, USA) baited with carbon dioxide and octenol. Midges were sorted manually into 91 specimens of the C. obsoletus group (comprising C. obsoletus, C. chiopterus, C. dewulfi, and C. scoticus) and 17 of the C. punctatus sensu stricto group, then stored at −20°C. Pools of culicoids were homogenized in water (100 µL) by using a 3-mm stainless steel bead (Dejay Distribution Ltd., Launceston, UK) in a TissueLyser II (QIAGEN, Hilden, Germany) for 1 min at 25 Hz (2). After homogenization, additional water (100 µL) was added to the samples, and then the mixture was centrifuged at 3,000 × g for 5 min. Nucleic acids were extracted from the supernatant (100 µL) by using a MagNA pure LC Total Nucleic Acid Isolation Kit on a MagNA pure LC (Roche Diagnostics, Basel, Switzerland) and eluted in water (50 µL). Two separate 1-step reverse transcription quantitative PCRs (RT-qPCRs), targeting the L segment and the S segment of SBV RNA, were performed according to protocols provided by the Friedrich Loeffler Institute in Germany (1) on the extracted nucleic acids by using a Mx3005p qPCR system (Agilent Technologies, Palo Alto, CA, USA). Another RT-qPCR targeting ruminant β-actin mRNA was performed as an internal endogenous control (3). Two of 22 pools tested strongly positive for the large (L) and small (S) segments of SBV RNA. Each positive sample was derived from 5 midges of the C. obsoletus group. One pool produced cycle threshold (Ct) values of 26.4 and 24.5 (in the L segment– and S segment–specific assays, respectively), whereas the second positive pool gave Ct values of 28.8 (L segment) and 27.6 (S segment). These pools were negative for the internal endogenous control that targeted the bovine/ovine β-actin mRNA. This result makes it unlikely that the detection of SBV RNA within the midges resulted from recent blood meals from infected animals remaining within the culicoids and suggests the virus has replicated within the midges. The PCR amplicons (145 bp; Figure) from the L segment–specific RT-qPCR were sequenced by using BigDye 1.1 chemistry on an ABI 3500 Genetic Analyzer (Applied Biosystems, Foster City, CA, USA). The sequences of 80 bp from the amplicons, excluding the primer sequences, had 100% identity with the expected region of the SBV segment L (1). Figure RNA extracted from pools of Culicoides obsoletus group midges was tested in 1-step reverse transcription quantitative PCRs (RT-qPCRs) for the Schmallenberg virus large segment, and the products were analyzed by agarose gel electrophoresis. Lanes 1–8, ... Reported Ct values generated by using the same assays from blood of naturally infected cattle were 24–35 (1). Usually, ≈100 µL of bovine/ovine blood is used for virus detection, whereas <1 µL of blood remains in a midge after a blood meal. This uptake of blood should therefore lead to a Ct value that is at least 6–7 units higher (≈100-fold lower level of RNA) when a single midge is tested by RT-qPCR (4). Thus, even if all 5 culicoids in a pool had recently taken a blood meal from a viremic animal, the Ct values observed here strongly suggest replication of SBV within the C. obsoletus group midges. However, in principle, other hosts of SBV could have a much higher level of viremia than cattle and could provide the levels of SBV RNA detected. C. punctatus s.s. midges cannot be ruled out as a possible vector of SBV because of the limited number of insects tested. Our study demonstrates the presence of SBV RNA in C. obsoletus group midges caught in Denmark during October 2011. The low Ct values (i.e., high SBV RNA levels) and the absence of ruminant β-actin mRNA in these samples strongly suggest that SBV replicates in these midges and hence that the C. obsoletus group midges are natural vectors for this virus.

Quantifying Dispersal of European Culicoides (Diptera: Ceratopogonidae) Vectors between Farms Using a Novel Mark-Release-Recapture Technique

Studying the dispersal of small flying insects such as Culicoides constitutes a great challenge due to huge population sizes and lack of a method to efficiently mark and objectively detect many specimens at a time. We here describe a novel mark-release-recapture method for Culicoides in the field using fluorescein isothiocyanate (FITC) as marking agent without anaesthesia. Using a plate scanner, this detection technique can be used to analyse thousands of individual Culicoides specimens per day at a reasonable cost. We marked and released an estimated 853 specimens of the Pulicaris group and 607 specimens of the Obsoletus group on a cattle farm in Denmark. An estimated 9,090 (8,918–9,260) Obsoletus group specimens and 14,272 (14,194–14,448) Pulicaris group specimens were captured in the surroundings and subsequently analysed. Two (0.3%) Obsoletus group specimens and 28 (4.6%) Pulicaris group specimens were recaptured. The two recaptured Obsoletus group specimens were caught at the release point on the night following release. Eight (29%) of the recaptured Pulicaris group specimens were caught at a pig farm 1,750 m upwind from the release point. Five of these were recaptured on the night following release and the three other were recaptured on the second night after release. This is the first time that movement of Culicoides vectors between farms in Europe has been directly quantified. The findings suggest an extensive and rapid exchange of disease vectors between farms. Rapid movement of vectors between neighboring farms may explain the the high rate of spatial spread of Schmallenberg and bluetongue virus (BTV) in northern Europe.

The range of attraction for light traps catching Culicoides biting midges (Diptera: Ceratopogonidae)

BackgroundCulicoides are vectors of e.g. bluetongue virus and Schmallenberg virus in northern Europe. Light trapping is an important tool for detecting the presence and quantifying the abundance of vectors in the field. Until now, few studies have investigated the range of attraction of light traps.MethodsHere we test a previously described mathematical model (Model I) and two novel models for the attraction of vectors to light traps (Model II and III). In Model I, Culicoides fly to the nearest trap from within a fixed range of attraction. In Model II Culicoides fly towards areas with greater light intensity, and in Model III Culicoides evaluate light sources in the field of view and fly towards the strongest. Model II and III incorporated the directionally dependent light field created around light traps with fluorescent light tubes. All three models were fitted to light trap collections obtained from two novel experimental setups in the field where traps were placed in different configurations.ResultsResults showed that overlapping ranges of attraction of neighboring traps extended the shared range of attraction. Model I did not fit data from any of the experimental setups. Model II could only fit data from one of the setups, while Model III fitted data from both experimental setups.ConclusionsThe model with the best fit, Model III, indicates that Culicoides continuously evaluate the light source direction and intensity. The maximum range of attraction of a single 4W CDC light trap was estimated to be approximately 15.25 meters. The attraction towards light traps is different from the attraction to host animals and thus light trap catches may not represent the vector species and numbers attracted to hosts.

Spatial abundance and clustering of Culicoides (Diptera: Ceratopogonidae) on a local scale

BackgroundBiting midges, Culicoides, of the Obsoletus group and the Pulicaris group have been involved in recent outbreaks of bluetongue virus and the former was also involved in the Schmallenberg virus outbreak in northern Europe.MethodsFor the first time, here we investigate the local abundance pattern of these two species groups in the field by intensive sampling with a grid of light traps on 16 catch nights. Neighboring trap catches can be spatially dependent on each other, hence we developed a conditional autoregressive (CAR) model framework to test a number of spatial and non-spatial covariates expected to affect Culicoides abundance.ResultsThe distance to sheep penned in the corner of the study field significantly increased the abundance level up to 200 meters away from the sheep. Spatial clustering was found to be significant but could not be explained by any known factors, and cluster locations shifted between catch nights. No significant temporal autocorrelation was detected. CAR models for both species groups identified a significant positive impact of humidity and significant negative impacts of precipitation and wind turbulence. Temperature was also found to be significant with a peak at just below 16 degrees Celcius. Surprisingly, there was a significant positive impact of wind speed. The CAR model for the Pulicaris group also identified a significant attraction to the smaller groups of sheep placed in the field. Furthermore, a large number of spatial covariates which were incorrectly found to be significant in ordinary regression models were not significant in the CAR models. The 95% C.I. on the prediction estimates ranged from 20.4% to 304.8%, underlining the difficulties of predicting the abundance of Culicoides.ConclusionsWe found that significant spatial clusters of Culicoides moved around in a dynamic pattern varying between catch nights. This conforms with the modeling but was not explained by any of the tested covariates. The mean abundance within these clusters was up to 11 times higher for the Obsoletus group and 4 times higher for the Pulicaris group compared to the rest of the field.

Observations of movement dynamics of flying insects using high resolution lidar

Insects are fundamental to ecosystem functioning and biodiversity, yet the study of insect movement, dispersal and activity patterns remains a challenge. Here we present results from a novel high resolution laser-radar (lidar) system for quantifying flying insect abundance recorded during one summer night in Sweden. We compare lidar recordings with data from a light trap deployed alongside the lidar. A total of 22808 insect were recorded, and the relative temporal quantities measured matched the quantities recorded with the light trap within a radius of 5 m. Lidar records showed that small insects (wing size <2.5 mm2 in cross-section) moved across the field and clustered near the light trap around 22:00 local time, while larger insects (wing size >2.5 mm2 in cross-section) were most abundant near the lidar beam before 22:00 and then moved towards the light trap between 22:00 and 23:30. We could distinguish three insect clusters based on morphology and found that two contained insects predominantly recorded above the field in the evening, whereas the third was formed by insects near the forest at around 21:30. Together our results demonstrate the capability of lidar for distinguishing different types of insect during flight and quantifying their movements.

Simulating the Epidemiological and Economic Impact of Paratuberculosis Control Actions in Dairy Cattle

We describe a new mechanistic bioeconomic model for simulating the spread of Mycobacterium avium subsp. paratuberculosis (MAP) within a dairy cattle herd. The model includes age-dependent susceptibility for infection; age-dependent sensitivity for detection; environmental MAP build up in five separate areas of the farm; in utero infection; infection via colostrum and waste milk, and it allows for realistic culling (i.e., due to other diseases) by including a ranking system. We calibrated the model using a unique dataset from Denmark, including 102 random farms with no control actions against spread of MAP. Likewise, four control actions recommended in the Danish MAP control program were implemented in the model based on reported management strategies in Danish dairy herds in a MAP control scheme. We tested the model parameterization in a sensitivity analysis. We show that a test-and-cull strategy is on average the most cost-effective solution to decrease the prevalence and increase the total net revenue on a farm with low hygiene, but not more profitable than no control strategy on a farm with average hygiene. Although it is possible to eradicate MAP from the farm by implementing all four control actions from the Danish MAP control program, it was not economically attractive since the expenses for the control actions outweigh the benefits. Furthermore, the three most popular control actions against the spread of MAP on the farm were found to be costly and inefficient in lowering the prevalence when used independently.

Culicoides (Avaritia) gornostaevae Mirzaeva, 1984 (Diptera: Ceratopogonidae) - a possible vector species of the Obsoletus group new to the European fauna

BackgroundCulicoides gornostaevae Mirzaeva, 1984, known previously only from Siberia, is a boreal species included into the Obsoletus group of Culicoides sg. Avaritia. Members of the subgenus can act as vectors of various diseases. In Europe they are involved in the transmission of the Schmallenberg virus and bluetongue virus.FindingsCulicoides gornostaevae Mirzaeva, 1984 is reported for the first time in Europe with new country records from Norway, Poland and Sweden.ConclusionsCulicoides gornostaevae Mirzaeva, 1984 has not been previously mentioned from Europe, even though there has been an extensive monitoring of Culicoides species during the last decades. Most probably this species has been notoriously overlooked in the materials, because of the problems with identification of the females of the subgenus Avaritia. Similar to other species of the Obsoletus group, C. gornostaevae should be regarded as a possible vector for Schmallenberg and bluetongue virus.

Models to Estimate Lactation Curves of Milk Yield and Somatic Cell Count in Dairy Cows at the Herd Level for the Use in Simulations and Predictive Models

Typically, central milk recording data from dairy herds are recorded less than monthly. Over-fitting early in lactation periods is a challenge, which we explored in different ways by reducing the number of parameters needed to describe the milk yield and somatic cell count of individual cows. Furthermore, we investigated how the parameters of lactation models correlate between parities and from dam to offspring. The aim of the study was to provide simple and robust models for cow level milk yield and somatic cell count for fitting to sparse data to parameterize herd- and cow-specific simulation of dairy herds. Data from 610 Danish Holstein herds were used to determine parity traits in milk production regarding milk yield and somatic cell count of individual cows. Parity was stratified in first, second, and third and higher for milk, and first to sixth and higher for somatic cell count. Fitting of herd level parameters allowed for cow level lactation curves with three, two, or one parameters per lactation. Correlations of milk yield and somatic cell count were estimated between lactations and between dam and offspring. The shape of the lactation curves varied markedly between farms. The correlation between lactations for milk yield and somatic cell count was 0.2–0.6 and significant on more than 95% of farms. The variation in the daily milk yield was observed to be a source of variation to the somatic cell count, and the total somatic cell count was less correlated with the milk production than somatic cells per milliliter. A positive correlation was found between relative levels of the total somatic cell count and the milk yield. The variation of lactation and somatic cell count curves between farms highlights the importance of a herd level approach. The one-parameter per cow model using a herd level curve allows for estimating the cow production level from first the recording in the parity, while a two-parameter model requires more recordings for a credible estimate, but may more precisely predict persistence, and given the independence of parameters, these can be easily drawn for use in simulation models. We also conclude that using total somatic cell count may stabilize models, and therefore, the dilution factor is of importance in Danish Holstein.

Vector-borne disease surveillance in livestock populations: A critical review of literature recommendations and implemented surveillance (BTV-8) in five European countries

Preparedness against vector-borne threats depends on the existence of a long-term, sustainable surveillance of vector-borne disease and their relevant vectors. This work reviewed the availability of such surveillance systems in five European countries (Denmark, France, The Netherlands, Sweden and United Kingdom, part of the CoVetLab network). A qualitative assessment was then performed focusing on surveillance directed particularly to BTV-8. Information regarding surveillance activities were reviewed for the years 2008 and 2012. The results were then complemented with a critical scoping review of the literature aimed at identifying disease surveillance strategies and methods that are currently suggested as best suited to target vector-borne diseases in order to guide future development of surveillance in the countries in question. Passive surveillance was found to be efficient for early detection of diseases during the early phase of introduction into a free country. However, its value diminished once the disease has been established in a territory. Detection of emerging diseases was found to be very context and area specific, and thus active surveillance designs need to take the available epidemiological, ecological and entomological information into account. This was demonstrated by the effectiveness of the bulk milk surveillance in detecting the first case in Sweden, highlighting the need for output based standards to allow the most effective, context dependent, surveillance strategies to be used. Preparedness was of fundamental importance in determining the timeliness of detection and control in each country and that this in turn was heavily influenced by knowledge of emerging diseases in neighboring countries. Therefore it is crucial to share information on outbreaks between researchers and decision-makers and across borders continuously in order to react timely in case of an outbreak. Furthermore, timely reaction to an outbreak was heavily influenced by availability of control measures (vaccines), which is also strengthened if knowledge is shared quickly between countries. The assessment of the bluetongue surveillance in the affected countries showed that the degree of voluntary engagement varied, and that it is important to engage the public by general awareness and dissemination of results. The degree of engagement will also aid in establishing a passive surveillance system.

Spatio-temporal optimization of sampling for bluetongue vectors (Culicoides) near grazing livestock

BackgroundEstimating the abundance of Culicoides using light traps is influenced by a large variation in abundance in time and place. This study investigates the optimal trapping strategy to estimate the abundance or presence/absence of Culicoides on a field with grazing animals. We used 45 light traps to sample specimens from the Culicoides obsoletus species complex on a 14 hectare field during 16 nights in 2009.FindingsThe large number of traps and catch nights enabled us to simulate a series of samples consisting of different numbers of traps (1-15) on each night. We also varied the number of catch nights when simulating the sampling, and sampled with increasing minimum distances between traps. We used resampling to generate a distribution of different mean and median abundance in each sample. Finally, we used the hypergeometric distribution to estimate the probability of falsely detecting absence of vectors on the field. The variation in the estimated abundance decreased steeply when using up to six traps, and was less pronounced when using more traps, although no clear cutoff was found.ConclusionsDespite spatial clustering in vector abundance, we found no effect of increasing the distance between traps. We found that 18 traps were generally required to reach 90% probability of a true positive catch when sampling just one night. But when sampling over two nights the same probability level was obtained with just three traps per night. The results are useful for the design of vector monitoring programmes on fields with grazing animals.