Background {#Sec1}
==========

*Shigella* is one of the most important etiological agents of diarrhea, in particular, dysentery commonly known as bacillary dysentery. It is a more severe form than gastroenteritis and is responsible for childhood morbidity and mortality, especially in developing countries \[[@CR1]--[@CR4]\]. Ingestion of even 100 microorganisms leads after 4--7 days to an acute diarrhea. Because of delay in humoral responses, complication and mortality rate due to shigellosis in children is higher than in other age groups \[[@CR5]\]. A literature review concluded that of the estimated 165 million cases of *Shigella* diarrhea that occur annually in the world, 99 % occur in the developing world and the remaining 1 % occurs in industrialized countries. In developing countries, 69 % of these episodes occur in children under 5 years of age. Moreover, of the 1.1 million deaths attributed to *Shigella* infections in developing countries, 60 % occur in the under 5 age group in children \[[@CR6]\].

The genus *Shigella* comprise of four different species which are referred to by a letter designation based on their serological antigen: serotype A-*Shigella dysenteriae*, Serotype B-*Shigella flexneri*, Serotype C-*Shigella boydii* and Serotype D-*Shigella sonnei*. The first three are common in developing countries while *S. sonnei* is common in developed countries and Iran \[[@CR7]\].

Shifts in the prevalent serogroups have been observed in many parts of the world and so is true for Nepal. Studies conducted during1999--2002 shows high incidence of *S. dysentriae* infection in Nepal \[[@CR8]\]. Studies conducted after 2002 showed prevalence of *S. flexneri* strains among *Shigella* isolates from all parts of Nepal \[[@CR9]--[@CR12]\].

According to WHO report, antimicrobial resistance pattern for *Shigella* varies in different parts of the world and with the time. Over the past decades, *Shigella* species have become progressively resistant to most widely used antimicrobials \[[@CR13], [@CR14]\]. Since the late 1990s, the fluoroquinolones, including ciprofloxacin, norfloxacin, and ofloxacin, have been the drug of choice for multidrug-resistant *Shigella* infections. However, strains resistant to the first-line drugs commonly used in the treatment of shigellosis like ampicillin, cotrimoxazole, and nalidixic and even ciprofloxacin and norfloxacin has also been reported from various studies in Nepal and worldwide \[[@CR9]--[@CR12], [@CR15]\]. Ceftriaxone has been used as a reserved antibiotic for treatment of multidrug-resistant (MDR) *Shigella* infection \[[@CR16]\]. However, some studies have also reported emerging ceftriaxone resistance in *Shigella*. This increasing level of antimicrobial resistance in *Shigella* have limited the therapeutic options and complicated the treatment of shigellosis.

The purpose of the present study is to determine the trend of change in *Shigella* species and their antimicrobial resistance patterns over a decade for the better management of shigellosis.

Methods {#Sec2}
=======

Surveillance sites {#Sec3}
------------------

Nepal is divided into five development regions. The study was carried out in 12 participating laboratories that are under Antimicrobial Resistance Surveillance network. The study was conducted in three development regions: central, eastern, and western development regions and covered six districts. In the central development region, there were five surveillance sites in Kathmandu district (National Public Health Laboratory (NPHL), Bir Hospital, Sukraraj Tropical Hospital, Kanti Children Hospital, Teaching Hospital) and two in Lalitpur district (Patan Hospital and KIST hospital). In the western development region, Western Regional Hospital and Manipal Teaching Hospital fall in Kaski district whereas Lumbini Zonal Hospital falls in Rupandehi district and United Mission Hospital falls in Palpa district. Only one site BPKIHS falls in Morang district of the eastern development region (Fig. [1](#Fig1){ref-type="fig"}).Fig. 1Map of Nepal showing location of participating laboratories in various development regions

Selection of surveillance sites {#Sec4}
-------------------------------

Surveillance sites were selected preferring large private/government hospital laboratories from various development regions where the diagnostic facilities were adequate and the probability of patient inflow during outbreaks was higher.

Sample population and data collection {#Sec5}
-------------------------------------

Diarrheal samples were collected from individual patients referring to the private or governmental laboratories. A questionnaire for each subject was filled at respective laboratories mentioning basic details of the patient, age (months), sex, location of residence, and travel. A history from all patients was also recorded, including duration of illness prior to admission to hospital (days), fever (defined as a prolonged temperature \>37.5 °C), abdominal discomfort, vomiting, watery diarrhea (defined as three or more loose bowel movements during a 24-h period), bloody or mucoid diarrhea (defined as \>3 loose stools with obvious blood or mucus), estimated number of episodes of diarrhea before attending hospital, and if there was any known pretreatment with antimicrobials.

Sources of *Shigellae* {#Sec6}
----------------------

Isolates recovered from all 12 sites of Antimicrobial Resistance surveillance network in Nepal during 2003--2015 were reported to NPHL and were included in the study.

Laboratory methods {#Sec7}
------------------

Stool specimens received at respective laboratories were cultured for enteric pathogens, including *Shigella* spp. by standard methods \[[@CR17]\]. Sample was streaked onto MacConkey agar and a selective medium (either deoxycholate citrate agar or xylose lysine deoxycholate agar or *Salmonella*-*Shigella* agar) (all media by HiMedia, India) by the quadrant isolation technique \[[@CR18]\]. The plates were incubated at 37 °C for 24 h. Isolates were identified by Gram's staining, colony characteristics, and biochemical tests.

After identification of isolates, they were subcultured onto nutrient agar slants in screw capped tube, incubated overnight and then transported to NPHL in a cold box. All reported isolates of *Shigella* were re-confirmed by serotyping using commercially available antisera from Denka Seiken, Japan, and were preserved in tryptic soy broth with 20 % glycerol at −75 °C for further use. Antibiotic susceptibility test was performed for all the identified *Shigella* isolates by Standard Kirby Bauer's disc diffusion technique. The antibiotics used for analysis were ampicillin (Amp-10 mcg), cotrimoxazole (Sxt-25 mcg), nalidixic acid (NA-30 mcg), ciprofloxacin (Cip-5 mcg), mecillinam (Mec), and ceftriaxone (CRO-5 mcg).

Statistical methods {#Sec8}
-------------------

Data obtained were analyzed using SPSS software for windows version 18. Data on variables like shigella serotypes, age group, and sex were calculated as percentages and compared using chi-square test. *P* \< 0.05 was considered to be statistically significant. Duplicate samples were counted as one during data entry.

Results {#Sec9}
=======

A total of 332 *Shigella* isolates were reported in a period of 2003--2015 from 12 participating laboratories of Antimicrobial Resistance surveillance network in Nepal. The highest number was reported in 2005. Peak isolation was seen in post monsoon season, i.e., September (Fig. [2](#Fig2){ref-type="fig"}).Fig. 2Region wise reported cases of total *Shigella* isolates (*n* = 232)

*Shigella dysenteriae* was the prevalent species in Nepal in the years 2003--2004 but since 2005, *S. flexneri* replaced it and became the most prevalent species. A marked increase in the number of *S. sonnei* (*N* = 6) was noticed in the year 2005 which is highest reported till date. In contrast, no isolate of *S. boydii* was reported since 2009 till 2014 where one isolate was reported. Since 2008, few cases of unidentified *Shigella* isolates were also reported (Table [1](#Tab1){ref-type="table"}).Table 1Annual distribution of different serotypes of *Shigella* (2003--2015)Year*S. dysentriae N* (%)*S. flexneri N* (%)*S. sonnei N* (%)*S. boydii N* (%)Unidentified *Shigella* spp.Total200313 (52 %)7 (28 %)2 (8 %)3 (12 %)--25200428 (63.7 %)12 (27.3 %)2 (4.5 %)2 (4.5 %)--44200510 (19.6 %)31 (60.8 %)6 (11.8 %)4 (7.8 %)--5120066 (26.1 %)14 (60.9 %)2 (8.7 %)1 (4.3 %)--23200710 (27 %)22 (59.5 %)3 (8.1 %)2 (5.4 %)--3720082 (11.8 %)8 (47.2 %)1 (5.8 %)1 (5.8 %)5 (29.4 %)1720098 (40 %)5 (25 %)1 (5 %)1 (5 %)5 (25 %)2020105 (50 %)4 (40 %)----1 (10 %)1020112 (18.2 %)5 (45.4 %)2 (18.2 %)--2 (18.2 %)1120122 (18.2 %)6 (54.5 %)3 (27.3 %)----1120135 (16.1 %)21 (67.7 %)3 (9.7 %)--2 (6.5 %)3120141 (4 %)12 (48 %)5 (20 %)1 (4 %)6 (24 %)2520153 (11.1 %)19 (70.4 %)4 (14.8 %)--1 (3.7 %)27Total95 (28.6 %)166 (50 %)34 (10.2 %)15 (4.5 %)22 (6.6 %)332

Region wise distribution of reported isolates showed that 18 isolates of *Shigella* were reported from one laboratory in the eastern development region, 188 isolates were reported from seven participating laboratories in the central development region, and 126 were reported from four participating laboratories in the western development region (Fig. [3](#Fig3){ref-type="fig"}).Fig. 3Seasonal distribution of total reported cases (*n* = 232)

Of the total, 24.5 % of cases were demonstrated in children below 10 years followed by 21--30 years age group (22.9 %) and 11--20 years age group (18.5 %). Significant difference was not noted (*p* = 0.059) in age wise distribution of cases over the years. *Shigella* infection was reported higher among male in age groups below 30 years and 71--80 years but higher cases were reported in females of age groups 31--70 years and 80+ years (Fig. [4](#Fig4){ref-type="fig"}).Fig. 4Percentage of laboratory-confirmed *Shigella* infection reported to NPHL, by age group and sex (*n* = 196 with age and sex information reported)

The antimicrobial susceptibility profile of *Shigella* isolates in the years 2003--2015 is shown in Table [2](#Tab2){ref-type="table"}. Different serotypes of *Shigella* exhibited various resistance patterns. *S. dysentriae* showed varying resistance to ampicillin (80--100 %), co-trimoxazole (20--100 %), and nalidixic acid (20--100 %). Likewise, *S. flexneri* also exhibited varying resistance to ampicillin (80--100 %), co-trimoxazole (58--100 %), and nalidixic acid (0--100 %). Until 2003, only *S. dysentriae* strains exhibited resistance to ciprofloxacin; however, 50 % of the *S. flexneri* isolates were resistant to ciprofloxacin by 2004. By 2012, *S. dysentriae* strains were found susceptible to ciprofloxacin in contrast to *S. flexneri*, which exhibited highest resistance. Among *S. sonnei* and *S. boydii* serotypes, resistance to ampicillin, cotrimoxazole, and nalidixic acid kept varying. Only three *S. sonnei* isolates were found resistant to ciprofloxacin while no ciprofloxacin-resistant *S. boydii* was reported. All the serotypes of *Shigella* were 100 % resistant to mecillinam since 2007; so by 2010, mecillinam was no longer included in antimicrobial testing of *Shigella.* Among the serotypes, resistance to ceftriaxone was first reported in *S. flexneri* (10--33 %) followed by non-typable *Shigella* spp. with around 40 % resistance; 22.2% of the total isolates tested were resistant to all three antimicrobials. Figure [5](#Fig5){ref-type="fig"} demonstrates the cumulative antibiogram of various species of *Shigella*.Table 2Antimicrobial susceptibility profile of various species of *Shigella* (2003--2015)*N* (%) resistance toYearSpeciesNumber testedAmpicillinCiprofloxacinNalidixic acidCotrimoxazoleMecillinamCeftriaxone2003*S. dysentriae*1311 (85)8 (62)10 (77)9 (69)5 (38)0*S. flexneri*77 (100)006 (86)2/6 (33)^a^0*S. boydii*32 (66)01 (33)2 (66)2 (66)0*S. sonnei*21 (50)02 (100)2 (100)1 (50)02004*S. dysentriae*2825 (89)13 (46)14 (50)18 (64)20/27 (74)^a^0*S. flexneri*129 (75)6 (50)6 (50)8 (67)4 (33)0*S. boydii*21 (50)001 (50)00*S. sonnei*2002 (100)2 (100)002005*S.dysentriae*108 (80)6/9 (66)^a^7 (70)8 (80)2/8 (25)^a^0*S. flexneri*3118/30 (60)^a^017/28 (59)^a^21 (68)13 (42)0*S. boydii*64 (66)03 (50)04 (66)0*S. sonnei*41 (25)02 (50)2 (50)002006*S. dysentriae*64/5 (80)^a^04/5 (80)^a^2/5 (40)^a^4 (66)0*S. flexneri*148/13 (62)^a^3/12 (25)^a^4/10 (40)^a^7/12 (58)^a^14 (100)0*S. boydii*2000000*S. sonnei*1SRRSSS2007*S. dysentriae*105/6 (83)^a^4 (40)2 (20)2 (20)10 (100)0*S. flexneri*2217/21 (81)^a^12/21 (57)^a^21 (95)19/21 (90)^a^22 (100)2/20 (10)^a^*S. boydii*33 (100)0003 (100)0*S. sonnei*200002 (100)02008*S. dysentriae*202 (100)2 (100)02 (100)0*S. flexneri*88 (100)3/5 (60)^a^8 (100)4/5 (80)^a^8 (100)0*S. sonnei*22 (100)02 (100)02 (100)0*Shigella* spp.53 (60)1/4 (25)^a^2/4 (50)^a^3/4 (75)^a^5 (100)1/3 (33)^a^2009*S. dysentriae*88 (100)2 (25)7 (88)7 (88)8 (100)0*S. flexneri*64/5 (80)^a^6 (100)4/5 (80)^a^6 (100)6 (100)2 (33)*S. sonnei*1RSRSRS*Shigella* spp.53 (60)3 (60)4 (80)4 (80)5 (100)2 (40)2010*S. dysentriae*55 (100)1 (20)5 (100)5 (100)--0*S. flexneri*44 (100)1 (25)4 (100)4 (100)--0*Shigella* spp.1RSRR--S2011*S. dysentriae*22 (100)1 (50)1 (50)2 (100)--0*S. flexneri*53 (60)003/4 (75)^a^--0*S. sonnei*21 (50)02 (100)0--0*Shigella* spp.22 (100)02 (100)2 (100)--02012*S. dysentriae*44 (100)03 (75)3 (75)--0*S. flexneri*77 (100)4/6 (66)^a^3/4 (75)^a^4/6 (67)^a^--2/6 (33)^a^*S. sonnei*3001 (33)3 (100)--02013*S. dysentriae*55 (100)04 (80)3 (60)--0*S. flexneri*1915/16 (93)^a^12/17 (71)^a^16/18 (89)^a^14 (74)--5/18 (27)^a^*S. sonnei*33 (100)2 (66)2 (66)2 (66)--0*Shigella* spp.33 (100)03 (100)3 (100)--02014*S. dysentriae*1SRSS--S*S. flexneri*127/11 (64)^a^8 (66)11 (92)7 (58)--0*S. sonnei*53 (60)5 (100)5 (100)5 (100)--0*S. boydii*1SRRR--S*Shigella* spp.66 (100)1/5 (20)^a^3/4 (75)^a^3 (50)--02015*S. dysentriae*3003 (100)2 (75)--0*S. flexneri*1916/18 (88)^a^9 (47)15 (79)17 (89)--3 (16)*S. sonnei*43 (75)4 (100)4 (100)4 (100)--0*Shigella* spp.1RSRR--S^a^The number in which antibiotics were incorporated, if less than the total number of isolates Fig. 5Year wise distribution of MDR isolates (in number and percentage)

The prevalence of multidrug-resistant (resistance to two or more classes of antibiotics) *Shigella* isolates is demonstrated in Fig. [5](#Fig5){ref-type="fig"}. The prevalence of MDR increased significantly in 2010 as compared to 2003. On the basis of their pattern of antimicrobial resistance to the six antibiotics used in the study, the reported *Shigellae* were grouped into seven phenotypes (Table [3](#Tab3){ref-type="table"}). Among the *Shigella* serotypes, high prevalence of MDR was observed in *S. dysentriae* isolates occupying 52 % of the total isolates.Table 3Multidrug-resistant phenotypes of *Shigella*MDR phenotypeNumber (%) of isolates exhibiting resistance*S. dysentriaeS. flexneriS. boydiiS. sonneiShigella* spp.Amp/Cip or NA/Mec1 (1.9 %)3 (5.2 %)2 (50 %)----Amp/Mec/Ts5 (9.4 %)3 (5.2 %)------Mec/Cip or NA/Ts2 (3.8 %)3 (5.2 %)--1 (10 %)--Amp/NA or Cip/Ts25 (47.2 %)36 (62.1 %)1 (25 %)9 (90 %)3 (50 %)Amp/NA or Cip/Ts/Mec20 (37.7 %)5 (8.6 %)----3 (50 %)Amp/NA or Cip/Ts/Mec/CRO--1 (1.7 %)1 (25 %)----Amp/NA or Cip/Ts/CRO--7 (11.9 %)------Total number MDR (%)53 (55.8 %)58 (34.9 %)4 (26.6 %)10 (29.4 %)6 (27.3 %)

Discussion {#Sec10}
==========

Shigellosis still accounts for a significant proportion of mortality and morbidity, especially in developing countries \[[@CR19]\]. The present study, covering 13 years surveillance (2003--2015), demonstrates the changing serotype and trend of antimicrobial resistance of *Shigella* isolates in Nepal.

An analysis of serotypes revealed that *S. dysenteriae* was the prevalent species in Nepal in the years 1999--2004 but a change in serotype in *Shigella* species was noticed during 2005. The changing patterns in the distribution of *Shigella* serogroups and serotypes have been reported from time to time \[[@CR20]--[@CR23]\]. In our study, *S. flexneri* remained the most prevalent serotype in Nepal since 2005 which is in accordance with many recent studies in Nepal and other developing countries \[[@CR9], [@CR15], [@CR17], [@CR24]--[@CR27]\] but dissimilar to some other studies \[[@CR28]--[@CR30]\]. It has been hypothesized that these changes may be due to varying behavourial and environmental factors such as reduction of overall malnutrition, use of safe drinking water instead from tubewells or pipes, and improved sanitary facilities for defecation \[[@CR31]\].

In our study, like in other studies, the seasonal tendency of shigellosis was summer-monsoon \[[@CR9], [@CR15]\]. This may be because during this period, there is a climatic variation due to which the sources of drinking water gets rapidly contaminated. Males were more commonly infected as compared to females as in other studies made by Khan et al. and Taneja \[[@CR15], [@CR32]\]. This male preponderance may be due to greater chances of exposure to male population who go out of the house much frequently than females for eating food from restaurants and street vendors. The susceptible age group in our study, like in other studies, was below 10 years \[[@CR12], [@CR15], [@CR25]\]. In our study, very few cases were reported from children below 5 years. This may be because the laboratories in the surveillance network are limited; many cases may be referred to private institutions which are not included in the surveillance network. Besides that, other factors like incomplete information on age (16 % missing) available from participating laboratories, inclusion of only one children's hospital in the surveillance and low *Shigella* isolation rate may have contributed to under reporting of true number of victims.

Treatment of shigellosis depends on appropriate antimicrobial therapy. A variety of antibiotics are effective for treatment of shigellosis, although options are becoming limited due to globally emerging drug resistance. The present study also shows high resistance among all *Shigella* species to the first line low cost drugs commonly used in the treatment of shigellosis like Ampicillin, Cotrimoxazole and Nalidixic acid. Considering resistance pattern in different groups, *S. flexneri* and *S. dysentriae*, the two most common serotypes revealed similar pattern in resistance with slightly higher range in *S. flexneri*. However, the difference was not statistically significant. Similar studies conducted by Khan et al. reveals that *S. flexneri* exhibited increased resistance to ampicillin and nalidixic acid as compared to *S. dysentriae* whereas high resistance was observed against ciprofloxacin and cotrimoxazole in *S. flexneri* as compared to *S. dysentriae* \[[@CR15]\]*.* This finding however does not agree with some previous studies done by Wilson et al. and Kansakar et al. \[[@CR9], [@CR12]\] in Nepal which demonstrates high resistance in *S. dysentriae*. Among remaining serotypes, untypable *Shigella* spp. exhibited high resistance to all antimicrobials except ciprofloxacin as compared to *S. sonnei* and *S. boydii. S. boydii* exhibited least resistance to all antimicrobials tested. This finding is also in accordance to studies made by Khan et al. \[[@CR15]\] but contradictory to findings by Kansakar et al. \[[@CR9]\] which reveals *S. sonnei* as the most susceptible serotype. Emerging resistance to ceftriaxone in *S. dysentriae*, *S. flexneri*, and *Shigella* spp. is noted which raises concern for the use of third-generation cephalosporins as last therapy of treatment.

This difference in antimicrobial resistance among various serotypes may be attributed to change in time period of the studies. Studies conducted during 2004/2005 showed contradictory result in terms of prevalent serotype and resistance profile among serotypes as compared to that of 2013/2014 which seems obvious.

MDR (resistance to three or more than three antibiotics) was found in 46.1 % of total isolates, and the most common MDR profile was ampicillin, nalidixic acid, and cotrimoxazole. This finding is very low as compared to 90.2 % reported by Reema et al. However, the MDR profile is similar \[[@CR33]\].

Limitations {#Sec11}
-----------

Data represented from our 12 surveillance sites may not be representative of the sick population at large in the community. Lack of total number of reported cases and inadequate data management systems also limited the results of our study.

Conclusions {#Sec12}
===========

With increasing resistance to first-line agents along with ciprofloxacin, only ceftriaxone remains the drug of choice; however, increasing MDR along with resistance to third-generation cephalosporins as last resort creates an alarming situation warranting proper usage of antimicrobial agents and its continuous monitoring.

Amp

:   Ampicillin

BPKIHS

:   B. P. Koirala Institute of Health Sciences

Cip

:   Ciprofloxacin

CRO

:   Ceftriaxone

KIST

:   KIST Medical College and Teaching Hospital

MDR

:   Multidrug resistant

Mec

:   Mecillinam

NA

:   Nalidixic acid

NPHL

:   National Public Health Laboratory

Ts

:   Cotrimoxazole

WHO

:   World Health Organization

We acknowledge all the laboratory staffs of the participating laboratories for their contribution in laboratory diagnosis of *Shigella* in their respective units. We are thankful for their invaluable support.

Funding {#FPar1}
=======

Not applicable.

Availability of data and materials {#FPar2}
==================================

Please contact author for data requests.

Authors' contributions {#FPar3}
======================

GS critically reviewed the manuscript and approved the final version of the manuscript. JA conceived and designed the research. SA analyzed and interpreted the data. NR finalized and drafted the manuscript. All authors have read and approved the final manuscript.

Competing interests {#FPar4}
===================

The authors declare that they have no competing interests.

Ethics approval and consent to participate {#FPar5}
==========================================

At the time of enrollment, verbal consent was taken from the caregivers or guardians on behalf of the patients. Parents or guardians were assured about the non-disclosure of information collected from them and were also informed about the use of data for analysis and using the results for improving patient care activities as well as publication without disclosing the name or identity of their children. The consent to utilize the data was obtained from respective participating laboratories who reported the cases, and the laboratories have no objection in its publication. The ethical review board of NPHL has approved the publication of this data.
