Skip to main content

Mosquitoes established in Lhasa city, Tibet, China

Abstract

Background

In 2009, residents of Lhasa city, Tibet Autonomous Region (TAR), China reported large numbers of mosquitoes and bites from these insects. It is unclear whether this was a new phenomenon, which species were involved, and whether these mosquitoes had established themselves in the local circumstances.

Methods

The present study was undertaken in six urban sites of Chengguan district Lhasa city, Tibet. Adult mosquitoes were collected by bed net trap, labor hour method and light trap in August 2009 and August 2012. The trapped adult mosquitoes were initially counted and identified according to morphological criteria, and a proportion of mosquitoes were examined more closely using a multiplex PCR assay.

Results

907 mosquitoes of the Culex pipiens complex were collected in this study. Among them, 595 were females and 312 were males. There was no significant difference in mosquito density monitored by bed net trap and labor hour method in 2009 and 2012. Of 105 mosquitoes identified by multiplex PCR, 36 were pure mosquitoes (34.29%) while 69 were hybrids (65.71%). The same subspecies of Culex pipiens complex were observed by bed net trap, labor hour method and light trap in 2009 and 2012.

Conclusion

The local Culex pipiens complex comprises the subspecies Cx. pipiens pipiens, Cx. pipiens pallens, Cx. pipiens quinquefasciatus and its hybrids. Mosquitoes in the Cx. pipiens complex, known to be, potentially, vectors of periodic filariasis and encephalitis, are now present from one season to the next, and appear to be established in Lhasa City, TAR.

Background

Once established in high altitude regions, some mosquito species may threaten the health of humans and vertebrates due to their ability to transmit numerous diseases [15]. With a permanent resident population of 559, 423 at the 2010 Chinese Census, Lhasa city, which is the administrative capital of the Tibet Autonomous Region (TAR), China is situated on the northern bank of the Lhasa River, a tributary of the Yarlung Zangbo, in the mid-south of TAR. To the east and southeast of Lhasa are the regions of Nyingchi and Sharman; Nagqu neighbours Lhasa on the north and west; Xigaze lies on its southwest. Among a total of 29,518 square kilometers, the urban area of Lhasa is just 50 square kilometers. Standing on a plain over 3,650 meters (13,000 feet) above sea level and surrounded by towering mountains, Lhasa is known as the “city of the sun.” With an annual average temperature of 7.5°C, its average temperature in January is 2.3°C and 15.4°C in July. The climate here is of the temperate plateau monsoon type. Lhasa has an annual precipitation of 426 millimeters with rain falling mainly in July, August and September.

The subspecies in the Culex pipiens complex, may transmit a range of pathogens including the West Nile [69] and St. Louis encephalitis viruses [10], avian malaria, and filarial worms [1113]. Members of the Cx. pipiens complex includes Cx. pipiens, Cx. quinquefasciatus, Cx. australicus, and Cx. globocoxitus[11]. Cx. pipiens has two recognized subspecies: Cx. pipiens pipiens and Cx. pipiens pallens. Cx. pipiens pipiens is an Old World taxa originally distributed from Northern Europe to the highlands of South Africa. Cx. pipiens pallens is found east of the Urals across temperate Asia [14]. Cx. quinquefasciatus typically thrives in tropical and sub-tropical regions, including the African lowlands, Americas, Asia, and Australia [15]. As to ecological preference, Cx. pipiens (Cx. pipiens pipiens & Cx. pipiens pallens) and Cx. pipiens quinquefasciatus[15] are found in most inhabited areas globally and are often closely associated with humans, earning them the names of northern and southern house mosquitoes, respectively [11]. In addition, Cx. p. pipiens also has two recognized forms “pipiens” and “molestus”, which differ dramatically in ecology. Though several members of the complex have limited geographic distributions [12]. Cx. pipiens (Cx. pipiens pipiens & Cx. pipiens pallens) and Cx. quinquefasciatus can hybridize extensively when their ranges overlap. Extensive introgression exists between populations of Cx. pipiens and Cx. quinquefasciatus in North America, Argentina, Madagascar, Japan and the Republic of South Korea [1619]. Cx. australicus and Cx. globocoxitus are restricted to Australia [12].

In China, Cx. pipiens complex consists of four subspecies [20], including Cx. pipiens pipiens, Cx. pipiens quinquefasciatus, Cx. pipiens pallens and Cx. pipiens molestus. Cx. pipiens quinquefasciatus cannot be considered as a separate species and Cx. pipiens pallens is not an intermediate form. Cx. pipiens molestus is present in the underground water system in Beijing and Shenyang, China [21]. The usual altitude of Cx. pipiens pipiens is lower than 3,000 m. In eastern Colorado, Cx. pipiens pipiens activity occurs primarily in the populated valleys at lower elevation, diminishing rapidly at higher levels (>3,000 m) [22]. Alvaro Diaz-Badillo et al. reported that Cx. pipiens pipiens, Cx. pipiens quinquefasciatus, and their hybrids were all present in Mexico City (2,200 m) [7]. In China, Cx. pipiens pipiens has been identified in Xinjiang Uygur Autonomous Region. Cx. pipiens quinquefasciatus occurred in areas south of 32°N. Xiaohong Sun et al. collected seventy-five Cx. pipiens quinquefasciatus in northeastern Yunnan Province (2,500-3,000 m) during 2005 and 2006 [23]. Cx. pipiens pallens distributed in areas north of the Yangtse River [24]. The highest elevation at which Cx. pipiens pallens has been observed in China is 2,900 metres, in Mainling County, Nyingchi area, Tibet [25].

Identifying members of the Cx. pipiens complex and other sibling species by morphologic methods is time-consuming and restricted to adult males [12, 26]. Other techniques, such as allozyme analyses [17], restriction fragment length polymorphism (RFLP) analysis of PCR products [27], only distinguish between the two major taxa of the complex: Cx. pipiens and Cx. quinquefasciatus. To solve this problem, Smith and Fonseca developed assays that use polymorphisms in the second intron of the acetylcholinesterase-2 (ace-2) locus to identify members of the Cx. pipiens complex and other sibling species. The same method may be used to detect introgression between Cx. pipiens and Cx. quinquefasciatus[12, 27, 28]. Extensive population level examination of most of the species shows they consistently generate unique fragments that may be easily resolved by electrophoresis on agarose gels. This method permits the rapid and reliable identification of local mosquitoes.

In recent years, there have been numerous changes that might assist mosquitoes to reach Lhasa and become established there. These include: global warming [2935], increasing international trade and tourism, population growth and mobility [36], transport improvements (such as completion of the Qinghai-Tibet Railway in 2006, the Qinghai-Tibet Highway, Sichuan-Tibet Highway and China-Nepal International Road, and the construction of the Gonggar Airport) [37], changing rainfall patterns [38], and developments in agriculture, urbanization and industrialization [39]. There are no official records to show whether mosquitoes existed in Lhasa city before 2009. In 2009, reports appeared in public media concerning the emergence of mosquitoes in Lhasa city. In addition, approximately 85.3 percent of local respondents said they were bitten by mosquitoes from the beginning of 2009 to the end of 2012, and almost one in 20 (4.5%) had to attend hospital for treatment for severe inflammation and local complications (Qiyong Liu et al., unpublished questionnaire survey in Lhasa in 2012). Therefore, this phenomenon is already perceived to be a serious public health problem. However, it is unclear which species of mosquitoes were involved, and whether these mosquitoes have indeed established themselves locally. This study was undertaken to test the media reports and to determine whether mosquitoes are now established in the city. The results provide the first scientific assessment of mosquitoes in Lhasa and provide a foundation for development of measures to control mosquito-borne diseases in Lhasa in the future.

Methods

Study sites

The present study was undertaken in six urban sites of Chengguan district during August 2009 and August 2012. The sites were selected to be broadly representative of the geographic conditions and socio-economic characteristics of urban Lhasa. They included Tibet Center for Disease Control and Prevention (Tibet CDC), Longwangtan Park, Tibet Post Hotel, Gamagongsang Community, Xiashasu Community and Jiacuo Community (Figure 1).

Figure 1
figure 1

The research sites of this study. These sites include Tibet Center for Disease Control and Prevention (Tibet CDC) (N29°39.682’, E91°07.433’), Longwangtan Park (N29°39.582′, E91°07.151′), Tibet Post Hotel (N29°39.164′, E91°07.398′), Gamagongsang Community (N29°29.307′, E91°08.872′), Xiashasu Community (N29°39.284′, E91°07.867′) and Jiacuo Community (N29°39.429′, E91°05.375′).

Tibet CDC lies to the northeast of the Potala Palace. The campus includes many family dormitory buildings and well-established trees (cypresses). Leaks from water pipes and the irrigation of lawns provide potential breeding sites for mosquitoes. Longwangtan Park lies to the northeast of the Potala Palace and features dense vegetation and a lake, with many fish and water birds. Tibet Post Hotel lies to the southeast of the Potala Palace, close to Longwangtan Park, it has many cypresses in its courtyard. Gamagongsang Community lies to the east of the Potala Palace with a population of 2,149 people within 837 households. There is limited infrastructure such as drainage systems and roads, and there are no parks or other urban green spaces. Xiashasu Community lies to the southeast of the Potala Palace with a population of 1,519 people within 907 households. Because of the famous “Dazhao Temple”, this community is the most crowded in urban Lhasa. Residents tend to be less educated, with lower incomes, restricted living spaces and poor dwelling conditions. Jiacuo Community lies to the northwest of the Potala Palace and adjacent to a large park. Residents are relatively wealthy and mainly live in self-built single family houses with small yards.

Mosquito collection and initial morphological identification

In this study, outdoor bed net traps (Tibet CDC, Longwangtan Park), light trap collection (Tibet Post Hotel, Tibet CDC, Gamagongsang Community, Xiashasu Community and Jiacuo Community) and labor hour method (Outpatient building and Residential area in Tibet CDC) were employed to collect adult mosquitoes in August 2009 and August 2012. Studies were carried out on the same days of the month, three years apart. The surveillance method was the standard method released by the General Administration of Quality Supervision, Inspection and Quarantine of the People’s Republic of China and Standardization Administration of the People’s Republic of China (surveillance methods for vector density-mosquito, GB/T 23797–2009).

The bed net traps were applied between 19:00 and 24:00 (the peak time for mosquitoes), taking account of the time for sunset in Lhasa city (generally 20:00 in August). Bed net traps were placed close to potential breeding habitats, at intervals of 100 m. The distance from the bed net traps to the nearest resident’s house was also about 100 m. The size of bed net traps was 1.5 m × 1.2 m × 1.5 m, with twenty-five centimeters between the floor and the bottom of the bed net traps. Some members of staff at Tibet CDC and China CDC were selected as human baits. These members (under double bed net traps to avoid mosquitoes bites) were also used repeatedly throughout the entire duration of the study [40]. Every hour, all mosquitoes inside the bed net traps were collected by an electrical aspirator for 15 minutes per hour throughout the 5 hour period. In 2009, the bed net traps were carried out on Aug.3rd - 4th in Tibet CDC (Lawn) (total of 3 bed net traps) and Aug.3rd in Longwangtan Park (total of 3 bed net traps). In 2012, the bed net traps were carried out on Aug. 7th - 8th in Tibet CDC (Lawn) (total of 4 bed net traps) and Aug. 10th - 11th in Longwangtan Park (total of 4 bed net traps).

An electric aspirator was employed for 15 minutes to collect mosquitoes inside an outpatient building and the residential area of Tibet CDC. In 2009, the labor hour method was carried out on Aug.4th in Tibet CDC (outpatient building) (total of 1 person) and Aug.3rd in Tibet CDC (residential area) (total of 1 person). In 2012, the labor hour method was carried out on Aug.8th in Tibet CDC (outpatient building) (total of 1 person) and Aug.7th in Tibet CDC (residential area) (total of 1 person).

Kung Fu Xiaoshuai miniature light traps (Photocatalytic Miewen Ying supply device; Wavelength: 2537Å; Power: 8W; Corporation: Wuhan Environmental Protection Technology Co., Ltd. Gemstar) were used to collect adult mosquitoes. The light traps were placed in the campus of Tibet Post Hotel, Tibet CDC and Gamagongsang Community, Xiashasu Community and Jiacuo Community. Traps were hung away from interference by light sources, 1.5 m above the floor. They were turned on 1 hour before sunset (20:00) and turned off 1 hour after sunrise (08:00). In 2012, the light traps were employed from Aug.5th - 12th in Tibet Post Hotel (total of 18 light traps), on Aug.7th - 9th in Tibet CDC (Lawn) (total of 9 light traps), on Aug.10th in Gamagongsang Community (total of 2 light traps), on Aug.12th in Xiashasu Community (total of 4 light traps) and on Aug.9th in Jiacuo Community (total of 4 light traps), respectively.

Information on temperature (°C) and relative humidity (%) was obtained from http://www.weather.com.cn. During collections, ambient outdoor air temperature and relative humidity was recorded hourly using a WS-1 Thermo-Hygrometer device.

Mosquito species identification

Each morning, the trapped adult mosquitoes were initially counted and identified according to morphological criteria using the key developed by Lu BL [24]. All collected mosquitoes were put into 1.5 ml centrifuge tubes individually and then transported to the laboratory of the Department of Vector Biology and Control in China CDC for further molecular identification. Genomic DNA was extracted from individual mosquitoes. A Qia Amp DNA Mini Kit (Qiagen Inc., CA) was adopted and DNA was extracted from the thorax of mosquitoes according to the manufacturer’s instructions.

To reveal the species composition of mosquitoes in Lhasa city, a multiplex PCR protocol was adopted using polymorphisms in the second intron of the acetylcholinesterase-2 (ace-2) locus, developed by Smith, J. L. & Fonseca, D. M [12]. Three forward primers (ACEquin, ACEpall and ACEpip) and one backward primer (B1246s) were adopted simultaneously. Each of the three primers was used in conjunction with the reverse primer B1246s [12, 26], (Table 1). Because of limited distribution of Cx. pipiens molestus in China [21], the primer of Cx. pipiens molestus was not included in this study. Approximately 105 (14.4%) mosquitoes that were selected from four sites (two institutions and two communities) in 2012, were further identified to sub-species level.

Table 1 Primer sequences for the acetylcholinesterase-2 locus-based polymerase chain reaction assay

The PCR assay was optimized for 25 ul volumes. Reactions contained 10 × PCR buffer, 250 uM of each dNTP, one unit of Taq polymerase, and genomic DNA. The amplification program consisted of one cycle at 94°C for five minutes, followed by 35 cycles at 94°C for 30 seconds, 55°C for 30 seconds, 72°C for one minute, and one cycle at 72°C for five minutes.

In addition, to further verify the subspecies of the Cx. pipiens complex, further sequence analysis of the Ace-2 gene for both some pure Cx. pipiens pipiens, Cx. pipiens quinquefasciatus, Cx. pipiens pallens and possible hybrids among them were conducted by Tsingke Company (Beijing, China) using the same mosquitoes as the multiplex PCR assay. Approximately three each of pure and possible hybrid mosquitoes were further sequenced in this study.

Statistical analysis

Information was recorded on the date of the collections, number of bed net traps, number of light traps, duration of mosquito catch (h), the presence and gender of Cx. pipiens complex mosquitoes. An independent-sample T test was adopted to compare the density of mosquitoes between 2009 and 2012 after a satisfactory check for normality of the distribution and homogeneity of variance of the data. Numbers of species identified by multiplex PCR at different sites in different years were recorded and calculated. Analysis was conducted by SPSS (Statistical Package for the Social Sciences) statistical software (version 17.0).

Ethics statement

We obtained ethical approval from the Ethical Review Committee of Chinese Center for Disease Control and Prevention for this study (No. 201214). Permission was also obtained from the Government, the Municipal Health Bureau and Tibet CDC in the Tibet Autonomous Region.

Results

Morphological identification

In this study, 907 mosquitoes in total were captured including 595 female and 312 male mosquitoes (Table 2). Preliminary morphological identification demonstrated that all these mosquitoes belonged to subspecies of the Cx. pipiens complex [24].

Table 2 The presence of the Culex pipiens complex in different collections during the mosquito season in Lhasa city, Tibet

Mosquitoes collected by bed net traps in different years

Using bed net traps, 132 mosquitoes (132 females) were collected in Tibet CDC (Lawn) and 5 mosquitoes (4 females, 1 male) were collected in Longwangan Park in 2009. Three years later, 25 mosquitoes (25 females) were collected in Tibet CDC (Lawn) and 24 mosquitoes (20 females, 4 males) were collected in Longwangan Park, (Table 2). The mean mosquito density was 23.72 (mosquitoes per hour per net) in 2009 while it was 3.06 (mosquitoes per hour per net) in 2012 (Table 3). There was no significant difference of mosquito density monitored by bed net traps in 2009 and 2012 (t = 1.299, df = 12, P = 0.218>0.05).

Table 3 The density of the Culex pipiens complex in different years during the mosquito season in Lhasa city, Tibet

Mosquitoes collected by labor hour method in different years

Using the labor hour method, 34 mosquitoes (34 females) were collected in Tibet CDC (Outpatient building) and 7 mosquitoes (4 females, 3 males) were collected in Tibet CDC (Residential area) in 2009. In 2012, 26 mosquitoes (17 females, 9 males) were collected in Tibet CDC (Residential area) (Table 2). The mean mosquito density was 62.10 (mosquitoes per hour per person) and 8.54 (mosquitoes per hour per person) in 2009 and 2012, respectively (Table 3). There was no significant difference of mosquito density monitored by labor hour method in 2009 and 2012 (t = 1.291, df = 2, P = 0.326>0.05).

Mosquitoes collected by light traps in 2012

In 2012, light traps collected 83 female and 58 male mosquitoes in Tibet Post Hotel and 12 females and 13 males in Tibet CDC (Lawn). 2 females were collected in Gamagongsang. 124 females and 53 males were collected in Xiashasu. 135 females and 171 males were collected in Jiacuo (Table 2). The mean mosquito density was 17.59 (mosquitoes per trap per night) in 2012 (Table 3).

Multiplex PCR assay for molecular identification

In this study, 105 mosquitoes from Lhasa city and another 17 mosquitoes used as positive controls from other provinces (positive control: Cx. pipiens pipiens from Urumchi, Xinjiang, 610 bp;Cx. pipiens quinquefasciatus from Dali, Yunnan, 274 bp; Cx. pipiens pallens from Beijing, 478 bp) were examined using the multiplex PCR assay developed by Smith, J. L. & Fonseca, D. M [12]. The results revealed that the size of the amplified product was 274 bp for Cx. pipiens quinquefasciatus, 478 bp for Cx. pipiens pallens and 610 bp for Cx. pipiens pipiens. Primers were successfully designed for the identification of these mosquitoes (Figure 2).

Figure 2
figure 2

Ethidium bromide-stained agarose gel showing multiplex PCR assay for molecular identification. M: Maker; Positive control products from Cx. pipiens complex are shown in lanes1-7. Lane 1 (Cx. pipiens pipiens from Urumchi, Xinjiang, 610 bp); Lane 2 (Cx. pipiens quinquefasciatus from Dali, Yunnan, 274 bp); Lane 3 (Cx. pipiens pallens from Beijing, 478 bp); Lane 4 (Cx. pipiens pipiens &Cx. pipiens pallens, 1 ul ofeach DNA, 610 bp &478 bp); Lane 5 (Cx. pipiens pipiens &Cx. pipiens quinquefasciatus, 1 ul of each DNA, 610 bp&274 bp); Lane 6 (Cx. pipiens quinquefasciatus&Cx. pipiens pallens, 1 ul of each DNA,274 bp&610 bp); Lane 7 (Cx. pipiens pipiens&Cx. pipiens pallens&Cx. pipiens quinquefasciatus, 0.67 ul of each DNA, 610 bp&274 bp&478 bp), respectively. Mosquitoes from Lhasa city and other provinces are shown in lanes 8–22, respectively. Lane 23, negative control. Outside lanes are 100 bp DNA ladders. bp = basepaires.

Multiplex PCR revealed that mosquitoes from Lhasa city include 36 pure mosquitoes (34.29%) and 69 hybrids (65.71%) (Table 4). The pure mosquitoes included 13 Cx. pipiens pipiens, 3 Cx. pipiens quinquefasciatus and 20 Cx. pipiens pallens. Possible hybrids consisted of 41 Cx. pipiens pipiens and Cx. pipiens pallens, 14 Cx. pipiens pallens and Cx. pipiens quinquefasciatus, only 1 Cx. pipiens pipiens and Cx. pipiens quinquefasciatus, and 13 among Cx. pipiens pallens & Cx. pipiens quinquefasciatus and Cx. pipiens pipiens.

Table 4 Multiplex PCR assay for molecular identification

Sequence analysis further confirmed the accuracy of multiplex PCR in this study. Based on the similarity analysis using blast, the similarity of the Cx. pipiens quinquefasciatus from Lhasa to the sequence of Cx. pipiens quinquefasciatus from central Bangladesh was 100% (273/273). The similarity of the Cx. pipiens pallens from Lhasa to the sequence of Cx. pipiens pallens from Iran was 99% (436/438). The similarity of the Cx. pipiens pipiens from Lhasa to the sequence of Cx. pipiens pipiens from Iran was 95% (530/560).

Discussion

This is the first investigation to verify media reports of mosquitoes in Lhasa. We observed subspecies of Cx. pipiens complex and its hybrids on two occasions, three years apart. Our findings were based on entomological investigations in the field and multiplex PCR methods in the laboratory. In this study, there was no significant difference of mosquito density monitored by bed net trap and labor hour method in 2009 and 2012. In urban Lhasa, we observed that the ecological and geographical factors did not change significantly three years later. However, it seemed that the mean mosquito density, both using bed net traps and labor hour method, were relatively higher in 2009 than in 2012 though no statistical significance was observed. We note that the summer of 2012 did not match the high temperatures of three years earlier: the maximum temperature in 2012 was 29.0°C compared with 30.4°C in 2009 (http://www.tianqi.com). In addition, the public health campaign may also play a major role in the relatively lower density in 2012. In recent years, Lhasa accelerated the process of establishing of the National Sanitary City in China. As a key indicator, mosquito density was controlled by local health authorities and related agencies using some insecticides and similar products. This campaign might have exerted some adverse impact on the density of mosquitoes in 2012. Furthermore, with the huge development of the economy and culture in Lhasa, local citizens focus more on their health status than ever before. A variety of insecticides have been adopted to protect themselves from mosquitoes bites.

Previous studies showed that the most northerly area of Cx. pipiens pipiens is about 45°N in the New World and the southerly area is about 39°S, and the usual altitude of this subspecies is lower than 3,000 m [41, 42]. In China, Cx. pipiens pipiens has been recorded only in Xinjiang Uygur Autonomous Region (northwestern China). However, previous data on the distribution of Cx. pipiens pipiens is limited: reliable identification depends on collection of males, which was not always the case. Furthermore, species and subspecies classification has been difficult because there were no populations of Cx. pipiens pipiens in Chinese laboratories. In the current study, pure Cx. pipiens pipiens (subspecies of Cx. pipiens complex) was definitively identified in urban Lhasa (an area of elevation higher than 3,600 m). This finding significantly extends present knowledge of the distribution of Cx. pipiens pipiens in China, and has important implications for the control of mosquitoes and mosquito-borne diseases in Lhasa city.

In Eastern Asia Cx. pipiens pallens transmits lymphatic filariasis and canine dogworm and may act as a vector for West Nile virus [4346]. Cx. p. pallens differs from hybrids of Cx. p. pipiens and Cx. quinquefasciatus[14]. In China, Cx. pipiens pallens has been found north of the Yangtse River [24], but not previously at an altitude of greater than 2,900 m [25].

This study has uncovered possibly extensive hybridization among subspecies of Cx. pipiens complex in Lhasa city. Natural hybridization is defined as “successful matings in nature between individuals from two populations, or groups of populations, that are distinguishable on the basis of one or more heritable characters” [47]. Combinations of this kind enhance rapid evolution and may lead to speciation [48]. According to the existing literature, recurring hybridization occurs in the Cx. pipiens complex mostly between the two most widespread species, Culex (Culex) pipiens and Cx. (Cx.) quinquefasciatus[14]. Hybrids have been reported in North America [17, 42], Argentina [19], as well as in Madagascar [49]. A multilocus genotype analysis revealed current hybridization between Cx. p. pallens and Cx. quinquefasciatus in southern Japan, Republic of Korea, and China [12, 50].

In the present study, primers specifically designed for East Asia by Smith & Fonseca were adopted. These primers were successfully designed for the identification of mosquitoes in Lhasa city. In this study, positive controls were also included, in other words, Cx. pipiens pipiens from Urumchi, Xinjiang, Cx. pipiens quinquefasciatus from Dali, Yunnan, and Cx. pipiens pallens from Beijing. The study of the indigenous populations of mosquitoes using molecular markers allowed us to confirm the occurrence of Cx. pipiens complex in Lhasa city (southwest of China). Of some interest is the discovery of hybrid populations including Cx. pipiens pipiens, Cx. pipiens quinquefasciatus, Cx. pipiens pallens and their hybrids (65.71%) in Lhasa city.

Climate change may have played a part in the arrival of mosquitoes in Lhasa. Average temperatures increased over the Tibetan plateau from 1955 to 1996 by about 0.16°C/decade [51], much more than in China generally. From 1961 to 2000, the greatest increase in daily mean temperatures in summer (June to August) in Tibet occurred in Lhasa city [52, 53]. In 2009, Tibet experienced unusually warm conditions and the maximum temperature in Lhasa reached 30.4°C, higher than the previously reported record (29.9°C in 1971). In other words, the first public reports of mosquitoes coincided with the warmest summer in Lhasa since records were first kept. It is possible that mosquitoes were introduced earlier, but numbers multiplied during the particularly hot summer of 2009. In the future, further warming is expected, and further economic development in Tibet will lead to even greater movement of freight and people. These conditions raise the risk of outbreak of mosquito-borne diseases in a population with no prior exposure to such infections [5456]. Therefore, it is urgent to strengthen the detection and monitoring of mosquito-borne diseases in the region.

Other factors, such as demographic and environmental factors, may also play a more important role in establishing the mosquito population in Lhasa. In the last 30 years, China has undergone enormous economic growth, largely due to greatly increased international trade. This burgeoning trade has triggered environmental threats from an expanding list of biological invaders and has already caused damage to China’s environment and economy. Huge construction projects, such as the Qinghai-Tibet Railway [57], could further spread invasive mosquitoes to Lhasa city [37, 58, 59]. As to urbanization, Tibet entered a stage of accelerated urbanization after 1995. The large floating population from outside Tibet has become the driving force for urban expansion and the rising urbanization rate. Tibet’s urbanization rate will be up to 43% by 2020 based on a website (http://chinatibet.people.com.cn/6828539.html). At present, Lhasa is claimed to be a modern city on the “roof of the world” with a forest of new buildings and luxury hotels, restaurants and stores. Previous research revealed that urbanization serves in the formation of appropriate habitat of culicines. In Macau, recent urbanization has provided optimal habitat for the population increase in culicines [60]. In Tanzania, urbanization resulted in some changes in mosquito populations [61]. Now, Lhasa, which has a large number of tourism resources, such as the Potala Palace, Jokhang Temple, Sara Monastery, and Barkhor Street, is a popular destination for both domestic and international travelers. By April, 2008, there were over 1,600 licensed tour guides in Lhasa according to The Chinese National Tourism Administration. Tourist aircraft or trains may carry mosquitoes to urban Lhasa and subsequently threaten the health and lives of local citizens. It was reported that labor flow and travelers are significant factors contributing to the spread of dengue virus infection and chikungunya fever [62].

In summary, our investigation provides insight into the new distribution of subspecies of Cx. pipiens complex and its hybrids in Lhasa, Tibet. The findings mentioned above have a significant implication in public health areas, both at policy making and practical levels. The multiplex PCR assay adopted in this study will be helpful to researchers and will aid vector control programs by facilitating the rapid and reliable identification of local Cx. pipiens complex and its hybrids. The future focus of the control and prevention of mosquito-borne diseases in Lhasa is West Nile virus, St. Louis encephalitis viruses, avian malaria, and filarial worms. Strengthened community health education and engagement should be conducted to better guarantee the health and life safety of local citizens. The results could provide a reference for development of varieties of strategies and measures to control mosquitoes and mosquito-borne diseases at high elevation regions in the world in future.

This study has limitations since it was planned and implemented, initially, in response to public concerns, and includes information from only two time points. However, the results indicate that mosquitoes are established in a high altitude urban setting in Tibet. Further studies are needed to confirm the continuing presence of mosquitoes, to clarify the patterns of hybridization, and to shed further light on likely origins and factors influencing their distribution and establishment in Lhasa city [63].

Conclusion

In summary, the results revealed subspecies of Cx. pipiens complex and its hybrids on two occasions, three years apart in urban Lhasa. There was no significant difference in mosquito density monitored by bed net trap and labor hour method in 2009 and 2012. Mosquitoes in the Cx. pipiens complex appear to be established in Lhasa City, TAR. Strengthened community health education and engagement should be conducted to better guarantee the health and life safety of local citizens.

Abbreviations

TAR:

Tibet autonomous region

RFLP:

Restriction fragment length polymorphism

ace-2:

Acetylcholinesterase-2.

References

  1. Lounibos LP: Invasions by insect vectors of human disease. Ann Rev Entomol. 2002, 47: 233-266. 10.1146/annurev.ento.47.091201.145206.

    Article  CAS  Google Scholar 

  2. Delatte H, Bagny L, Brengue C, Bouetard A, Paupy C, Fontenille D: The invaders: phylogeography of dengue and chikungunya viruses Aedes vectors, on the South West islands of the Indian Ocean. Infect Genet Evol. 2011, 11 (7): 1769-1781. 10.1016/j.meegid.2011.07.016.

    Article  CAS  PubMed  Google Scholar 

  3. Tolle MA: Mosquito-borne diseases. Curr Probl Pediatr Adolesc Health Care. 2009, 39 (4): 97-140. 10.1016/j.cppeds.2009.01.001.

    Article  PubMed  Google Scholar 

  4. Iturbe-Ormaetxe I, Walker T, O’ Neill SL: Wolbachia and the biological control of mosquito-borne disease. EMBO reports. 2011, 12 (6): 508-518. 10.1038/embor.2011.84.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. Benitez MA: Climate change could affect mosquito-borne diseases in Asia. Lancet. 2009, 373 (9669): 1070-10.1016/S0140-6736(09)60634-6.

    Article  PubMed  Google Scholar 

  6. Fonseca DM, Keyghobadi N, Malcolm CA, Mehmet C, Schaffner F, Mogi M, Fleischer RC, Wilkerson RC: Emerging vectors in the Culex pipiens complex. Science. 2004, 303 (5663): 1535-1538. 10.1126/science.1094247.

    Article  CAS  PubMed  Google Scholar 

  7. Diaz-Badillo A, Bolling BG, Perez-Ramirez G, Moore CG, Martinez-Munoz JP, Padilla-Viveros AA, Camacho-Nuez M, Diaz-Perez A, Beaty BJ, Munoz Mde L: The distribution of potential West Nile virus vectors, Culex pipiens pipiens and Culex pipiens quinquefasciatus (Diptera: Culicidae), in Mexico City. Parasit Vectors. 2011, 4: 70-10.1186/1756-3305-4-70.

    Article  PubMed Central  PubMed  Google Scholar 

  8. Amraoui F, Krida G, Bouattour A, Rhim A, Daaboub J, Harrat Z, Boubidi SC, Tijane M, Sarih M, Failloux AB: Culex pipiens, an experimental efficient vector of West Nile and Rift Valley fever viruses in the Maghreb region. PloS one. 2012, 7 (5): e36757-10.1371/journal.pone.0036757.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Hamer GL, Kitron UD, Brawn JD, Loss SR, Ruiz MO, Goldberg TL, Walker ED: Culex pipiens (Diptera: Culicidae): a bridge vector of West Nile virus to humans. J Med Entomol. 2008, 45 (1): 125-128. 10.1603/0022-2585(2008)45[125:CPDCAB]2.0.CO;2.

    Article  PubMed  Google Scholar 

  10. Bailey CL, Eldridge BF, Hayes DE, Watts DM, Tammariello RF, Dalrymple JM: Isolation of St. Louis encephalitis virus from overwintering Culex pipiens mosquitoes. Science. 1978, 199 (4335): 1346-1349. 10.1126/science.628843.

    Article  CAS  PubMed  Google Scholar 

  11. Farajollahi A, Fonseca DM, Kramer LD, Marm Kilpatrick A: “Bird biting” mosquitoes and human disease: a review of the role of Culex pipiens complex mosquitoes in epidemiology. Infect Genet Evol. 2011, 11 (7): 1577-1585. 10.1016/j.meegid.2011.08.013.

    Article  PubMed Central  PubMed  Google Scholar 

  12. Smith JL, Fonseca DM: Rapid assays for identification of members of the Culex (Culex) pipiens complex, their hybrids, and other sibling species (Diptera: culicidae). Am J Trop Med Hyg. 2004, 70 (4): 339-345.

    CAS  PubMed  Google Scholar 

  13. Solomon T: Flavivirus encephalitis. N Engl J Med. 2004, 351 (4): 370-378. 10.1056/NEJMra030476.

    Article  CAS  PubMed  Google Scholar 

  14. Fonseca DM, Smith JL, Kim HC, Mogi M: Population genetics of the mosquito Culex pipiens pallens reveals sex-linked asymmetric introgression by Culex quinquefasciatus. Infect Genet Evol. 2009, 9 (6): 1197-1203. 10.1016/j.meegid.2009.06.023.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Fonseca DM, Smith JL, Wilkerson RC, Fleischer RC: Pathways of expansion and multiple introductions illustrated by large genetic differentiation among worldwide populations of the southern house mosquito. Am J Trop Med Hyg. 2006, 74 (2): 284-289.

    PubMed  Google Scholar 

  16. Wang E, Ni H, Xu R, Barrett AD, Watowich SJ, Gubler DJ, Weaver SC: Evolutionary relationships of endemic/epidemic and sylvatic dengue viruses. J Virol. 2000, 74 (7): 3227-3234. 10.1128/JVI.74.7.3227-3234.2000.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Urbanelli S, Silvestrini F, Reisen WK, De Vito E, Bullini L: Californian hybrid zone between Culex pipiens pipiens and Cx. p. quinquefasciatus revisited (Diptera: Culicidae). J Med Entomol. 1997, 34 (2): 116-127.

    Article  CAS  PubMed  Google Scholar 

  18. Kothera L, Zimmerman EM, Richards CM, Savage HM: Microsatellite characterization of subspecies and their hybrids in Culex pipiens complex (Diptera: Culicidae) mosquitoes along a north–south transect in the central United States. J Med Entomol. 2009, 46 (2): 236-248. 10.1603/033.046.0208.

    Article  CAS  PubMed  Google Scholar 

  19. Humeres SG, Almiron WR, Sabattini MS, Gardenal CN: Estimation of genetic divergence and gene flow between Culex pipiens and Culex quinquefasciatus (Diptera: Culicidae) in Argentina. Mem Inst Oswaldo Cruz. 1998, 93 (1): 57-62. 10.1590/S0074-02761998000100011.

    Article  CAS  PubMed  Google Scholar 

  20. Zhao T, Baolin L: Biosystematics of Culex pipiens complex in China. Insect Sci. 1995, 2 (1): 1-8. 10.1111/j.1744-7917.1995.tb00016.x.

    Article  Google Scholar 

  21. Zhao T, Baolin L: A new record of Culex Pipiens Molestus in China and Studies of autogeny and systematics. China J Vector Bio Control. 1993, 4 (4): 241-243.

    Google Scholar 

  22. Shaman J, Day JF, Komar N: Hydrologic conditions describe West Nile virus risk in Colorado. Int J Environ Res Public Health. 2010, 7 (2): 494-508. 10.3390/ijerph7020494.

    Article  PubMed Central  PubMed  Google Scholar 

  23. Sun X, Fu S, Gong Z, Ge J, Meng W, Feng Y, Wang J, Zhai Y, Wang H, Nasci R: Distribution of arboviruses and mosquitoes in northwestern Yunnan Province, China. Vector Borne Zoonotic Dis. 2009, 9 (6): 623-630. 10.1089/vbz.2008.0145.

    Article  PubMed  Google Scholar 

  24. Lu B: Fauna Sinaca, Insecta, Diptera: Culicidae II, Volume 9. 1997, Beijing, China: Science Press

    Google Scholar 

  25. Yixing L, Minghua L, Shihong F, Weixin C, Qiyong L, Hailin Z, Wa D, Songlin H, Sang Dan La M, Ju B: Japanese Encephalitis, Tibet, China. Emerg Infect Diss. 2011, 17 (5): 934-936. 10.3201/eid1705.101417.

    Article  Google Scholar 

  26. Brownie J, Shawcross S, Theaker J, Whitcombe D, Ferrie R, Newton C, Little S: The elimination of primer-dimer accumulation in PCR. Nucleic Acids Res. 1997, 25 (16): 3235-3241. 10.1093/nar/25.16.3235.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  27. Bourguet D, Fonseca D, Vourch G, Dubois MP, Chandre F, Severini C, Raymond M: The acetylcholinesterase gene ace: a diagnostic marker for the Pipiens and Quinquefasciatus forms of the Culex pipiens complex. J Am Mosq Control Assoc. 1998, 14 (4): 390-396.

    CAS  PubMed  Google Scholar 

  28. Kasai S, Komagata O, Tomita T, Sawabe K, Tsuda Y, Kurahashi H, Ishikawa T, Motoki M, Takahashi T, Tanikawa T: PCR-based identification of Culex pipiens complex collected in Japan. Jpn J Infect Dis. 2008, 61 (3): 184-191.

    CAS  PubMed  Google Scholar 

  29. Reiter P: Climate change and mosquito-borne disease. Environ Health Perspect. 2001, 109 (Suppl 1): 141-161.

    Article  PubMed Central  PubMed  Google Scholar 

  30. Zell R: Global climate change and the emergence/re-emergence of infectious diseases. Int J Med Microbiol. 2004, 293 (Suppl 37): 16-26.

    PubMed  Google Scholar 

  31. Patz JA, Epstein PR, Burke TA, Balbus JM: Global climate change and emerging infectious diseases. JAMA. 1996, 275 (3): 217-223. 10.1001/jama.1996.03530270057032.

    Article  CAS  PubMed  Google Scholar 

  32. Rydzanicz K, Kiewra D, Lonc E: Changes in range of mosquito-borne diseases affected by global climatic fluctuations. Wiad Parazytol. 2006, 52 (2): 73-83.

    PubMed  Google Scholar 

  33. Stachowicz JJ, Terwin JR, Whitlatch RB, Osman RW: Linking climate change and biological invasions: ocean warming facilitates nonindigenous species invasions. Proc Natl Acad Sci U S A. 2002, 99 (24): 15497-15500. 10.1073/pnas.242437499.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  34. Liu XB, Liu QY, Guo YH, Jiang JY, Ren DS, Zhou GC, Zheng CJ, Liu JL, Chen Y, Li HS: Random repeated cross sectional study on breeding site characterization of Anopheles sinensis larvae in distinct villages of Yongcheng City, People’s Republic of China. Parasit Vectors. 2012, 5: 58-10.1186/1756-3305-5-58.

    Article  PubMed Central  PubMed  Google Scholar 

  35. Liu Q, Liu X, Zhou G, Jiang J, Guo Y, Ren D, Zheng C, Wu H, Yang S, Liu J: Dispersal range of Anopheles sinensis in Yongcheng City, China by mark-release-recapture methods. PloS one. 2012, 7 (11): e51209-10.1371/journal.pone.0051209.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  36. Yu C, Zhang Y, Claus H, Zeng R, Zhang X, Wang J: Ecological and environmental issues faced by a developing Tibet. Environ Sci Technol. 2012, 46 (4): 1979-1980. 10.1021/es2047188.

    Article  CAS  PubMed  Google Scholar 

  37. Li WJ, Wang JL, Li MH, Fu SH, Wang HY, Wang ZY, Jiang SY, Wang XW, Guo P, Zhao SC: Mosquitoes and mosquito-borne arboviruses in the Qinghai-Tibet Plateau–focused on the Qinghai area, China. Am J Trop Med Hyg. 2010, 82 (4): 705-711. 10.4269/ajtmh.2010.09-0649.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  38. Ruiz MO, Chaves LF, Hamer GL, Sun T, Brown WM, Walker ED, Haramis L, Goldberg TL, Kitron UD: Local impact of temperature and precipitation on West Nile virus infection in Culex species mosquitoes in northeast Illinois, USA. Parasit Vectors. 2010, 3 (1): 19-10.1186/1756-3305-3-19.

    Article  PubMed Central  PubMed  Google Scholar 

  39. Jenifer Huang MB, McBeath J: Invasive species and food security: invasive species and food security. Adv in Global Change Res. 2010, 35: 157-176. 10.1007/978-1-4020-9180-3_6.

    Article  Google Scholar 

  40. Liu XB, Liu QY, Guo YH, Jiang JY, Ren DS, Zhou GC, Zheng CJ, Zhang Y, Liu JL, Li ZF: The abundance and host-seeking behavior of culicine species (Diptera: Culicidae) and Anopheles sinensis in Yongcheng city, People’s Republic of China. Parasit Vectors. 2011, 4: 221-10.1186/1756-3305-4-221.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  41. Barr AR: Occurrence and distribution of the Culex pipiens complex. Bull World Health Organ. 1967, 37 (2): 293-296.

    PubMed Central  CAS  PubMed  Google Scholar 

  42. Barr AR: The distribution of Culex p. pipiens and C.P. quinquefasciatus in North America. Am J Trop Med Hyg. 1957, 6 (1): 153-165.

    CAS  PubMed  Google Scholar 

  43. Oda T, Eshita Y, Uchida K, Mine M, Kurokawa K, Ogawa Y, Kato K, Tahara H: Reproductive activity and survival of Culex pipiens pallens and Culex quinquefasciatus (Diptera: Culicidae) in Japan at high temperature. J Med Entomol. 2002, 39 (1): 185-190. 10.1603/0022-2585-39.1.185.

    Article  PubMed  Google Scholar 

  44. Chen L, Zhu C, Zhang D: Naturally occurring incompatibilities between different Culex pipiens pallens populations as the basis of potential mosquito control measures. PLoS Negl Trop Dis. 2013, 7 (1): e2030-10.1371/journal.pntd.0002030.

    Article  PubMed Central  PubMed  Google Scholar 

  45. Shirafuji H, Kanehira K, Nishiguchi A, Kamio T: Nationwide surveillance of West Nile virus targeting mosquitoes and dead birds from April 2004 through March 2007 in Japan. Zoonoses Public Health. 2011, 58 (3): 153-157. 10.1111/j.1863-2378.2009.01321.x.

    Article  CAS  PubMed  Google Scholar 

  46. Tsuda Y, Komagata O, Kasai S, Hayashi T, Nihei N, Saito K, Mizutani M, Kunida M, Yoshida M, Kobayashi M: A mark-release-recapture study on dispersal and flight distance of Culex pipiens pallens in an urban area of Japan. J Am Mosq Control Assoc. 2008, 24 (3): 339-343. 10.2987/5754.1.

    Article  PubMed  Google Scholar 

  47. Arnold ML: Natural hybridization and the evolution of domesticated, pest and disease organisms. Mol Ecol. 2004, 13 (5): 997-1007. 10.1111/j.1365-294X.2004.02145.x.

    Article  CAS  PubMed  Google Scholar 

  48. Jiggins CD, Salazar C, Linares M, Mavarez J: Review. Hybrid trait speciation and Heliconius butterflies. Philos Trans R Soc Lond B Biol Sci. 2008, 363 (1506): 3047-3054. 10.1098/rstb.2008.0065.

    Article  PubMed Central  PubMed  Google Scholar 

  49. Urbanelli S, Silvestrini F, Sabatinelli G, Raveloarifera F, Petrarca V, Bullini L: Characterization of the Culex pipiens complex (Diptera: Culicidae) in Madagascar. J Med Entomol. 1995, 32 (6): 778-786.

    Article  CAS  PubMed  Google Scholar 

  50. Cui F, Qiao CL, Shen BC, Marquine M, Weill M, Raymond M: Genetic differentiation of Culex pipiens (Diptera: Culicidae) in China. Bull Entomol Res. 2007, 97 (3): 291-297. 10.1017/S0007485307004968.

    Article  CAS  PubMed  Google Scholar 

  51. Xiaodong L, Baode C: Climatic warming in the Tibetan Plateau during recent decades. INT J CLIMATOL. 2000, 20 (14): 1729-1742. 10.1002/1097-0088(20001130)20:14<1729::AID-JOC556>3.0.CO;2-Y.

    Article  Google Scholar 

  52. Jun D: Change of temperature in Tibetan Plateau From 1961 to 2000. Acta Geographica Sinica. 2001, 56 (6): 682-690. (in Chinese)

    Google Scholar 

  53. Huang X: Climatic change of maximum and minimum temperature in Lasa. Meteorol Monthly. 2000, 26 (3): 47-50. in Chinese

    Google Scholar 

  54. Loevinsohn ME: Climatic warming and increased malaria incidence in Rwanda. Lancet. 1994, 343 (8899): 714-718. 10.1016/S0140-6736(94)91586-5.

    Article  CAS  PubMed  Google Scholar 

  55. Reiter P: Global warming and mosquito-borne disease in USA. Lancet. 1996, 348 (9027): 622-

    Article  CAS  PubMed  Google Scholar 

  56. Patz JA, Campbell-Lendrum D, Holloway T, Foley JA: Impact of regional climate change on human health. Nature. 2005, 438 (7066): 310-317. 10.1038/nature04188.

    Article  CAS  PubMed  Google Scholar 

  57. West JB: Completion of the railroad to Lhasa. High Alt Med Biol. 2006, 7 (4): 263-264. 10.1089/ham.2006.7.263.

    Article  PubMed  Google Scholar 

  58. West JB: A new approach to very-high-altitude land travel: the train to Lhasa, Tibet. Ann Intern Med. 2008, 149 (12): 898-900. 10.7326/0003-4819-149-12-200812160-00009.

    Article  PubMed  Google Scholar 

  59. Juliano SA, Lounibos LP: Ecology of invasive mosquitoes: effects on resident species and on human health. Ecol Lett. 2005, 8 (5): 558-574. 10.1111/j.1461-0248.2005.00755.x.

    Article  PubMed Central  PubMed  Google Scholar 

  60. Easton ER: Urbanization and its effects on the ecology of mosquitoes in Macau, Southeast Asia. J Am Mosq Control Assoc. 1994, 10 (4): 540-544.

    CAS  PubMed  Google Scholar 

  61. Bang YH, Mrope FM, Sabuni IB: Changes in mosquito populations associated with urbanization in Tanzania. East Afr Med J. 1977, 54 (7): 403-410.

    CAS  PubMed  Google Scholar 

  62. Yoshikawa MJ, Kusriastuti R: Surge of dengue virus infection and chikungunya Fever in bali in 2010: the burden of mosquito-borne infectious diseases in a tourist destination. Tropical med health. 2013, 41 (2): 67-78. 10.2149/tmh.2011-05.

    Article  Google Scholar 

  63. Trawinski PR, Mackay DS: Identification of environmental covariates of West Nile virus vector mosquito population abundance. Vector Borne Zoonotic Dis. 2010, 10 (5): 515-526. 10.1089/vbz.2008.0063.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

We thank Professor Peng Bi and Abdelrafie M. Makhawi for comments on the manuscript. This study was funded by the National Basic Research Program of China (973 Program) (Grant No. 2012CB955504) and the National Natural Science Foundation of China (NSFC) (Grant No.81273139).

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Qiyong Liu or Cordia Chu.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

QL, XL, AW and CC planned the project and wrote the paper. QL, XL, C, P, FW, B, LB, YG, D, GL, JW, SS, D and X conducted the field survey. XL, LL, L and HW carried out the multiplex PCR assay in the lab. XL, QL, CC, and AW contributed to data analysis. All authors read and approved the final manuscript.

Qiyong Liu, Xiaobo Liu, Cirendunzhu contributed equally to this work.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

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.

Reprints and permissions

About this article

Cite this article

Liu, Q., Liu, X., Cirendunzhu et al. Mosquitoes established in Lhasa city, Tibet, China. Parasites Vectors 6, 224 (2013). https://doi.org/10.1186/1756-3305-6-224

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1756-3305-6-224

Keywords