Estimating the temporal and spatial risk of bluetongue related to the incursion of infected vectors into Switzerland
© Racloz et al; licensee BioMed Central Ltd. 2008
Received: 05 February 2008
Accepted: 15 October 2008
Published: 15 October 2008
The design of veterinary and public health surveillance systems has been improved by the ability to combine Geographical Information Systems (GIS), mathematical models and up to date epidemiological knowledge. In Switzerland, an early warning system was developed for detecting the incursion of the bluetongue disease virus (BT) and to monitor the frequency of its vectors. Based on data generated by this surveillance system, GIS and transmission models were used in order to determine suitable seasonal vector habitat locations and risk periods for a larger and more targeted surveillance program.
Combined thematic maps of temperature, humidity and altitude were created to visualize the association with Culicoides vector habitat locations. Additional monthly maps of estimated basic reproduction number transmission rates (R0) were created in order to highlight areas of Switzerland prone to higher BT outbreaks in relation to both vector activity and transmission levels. The maps revealed several foci of higher risk areas, especially in northern parts of Switzerland, suitable for both vector presence and vector activity for 2006.
Results showed a variation of R0 values comparing 2005 and 2006 yet suggested that Switzerland was at risk of an outbreak of BT, especially if the incursion arrived in a suitable vector activity period. Since the time of conducting these analyses, this suitability has proved to be the case with the recent outbreaks of BT in northern Switzerland.
Our results stress the importance of environmental factors and their effect on the dynamics of a vector-borne disease. In this case, results of this model were used as input parameters for creating a national targeted surveillance program tailored to both the spatial and the temporal aspect of the disease and its vectors. In this manner, financial and logistic resources can be used in an optimal way through seasonally and geographically adjusted surveillance efforts. This model can serve as a tool for other vector-borne diseases including human zoonotic vectors which are likely to spread into Europe.
Bluetongue disease virus (BT) is a vector-borne, infectious but non-contagious animal pathogen. This emerging disease affects all ruminants and has been responsible for an unprecedented continuing European epidemic which has been occurring for the past decade . Belonging to the Orbivirus genus and Reoviridae family, there are currently 24 recognized serotypes transmitted globally by a multitude of Culicoides midge species, each with their own habitat preferences and geographical distribution albeit with some overlapping occurrence. Several serotypes, mostly affecting sheep, have been circulating in the Balkan and Mediterranean areas since the late 90's, which could be predicted by the advance of its Old World vector C. imicola. Yet, an outbreak of BT serotype 8 (BTV-8) in 2006, which was last recorded in the African and the Caribbean region , suddenly occurred in Northern Europe, an area previously free of BT infection . Preceding this event, outbreaks had been reported on a regular seasonal basis in southern Europe, mainly the Mediterranean region involving several serotypes namely BTV -1, -2, -4, -6, -9 and -16 .
At the time of writing, the first occurrence of BTV-8 in Switzerland was reported in northern Switzerland in the canton of Basel-City in late October 2006. This was shortly followed by cases in the canton of Solothurn, and Basel-Land respectively. Further cases of BTV-8 were consecutively detected in Basel-Land and Solothurn. In January 2007, a surprising case was discovered in the southern canton of Valais. In the first outbreaks of BT ever recorded in Switzerland dating from October 2006 to February 2007, a total of 12 cattle and two goats tested positive for BTV-8 originating from seven different farms. Although no firm conclusion have yet arisen as to the cause of these cases, wind direction patterns along with temperature records of the affected areas suggest likely intrusion of infected vectors originating from the surrounding BT-reporting areas.
Due to the dynamics of the pathogen, combined with the fact that the geographical limits of other vector-borne diseases are also expanding, Switzerland conducted a nation-wide survey in 2003 to determine the status of BT disease and the presence of its vectors . Additionally since 2005 , studies conducted have recorded the presence of the following potential BT vectors: C. obsoletus (Meigen), C. pulicaris, C. scoticus and C. dewulfi, species which have all been implicated in BT outbreaks in several countries located in Europe and the Balkan areas [3, 6, 7]
The aim of this study was to combine GIS maps of data collected in the field with maps of predicted basic reproduction number estimates for potential BT outbreaks in Switzerland, in order to explore the spatial and temporal regions more prone to 1) the establishment of important vector populations, and 2) enabling the spread of the disease due to the nature of geographical and climatic features.
Temperature variability and R0calculations
The combined results of the suitability maps and the R0 maps for the months of January, July, September and October of 2006 are shown in Figure 4. The suitability maps (Fig 4, left) highlight various localized microclimates occurring in Switzerland and confirm the role of the Alps in separating the northern and southern parts of the country. The maps (Fig. 4, left) showed that some areas remain suitable for Culicoides survival throughout the year despite colder temperatures, agreeing with virus over-wintering theories. In contrast though, the R0 maps showed negligible risk of BT spread in the colder months of 2006 (Fig. 4, middle), due to low vector activity rates. However when combined, the risk maps (Fig 4, right) showed that the risk of an outbreak and spread of BT in Switzerland was not negligible even in the colder seasons due to the continued suitability for vector survival in certain areas, noting that Culicoides vectors have also been recorded within farms where local temperatures (personal communication) would be efficient for virus transmission.
The difference between the vector suitability maps and the calculated transmission values are most distinct when comparing the maps for July and September of 2006. July is less suitable for vector reproduction due to the high temperatures recorded in that month which could, among other factors, affect larval death rates through desiccation, yet R0 values are higher than in the month of September, due to the shortened extrinsic incubation period. On the other hand, the maps for September show more suitable areas for vector activity through the effect of less extreme temperatures.
Our results showed that the months of July and September 2006 were the most suitable period in Switzerland for vector presence in terms of climatic conditions (suitability maps), while the highest R0 value indicating vector activity and therefore spread, occurred in June. This is an interesting finding in relation to the BT incursion and subsequent outbreak in Northern Europe, which was first detected in August 2006, possibly indicating similar conditions in the affected countries, such as The Netherlands, Belgium and Germany. The R0 values estimated in this study were used as a risk indicator levels for vector activity in different geographic areas (Fig. 4 middle) in terms of BT transmission rates. They were then combined with maps indicating areas with predicted higher vector suitability likelihood (Fig. 4 left). The actual cases of BT correspond to the maps generated for 2006, yet due to the data originating from this year, cannot be replicated for future outbreaks of BT. Similarly to the mapping methods for climate, R0 values were classified into 4 risk levels, signifying most rapid spread of disease in areas presenting higher range of R0 values.
The basic reproduction number (R0) for vector-borne diseases is a more complex number to calculate due to the influence of seasonal fluctuations , local climate and environmental features as well as the abundance of breeding sites available near hosts which affect vector dynamics . The transmission rates of the disease will also change depending on temperature factors affecting vector to host or host to vector interaction, along with the extrinsic incubation time, biting rates and vector mortality rates . In this study, vector density numbers along with temperature values were used from field data collected throughout the project in order to produce R0 values specific to the areas and time frame studied. The values which were most affected in terms of location and time period were those of vector density, vector death rates and extrinsic incubation rates. These were mainly affected by differences in temperature and humidity levels which in turn are very specific to the various microclimates present in Switzerland targeted in the sampling. Concerning the effect of host availability on transmission rates, since data was taken from sentinel farms where minimal transhumance activity took place, cattle density was considered as stable, with a Swiss average size of 30 cattle per farm.
A recent study by Gubbins et al.  assessed the risk of BT in the United Kingdom using the basic reproduction number R0. Similarly to their findings, the R0 values were highest when temperatures were between 15°C–25°C. Whilst the mentioned study took into account both the ratio of vectors to cattle and sheep in the formula, our study only considered the relation of BTV-8 in cattle farms and not sheep due to much larger cattle population of larger ruminants present in Switzerland. Additionally in financial terms, the cattle industry and more specifically the milk sector is an important part of the Swiss economy http://www.blw.admin.ch.
For Switzerland, the peak of R0 can be explained by the fact that the largest amount of Culicoides midges in 2006 were recorded in June, probably due to optimal breeding and hatching conditions as May was a mild and humid month. July held the record for maximum temperature, which leads to higher vector activity and successful virus transmission. However high temperatures also increase vector mortality rate and thereby lower R0 values for this month. Due to the very different meteorological patterns in the past three years, we would therefore expect that maps and R0 for 2007 (comparable to conditions in 2003) to be quite different compared to 2005 and 2006 in terms of increased vector activity and BT transmission rates. http://www.meteosuisse.admin.ch/web/fr/climat/climat_des_derniers_mois/Jahresbilanzen.html.
Affected by the initial high temperatures recorded in July 2006, transmission values as well as vector density decreased significantly following a notably milder August. In terms of the Northern European BTV8 outbreak data, this would match reports stating that the maximum number of cases occurred in October which would originate from a high number of vectors present a few weeks previously. Due to the nature of Culicoides development, it has been suggested that cases occur at a time lag of circa four weeks from peak vector density periods, which corroborate the evidence from the trapping data and the transmission values . As seen in The Netherlands which recorded high temperatures in July, this climatic event could have primed the vectors in terms of their competence and capacity levels, translating in an effective transmission period, as observed in many other affected countries . As mentioned in a study conducted by Murray , rapid changes in climatic conditions can affect vector population age structure, along with vector density and alter transmission rates of disease.
A limiting factor in our model is that vector density data resulted from trapping sites representing high risk areas for entomological surveillance, located in southern parts of Switzerland. In these locations, high numbers of vectors were expected. Therefore, the transmission values may be overestimated in some parts of the country, yet due to the targeted nature of setting up sentinel herds in high risk areas this limitation was acceptable in this scenario.
Another factor to consider is that only the dynamics between BTV-8 and the vectors belonging to the Culicoides obsoletus group were studied, and considering that a majority of the biological parameters used in calculating values for R0 were derived from the available literature, this should be taken into account when applying such methods to different countries, especially when micro-climates or other BT vectors are present. Outbreaks involving other BT serotypes and their dynamics in Culicoides vectors produce different R0 values . They may also have a different preference for distinct geographic and climatic conditions. Such differences have been described for the behaviour of BTV-2 and Culicoides imicola in southern France , and the role of C. imicola in South Africa .
The recent cases of BT in Switzerland in the month of October 2007 occurred in the northern part of the country, where temperatures for the affected region were similar to those of 2006, and the number of cattle affected on the seven farms correspond to the R0 figures calculated for that area and time period. Out of a total of 608 susceptible cattle, 16 sheep and six goats, originating from the farms located in Basel-City, Solothurn, Basel-Land, Valais and Jura, BT was detected in a total of 12 cattle and 2 goats, with prevalence rates ranging in chronological terms of 10.3% to 0.6% in cattle and 33% in goats http://www.bluetongue.ch. The R0 values provided by the model were similar to the actual prevalence rates for the time of the year, with the exception of goats, which were not included in the study.
Our findings highlight the potential for establishing a flexible surveillance system taking into account environmental factors. In a targeted surveillance system, this could mean increased serological testing during a specific warmer period or in specific geographical areas. Given the scarcity of epidemiological data available for Switzerland concerning BT cases prior to this study, the creation of thematic and risk maps on a monthly and annual basis, illustrated the variability in the behaviour of vector borne diseases and the possible consequences of virus introduction. It also provided the basis for creating a surveillance system targeting higher BT-risk regions and months. In the case of Switzerland, the maps and R0 values were used as input parameters for the creation of a BT surveillance Scenario Tree model , with the aim of comparing the efficacy of alternative surveillance system designs. A risk-based surveillance program was implemented in July 2007 consisting of three surveillance system components; serological bulk milk testing of 200 sentinel herds located in areas considered of higher risk to BT occurrence, as well as clinical surveillance programs for cattle and sheep farmers throughout the country.
GIS mapping techniques combined with statistical and mathematical models can help improve disease surveillance and control methods by providing a basis for targeting monitoring efforts. Targeted surveillance in this scenario meant focusing financial and monitoring efforts to pre-defined areas originating from the maps created in this study.
An advantage of GIS methods includes the ability to improve prediction maps once more comprehensive field data has been collected, and adjust surveillance efforts in a timely and accurate manner. Flexible surveillance programs should be used in order to attribute financial and human resources to high-risk areas considering temporal and spatial factors.
In this study, maps were used to highlight areas in Switzerland which presented higher likely risk for vector habitat, which incorporated climatic elements, (Figure 4. left margin), combined with maps highlighting different vector activity levels based on R0 values for BT (Figure 4. middle column), in order to produce a set of final combination maps representing areas in Switzerland most likely to have outbreaks of BT due to vector presence in relation to vector activity levels on a monthly basis for 2006. As mentioned, the different risk zone distribution originating from the final combination maps were used as input parameters in a Scenario Tree model to aid in decision making processes concerning BT surveillance in Switzerland http://www.bluetongue.ch.
Approaches using GIS maps, and the combination of spatial analysis and mathematical modelling as predictive tools have been used in other countries such as Italy, Spain, and France concerning BT and its vectors [22–24], and stress the advantages of using different technological methods in supporting surveillance efforts, as well as for other purposes aiding in the improvement of health and prevention of disease both in public and veterinary terms.
Parameters and values used for calculation of R0
Vector density (average)
Vector biting rate
Transmission from vector to cattle
%infectious bite/infected host
Transmission from cattle to vector
Recovery rate of cattle
Cattle death rate
Vector death rate
Extrinsic incubation period
In terms of vector-borne diseases, the basic reproduction number (R0), is defined as the number of new infections that would result from the introduction of a single infectious vector specimen into a completely susceptible/naive population of hosts . The R0equation is made of three components (Fig. 2). Component (I) consists of the following symbols: vector density (m), derived from Culicoides catches in the national entomological surveillance, vector biting rates (a), transmission rates from vector to cattle (b1), cattle recovery rates (r) and cattle death rates (λ). This part of the equation involves the stage of vectors infecting their hosts. The second section (II) includes the extrinsic incubation period (τ), vector death rates (μ) and signifies the amount of time the virus is developing inside the vector taking into account the lifespan of the vector. Finally, the third component (III), which relates to the percentage of successful infectious bites per infected host, includes the vector biting rates (a), the transmission rates from cattle to vector (b2) and the vector death rates (μ).
Through the variation of extrinsic incubation periods relating to temperature and humidity levels, as well as the collected vector abundance data derived from field data, R0 values were estimated for each month for the years 2005 and 2006. The R0 values calculated were then categorised into 4 levels (high, medium-high, medium-low and low), symbolising vector activity levels in terms of BT virus transmission, for incorporation into GIS maps (Fig. 4. middle). Additionally, the R0values were plotted against monthly minimum, mean and maximum temperatures for the area where entomological trapping occurred (Fig. 3a and 3b).
Separate thematic maps were created using ArcGis (Version 8.3, Environmental Systems Research Institute, Inc.) for monthly mean temperature, altitude and humidity for the years of 2005 and 2006 using data from 50 meteorological measuring stations provided by the Swiss Meteorological Office, as previously described . The aim was to create combined monthly vector suitability maps using these parameters to visualize the variation in potential risk areas during each season. Once monthly datasets were incorporated into the map, smoothing out was performed through kriging, apart from the altitude map which was derived from an elevation model. Suitability categories, based on Culicoides obsoletus group biology and habitat data [1, 6, 7, 26–28], were used to reclassify the output values, in order to grade all monthly maps on a standard scale, as done for the R0 values mentioned above. The 'environmental envelope' of the Obsoletus group of Culicoides was concentrated upon in contrast to the classical Old World vector C. imicola, due the fact that the former is the most abundant group caught in the Swiss entomology surveillance program  and has been shown to transmit BT virus in other countries . The maps were then layered together using the addition function in the raster map calculator which created a single combined suitability map for each month.
In a similar fashion, R0 values for each month were assigned to the 50 geographical coordinates based on the location of the meteorological stations, and kriging was performed in order to smooth out the data. The R0 values were divided into four categories, and reclassified to share a standardized scale for each month, as well as being comparable to the suitability map scales. The two sets of maps were then combined by adding the respective layers for each month through the raster map function to produce final combination maps incorporating both spatial and temporal factors (Fig. 4 right).
VR would like to acknowledge Dr. Patrick Presi and Dr. Heinzpeter Schwermer for their help in GIS matters, as well as Ms Monika Kuhn and Dr. Simona Casati for their contribution to the entomological field work.
- Purse BV, Mellor PS, Rogers DJ, Samuel AR, Mertens PP, Baylis M: Climate change and the recent emergence of bluetongue in Europe. Nat Rev Microbiol. 2005, 3: 171-181. 10.1038/nrmicro1090.PubMedView ArticleGoogle Scholar
- Mo C, Thompson LH, Homan EJ, Oveido MT, Greiner EC GJ, MR S: Bluetongue virus isolations from vectors and ruminants in Central America and the Caribbean. Interamerican Bluetongue Team. Am J Vet Res. 1994, 55: 211-215.PubMedGoogle Scholar
- Mehlhorn H, Walldorf V, Klimpel S, Jahn B, Jaeger F, Eschweiler J, Hoffmann B, Beer M: First occurrence of Culicoides obsoletus-transmitted Bluetongue virus epidemic in Central Europe. Parasitol Res. 2007, 101: 219-228. 10.1007/s00436-007-0519-6.PubMedView ArticleGoogle Scholar
- Cagienard A, Griot C, Mellor PS, Denison E, Stark KD: Bluetongue vector species of Culicoides in Switzerland. Med Vet Entomol. 2006, 20: 239-247. 10.1111/j.1365-2915.2006.00621.x.PubMedView ArticleGoogle Scholar
- Casati S, Racloz V, Delécolle JC, Kuhn M, Mathis A, Griot C, Stärk KDC, Vanzetti T: Entomological monitoring of bluetongue vectors at sentinel sites in southern Switzerland: an investigation on the Culicoides population composition. Medical and Veterinary Entomology.Google Scholar
- De Liberato C, Scavia G, Lorenzetti R, Scaramozzino P, Amaddeo D, Cardeti G, Scicluna M, Ferrari G, Autorino GL: Identification of Culicoides obsoletus (Diptera: Ceratopogonidae) as a vector of bluetongue virus in central Italy. Vet Rec. 2005, 156: 301-304.Google Scholar
- Purse BV, Nedelchev N, Georgiev G, Veleva E, Boorman J, Denison E, Veronesi E, Carpenter S, Baylis M, Mellor PS: Spatial and temporal distribution of bluetongue and its Culicoides vectors in Bulgaria. Med Vet Entomol. 2006, 20: 335-344. 10.1111/j.1365-2915.2006.00636.x.PubMedView ArticleGoogle Scholar
- Racloz V, Griot C, Stark KD: Sentinel surveillance systems with special focus on vector-borne diseases. Anim Health Res Rev. 2006, 7: 71-79. 10.1017/S1466252307001120.PubMedView ArticleGoogle Scholar
- Racloz V, Straver R, Kuhn M, Thur B, Vanzetti T, Stark KD, Griot C, Cagienard A: Establishment of an early warning system against bluetongue virus in Switzerland. Schweiz Arch Tierheilkd. 2006, 148: 593-598. 10.1024/0036-72220.127.116.113.PubMedView ArticleGoogle Scholar
- Racloz V, Presi P, Vounatsou P, Schwermer H, Casati S, Vanzetti T, Griot C, Stärk KDC: Use of mapping and statistical modeling for the prediction of Bluetongue occurrence in Switzerland based on vector biology. Veterinaria Italiana. 2007, 43 (3): 513-518.PubMedGoogle Scholar
- Roberts MG, Heesterbeek JA: A new method for estimating the effort required to control an infectious disease. Proc Biol Sci. 2003, 270: 1359-1364. 10.1098/rspb.2003.2339.PubMed CentralPubMedView ArticleGoogle Scholar
- Smith DL, McKenzie FE, Snow RW, Hay SI: Revisiting the basic reproductive number for malaria and its implications for malaria control. PLoS Biol. 2007, 5: e42-10.1371/journal.pbio.0050042.PubMed CentralPubMedView ArticleGoogle Scholar
- Wonham MJ, de-Camino-Beck T, Lewis MA: An epidemiological model for West Nile virus: invasion analysis and control applications. Proc Biol Sci. 2004, 271: 501-507. 10.1098/rspb.2003.2608.PubMed CentralPubMedView ArticleGoogle Scholar
- Lord CC, Woolhouse ME, Heesterbeek JA, Mellor PS: Vector-borne diseases and the basic reproduction number: a case study of African horse sickness. Med Vet Entomol. 1996, 10: 19-28. 10.1111/j.1365-2915.1996.tb00077.x.PubMedView ArticleGoogle Scholar
- Gubbins S, Carpenter S, Baylis M, Wood JL, Mellor PS: Assessing the risk of bluetongue to UK livestock: uncertainty and sensitivity analyses of a temperature-dependent model for the basic reproduction number. J R Soc Interface. 2008, 5 (20): 363-71. 10.1098/rsif.2007.1110.PubMed CentralPubMedView ArticleGoogle Scholar
- Bacaër N: Approximation of the basic reproduction number R0 for vector-borne diseases with a periodic vector population. Bull Math Biol. 2007, 69: 1067-1091. 10.1007/s11538-006-9166-9.PubMedView ArticleGoogle Scholar
- Anonymous: Report on Epidemiological analysis of the 2006 bluetongue virus serotype 8 epidemic in north-western Europe: provisional findings through 31 January 2007. EFSA report. 2007Google Scholar
- Murray MD: Influences of vector biology on transmission of arboviruses and outbreaks of disease: the Culicoides brevitarsis model. Vet Microbiol. 1995, 46: 91-99. 10.1016/0378-1135(95)00075-L.PubMedView ArticleGoogle Scholar
- Venter GJ, Paweska JT, Van Dijk AA, Mellor PS, Tabachnick WJ: Vector competence of Culicoides bolitinos and C. imicola for South African bluetongue virus serotypes 1, 3 and 4. Med Vet Entomol. 1998, 12: 378-385. 10.1046/j.1365-2915.1998.00116.x.PubMedView ArticleGoogle Scholar
- Breard E, Hamblin C, Hammoumi S, Sailleau C, Dauphin G, Zientara S: The epidemiology and diagnosis of bluetongue with particular reference to Corsica. Res Vet Sci. 2004, 77: 1-8. 10.1016/j.rvsc.2003.08.004.PubMedView ArticleGoogle Scholar
- Venter GJ, Nevill EM, Linde Van der TC: Geographical distribution and relative abundance of stock-associated Culicoides species (Diptera: Ceratopogonidae) in southern Africa in relation to their potential as viral vectors. Onderstepoort J Vet Res. 1996, 63: 25-38.PubMedGoogle Scholar
- Purse BV, Tatem AJ, Caracappa S, Rogers DJ, Mellor PS, Baylis M, Torina A: Modelling the distributions of Culicoides bluetongue virus vectors in Sicily in relation to satellite-derived climate variables. Med Vet Entomol. 2004, 18: 90-101. 10.1111/j.0269-283X.2004.00492.x.PubMedView ArticleGoogle Scholar
- Conte A, Colangeli P, Ippoliti C, Paladini C, Ambrosini M, Savini L, Dall'Acqua F, Calistri P: The use of a web-based interactive geographical information system for the surveillance of bluetongue in Italy. Rev Sci Tech. 2005, 24: 857-868.PubMedGoogle Scholar
- Guis H, Tran A, de La Rocque S, Baldet T, Gerbier G, Barragué B, Biteau-Coroller F, Roger F, Viel JF, Mauny F: Use of high spatial resolution satellite imagery to characterize landscapes at risk for bluetongue. Vet Res. 2007, 38 (5): 669-683. 10.1051/vetres:2007025.PubMedView ArticleGoogle Scholar
- Martin PA, Cameron AR, Greiner M: Demonstrating freedom from disease using multiple complex data sources 1: a new methodology based on scenario trees. Prev Vet Med. 2007, 79: 71-97. 10.1016/j.prevetmed.2006.09.008.PubMedView ArticleGoogle Scholar
- Carpenter S, Lunt HL, Arav D, Venter GJ, Mellor PS: Oral susceptibility to bluetongue virus of Culicoides (Diptera: Ceratopogonidae) from the United Kingdom. J Med Entomol. 2006, 43: 73-78. 10.1603/0022-2585(2006)043[0073:OSTBVO]2.0.CO;2.PubMedView ArticleGoogle Scholar
- Mellor PS, Boorman J, Baylis M: Culicoides biting midges: their role as arbovirus vectors. Annu Rev Entomol. 2000, 45: 307-340. 10.1146/annurev.ento.45.1.307.PubMedView ArticleGoogle Scholar
- Osmani A, Murati B, Kabashi Q, Goga I, Berisha B, Wilsmore AJ, Hamblin C: Evidence for the presence of bluetongue virus in Kosovo between 2001 and 2004. Vet Rec. 2006, 158: 393-396.PubMedView ArticleGoogle Scholar
- Baylis M, Touti J, Bouayoune H, Moudni L, Taoufiq B, el Hasnaoui H: Studies of the mortality rate of Culicoides imicola in Morocco. Arch Virol Suppl. 1998, 14: 127-136.PubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.