- Open Access
Potential vectors of bluetongue virus in high altitude areas of Yunnan Province, China
Parasites & Vectors volume 12, Article number: 464 (2019)
Bluetongue disease of ruminants is a typical insect-borne disease caused by bluetongue virus (BTV) of the genus Orbivirus (family Reoviridae) and transmitted by some species of Culicoides (Diptera: Ceratopogonidae). Recently, the detection of BTV in yaks in high altitude meadows of the Shangri-La district of Yunnan Province, China, prompted an investigation of the Culicoides fauna as potential vectors of BTV.
A total of 806 Culicoides midges were collected by light trapping at three sites at altitudes ranging from 1800 to 3300 m. The species were identified based on morphology and the DNA sequences of cytochrome c oxidase subunit 1 (cox1). PCR and quantitative PCR following reverse transcription were used to test for the presence of BTV RNA in Culicoides spp. A phylogenetic analysis was used to analyze the cox1 sequences of some specimens.
Four species dominated these collections and cox1 barcoding revealed that at least two of these appear to belong to species new to science. Culicoides tainanus and a cryptic species morphologically similar to C. tainanus dominated low altitude valley collections while C. nielamensis was the most abundant species in the high-altitude meadow. A species related to C. obsoletus occurred at all altitudes but did not dominate any of the collections. BTV RT-qPCR analysis detected BTV RNA in two specimens of C. tainanus, in one specimen closely related to C. tainanus and in one specimen closely related to C. obsoletus by barcode sequencing.
This study suggests that BTV in high altitude areas of Yunnan is being transmitted by three species of Culicoides, two of which appear to be new to science. This research may be useful in improving understanding of the effects of global warming on arboviral disease epidemiology and further study is important in research into disease control and prevention.
Bluetongue disease of ruminants is caused by bluetongue virus (BTV) of the genus Orbivirus (family Reoviridae), which is a notifiable disease to the World Organization for Animal Health (OIE, Office International des Epizooties). The virus is transmitted by species of Culicoides . There are 1368 extant species of Culicoides  but only about 30 of these have been associated with transmission of BTV [1, 3,4,5,6]. Bluetongue virus has mainly been reported in tropical and subtropical areas of Africa, Australia, North and South America, Asia and the Middle East with seasonal outbreaks reported from Mediterranean countries and Europe . The unprecedented northward spread of BTV into Europe in the early 2000s has been linked to the effects of climate change and its influence on both the distribution and vector capacity of local Culicoides populations [3, 8].
BTV has been isolated from many parts of China since 1979 [9,10,11] with at least 14 BTV serotypes identified by serological investigation and virus isolation [10, 12]. Some species of Culicoides have been considered as vectors of BTV, including C. oxystoma Kieffer (as C. schultzei Enderlein), C. peregrinus Kieffer, C. arakawai (Arakawa) and C. circumscriptus Kieffer ; however, solid evidence to support these assertions is lacking. A number of proven and suspected vector species are present in China including C. actoni Smith, C. brevitarsis Kieffer, C. fulvus Sen & Das Gupta, C. obsoletus (Meigen) and C. wadai Kitaoka, although some of these species are actually a complex of species and therefore require confirmation . Recently, Ma et al.  reported the presence of BTV in yaks and sheep on the Tibetan Plateau, western China, at an altitude of over 3000 m but little is known of the species of midges living at such high altitudes or their potential to transmit BTV.
Yunnan Province is located in the southwest of China and shares a border with Laos, Myanmar and Vietnam. While the southern tropical and sub tropical districts of Yunnan are endemic for BTV , the north of Yunnan is mostly cold, high altitude mountains and river valleys, and there is little information on the prevalence of BTV or other arboviruses in this area.
Shangri-La district is located in the northwest of Yunnan, next to the Tibetan Plateau, with an area of 11,613 km2. The Shangri-La district consists of relatively warm river valleys and cold highland meadows. The river valleys are located about 1900 m above sea level along the Changjiang River and contain many small tributary streams. The highland meadows average about 3450 m above sea level and are covered by snow for almost four months of the year. There is, however, some evidence that this extreme climate is becoming milder as historical weather records indicate that the average monthly minimum temperatures between 2011 and 2018 are about 1 °C higher than that of the years 1971 to 2000 (Fig. 1a). Furthermore, the monthly minimum temperatures between 2011 and 2018 were milder than between 1971 and 2000 (Fig. 1b) . The impact of such changes is difficult to predict but could result in increased activity of Culicoides spp. and associated arboviruses in this region.
During an investigation between 2014 and 2017, an average of 5% of yaks living in a 3290 m altitude highland meadow of Shangri-La were found to have BTV antibodies . In a nearby river valley at a lower altitude, a 35–65% seroprevalence of BTV and a BTV-21 strain was isolated from some goat farms  indicating that BTV is prevalent at these altitudes. There are scarce data on the species diversity and abundance of Culicoides spp. at high altitudes in Asia so it is unclear if any local species can survive the cold and long winters experienced in these areas. Some midge species are, however, able to survive at high altitudes in South America  and Europe  so it is likely that at least some Asian species are also able to survive in the Shangri-La area. Also unclear is whether the BTV found in yaks in the highland meadows is endemic or introduced seasonally by the immigration of infected midges from the warmer river valleys.
This paper provides a brief account of the diversity and abundance of Culicoides spp. on livestock farms in the Shangri-La district of Yunnan Province, China and the screening of the most common species for the presence of BTV RNA using real-time quantitative PCR.
Culicoides spp. collection
Culicoides were collected from farms at three different sites. The first site was on a farm (27°03′36″N, 100°04′12″E) with 4 cattle and 24 goats near Jinjiang village along the Changjiang River, at an altitude of 1800 m. The second site was on a farm (28°06′00″N, 99°25′12″E) with 2 cattle and 85 goats near Nixi village, located in a small tributary river valley named Longchi at an altitude of 2200 m. The third site was on a farm (27°30′36″N, 99°48′36″E) with 15 yaks and 41 goats near Xiaozhongdian village, in the middle of a highland meadow area at an altitude of 3300 m. Mains-powered UV light traps (LTS-M02, Wuhan Lucky Star Medical Treatment Technology Co., Wuhan, China) were used for trapping Culicoides. Two light traps were placed close (10 m) together at each site for one night and ran from 16:00 h to 9:00 h the following day. The insects were collected directly into 70% ethanol and stored at 4 °C. All collections were made in the period 8th to 15th August 2017.
Insect specimens were sorted into morphospecies based on gross morphology and wing pattern using keys and illustrations of Yu et al. . Parous and gravid females free from visible blood were selected based on abdominal pigmentation  and contents and submitted individually for nucleic acid extraction. Following digestion, representative specimens were mounted onto microscope slides using the technique of Bellis et al.  and identified morphologically using the keys and descriptions of Yu et al. .
DNA and RNA extraction
Whole midges were digested non-destructively following the methods of Bellis et al. . Briefly, insects were immersed in 300 µl of DXT tissue digestion buffer (Qiagen, Hilden, Germany) incorporating 1% DX digestion buffer additive (Qiagen) and incubated at 40 °C overnight. Following digestion, the insect cadaver was removed from the digestion lysate and preserved in 70% ethanol for subsequent morphological examination as described above. From the lysate, 220 µl was used for DNA extraction and the remaining 50 µl was used for RNA extraction. The RNA was extracted using a MagMAX™-96 viral RNA Isolation kit (Am1836; Ambion, Austin, TX, USA) following the manufacturer’s directions and using a MagMAX™ Express-96 machine (Ambion). The DNA was extracted using a DNeasy Blood & Tissue kit (#69560; Qiagen) following the manufacturer’s directions. The final DNA elution of 120 µl and RNA of 50 µl were stored at − 20 °C.
Pan BTV serotype RT-qPCR
The primers BTVF-MH, BTVR-MH and probe BTVP-MH for real-time quantitative polymerase chain reaction (qPCR) described by Hofmann et al.  were used to detect the NS3 gene of BTV in insect specimens. Reverse transcription quantitative PCR (RT-qPCR) was performed on a Fast7500 Realtime PCR machine (Applied Biosystems, Carlsbad, CA, USA) using an AgPath-ID™ one step RT-PCR Kit (Ambion). The reaction was performed as described by the manufacturer on a 2 µl RNA sample in a total volume of 20 µl, using PCR master mix with 6 mM MgCl2, 1 μM of each PCR primer and 0.2 μM of each probe. The RT-qPCR program consisted of the following: 10 min reverse transcription at 45 °C; 10 min denaturation at 95 °C; 45 cycles of 15 s denaturation at 95 °C, 45 s annealing and extension at 65 °C. Fluorescence was measured at the end of each annealing step. Only reactions with Cq values of 25 or lower were regarded as indicative of replication of BTV in the midge [4, 24, 25].
cox1 gene amplification
Two specimens morphologically referable to C. nielamensis Liu & Deng, 19 specimens morphologically referable to C. tainanus Kieffer and 9 specimens morphologically referable to C. obsoletus (Meigen) including the 4 specimens found to be infected with BTV, were processed and sequenced for cytochrome c oxidase subunit 1 (cox1). Primers BC1culicFm and JerR2m, and techniques described by Bellis et al. , were used to amplify the cox1 barcode region. Resultant sequences were compared with published data on Barcode of Life Data System (BOLD) and the National Center for Biotechnology Information (NCBI).
The PCR products of cox1 (692 bp) genes from Culicoides spp. were sent to a local company for sequencing using the Sanger method on an ABI3739XL machine (Kunming Shuoqing Biological Technology Company, Kunming, China).
cox1 sequences of representative Culicoides specimens and the most similar Culicoides spp. from NCBI data were analyzed by alignment using MEGA-X. Pairwise Kimura 2-parameter (K2P) distances among sequences were calculated and compared as a neighbor-joining (NJ) tree (bootstrap = 1000) using MEGA-X .
Weather conditions at the three sites were typical for the season. The temperatures in the afternoon at the collection sites were 28 °C (Jinjiang), 25 °C (Nixi) and 20 °C (Xiaozhongdian). The total number of midges collected at Jinjiang, Xiaozhongdian and Nixi was 591, 108 and 107, respectively. Specimens with morphology consistent with C. tainanus, C. nielamensis and C. obsoletus dominated collections, comprising 98% of all specimens collected (Fig. 2 and Table 1). Other species collected included a species similar to C. punctatus Latreille, another close to C. sinanoensis Tokunaga, and five species which remain unidentified due to lack of description in reference books (Table 1).
cox1 sequences of 8 specimens morphologically identified as C. obsoletus (BOLD accession numbers SGRL234-SGRL240 and SGRL242) and 16 specimens morphologically identified as C. tainanus (BOLD accession numbers SGRL089-SGRL094, SGRL096-SGRL098 and SGRL101-SGRL107) were registered on BOLD. A neighbor-joining tree of these specimens (Fig. 3) revealed the presence of two clades of C. tainanus. One of these clades was 99.47% similar to C. tainanus (= C. maculatus Shiraki) sequences on BOLD while the other was only 90.12% similar to those specimens (Table 2). So far, we have been unable to distinguish specimens from these two clades morphologically, so we are tentatively calling this second species Culicoides sp. close to C. tainanus.
Similarly, we found no close matches for our specimens resembling C. obsoletus (Figs. 3, 4). The closest match available on NCBI or BOLD was for C. obsoletus from Canada, suggesting that our specimens are close to, but different from, true C. obsoletus and possibly represent a cryptic species.
No reference sequences in BOLD nor NCBI were close to our cox1 sequences from C. nielamensis; the nearest sequence was 84.14% similar and belonged to C. minutissimus Zett (Table 2).
The relative abundance of these four dominant species varied with altitude, C. tainanus and Culicoides sp. close to C. tainanus were only found at low and mid-altitudes, C. nielamensis was most abundant at high altitude and Culicoides sp. close to C. obsoletus was found at all altitudes. From high to low altitude, the relative abundance of C. nielamensis was 97.22, 49.53 and 3.55%, C. tainanus and Culicoides sp. close to C. tainanus was 0, 35.51 and 91.70% and Culicoides sp. close to C. obsoletus was 2.78, 8.4% and 3.89% (Fig. 2 and Table 1).
BTV RT-qPCR of Culicoides
A total of 150 specimens of C. tainanus and Culicoides sp. close to C. tainanus, 35 Culicoides sp. close to C. obsoletus and 80 C. nielamensis were analyzed for the presence of BTV. Many specimens of all species had detectable Cq values, including 70 C. tainanus and Culicoides sp. close to C. tainanus, 32 Culicoides sp. close to C. obsoletus and 42 C. nielamensis. Cq values exhibited a bimodal distribution with 4 specimens having values of 25 or less and the remainder of specimens with values between 30 and 35 (Fig. 5).
Cq values below 25 have been shown in both colony  and field-collected [4, 25] specimens to indicate that viral replication has occurred in the insect so the 4 specimens with Cq values less than 25 were considered infected with BTV. These 4 specimens were all collected from Nixi and were all parous females; 2 were identified as C. tainanus, 1 was Culicoides sp. close to C. tainanus and 1 was Culicoides sp. close to C. obsoletus (Fig. 3 and Table 3).
To our knowledge, this is the first investigation of BTV vectors at high altitude areas of Yunnan and follows the detection of an average 5% BTV seroprevalence in yaks at 3290 m altitude of Jiantang and 35–65% BTV seroprevalence in goats at 2000–2800 m altitude of Nixi in Shangri-La . Based on our limited data, the abundance of midges in this area is low compared to tropical or sub-tropical areas of Yunnan where over 10,000 midges can be collected after 12 hours trapping (our unpublished data).
The low prevalence of BTV observed in yaks at 3300 m could be due to a variety of factors, many of which are related to the impact of the harsh environment on vector biology, for example low vector abundance, low vector competence or short season of vector activity. The data from this brief study are not sufficient to allow firm conclusions about the absolute abundance of vector populations at different altitudes but some differences in midge diversity were evident between the three sites studied. The change in dominance from C. tainanus and Culicoides sp. close to C. tainanus at low altitude to C. nielamensis at high altitude is evident and coupled with the detection of BTV in C. tainanus and Culicoides sp. close to C. tainanus but not in C. nielamensis may be a contributing factor to the lowered BTV prevalence at high altitude.
Culicoides sp. close to C. obsoletus was present only in low numbers at all sites which may reduce its potential to act as a significant vector of BTV. Specimens from Shangri-La are morphologically similar to the European C. obsoletus, but significantly different in cox1 (about 9%). The Obsoletus group of Culicoides subgenus Avaritia has been shown to contain several species with very similar morphology and it is likely the species present in Shangri-La is a further cryptic species within this complex. Species from this complex are Holarctic so it is not surprising that populations are able to survive in high altitude sites such as Shangri-La. Further study on this group in Asia may reveal that the species present in Shangri-La is more widely distributed than currently known. Some species of the Obsoletus group have also been implicated in the transmission of BTV [6, 27, 28], so it is not surprising that a further member of this group may also be acting as a vector of BTV.
Recent work in Japan  found evidence that C. tainanus (= C. maculatus Shiraki) may be acting as a vector of BTV and this is supported by the detection of BTV in specimens from Shangri-La. The presence of a morphologically similar species in Asia that may also be a vector of BTV highlights the need for ongoing studies of the taxonomy of species in this region, particularly species belonging to the subgenus Avaritia of Culicoides. Integrative taxonomy using a combination of morphological and molecular analyses is currently the best method of establishing the identity of specimens  and should be applied to any studies of the vector status of species.
Specimens of Culicoides sp. close to C. obsoletus, C. tainanus and Culicoides sp. close to C. tainanus were found to contain BTV RNA in concentrations consistent to indicate virus had replicated in the insect thus fulfilling one of the criteria of these species being proven vectors of BTV [30, 31]. Consistent with the findings of previous studies [4, 24, 25], many specimens of midges from Shangri-La exhibited a positive response in the RT-PCR but the high Cq values from these specimens suggest either a non-specific response or a response due to remnant virus in the midge following a viraemic blood meal. The absence of specimens with Cq values between 25 and 30 supports the conclusion of Veronesi et al.  and van der Saag et al. [4, 25] that the virus has indeed multiplied in those midges with Cq values of less than 25.
None of the 80 specimens of C. nielamensis tested showed evidence of infection with BTV suggesting a prevalence of less than 1 in 80 (< 1.25%) compared to a prevalence of > 2.0% for C. obsoletus, C. tainanus and Culicoides sp. close to C. tainanus. More specimens would need to be tested to discount the possibility of C. nielamensis being capable of transmitting BTV. Few studies have reported on the prevalence of BTV in field populations of vectors [32,33,34], so it is difficult to establish when a species is not acting as a vector of BTV. However, the lower prevalence of BTV at high altitude, despite the higher proportion of C. nielamensis at these altitudes, suggests that C. nielamensis may not be acting as a vector, or at least not an efficient vector.
This study has shown that biting midge species capable of carrying BTV are present at 3300 m altitude sites in Shangri-La, which may explain the occurrence of BTV in yaks at that altitude. It is unknown if BTV moves seasonally to these highland meadows from nearby low altitude areas but as the seasonal movement of midges and viruses has been documented elsewhere , it is possible that this is also occurring in Yunnan or on the neighboring Tibetan Plateau. Based on the available evidence, any such movement could possibly be due to the dispersal of Culicoides sp. close to C. obsoletus which occurs at all altitudes studied here and has been shown to be infected with BTV; however, further work is needed to confirm that C. nielamensis or some other species is not also acting as a vector capable of moving virus between mid and high altitude sites.
Yaks are an important economic source for the Zang people who live in and around the Tibetan Plateau, and a recent study  indicated that the abortion rate of yaks in Qinghai Province, China, of 21.39% could be caused by BTV or other pathogens. Further study is needed to confirm the ability of BTV to cause clinical disease in yaks and subsequently to assist with disease control and prevention. Knowledge of the ecology of BTV in the extreme climate of Shangri-La may improve understanding of the effects of global warming on vector biology and disease epidemiology, particularly given the evidence that temperatures in the region are slowly rising.
This study has identified three species of Culicoides, two of which appear to be new to science, which may be acting as vectors of BTV in high altitude areas of Yunnan. Further studies including laboratory infection and transmission of BTV is required to confirm the vector status of these species. The seroprevalence of BTV at high altitudes of Yunnan and the neighboring Tibetan Plateau could be caused by the seasonal movement of BTV vectors, such as Culicoides sp. close to C. obsoletus, from lower altitude areas where BTV is endemic.
Availability of data and materials
Data supporting the conclusions of this article are provided within the article. Raw data are available from the corresponding author upon request. The newly generated sequences were submitted to the GenBank database under the accession numbers MN513299–MN513328.
basic local alignment search tool
Barcode of Life Data system
cytochrome c oxidase subunit 1
pairwise Kimura 2-parameter
National Center for Biotechnology Information
Office International des Epizooties
polymerase chain reaction
reverse transcription PCR
reverse transcription quantitative PCR
Meiswinkel R, Gomulski LM, Delecolle JC, Goffredo M, Gasperi G. The taxonomy of Culicoides vector complexes—unfinished business. Vet Ital. 2004;40:151–9.
Borkent A. World species of biting midges (Diptera: Ceratopogonidae). 2016. http://www.inhs.illinois.edu/research/FLYTREE/Borkent.html. Accessed 10 Jul 2018.
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–81.
Van der Saag M, Nicholas A, Ward M, Kirkland P. Evaluation of in vitro methods for assessment of infection of Australian Culicoides spp. with bluetongue viruses. Vet Ital. 2015;51:335–43.
Wilson A, Mellor P. Bluetongue in Europe: vectors, epidemiology and climate change. Parasitol Res. 2008;103(Suppl. 1):S69–77.
Caracappa S, Torina A, Guercio A, Vitale F, Calabro A, Purpari G, et al. Identification of a novel bluetongue virus vector species of Culicoides in Sicily. Vet Rec. 2003;153:71–4.
Purse BV, Carpenter S, Venter GJ, Bellis G, Mullens BA. Bionomics of temperate and tropical Culicoides midges: knowledge gaps and consequences for transmission of Culicoides-borne viruses. Annu Rev Entomol. 2015;60:373–92.
Conte A, Giovannini A, Savini L, Goffredo M, Calistri P, Meiswinkel R. The effect of climate on the presence of Culicoides imicola in Italy. J Vet Med B Infect Dis Vet Public Health. 2003;50:139–47.
Zhang N, Li Z, Zhang K. Bluetongue history, serology and virus isolation in China. In: St George TD, Kegao P, editors. Bluetongue disease in Southeast Asia and the Pacific. Canberra: Australian Centre for International Agricultural Research; 1996. p. 43–55.
Yang H, Xiao L, Wang J, Meng J, Lv M, Liao D, et al. Phylogenetic characterization genome segment 2 of bluetongue virus strains belonging to serotypes 5, 7 and 24 isolated for the first time in China during 2012 to 2014. Transbound Emerg Dis. 2016;64:1317–21.
Qin S, Yang H, Zhang Y, Li Z, Lin J, Gao L, et al. Full genome sequence of the first bluetongue virus serotype 21 (BTV-21) isolated from China: evidence for genetic reassortment between BTV-21 and bluetongue virus serotype 16 (BTV-16). Arch Virol. 2018;163:1379–82.
Kirkland PD, Zhang N, Hawkes RA, Li Z, Zhang F, Davis RJ, et al. Studies on the epidemiology of bluetongue virus in China. Epidemiol Infect. 2002;128:257–63.
Bi YL, Li CH, Li SY, Qing B, Zhong N, Hu LF, et al. An epidemiological survey of bluetongue in Yunnan Province, China. In: St George TD, Kegao P, editors. Bluetongue disease in Southeast Asia and the Pacific. Canberra: Australian Centre for International Agricultural Research; 1996. p. 51–6.
Gopurenko D, Bellis GA, Yanase T, Wardhana AH, Thepparat A, Wang J, et al. Integrative taxonomy to investigate species boundaries within Culicoides (Diptera: Ceratopogonidae): a case study using subgenus Avaritia from Australasia and Eastern Asia. Vet Ital. 2015;51:345–78.
Ma JG, Zhang XX, Zheng WB, Xu YT, Zhu XQ, Hu GX, et al. Seroprevalence and risk factors of bluetongue virus infection in Tibetan sheep and yaks in Tibetan Plateau, China. Biomed Res Int. 2017;2017:5139703.
Historic weather of Zhongdian (Shangri-La). http://php.weather.sina.com.cn/whm.php?city=%D6%D0%B5%E9&dpc=1. Accessed 11 Jan 2019.
Duan YL, Miao HS, Liao DF, Kou ML, Li ZH, Wang Z, et al. The serologic investigation and viral isolation of bluetongue virus in Shangri-La in southwest China. Transbound Emerg Dis. 2019. https://doi.org/10.1111/tbed.13292.
Wirth WW, Lee DJ. New species of Culicoides from high altitudes in the Colombian Andes (Diptera: Ceratopogonidae). Proc US Natl Mus. 1967;124:1–22.
Kaufmann C, Steinmann IC, Hegglin D, Schaffner F, Mathis A. Spatio-temporal occurrence of Culicoides biting midges in the climatic regions of Switzerland, along with large scale species identification by MALDI-TOF mass spectrometry. Parasit Vectors. 2012;5:246.
Liu JH, Liu GP, Liu ZJ, Yan G, Hao BS, Zhao TS, et al. Ceratopogonidae: Culicoides. In: Yu YX, editor. Ceratopogonidae of China, vol. 2. 1st ed. Beijing: Military Medical Science Press; 2005. p. 816–1323.
Dyce AL. The recognition of nulliparous and parous Culicoides (Diptera: Ceratopogonidae) without dissection. Aust J Entomol. 1969;8:11–5.
Bellis GA, Dyce AL, Gopurenko D, Mitchell A. Revision of the Immaculatus group of Culicoides Latreille (Diptera: Ceratopogonidae) from the Australasian Region with descriptions of two new species. Zootaxa. 2013;3680:15–37.
Hofmann M, Griot C, Chaignat V, Perler L, Thur B. Bluetongue disease reaches Switzerland. Schweiz Arch Tierheilkd. 2008;150:49–56.
Veronesi E, Antony F, Gubbins S, Golding N, Blackwell A, Mertens PP, et al. Measurement of the infection and dissemination of bluetongue virus in Culicoides biting midges using a semi-quantitative rt-PCR assay and isolation of infectious virus. PLoS ONE. 2013;8:e70800.
Van Der Saag MR, Ward MP, Kirkland PD. Application of an embryonated chicken egg model to assess the vector competence of Australian Culicoides midges for bluetongue viruses. Med Vet Entomol. 2017;31:263–71.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: Molecular Evolutionary Genetics Analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:2731–9.
Carpenter S, McArthur C, Selby R, Ward R, Nolan DV, Luntz AJ, et al. Experimental infection studies of UK Culicoides species midges with bluetongue virus serotypes 8 and 9. Vet Rec. 2008;163:589–92.
De Liberato C, Scavia G, Lorenzetti R, Scaramozzino P, Amaddeo D, Cardeti G, et al. Identification of Culicoides obsoletus (Diptera: Ceratopogonidae) as a vector of bluetongue virus in central Italy. Vet Rec. 2005;156:301–4.
Kato T, Shirafuji H, Tanaka S, Sato M, Yamakawa M, Tsuda T, et al. Bovine arboviruses in Culicoides biting midges and sentinel cattle in Southern Japan from 2003 to 2013. Transbound Emerg Dis. 2016;63:e160–72.
WHO. Arboviruses and human disease. World Health Organisation technical report series 369. Geneva: World Health Organisation; 1967. p. 460–2.
WHO. Arthropod-borne viruses. In: Anonymous, editors. World Health Organisation technical report series 219. Geneva: World Health Organisation; 1961. p. 68.
Bellis GA, Melville L, Walsh S. Preliminary assessment of the vector capacity of Culicoides nudipalpis and C. orientalis for bluetongue virus. Arbovirus Res Aust. 2005;9:28–32.
Melville LF, Hunt NT, Daniels P. Application of the polymerase chain reaction (PCR) test with insects in studying bluetongue virus activity. In: St George TD, Kegao P, editors. Bluetongue disease in Southeast Asia and the Pacific. Canberra: Australian Centre for International Agricultural Research; 1996. p. 141–5.
Foxi C, Delrio G, Falchi G, Marche MG, Satta G, Ruiu L. Role of different Culicoides vectors (Diptera: Ceratopogonidae) in bluetongue virus transmission and overwintering in Sardinia (Italy). Parasit Vectors. 2016;9:440.
Bishop AL, Kirkland PD, Mckenzie HJ, Spohr LJ, Barchia IM, Muller MJ. Distribution and seasonal movements of Culicoides brevitarsis Kieffer (Diptera: Ceratopogonidae) at the Southern limits of its distribution in New South Wales and their correlation with arboviruses affecting livestock. Aust J Entomol. 1995;34:289–98.
Hu GW, Shen YL, Zhao QB, Chao YL, Kan W, Zhou JZ, et al. Preliminary investigation on abortion of yaks in part of Qinghai Province. Chinese J Vet Med. 2016;52:3–5.
We are thankful to Peter Kirkland (Elizabeth Macarthur Agriculture Institute, Australia) for advice on BTV PCR and to David Gopurenko (Wagga Wagga Agricultural Institute, Australia) for advice on cox1 analysis.
The present study was funded by The Key Project of Yunnan Science and Technique Department (2016FA005).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.