ORIGINAL RESEARCH article

Front. Vet. Sci., 04 July 2019

Sec. Parasitology

Volume 6 - 2019 | https://doi.org/10.3389/fvets.2019.00207

Marmots and Yersinia pestis Strains in Two Plague Endemic Areas of Tien Shan Mountains

  • 1. Department of Natural Sciences, Issyk-Kul State University, Karakol, Kyrgyzstan

  • 2. Karakol Anti-plague Department, Republic Center of Quarantine and Dangerous Infections, Karakol, Kyrgyzstan

  • 3. Republic Center of Quarantine and Dangerous Infections, Bishkek, Kyrgyzstan

  • 4. Reference-Laboratory, Kazakh Scientific Center of Quarantine and Zoonotic Diseases, Almaty, Kazakhstan

  • 5. Institute of Evolutionary Sciences, University of Montpellier, Montpellier, France

  • 6. Department of Microbiology and Immunology, Medical Branch, University of Texas, Galveston, TX, United States

  • 7. Division of Vector-Borne Diseases, Centers for Disease Control and Prevention, Fort Collins, CO, United States

Abstract

The main purpose of this study was to clarify the role of gray marmots (Marmota baibacina) in the long-term maintenance of highly virulent strains of Yersinia pestis in two plague endemic foci of the Tien Shan Mountains in Kyrgyzstan. We present data from regular observations of populations of M. baibacina and small rodents cohabiting with marmots in the mountainous grasslands of the Sari-Dzhas (east of Issyk-Kul Lake) and the Upper-Naryn (south of Issyk-Kul Lake) natural foci. During 2012–2017, an abundance of marmots and their ectoparasites (fleas and ticks) was significantly higher in Upper-Naryn comparing to Sari-Dzhas, although there were no differences in a number and diversity of small rodents cohabiting with marmots. The plague bacterium was detected either in marmots or in their ectoparasites collected during 4 of 6 years of observation in Sari-Dzhas and during 2 of 4 years of observation in Upper-Naryn. Plague was found in three sectors situated closely to each other in Sari-Dzhas and in 1 of 8 repeatedly surveyed sectors in Upper-Naryn. During 6 years, we isolated 9 strains of Y. pestis from marmots, two from their fleas Oropsylla silantiewi, one from an unidentified tick, and one from the gray hamster (Cricetulus migratorius). All plague strains isolated from the rodents and their ectoparasites in this study were similar to Antiqua biovar specific for marmots. The results indicate that plague can circulate continuously in the Tien Shan Mountains in populations of gray marmots and their ectoparasites with a facultative involvement of other rodent species after significant changes in rodent communities that happened in Kyrgyzstan during the previous two decades. The simultaneous field survey of two natural foci of plague, Sari-Dzhas, and Upper-Naryn, would be important for further analysis of circulation of Y. pestis strains belonging to Antiqua biovar in the Tien Shan Mountains.

Introduction

One of the most active plague endemic areas in the world is located in Kyrgyzstan, Central Asia. There are two independent foci of plague in the Issyk-Kul province: the Sari-Dzhas (east of Issyk-Kul Lake) and the Upper-Naryn (south of Issyk-Kul Lake) sub-regions of the Tien Shan natural focus of plague (1, 2). Several plague outbreaks involving natural vectors, including the gray marmots (Marmota baibacina) and other rodents, have been registered there since 1907 (3). Most of the 5,000 km2 Sari-Dzhas natural plague focus is located in Kyrgyzstan (4,250 km2), with a part extending to Kazakhstan (750 km2). The 8,000 km2 Upper-Naryn area is located in the Naryn and Issyk-Kul regions of Kyrgyzstan.

The first plague outbreak with high human mortality was recorded in the Kyrgyzstan part of Tien Shan focus in 1907, and the last big outbreak of plague occurred in 1928 in the Bash-Kaindi settlement of Atbashi district resulting in 54 deaths (Upper-Naryn mezofocus) (4). Sporadic human cases were also seen in 1965 and 1982. In 2013, a 15-year-old boy who ate marmot meat died from plague in Kyrgyzstan. The Kyrgyz Ministry of Health established a temporal quarantine in parts of the country's mountainous northeast (5), where the risk of plague is well-known for a century. The endemic areas were investigated previously starting from 1942, and the last activity of plague was reported there in 1983 (4, 6, 7).

Gray marmots (M. baibacina) are well-known hosts of multiple zoonotic diseases, out of which plague presents the greatest danger to people (810). The gray marmots prefer grasslands and shrubs and avoid woodlands, including even the forest edges and forest-steppe areas with a tree cover <10% (11). Marmot fleas are actively involved in transmission of plague pathogen between animals (1214). Marmots can host the following species of fleas: Oropsylla silantiewi, Rhadinopsylla li ventricosa, Ceratophyllus lebedewi, and Pulex irritans; the first two flea species are host-specific for marmots (15, 16).

Human activities, such as agriculture, hunting, recreation, and degradation of natural habitats, may dramatically influence plague manifestations (17). Changes in landscape potentially can lead to altering the rodent and flea communities that in turn affect plague transmission cycle (18). Plague in Asian natural systems is commonly spatially stable and corresponds with distribution of primary rodent hosts that can support a strong association of specific strains of Yersinia pestis in mammalian communities (19). Marmots in mountainous plague foci of Central Asia carry a specific strain of Y. pestis (20). Disturbance of evolutionary relations between plague pathogen and their rodent hosts, such as marmots, may result in unpredictable dynamics of the infection.

The collapse of the Soviet Union and subsequent separation of anti-plague station in Kyrgyzstan from the federal anti-plague system led to dramatic reduction of plague investigations and plague control measures. During the last couple of decades, the density of marmots was significantly reduced as a result of chemical suppression in the 1960s, increased hunting of marmots, habitat destruction, and climate changes (68, 21). This raises the question whether the changes in rodent communities in these areas affected marmots playing a leading role in circulation of plague pathogen. Thus, the present study aimed to clarify the role of gray marmots in the long-term maintenance of the highly virulent strains of Y. pestis in two plague foci of the Tien Shan Mountains in Kyrgyzstan. The human case in 2013 attracted attention to the situation that plague is endemic to the territory of Kyrgyzstan and the urgent need to obtain new information about the status of the epidemiological situation in these plague foci.

Materials and Methods

Fieldwork

In June–August 2012–2017, the Karakol Anti-Plague Department (KAPD) of the Republic Center for Quarantine and Dangerous Infections (RCQDI) organized field surveys within the Sari-Dzhas and Upper-Naryn foci. The fieldwork was organized as described by Weaver et al. (22) and Sariyeva et al. (6, 7). All work with wild and laboratory animals and plague strains was conducted in accordance with the regulations and protocols approved by the Ministry of Health of Kyrgyzstan (23) in 2015. The procedures were similar to those described by Aytkuluyev (24) and Ezhlova et al. (25). The animal work in the field was performed according to the Regulations approved by the State Agency for Environmental and Forest Protection of the Kyrgyz Republic (details of permits are presented in Table 1). Each permit was issued for collecting a certain number of animals (400 marmots, 300–400 rodents of other species, and excavating 10 nests of marmots and rodents to look for nest parasites) and used during a fixed period of time (from June 1 to August 30) in specific areas within plague foci.

Table 1

YearSari-DzhasUpper Naryn
Number of licenseData of issuePermitted number of animalsNumber of licenseData of issuePermitted number of animals
2012000078-KC26.04.2012400 marmots, 400 small rodents
2013000110-KC06.05.2013400 marmots, 300 s.r.000111-KC06.05.2013300 marmots, 300 small rodents
2014000144-KC06.05.2014400 marmots, 300 s.r.000145-KC06.05.2014400 marmots, 300 s.r.
2015000168-KC22.05.2015400 marmots, 300 s.r.000169-KC22.05.2015400 marmots, 300 s.r.
2016000173-KC17.05.2016400 marmots, 300 s.r.000174-KC17.05.2016400 marmots, 300 s.r.
2017000003-KC05.03.2017400 marmots, 300 s.r.000004-KC05.03.2017400 marmots, 300 s.r.

Numbers and dates of permits for trapping of animals for epizootological research in 2012–2017.

Trapping Rodents

At the beginning of the annual anti-plague field work, zoologists visually estimated a number of marmots within study sites and accordingly planned trapping efforts in each sector (10 km2). The number of installed traps depends on the average density of marmot population in each sector. Around 8–10 traps were set next to marmot burrows each day at early morning. As the weather in May–June in high-altitude areas of Tien Shan is highly variable with snow, rains, and cold weather, the trapping was conducted during each sunny day. The marmot traps were set at the entrance of burrows and removed by evening. One trap remained at each burrow for half a day. The burrows for trapping were chosen based on external signs such as presence of fresh litter at the entrance, presence of freshly excavated and well-rammed entrances, grass coverage, presence of marmot footprints, and paths around the burrow. The traps used for capturing marmots were 20 cm wide and 8 cm high (Figure 5). The traps were attached to 50-cm-long metal rods with a strong metal chain. The rods are driven into the rocky soil to their full depth to prevent marmots or their predators from dragging the traps away. The metal parts were masked by horse manure to reduce smell of metal.

A number of small rodents, such as mice, voles, and hamsters, were estimated by counting animals captured by handmade snap traps set in transects (1, 25). The snap traps were also sprinkled with soil and baited with dry bread soaked in vegetable oil, vegetables, or fruit. We placed traps in the evening and checked them in the early morning. If there was a captured rodent, we transferred the trap to another spot. The predators were estimated by trapping and visual observations. Each trapped animal was morphologically identified by species. The captured marmots were euthanized by cervical dislocation (25). Then, the animal was wrapped in two linen sacs with label (date, place, number of sector, and name of catcher), the inlet was tightly wrapped so that the ectoparasites could not escape, and then the animal was placed in a plastic bag, in a canvas bag, and finally brought to camp.

Parasitological Analysis

Ectoparasites of marmots were collected by combing the captured animals (“body” ectoparasites) and by collecting “off-host” ectoparasites from rodent burrows and nests (25). Ectoparasites were then identified using entomological keys (26) and placed in labeled glass tubes with ether. Inside burrows, dry grass bedding from the nesting chamber was carefully removed with a wire and was searched for fleas and other ectoparasites. Fleas of one species from a single animal were pooled and triturated in saline solution.

Bacteriological Analysis

The triturated flea suspension was inoculated on Hottinger agar with pH 7.2 (27) (bacteriological approach). The inoculation pool consisted of 20 ectoparasites of one species, collected from the same sector during several days and stored before inoculation in glass tubes without ether. If the animal exhibited any pathological manifestation, the ectoparasites from the pool were inoculated individually. Identification of the isolated bacteria was carried out by standard microbiological microscopy of smears and by Pokrovsky's test for Y. pestis using both pseudo-tuberculosis bacterial phage and bacteriophage for Y. pestis (L-413-C) produced by the Kazakh Scientific Center for Quarantine and Zoonotic Diseases (KSCQZD). Additional diagnostic methods included serological indirect hemagglutination assay with both erythrocyte–immunoglobulin and erythrocyte–antigen diagnostic reagents produced by KSCQZD. The tissues of internal organs of sampled marmots (liver, lung, spleen, lymph nodes, and blood from the heart) were sterilely taken with forceps and immediately placed on an agar plate by directly touch. A Petri dish was divided into three segments to separately plate lungs, livers, and spleens. Each plate was marked and placed in an incubator at 37°C. We used individual and pool methods of inoculation. The individual method was used for ectoparasites collected from carcass of marmot or marmot with visual pathological abnormalities relevant to plague (pathology of inner organs, sick marmots). For inoculation of pooled tissue suspensions, we combined tissue pieces from five marmots collected on the same day from the same sector. Then, the organ suspension was inoculated on agar plate with a bacterial inoculating loop. After inoculation, tissue samples were placed in liquid nitrogen (Dewar flask) and transported to the laboratory in Karakol for further confirmation.

Genetic Analysis

The genetic analysis was performed using two types of material: whole-genome DNA samples and MLVA fragments (Table 2). The 14 studied Y. pestis strains were obtained from the RCQDI collection. The strains were isolated during the last several years in the territory of Sari-Dzhas and Upper-Naryn plague foci (Table 2). The suspension of Y. pestis cultures was heated at 100°C for 20 min and centrifuged at 12,000 rpm for 2 min. This allowed the inactivation of the pathogens and the release of the DNA in the cells. The DNA samples were further used for genotyping. The control MLVA amplification fragments of the reference strains CO92, Pestoides F, KIM10+, and Nepal 516 of Y. pestis were obtained from the University of Texas Medical Branch, Texas, USA. In addition, nucleotide sequences of nine Y. pestis species (GenBank accession numbers: CP010023, CP010247, CP006751, CP009935, AE017042, CP006806, CP000308, CP006794, and CP002956) and three Yersinia pseudotuberculosis strains (GenBank accession numbers: CP009712, CP008943, and CP001048) were used for phylogenetic analysis. The studied strains belonging to the Y. pestis species was confirmed using the “Pest-Quest” PCR assay (“Master-Gene” Limited Liability Partnership, Almaty, Kazakhstan). The genotyping of several Y. pestis strains was performed by the methods of Multi-Locus VNTR Analysis (MLVA) and Melt Analysis of Mismatch Amplification Mutation Assay (Melt-MAMA). For MLVA assay, seven VNTR (Variable Number Tandem Repeats) loci were studied using conventional PCR and agarose gel electrophoresis as described by Le Flèche et al. (28). For Melt-MAMA assay, three sets of PCR primers designed and produced by KSCQZD were used to identify three SNP loci, described by Morelli et al. (29). The selected SNP loci allowed differentiation of four branches of Y. pestis bv. Antiqua (0.ANT1, 0.ANT2, 0.ANT3, and 3.ANT). The Melt-MAMA analysis was performed as described by Birdsell et al. (30). Each primer set consisted of two forward primers (“ancestral” and “derived”) and a universal reverse primer. One of the forward primers was linked to the so-called GC-clamp, three to four repeats of a GGGGC motif, which increased the melting temperature of the corresponding amplicons and made it possible to differentiate the alleles. Phylogenetic analysis of the studied strains was carried out using the PAUP 4.0 software and the UPGMA (Unweighted Pair Group Method with Arithmetic Mean) algorithm.

Table 2

Strain codeObject, site, and period of strain isolation
KG-1Human outbreak, the Sari-Dzhas, 2013
KG-2Human outbreak, the Sari-Dzhas, 2013
KG-3Gray marmot, the Sari-Dzhas, 2014
KG-4Gray marmot, the Upper-Naryn natural plague focus, Ishtyk-Akshiyrak site, 2015
KG-5Gray marmot, the Upper-Naryn natural plague focus, Ishtyk-Akshiyrak site, 1963
KG-6Ectoparasites, there was a suspicion that this strain was Y. pseudotuberculosis, 1962
KG-7Gray marmot carcass, the Sari-Dzhas natural plague focus, Enylchek-Kaindy site, 2016
KG-8Ectoparasites—ticks from the carcass of gray marmot, the Sari-Dzhas natural plague focus, Enylchek-Kaindy site, 2016
KG-9Ectoparasites—fleas from the carcass of gray marmot, the Sari-Dzhas natural plague focus, Enylchek-Kaindy site, 2016
KG-10Gray marmot, the Sari-Dzhas natural plague focus, Enylchek-Kaindy site, 2016
KG-11Gray marmot, the Upper-Naryn natural plague focus, Ishtyk-Akshiyrak site, 1957
KG-12Gray marmot, the Upper-Naryn natural plague area, Ishtyk-Akshiyrak site, 1962
KG-13Gray marmot, the Upper-Naryn natural plague area, Ishtyk-Akshiyrak site, 1960
KG-14Gray marmot, the Upper-Naryn natural plague area, Ishtyk-Akshiyrak site, 1983

Strains of Y. pestis isolated in the Sari-Dzhas and Upper-Naryn plague foci used for DNA analysis.

Mapping

The map of Sari-Dzhas and Upper-Naryn areas of plague was prepared using GIS MapInfo Professional 7.8. For their construction, we used the geographic objects digitized from topographic maps of a scale of 1:100,000, as well as the Digital Elevation Model GDEM2. The maps are constructed in UTM—the Mercater projection (WGS 84). The boundaries of the foci are plotted along the boundaries of the sectors (21).

Results

Sari-Dzhas Plague Focus

Rodents

Overall, in the Sari-Dzhas plague focus for the period 2012 to 2017, we surveyed a total of 3,700 km2, out of which 500 km2 were surveyed repeatedly in 2012, 2013, 2014, and 2016 (Figure 1; Table 4). In total, we trapped and analyzed 972 marmots (M. baibacina), 442 narrow-headed voles (Microtus gregalis), 61 wood mouse (Apodemus uralensis), 40 gray dwarf hamsters (Cricetulus migratorius), and a small number of other mammals: the large-eared pika (Ochotona macrotis), the hare-tolai (Lepus tolai), and the beech marten (Martes foina) (Table 3). Additionally, we found carcasses of six dead rodents.

Figure 1

Table 3

FocusYearMonthArea and number of sectorsNumber of caught marmotsNumber of other mammals trappedNumber of collected ectoparasitesY. pestis isolated from marmotsY. pestis isolated from ectoparasitesY. pestis isolated from other rodents
Sari-Dzhas2012June–July800 km2, 8 sectors218Microtus gregalis−42,
Apodemus uralensis−47,
Cricetulus migratorius−12 + carcass,
Lepus tolai−2 (total 103)
240 fleas, 459 ticks51 (gray hamster)
2013June400, 4 sectors70Microtus gregalis−30,
Apodemus uralensis−6,
Cricetulus migratorius−10 (total 46)
10 fleas, 137 ticks
2014July400, 4 sectors197Microtus gregalis−42,
Ochotona macrotis−1,
Apodemus uralensis−1,
Cricetulus migratorius−11 (total 55)
84 fleas, 155 ticks, 6 lice1 from fleas
2015July–August600, 6 sectors177Microtus gregalis−129,
Martes foina−5 (total 134)
175 fleas, 42 ticks, 16 lice
2016June–July800, 8 sectors180Microtus gregalis−124,
Apodemus uralensis−3,
Cricetulus migratorius−1 (total 128)
145 fleas, 344 ticks, 9 lice21 from ticks
2017June–July700, 7 sectors130Microtus gregalis−75,
Apodemus uralensis−4 (total 79)
124 fleas, 15 ticks, 125 lice
Upper Naryn2013June–August800, 8 sectors154Microtus gregalis−47,
Cricetulus migratorius−13,
Mustela eversmanni−4 (total 64)
268 fleas, 1045 ticks, 127 lice
2015June–July800, 8 sectors260Microtus gregalis−107,
Lepus tolai−2,
Cuon alpinus−1,
Mustela eversmanni−2 (total 112)
489 fleas, 113 ticks, 52 lice2
2016June–August700, 7 sectors190Microtus gregalis−150,
Alticola argentatus−1,
Vulpes vulpes−1 (total 152)
261 fleas, 233 ticks, 2 lice1 from fleas
2017June–July600, 6 sectors224Microtus gregalis−38,
Lepus tolai−2,
Vulpes vulpes−2 (total 42)
311 fleas, 593 ticks, 38 lice

Field survey of the Sari-Dzhas and Upper Naryn sites of the Tien Shan high-altitude focus of plague in 2012–2017.

From all sampled marmots, we collected and tested 778 fleas, 1,152 ticks, and 156 lice (Table 3). The average number of sampled marmots ranged from 17 to 27 animals per sector in 2012 and 2013 (Table 4). The highest number of trapped marmots (49 ± 7.62) was recorded in 2014 in a forested area with preferable conditions such as vegetation and humidity. In 2015–2017, the number of collected marmots was similar to 2012. Among small rodents, the most common species was M. gregalis, which is widely distributed across all observed territories (Tables 3, 4). The density of population of M. gregalis and A. uralensis varied significantly. A reverse correlation between abundance of these two species was observed. A high number of M. gregalis was associated with limited or no A. uralensis (Table 4).

Table 4

Years of observation#sectorTime period of observationMarmota baibacina
(trapped)
Marmota baibacina
(carcass)
MicrotusgregalisApodemus uralensisOchotona macrotisCricetulus migratoriusLepus tolaiMartes foinaIsolated strain of Y. pestis
2012322440633429 May−3 June28111
322440751229 May−2 June18131
31244075216–11 June43143
31244063436–11 June26164 + carcass
312440634413–15 June1863
312440643314–15 June7
312440634218–29 June532*33121*15
312440641319–26 June2575
Average:27.25±10Total 4247Total 12+carcass2
2013 (human case)312440622222 August1
201331244063424–11 June28214
31244064138–10 June1396
312440634411–18 June14
312440634310–17 June156
Average:17.5±5Total 306Total 10
201431244063119–20 July582*231163
312440631216–21 July88195
312440513319–21 July29
312440514319–21 July22
Average:49.25±23Total 4211Total 11
2015312440523327 July−5 August54203
312440523427 July−16 August52761
31244052438–11 August30131
31244052419–11 August120
31244051445–6 August611
312440632113–16 August239
Average:29.5±15Total 1295
2016312440624228 June−6 July398231
312440624130 June−2 July160
312440631330 June−5 July2526
31244062232–9 July250
31244062245–6 July60
31244062328–9 July291*163
312440631115–19 July31
312440504318–22 July9
Average:22.5±9Total 12431
2017312420512315–29 June197
312440514116–19 June186
312440505217–21 June13112
312440503423–27 June27231
312440512127 June−6 July34211
31244051422–4 July197
Average:21.67±5Total 754

The number of trapped mammals in open stations of the Sari-Dzhas natural plague focus in 2012–2017.

*

Object of isolation of Y. pestis. Bold values indicate average or total meanings.

Comparing the total number of collected marmots in Sari-Dzhas for 2013 and 2015–2017 with that in Upper-Naryn, the last one has much higher value on the whole area observed in this period (Figure 2). At the same time, there were no differences in the total number of trapped rodents and other mammals in these periods (2013, 2015–2017) between the two plague foci (Figure 2).

Figure 2

Ectoparasites

The number of ticks collected from marmots was three to five times higher than fleas collected within all sectors of both investigated plague foci (Tables 5, 7; Figure 2). In the Upper-Naryn natural plague area, the average number of fleas and ticks was significantly higher than in the Sari-Dzhas in all years, except ticks in 2016 (Figure 2).

Table 5

Years of observation#sectorNumber of analyzed animal (M. baibacina)Number of fleasO. silantiewiRh. li ventricosaNumber of ticksNumber of lice
2012322440633439191995
322440751218434360
312440752143353592
312440634326444462
31244063441823237
3124406433722221
312440634253121295
3124406413252119236
Total:2292192172459
20133124406342285570
3124406413134410
3124406344141146
31244063431511
Total701010137
201431244063115815132622
31244063128818162814
31244051332911
3124405143228812
Total:197423841556
2015312440523354434217
312440523452363516
31244052433077
3124405241121111
312440514466425
312440632123282824
Total:177131127442
2016312440624239909091
3124406241161113
312440631325422105
3124406223255589
31244062246888
3124406232295537
312440631131
31244050439
Total:18011310310343
20173124205123191212
31244051411866
31244050521366
31244050342777
31244051213488
31244051421945351015
Total:13084741015

The number of “body-ectoparasites” collected from trapped marmots in open stations of the Sari-Dzhas natural plague focus in 2012–2017.

Isolation of Y. pestis

In 2012 after the 29-year period between outbreaks, one observed an acute epizootic of plague in its primary (M. baibacina) and secondary carriers (C. migratorius). In total, five strains of Y. pestis were isolated during the current study: two strains were from found carcasses of marmots, two were from trapped marmots, and one strain was from C. migratorius (Tables 3, 4). All strains were isolated from marmots trapped or found in one sector from eight studied in 2012 within the southeast part of the Sari-Dzhas focus (Figure 4).

In 2013, the human case of plague was registered on the north-west part of the Sari-Dzhas focus, sector #3124406222 (Table 4; Figure 4). In 2014, the plague epizootic in marmots with involved host-specific fleas (O. silantiewi) was in the neighbor sector #3124406311. Three strains of Y. pestis were isolated from found carcasses of marmots and fleas, collected from them, and confirmed serologically.

In 2016, the field study was conducted near the area observed in 2014–2015, with one sector overlapping between surveys (Figures 1, 4). As a result of this investigation, three strains of plague pathogen were isolated—two from M. baibacina (one strain from the pool inoculation of marmots captured in this area and second from a carcass of marmot found in another sector, Table 4). The third strain was isolated from pooled unidentified mites collected from a plague-positive carcass.

In 2013 and 2015, different areas of the Sari-Dzhas plague focus were surveyed with plague-negative result (Figure 1).

Three strains were isolated in the Upper-Naryn: two from marmots in 2015 and one from ectoparasites collected from marmots in 2016. All three strains were isolated in the same sector #3224407344 (Figure 4; Table 6).

Table 6

Years of observation# sectorTime period of observationMarmota baibacinaMicrotusgregalisApodemus uralensisOchotona roleyCricetulus migratoriusLepus tolaiCuon alpinusMustela eversmanniAlticola argentatusVulpes vulpesIsolated strain of Y. pestis
2013322430844310–15 June1012
322430843412–17 June28111
322430844414–18 June14
322430862119–24 June164
322430961222–29 June1831
322430843328 June−2 July165
32243096113–6 July2964
32243084317–12 July23822
Total:15442134
2015322440734315–21 June38752
322440851322–28 June629
322440851229 June−4 July4941
32244073445–9 July8819112 from marmots
322440743110–13 July90
322440734211–13 July90
322440743312–13 July50
Total:260107212
201632244073447–10 June9911211 from fleas
322440743111–15 June309
322440734212–17 June1081 (carcass)
322440734316–20 June2310
322440743320–25 June9
322440851326–30 June1111
322440851230 June−5 July8
Total:1901501
2017322440851428 May−6 June5817
322440852230 May−6 June521511
32244085122–7 June363
32244085233–10 June5121
322440852210–15 June1611
322440852416–25 June11
Total:2243822

The number of trapped mammals in open stations of the Upper Naryn natural plague focus in 2013–2017.

Table 7

Years of observation#sectorNumber of analyzed animalNumber of fleasO. silantiewiRh. li ventricosaNumber of ticksNumber of lice
20133224308443108444540
3224308434288823843
32243084441476
3224308621164416144
322430961218503911159
322430843316595
32243096112939354167
32243084312342192399
Total:154151114421040127
20153224407343382121119
322440851362999913
3224408512495351260
322440734488625483313
3224407431955
32244073429221
32244074335334
Total:2602452351011232
201632244073449940364572
322440743130221210128
3224407342102212101
322440734323137623
32244074339
3224408513111991013
3224408512810739
Total:19012683432312
201732244085145818218214617
32244085225225252475
32244085123635351043
322440852351121218
3224408522161313514
3224408524113315
Total:22427027058129

The number of “body-ectoparasites” collected from trapped marmots in open stations of the Upper-Naryn natural plague focus in 2013–2017.

The Upper-Naryn Plague Focus

During 2013–2017, a total of 2,900 km2 were investigated in the Upper-Naryn focus, out of which we surveyed 700 km2 repeatedly in 2015, 2016, and 2017 (Figure 1; Table 3). In total, 828 marmots, 342 narrow-headed voles (M. gregalis), 13 gray dwarf hamsters (C. migratorius), and a small number of other mammals: six steppe polecats (Mustela eversmanni), four tolai hares (L. tolai), one red fox (Vulpes vulpes) were trapped and screened (Table 6). Additionally, we found the carcass of a dead fox.

From all sampled marmots, we collected and tested 1,329 fleas, 1,984 ticks, and 219 lice (Table 3). As a result, in 2015–2016, a plague epizootic was recorded in the same sector (Figure 4; Table 6). Y. pestis strains were isolated from M. baibacina and their specific fleas—O. silantiewi.

Results of the Genetic Analysis

The application of the “Pest-Quest” PCR assay confirmed that all the studied strains from Sari-Dzhas plague focus belonged to the Y. pestis species. Based on the results of the MLVA-7 genotyping, all the studied Y. pestis isolates were presented by nine genotypes differing in one or several VNTR loci. Five strains (KG-1–KG-5) shared the same MLVA-7 profile (probably clones of the same strain). The software we used for phylogenetic analysis showed them not as one branch, but as a bunch of sub-branches (Table 8). All 14 strains were assigned to biovar Antiqua that was considered to be typical for marmot strains (Figure 3). According to the Melt-MAMA analysis, most of the strains belonged to the branch 0.ANT2, while one strain (KG-14) was apparently a member of the branch 0.ANT3 (Table 9). These two branches were among those reported by Eroshenko et al. (31) in Kyrgyzstan (along with branches 0.ANT5, 0.PE4t, and 2.MED1).

Table 8

Strain of Y. pestisMLVA−7Genotypes
ms01ms04ms06ms07ms46ms62ms70
KG_1264213547174259276164Genotype 1
KG_2264213547174259276164Genotype 1
KG_3264213547174259276164Genotype 1
KG_4264213547174259276164Genotype 1
KG_5264213547174259276164Genotype 1
KG_6246230547174259294146Genotype 2
KG_7246196546174252294146Genotype 3
KG_8246213546174259294146Genotype 4
KG_9246213546174259303146Genotype 5
KG_10246213546174259294146Genotype 4
KG_112461790174259276137Genotype 6
KG_12246196546174259294146Genotype 7
KG_13246179546184259276137Genotype 8
KG_14228213546174293240155Genotype 9
YP_EV_76210230606184252222155Genotype 10
YP_3770210179547174245285119Genotype 11
YP_Angola210162362134259276137Genotype 12
YP_91001210196303184252258128Genotype 13
YP_Antiqua192179487164238231209Genotype 14
YP_1045210196305174322303155Genotype 15
YP_790192213545174322249137Genotype 16
YP_A1122210230606184252294146Genotype 17
YP_Pest_G210179303174245312119Genotype 18
YP_Pest_B192179547174252240137Genotype 19
YP_KIM10192196305164378276146Genotype 20
YP_Nepal516192196185174259222137Genotype 21
YP_Pest_F210179487194245294119Genotype 22
YP_CO92228230606184252240146Genotype 23
Pt_IP329531921451088144266330155Genotype 24
Pt_69041921457861440285155Genotype 25
Pt_PB1192164606154301267200Genotype 26
Pt_28411742130154252240128Genotype 27
Pt_4331742130154259231137Genotype 28

Amplicon sizes of Y. pestis strains on 7 VNTR loci.

Figure 3

Figure 4

Figure 5

Table 9

StrainsLocus s87 (T/G)Locus s332 (G/T)Locus s645 (G/T)BiovarEstimated branch
KG-1–KG-13GGGAntiqua0.ANT2
KG-14GTGAntiqua0.ANT3

Results of Melt-MAMA analysis.

Discussion

The greater part of both Sari-Dzhas and Upper-Naryn plague foci were studied simultaneously in 2013 and 2015–2017. If an epizootic of plague was observed once, the area was repeatedly examined in the succeeding years. This allows us to compare two independent plague origin areas for its epizootological status during the same time period. In the Tien Shan mountain focus of plague, the gray marmot is known as the main carrier of plague pathogen, with other rodents serving as secondary hosts (2, 8, 32, 33). The spatial distribution of marmots within plague focus is related mainly to local landscape and climatic conditions, as well as to different human activity–animal grazing, hunting, and tourism that are intensively developing during the last decade, but unequally presented in different sectors (21). Actually, the high-altitude pastures for horses and sheep (so called “syrts”) are more distant from human settlements in the Upper-Naryn area than in Sari-Dzhas. A higher distance to humans may favor the increase in density of marmots in Upper-Naryn compared to Sari-Dzhas where most of the syrts are tightly used for summer seasonal pastures and for international tourism. The disturbance of the environment may affect the marmots leading to a decrease in their population sizes. The number of other small rodents and their variability and density are approximately equal in simultaneously studied years in both plague foci. Prevalence of M. gregalis confirms the role of this species as a potential secondary host of plague pathogen in mixed ecosystems of high-altitude Tien Shan focus. Isolation of one Y. pestis strain from C. migratorius in the Sari-Dzhas focus confirms the role of small mice-like rodents in the epidemiology of plague. Previously, plague strains were isolated in the neighboring Aksay focus from M. gregalis (1968) and C. migratorius (1983–1984) and in the Upper-Naryn focus from Alticola argentatus and C. migratorius (3). Similar observations were reported in the high altitude of Altay and Tuva plague foci of Russia (34) and the North Aral sandy plague focus of Kazakhstan (35). Increasing density of rodents co-inhabiting with the marmots altogether with an increased number of their specific ectoparasites as observed in the Upper-Naryn focus in 2013–2017 could be a sign of potential activity of plague epizootics in this area of Kyrgyzstan. The isolation of 9 of 13 strains from marmots in 2012–2017, 1 strain from other rodents, and 3 strains from their ectoparasites is in favor of this hypothesis. Both independent plague foci (Sari-Dzhas and Upper-Naryn), located on the border of Kyrgyzstan, Kazakhstan, and China, are among the most active natural high-altitude foci of Central Asia. This territory is characterized by special climatic conditions, different relief with many heavily rugged canyons, mountain river valleys, and patches with specific flora and fauna with a high degree of biological diversity (36). Apparently, these factors give optimal conditions for the long-term circulation of Y. pestis biovar Antiqua (0.ANT2, 0.ANT3) in the populations of its natural host (31, 37). In total, 462 strains of Y. pestis were registered in the Sari-Dzhas area from 1944 to 1976 (4, 6, 7). Ectoparasites are actively involved in the epizootic process, in particular O. silantiewi and Rh. li ventricosa fleas specific to M. baibacina, as well as ixodes mites. The Sari-Dzhas area was very active in the 1940s to 1980s, and then mass disinsection of marmot populations with dichlorodiphenyltrichloroethane (DDT) insecticide and the decrease in the number of animals as a result of human activities (special extermination conducted in 1960s, hunting) led to a significant decrease in epizootic activity. Whereas, human activity (agriculture development, tourism, and hunting) changed environment during this period significantly, an obvious activation of epizootic plague activity within the area was observed after 30 years. Overall, the obtained results confirm ongoing epidemiological risk and vulnerability of the territory to plague. In such situation, it is necessary to strengthen ecological and epidemiological monitoring and control over the entire endemic area in order to preserve safety of local populations.

Statements

Ethics statement

All manipulations with wild and laboratory animals and handling plague strains were conducted according to the protocols approved by the Kyrgyz government regulations that are related to guidelines of Russia, Kazakhstan (1, 25). The animal work was performed according to the Regulation of the State Agency for Environmental and Forest Protection Facilities at the Government of the Kyrgyz Republic. Animal work in the field was conducted according to the following protocols: protocol 000144-KC for trapping marmots, other rodents, and analysis of their nests approved on May 6, 2014; protocol # 000168-KC approved on May 22, 2015; protocol #000173-KC approved on May 17, 2016; and protocol #000003-KC approved on March 5, 2017. Each permit to collecting a certain number of animals (400 marmots, 400 rodents of other species, and 10 nests of marmots and rodents) was issued and used during the fixed period of time (from June 1 to August 30).

Author contributions

GS initiated the study and wrote the first draft. GB organized and conducted the field studies and performed zoological and entomological analysis. RMu performed agreement with the state agencies for animal catching and analysis. SA provided strains for DNA analysis and genotyping. BK and AA performed molecular-genetic analysis of strains. AS performed GIS-mapping and database design. ZS performed bacteriological analysis and spatial epizootical analysis. AD performed the bacteriological analysis. ZA analyzed the environmental effect on population structure. RMa provided consultative assistance. SM provided collaborative efforts. VM edited the draft and provided reference strains. MK edited the draft and provided collaborative efforts. All authors read and approved the final manuscript.

Funding

This study was supported by the International Science and Technology Center (ISTC). The field epizootical surveys and other analyses were funded through regular projects #KR-1784 and #KR-2111 with financial support from the Government of Canada, EU, and USA. The authors are very grateful for the assistance given by both ISTC and international collaborators of the projects.

Acknowledgments

We thank Kenneth Gage and Irina Goodrich for the careful reading of the manuscript and constructive comments.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  • 1.

    AubakirovSASerzhanovOSFomushkinVMKhotkoIIRudenchikUVPoleSBet al. Guidelines for the Landscape-Epizootic Zoning of Natural Foci of Plague in Central Asia and Kazakhstan.Kazakh Scientific Center of Qurantine and Zoonotic Diseases, Alma-Ata (1990).

  • 2.

    AtshabarBBBurdelovLAAgeyevVSAubakirovSADubyanskyVMGrazhdanovAKet al. Atlas of the Spread of Bacterial and Viral Zoonotic Infections in Kazakhstan. Almaty: Reforma (2010).

  • 3.

    Ibragimov AbdikarimovSTUsenbaevNTAltybaevKI. Some of the results of non-specific prevention of plague foci in the highland of the Kyrgyz Republic. Health Care Kyrgyzstan. (2017) 4:303. Available online at: https://elibrary.ru/item.asp?id=32741738

  • 4.

    AbdikarimovSSariyevaGSagiyevZBazarkanovaGShabuninAMaimulovRet al. The Database “Epidemiology and Epizootology of the Plague in the Sari-Dzhas Natural Focus of Kyrgyzstan.”Intellectual property—Official Bulletin of the State Service of Intellectual Property and Innovations under the Government of the Kyrgyz Republic (2017) 801.

  • 5.

    PopovNVKedrovaOVToporkovVP. Information Letter “On the Situation of the Plague in the Kyrgyz Republic in August 2013”. Available online at: http://www.microbe.ru/files/infletter_2013_1.pdf

  • 6.

    SariyevaGAbdelZShabuninASagiyevZAbdikarimovSBazarkanovaGet al. Current status of the Sari-Dzhas natural focus of plague, Kyrgyzstan: epizootic activity and marmot population. Vector Borne Zoonotic Dis. (2018) 18:52432. 10.1089/vbz.2017.2200

  • 7.

    SariyevaGEBazarkanovaGDSagiyevZAShabuninAGAbdikarimovSTMaimulovRKet al. Atlas of the Sari-Dzhas natural focus of plague in Kyrgyzstan: Spatial and temporal characteristics.Bishkek: Global Print (2018), p. 172.

  • 8.

    PoleSB. Marmots and zoonotic infections in CIS. In: RamousseRAllainéDLe BerreM editors. Adaptive Strategies and Diversity in Marmots/Stratégies Adaptatives et Diversité chez les Marmottes 2003.Lyon: International Network on Marmots (2003). p. 138.

  • 9.

    ArmitageKB. Climate change and the conservation of marmots. Nat Sci. (2013) 5:3643. 10.4236/ns.2013.55A005

  • 10.

    DaiRWeiBXiongHYangXPengYHeJet al. Human plague associated with Tibetan sheep originates in marmots. PLoS Negl Trop Dis. (2018) 12:e0006635. 10.1371/journal.pntd.0006635

  • 11.

    RicankovaVRiegertJSemancíkovaEHaisMCejkovaAPrachK. Habitat preferences in gray marmots (Marmota baibacina). Acta Theriol. (2014) 59:31724. 10.1007/s13364-013-0161-x

  • 12.

    ReijniersJBegonMAgeyevVSLeirsH. Plague epizootic cycles in Central Asia. Biol Lett. (2014) 10:20140302. 10.1098/rsbl.2014.0302

  • 13.

    EisenRJBorchertJNMpangaJTAtikuLAMacMillanKBoeglerKAet al. Flea diversity as an element for persistence of plague bacteria in an East African plague focus. PLoS ONE. (2012) 7:e35598. 10.1371/journal.pone.0035598

  • 14.

    SuleymenovBM. Enzootic and Epizootics of Plague.Almaty: Evero (2015). p. 514.

  • 15.

    VashchenokVS. The role of fleas (Siphonaptera) in the epizootology of plague. Parazitology. (1999) 33:198209.

  • 16.

    SuntsovVVSuntsovaNI. Plague. Origin and Evolution of the Epizootic System (Ecological, Geographical and Social Aspects).Moscow: KMK Publishing House (2006). p. 247.

  • 17.

    CollingeSKJohnsonWCRayCMatchettRGrenstenJCullyJFet al. Landscape structure and plague occurrence in black-tailed prairie dogs on grasslands of the Western USA. Landscape Ecol. (2005) 20:94155. 10.1007/s10980-005-4617-5

  • 18.

    GageKLKosoyMY. Natural history of plague: perspectives from more than a century of research. Annu Rev Entomol. (2005) 50:50528. 10.1146/annurev.ento.50.071803.130337

  • 19.

    BigginsDEKosoyM. Influences of introduced plague on North American mammals: implications from ecology of plague in Asia. J Mamm. (2001) 82:90616. 10.1644/1545-1542(2001)082<0906:IOIPON>2.0.CO;2

  • 20.

    LowellJLZhansarinaAYockeyBMeka-MechenkoTStybayevaGAtshabarBet al. Phenotypic and molecular characterizations of Yersinia pestis isolates from Kazakhstan and adjacent regions. Microbiology. (2007) 153:16977. 10.1099/mic.0.29059-0

  • 21.

    ShabuninAGSarievaGEAbdikarimovSTMaimulovRKBazarkanovaGDDzhaparovaAKet al. Assessment of the vulnerability of the population on Sari-Dzhas autonomous focus of plague, Kyrgyzstan. Acta Biomed Scientifica. (2017) 2:10714. 10.12737/article_59fad5258950d0.72474236

  • 22.

    WeaverMLSariyevaGEKendirbaevDU. Biosafety in the field: an account of a Kyrgyz anti-plague expedition. Appl Biosafety. (2014) 19:7486. 10.1177/153567601401900203

  • 23.

    OrderNo. 207 of 04/27/2015. “On the Improvement of Epidemiological and Epidemiological Surveillance of the Plague in the Kyrgyz Republic”. p. 8792.

  • 24.

    AytkuluyevNS. Instructions on the conducting of anti-epidemic and preventive measures for plague, approved by the order of the ministry of health care of Kyrgyzstan 102 from 24.02.2010. (2010) 1213.

  • 25.

    EzhlovaEBPakskinaNDPopovVPKnyazevaTBKuznetzovAAPorshakovAMet al. Collection, Recording and Preparation for Laboratory Investigation of Blood-Sucking Arthropods in Natural Foci of Dangerous Infectious Diseases: Methodological guidelines. MU 3.1.3012-12. Federal Center for Hygiene and Epidemiology of Russian Surveillance on Consumer Rights Protection and Human Wellbeing (2011) 55:14.

  • 26.

    IoffIGMichelinMASkalonOI. Determinant of Fleas of Central Asia and Kazakhstan.Moscow: Medicine (1965). p. 1121.

  • 27.

    SuleymenovBMAbdirasilovaAA. The Basis of Individual Biosafety Under Work With Plague Pathogen.Almaty: Evero (2014). p. 200.

  • 28.

    LeFlèche PHauckYOntenienteLPrieurADenoeudFRamisseVet al. A tandem repeats database for bacterial genomes: application to the genotyping of Yersinia pestis and Bacillus anthracis. BMC Microbiol. (2001) 1:211. 10.1186/1471-2180-1-2

  • 29.

    MorelliGSongYMazzoniCJEppingerMRoumagnacPWagnerDMet al. Phylogenetic diversity and historical patterns of pandemic spread of Yersinia pestis. Nat Genet. (2010) 42:114043. 10.1038/ng.705

  • 30.

    BirdsellDNPearsonTPriceEPHornstraHMNeraRDStoneNet al. Melt analysis of mismatch amplification mutation assays (Melt-MAMA): a functional study of a cost-effective SNP genotyping assay in bacterial models. PLoS ONE. (2012) 7:e32866. 10.1371/journal.pone.0032866

  • 31.

    EroshenkoGANosovNYKrasnovYMOglodinYGKuklevaLMGusevaNPet al. Yersinia pestis strains of ancient phylogenetic branch 0.ANT are widely spread in the high mountain plague foci of Kyrgyzstan. PLoS ONE. (2017) 12:e0187230. 10.1371/journal.pone.0187230

  • 32.

    BurdelovLA. Modern Ideas About State of the Art of Natural Plague Foci: An Overview. Problems of Quarantine and Zoonotic Infections in Kazakhstan, Alma-Ata: Kazakh Scientific Center of Quarantine and Zoonotic Diseases. (1999). p. 1823.

  • 33.

    BibikovDI. Mountain Marmots of Central Asia and Kazakhstan.Moscow: Moscow University Press (1967). p. 202.

  • 34.

    PopovNVBezsmertnyVEMatrosovANKuznetsovAAKnyzevaTVPopovVPet al. Epizootic activity of natural plague foci in the territory of the Russian Federation in 2011, and prognosis for 2012. Prob Danger Infect. (2012) 111:510. 10.21055/0370-1069-2012-1(111)-5-10

  • 35.

    VarshavskiyBS. Spatial and biocenotic structure of the plague focus on the northern edge of the area of the great gerbil (Northern Aral Sea). (Dissertation thesis) : Astrakhan anti-plague station, Astrakhan, Russia (1985).

  • 36.

    MamatkanovDMBazhanovaLBKuzmichenokVARomanovskyVV. Influence of Climate Changes on the Mountain Ecosystem of the Tien-Shan (on the Example of the Issyk-Kul-Chui basin).Bishkek: Nur-As Press (2014). p. 273.

  • 37.

    KurmanovBKAtshabarBBSariyevaGEAbdirassilovaAAAbdelZBakhtybekkyzySHet al. Genotyping of Yersinia pestis strains from the Central Asian deserts and Tien-Shan high mountain plague foci. Medicine. (2016) 12:807.

Summary

Keywords

grey marmot, ectoparasites, plague, rodent, Kyrgyzstan, Yersinia pestis

Citation

Sariyeva G, Bazarkanova G, Maimulov R, Abdikarimov S, Kurmanov B, Abdirassilova A, Shabunin A, Sagiyev Z, Dzhaparova A, Abdel Z, Mussagaliyeva R, Morand S, Motin V and Kosoy M (2019) Marmots and Yersinia pestis Strains in Two Plague Endemic Areas of Tien Shan Mountains. Front. Vet. Sci. 6:207. doi: 10.3389/fvets.2019.00207

Received

19 December 2018

Accepted

11 June 2019

Published

04 July 2019

Volume

6 - 2019

Edited by

David Modrý, University of Veterinary and Pharmaceutical Sciences Brno, Czechia

Reviewed by

Jan Votypka, Charles University, Czechia; Ettore Napoli, University of Messina, Italy

Updates

Copyright

*Correspondence: Gulmira Sariyeva

This article was submitted to Parasitology, a section of the journal Frontiers in Veterinary Science

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics