key: cord-0020530-u36iuuxh authors: Sanikhani, Rahimeh; Moeinirad, Mohammad; Solgi, Hamid; Hadadi, Azar; Shahcheraghi, Fereshteh; Badmasti, Farzad title: The face of hypervirulent Klebsiella pneumoniae isolated from clinical samples of two Iranian teaching hospitals date: 2021-08-31 journal: Ann Clin Microbiol Antimicrob DOI: 10.1186/s12941-021-00467-2 sha: 4b17ac70ffbd829a8fca161e3363d7771f7a2175 doc_id: 20530 cord_uid: u36iuuxh Hypervirulent Klebsiella pneumoniae (hvKp) has emerged as a pathogen of global concern. In this study, both phenotypic and genotypic tests were used to detect hvKp. Antimicrobial resistance profiles and clonal relatedness of clinical isolates were also determined. We found that 34.2% (163/477) of the isolates were tellurite resistant, and among them 102 hvKp isolates detected with iucA or iutA or peg-344 as molecular markers. The bla(SHV) (80.4%), followed by bla(CTX-M-15) (76.5%) and bla(TEM) (67.6%), bla(OXA-48) (53.9%), and bla(NDM-1) (32.3%) were detected, while bla(KPC-1) was not present in any hvKp isolates. It was found that the majority of hvKp isolates belonged to capsular serotype K20 and ompK36 group C, which is related to clonal group (CG) 23 (e.g. ST23). A high percentage of multidrug-resistant hvKp (76.6%) and high resistance to imipenem (67%) indicated a serious problem that should be addressed in the clinical setting. Hypervirulent Klebsiella pneumoniae (hvKp), an emerging pathotype of K. pneumoniae was first reported from Taiwan. It was identified as an important cause of pyogenic liver abscess [1, 2] . In hvKp isolates, pLVPK-like plasmids (Large Virulence Plasmid of K. pneumoniae) encoding virulence factor genes including capsular polysaccharide synthesis regulators (rmpA and rmpA2) and iron acquisition systems (iucA, iutA, and iro siderophore gene cluster), a metabolic transporter (peg-344) and also heavy metal resistance genes (copper, silver, lead, and tellurite), have been identified [3, 4] . Therefore, most hvKp isolates are able to reduce tellurite and form a black colony due to the presence of the major virulence plasmids containing a tellurite resistance gene [5] . The pLVPK-like plasmids may carry all virulence factor genes or have lost some of them [6, 7] . On the other hand, the acquisition of antibiotic resistance plasmids or insertion of resistant mobile genetic elements into the hvKp plasmid turns them into superbugs that can be termed hyper-resistant hvKp strains [8] [9] [10] . Some K. pneumoniae clones are characterized as high-risk clones that play an important role in the spread of antibiotic-resistant strains [11, 12] . The association of the porin ompK36 with clonal relatedness of K. pneumoniae isolates has been described in several studies [13, 14] . This typing method can be considered as a rapid method for characterizing the clonal relatedness of K. pneumoniae isolates. Four different genotypes for ompK36 porin (A to D) in K. pneumoniae were defined and the correlation of different variants of ompK36 with specific sequence types (STs) was illustrated [13] [14] [15] . It is important to distinguish hvKps from classical K. pneumoniae (cKp) isolates. To date, several methods have been used to identify hvKp isolates. Detection of hypermucoviscous phenotype on agar (string test), use of Galleria mellonella infection model, serum killing assay and mouse infection models are some of the methods used for phenotypic identification of these strains [16, 17] . However, there is no consensus phenotypic test for the early diagnosis of hvKps, and it appears that more accurate phenotypic tests are needed to rapidly identify these pathotypes of K. pneumoniae. In addition to phenotypic methods, the presence of different virulence genes has been studied to increase the sensitivity and accuracy of hvKps identification. Among these, the genes of iucA, iutA, and peg-344 have been introduced as the best genetic diagnostic markers with the highest accuracy [18] . Therefore, the aim of this study was to develop a rapid identification method for hypervirulent K. pneumoniae. We also investigated the genotypic characteristics, prevalence of virulence factors and antibiotic resistance of hvKps in clinical samples isolated from Iran. In this cross-sectional study, we collected a total of 477 non-repetitive K. pneumoniae as clinical isolates from two educational hospitals in Tehran over a period of time from June 2019 to December 2020. All bacterial isolates were identified using standard biochemical laboratory methods and then the isolates were stored in a freezer at −70 °C in nutrient broth containing 20% glycerol until further studies. We used tellurite agar culture as a rapid screening test in this study. The isolates that formed black colonies on this tellurite-containing selective medium were considered as presumptive hypervirulent strains for further study. For this purpose, 0.1 g of potassium tellurite powder was first dissolved in 10 ml of sterile distilled water and filtered using membrane filters of pore size 0.45 µm. Then we added 300 μl of the potassium tellurite solution to 100 ml of Mueller-Hinton agar medium, which was autoclaved and cooled to 45-50 °C. Finally, we dispensed into sterile plates. Colonies were examined after overnight incubation at 37 °C. Hypermucoviscous phenotype of the hvKp isolates was examined by the string test, and the positive result was confirmed via the formation of a 5-mm viscous filament by stretching of bacterial colonies on a blood agar after 24 h of incubation at 37 °C [19] . Plasmid DNA extraction Mini Kit (FAVORGEN Biotech Corporation, Taiwan) has been used for the detection of genes carried on plasmids. In addition, the boiling method was used for isolation of genomic DNA [20] . All amplification reactions for PCR assays were prepared in a total volume of 25 μl. The list of primer sequences, PCR product sizes, and PCR conditions is shown in Table 1 . Finally, all PCR amplification products were sequenced and then searched in the GenBank database using BLAST tool (http:// www. ncbi. nlm. nih. gov/ blast/). All tellurite-resistant K. pneumoniae were screened for the presence of the aerobactin (iucA), its receptor (iutA) genes and peg-344. The isolates containing the iucA or iutA or peg-344 genes were considered as hvKps [18] . Antimicrobial susceptibility testing was performed using the disc diffusion method according to the clinical and laboratory standards institute (CLSI) guidelines (CLSI 2018-M100-S28) by the following antibiotic discs including amikacin (AK), gentamicin (GN), cefotaxime (CTX), ceftazidime (CAZ), ceftriaxone (CRO), imipenem (IMI), meropenem (MRP), cefepime (FEP), ciprofloxacin (CIP), ampicillin (AMP) and aztreonam (AZM). Minimum inhibitory concentrations (MICs) of imipenem and ceftazidime were determined by broth dilution method. Escherichia coli ATCC 25922 was used as the quality control strain for antimicrobial susceptibility testing. The hvKp capsular serotypes K1, K2, K5, K20, K54, and K57 were identified using PCR method [21] . The hvKp virulence genes including salmochelin siderophore (iroB), mucoviscosity-associated gene (magA), Klebsiella ferric uptake (kfu), yersiniabactin (ybt), allantoin metabolism gene (allS), and rmpA were detected by specific primers listed in Table 1 . In addition, PCR assays were carried out for detection of bla TEM , bla SHV , and bla CTX-M-15 , bla KPC-1 , bla NDM-1 , and bla OXA-48 genes in all hvKp isolates. All hvKp isolates were subjected to ompK36 typing by the PCR-based method described by Yan et al., using four pairs of primers [13] . The statistical analyses of data were performed using SPSS software, version 16.0 (IBM, Armonk, NY, USA) and Chi-square tests (2 × 2 contingency table) were used to compare the data associated with hvKp and cKp strains. Finally, the P values < 0.05 was considered statistically significant. The accession numbers of bla OXA-48 , bla NDM-1 , iutA, iucA and peg-344 are MZ245618, MZ245619, MZ245620, MZ245621 and MZ245622, respectively in GenBank database. In this study, 163 (34.2%) out of 477K. pneumoniae isolates were able to grow on tellurite-containing MH medium and were considered tellurite-resistant strains, so they were selected for the molecular identification test. In addition, 62 out of 477K. pneumoniae isolates (13%) were reported with positive string test and hypermucoviscous phenotype. Table 2 . In this study, we investigated the antimicrobial susceptibility profile in 90 hvKp isolates. Susceptibility profiles against antimicrobials agents are shown in Table 2 . The highest rate of antibiotic resistance was related to ampicillin (100%), followed by cefotaxime and ceftazidime (91% Capsular genotyping (K genotyping) of hvKp isolates showed that capsular serotype K20 was detected in more than half of the hvKp strains (54.9%). K2 and K1 were identified in only 3 (2.9%) and one isolate (1%), respectively, while K5, K54 and K57 were not detected in any of the hvKp isolates. In addition, 42 isolates (41.2%) did not belong to serotypes K1, K2, K5, K20, K54 and K57 as shown in Table 2 . PCR for virulence-associated genes revealed that ybt (77.5%) was the most common virulence factor gene after iucA. The other virulence factor genes including rmpA, iroB, magA, kfu and allS were detected in 48%, 3.9%, 0.98, 21.6% and 1.96% hvKp isolates, respectively. See Table 2 . The distribution of ESBLs and carbapenemase genes among hvKp isolates are shown in Table 2 . The results showed that 92.2% (94/102) of hvKp isolates carried at least one antibiotic resistance gene and only 7.8% (8/102) had no resistance gene. The bla SHV was the most common Beta-lactamase gene (80.4%), followed by bla CTX-M-15 (76.5%) and bla TEM (67.6%). Also, PCR amplification of carbapenemase genes showed that bla OXA-48 (53.9%) was the dominant genotype of carbapenem-resistant strains, followed by bla NDM-1 (32.3%), while bla KPC-1 was not detected in any hvKp isolate. Thus, 56.8% (58/102) of the hvKp isolates that co-carried bla TEM , bla SHV , bla CTX-M-15 , bla OXA-48 and bla NDM-1 genes was the predominant MDR-hvKp genotype. On the other hand, as shown in Table 2 , more than 2 virulence factor genes were detected simultaneously in the majority of resistant hvKp strains and virulence profiles including: iucA, iutA, ybt and rmpA genes have been reported in 40.4% (38/94) resistant hvKp isolates. PCR-based ompK36 typing revealed that ompK36 group C was the most common type with 70.6% (72/102) frequency. The prevalence of the other types, including ompK36 groups A, B, and C, was (14/102) 13.7%, (5/102) 4.9%, and 10.8% (11/102), respectively. See Table 2 . In this study, 78 cKps and 102 hvKps isolates were examined. Demographic data and antimicrobial resistance profile were compared using chi-square tests. The data showed that there were no significant differences in demographic data between two groups. However, significant differences were found in antimicrobial resistance (e.g. amikacin, cefotaxime and gentamicin) and the presence of carbapenemases (bla OXA-48 , bla NDM-1 ). See Table 3 . In this study, to identify these K. pneumoniae superbugs, we used a combination of phenotypic and genotypic methods including tellurite resistance, and preferential gene markers. Previously, only the string test was used as a phenotypic method to identify hvKp isolates, but the string test is not a reliable rapid test for hvKp detection [1, 22, 23] . MacConkey inositol potassium tellurite agar (MCIK) has been used as a selective medium for the detection of K. pneumoniae from environmental sources or animal and human fecal samples [24] . This study showed that the trait of tellurite resistance is strongly associated with CG23, CG65 and CG86, which are mostly invasive community-acquired strains of K. pneumoniae. It appears that the large virulence plasmids of hvKp harbor tellurium resistance genes [24] . In silico analysis revealed that the tellurium gene cluster is highly prevalent among hypervirulent plasmids which is present in different sequence types (data not yet published). This prompted us to use Mueller-Hinton agar containing potassium tellurite as a selective medium for the rapid detection of hypervirulent K. pneumoniae strains. In this study, out of 163 tellurite-resistant isolates, 102 strains were genetically confirmed as hypervirulent K. pneumoniae, so making this method superior than string test for rapid phenotypic hvKp identification. We also used three key virulence genes as molecular biomarkers previously introduced by Russo et al. to increase the accuracy and sensitivity of hvKp detection [18] . In addition, all hvKp isolates were examined for the presence of other virulence factor genes. In general, the frequencies of virulence factor genes, from highest to lowest, iucA, ybt, iutA, rmpA, kfu, iroB, peg-344, allS, and magA, respectively, were reported. Other studies have also shown that the aerobactin is produced by more than 90% of hvKp, whereas only 6% of cKp strains can express it [18, 25] . In a study by XU et al. the prevalence of iucA, iutA, rmpA and iro was reported to be 56.8%, 56.8%, 43.2% and 40.9%, respectively. The prevalence of iutA and rmpA was similar to our study, but in the present study, the prevalence of iucA was higher and iro was lower than the results of the study by Xu et al. [26] . The ybt was the second most prevalent virulence factor gene among hvKp isolates in this study. The yersiniabactin gene and its receptor, which is an important virulence factor for the survival of Klebsiella strains under severe conditions, can transmit both an integrative conjugative element (ICEKp) and a plasmid (recently reported) [27] . Some studies have described the correlation between yersiniabactinproducing hvKps and pulmonary infectious diseases [28, 29] . In Iran, a study conducted by Tabrizi et al. reported that 5 of 53 K. pneumoniae strains isolated from ventilator-associated pneumonia were hvKp [30] . In the current study, of 33 hvKps isolated from lung-related samples, 27 isolates were ybt-positive, confirming the results of previous studies. The rmpA was identified as the fourth most virulence factor. Because rmpA increases the expression of capsular polysaccharide (CPS), we expected that the rmpA-producing hvKp that were isolated would be string test-positive, but this hypothesis was refuted by our results, such that only 36.7% of the rmpA-positive isolates were reported as hypercystic phenotype. Studies have shown that other genes besides rmpA are involved in capsular gene expression, such as regulation of capsular synthesis B (rcsB). Both rmpA and rcsB genes have been shown to co-occur [31] . In addition, the data show that the co-presence of four genes (iucA, ybt, iutA and rmpA) was more frequent in hvKp. In addition, other plasmid-born genes such as iro, peg-344 were less frequent and were reported only sporadically. Sequencing and analysis of large virulence plasmids from hvKp strains revealed that virulence-associated genes were mainly found in two regions. The rmpA2, iucABCD and iutA genes are located close to each other, followed by the rmpA, peg-344 and iroBCDN genes in the second region. Some virulence plasmids carry all virulence genes (e.g. pLVPK, GenBank accession number: AY378100), but others have lost one or more virulenceassociated loci, confirming our result (e.g. pVir, GenBank accession number: CP029383.2) [6, 32] . Despite most Asian countries having introduced K1 and K2 as the most common capsular serotypes [21, [33] [34] [35] , we identified K20 as the most common capsular type in Iran. This phenomenon suggests that the prevalence of the different serotypes may vary depending on the geographical area. Although there has been no comprehensive study on the hvKp isolates in Iran and little information is available on them, no K1 and K2 were found among the K. pneumoniae isolates in the study conducted by Aghamohammad et al. [36] . Also, in another study, one K1 and 15 K2 were identified among 122 K. pneumoniae isolates from Semnan, Iran, which are in agreement with our results (we detected only one K1 and three K2) [37] . Another study from Iran conducted by Solgi et al. reported that the prevalence of K1 and K2 was 45.9% and 13.5% respectively, which was more than the present study [38] . Also, in the present study, according to the Table 3 , in terms of frequency in the type of sample, and the hospital wards did not show a significant difference between these two variants. In both groups, almost half of the samples were isolated from the ICU. However, patients admitted to the ICU due to prolonged hospitalization were more susceptible to hvKp infection. In addition, it may increase the probability of horizontal gene transfer in clinical settings [39, 40] . In this study, it was also found that inhospital B, most hvKps had the same genotypic characteristics such as capsular serotype K20, ompK36 type C, and similar antibiotic resistance profiles. Therefore, the hvKp regional expansion hypothesis seems reasonable. Most hvKps are sensitive to most antibiotics except for intrinsic resistance to ampicillin, similarly in this study all hvKp isolates were ampicillin resistant [40] . Studies have shown that hvKps are unlikely to take up DNA from other resistant bacteria due to the large size of the capsule and increased expression of capsule-related genes, therefore antibiotic resistance is less common in hvKps than in cKp isolates [41] . But contrary to expectations, the rate of resistance was not much different between the two variants hvKps and cKps. Our results are consistent with other studies that have shown that the rate of hvKp resistance is increasing worldwide [42, 43] . The current study revealed the high prevalence of MDR-hvKp and high resistance to imipenem (66%). Moreover, the presence of bla TEM , bla SHV , bla CTX-M-15 , bla OXA-48 and bla NDM-1 was detected simultaneously in 56.8% of hvKp isolates. There was a significant difference between hvKps and cKps about carbapenem-resistance genes so that bla OXA-48 was more frequent in carbapenem-resistant hvKps, and in contrast, bla NDM-1 was more detected in carbapenem-resistant cKp isolates. The reason for this difference is not clear. However, two pathways have been introduced for the emergence of MDR-hvKp strains, the horizontal acquisition of resistance genes by plasmids and mobile genetic elements (MGEs) by hvKp isolates (type I), and another pathway is the acquisition of the virulence-associated plasmid (e.g., pLVPK and pVir) by MDR-cKp (type II) [42, 44] . Ultimately, both mechanisms lead to the development of MDR-hvKp strains that are resistant to antibiotic treatment in addition to having a very high pathogenicity that poses a serious threat to public health. The correlation of ompK36 porin variants with specific sequence types (STs) of K. pneumoniae was first described by Papagiannitsis et al. K. pneumoniae isolates can be classified into four groups (designated groups A to D) by ompK36 genotyping [13, 45] . There is a relationship between ompK36 type and clonal group (CG). Also, some STs were reported to be associated with hvKp isolates, e.g. ST11 and ST23, that ST11 (CG258) belonged to ompK36 group A and ST23 (CG23) belonged to ompK36 group C [13, 46] . In this study, clonal relatedness by ompK36 typing revealed that group C (70.6%) was the most common ompK36 porin type among hvKp isolates. A study in Taiwan showed that ompK36 group C was significantly more abundant among K. pneumoniae isolates [46] . This study was in agreement with our study in Iran. This study presented a new rapid screening method based on the resistance of hvKp to tellurite, which was superior than string test in phenotypic identification of hvKp isolates. The consideration of phenotypic detection along with genotyping of hvKp render a reliable identification of hypervirulent strains. In this study, a high prevalence of MDR-hvKp and a high level of resistance to imipenem (66%) were detected. In addition, co-existence of bla TEM , bla SHV , bla CTX-M-15 , bla OXA-48 and bla NDM-1 was identified in 56.8% of hvKp isolates. Using the PCR-based ompK36 typing method, which was simpler and less expensive than MLST, we were also able to investigate the clonal relatedness of the strains. It was found that the majority of hvKp isolates belonged to capsular serotype K20 and ompK36 group C, which is related to CG23 (e.g. ST23). It seems that the expansion of MDR-hvKp in clinical settings is an inevitable event and this needs an urgent infection control program in the healthcare setting. The plasmids harboring virulence and antimicrobial resistance factors will change the clinical face of K. pneumoniae soon. Hypervirulent Klebsiella pneumoniae (hypermucoviscous and aerobactin positive) infection over 6 years in the elderly in China: antimicrobial resistance patterns, molecular epidemiology and risk factor Hypervirulent Klebsiella pneumoniae in cryptogenic liver abscesses Tracking microevolution events among ST11 carbapenemase-producing hypervirulent Klebsiella pneumoniae outbreak strains Molecular epidemiology of hypervirulent Klebsiella pneumoniae: a systematic review and meta-analysis Klebsiella species (K. pneumoniae, K. oxytoca, K. ozaenae and K. rhinoscleromatis) Hypervirulent Klebsiella pneumoniae A hybrid plasmid formed by recombination of a virulence plasmid and a resistance plasmid in Klebsiella pneumoniae High Prevalence of 16s rRNA methylase genes among carbapenem-resistant hypervirulent Klebsiella pneumoniae isolates in a Chinese Tertiary Hospital Distinct evolutionary dynamics of horizontal gene transfer in drug resistant and virulent clones of Klebsiella pneumoniae Carbapenem-resistant hypervirulent Klebsiella pneumoniae of sequence type 36 Emerging antimicrobial-resistant high-risk Klebsiella pneumoniae clones ST307 and ST147 Betalactamasas: la evolución del problema. Rev peru investig salud Associations of the major international high-risk resistant clones and virulent clones with specific ompK36 allele groups in Klebsiella pneumoniae in Taiwan Role of association of OmpK35 and OmpK36 alteration and blaESBL and/or bla ampC genes in conferring carbapenem resistance among non-carbapenemase-producing Klebsiella pneumoniae Multiresistant Gram-negative bacteria: the role of high-risk clones in the dissemination of antibiotic resistance The Galleria mellonella infection model does not accurately differentiate between hypervirulent and classical Klebsiella pneumoniae Identification of hypervirulent Klebsiella pneumoniae isolates using the string test in combination with Galleria mellonella infectivity Identification of biomarkers for differentiation of hypervirulent Klebsiella pneumoniae from classical K. pneumoniae Clinical and molecular characteristics of multi-clone carbapenem-resistant hypervirulent (hypermucoviscous) Klebsiella pneumoniae isolates in a tertiary hospital in Beijing, China -extended spectrum β-lactamases and class 1 integron among Enterobacter cloacae isolates collected from hospitals of Tehran and Qazvin Klebsiella pneumoniae genotype K1: an emerging pathogen that causes septic ocular or central nervous system complications from pyogenic liver abscess Hypervirulence and hypermucoviscosity: two different but complementary Klebsiella spp. phenotypes? Virulence Detection of multiple hypervirulent Klebsiella pneumoniae strains in a New York City hospital through screening of virulence genes Association of tellurite resistance with hypervirulent clonal groups of Klebsiella pneumoniae Klebsiella pneumoniae: going on the offense with a strong defense High prevalence of KPC-2-producing hypervirulent Klebsiella pneumoniae causing meningitis in Eastern China Tracking key virulence loci encoding aerobactin and salmochelin siderophore synthesis in Klebsiella pneumoniae Klebsiella pneumoniae yersiniabactin promotes respiratory tract infection through evasion of lipocalin 2 Yersiniabactin is a virulence factor for Klebsiella pneumoniae during pulmonary infection Outbreak of hypervirulent Klebsiella pneumoniae harbouring bla VIM-2 among mechanicallyventilated drug-poisoning patients with high mortality rate in Iran RmpA regulation of capsular polysaccharide biosynthesis in Klebsiella pneumoniae CG43 Mapping the evolution of hypervirulent Klebsiella pneumoniae Virulence characteristics of Klebsiella and clinical manifestations of K. pneumoniae bloodstream infections Clinical and molecular characteristics of emerging hypervirulent Klebsiella pneumoniae bloodstream infections in mainland China Capsular serotype K1 or K2, rather than magA and rmpA, is a major virulence determinant for Klebsiella pneumoniae liver abscess in Singapore and Taiwan First report of extended-spectrum betalactamase-producing Klebsiella pneumoniae among fecal carriage in Iran: high diversity of clonal relatedness and virulence factor profiles The emergence of the hypervirulent Klebsiella pneumoniae (hvKp) strains among circulating clonal complex 147 (CC147) harbouring bla NDM/OXA-48 carbapenemases in a tertiary care center of Iran Molecular characterization of carbapenem-resistant serotype K1 hypervirulent Klebsiella pneumoniae ST11 harbouring bla NDM-1 and bla OXA-48 carbapenemases in Iran An outbreak of carbapenem-resistant and hypervirulent Klebsiella pneumoniae in an intensive care unit of a major teaching hospital in Wenzhou Epidemiological characteristics and formation mechanisms of multidrug-resistant hypervirulent Klebsiella pneumoniae Virulence factors in hypervirulent Klebsiella pneumoniae. Front microbiol A fatal outbreak of ST11 carbapenem-resistant hypervirulent Klebsiella pneumoniae in a Chinese hospital: a molecular epidemiological study Increasing occurrence of antimicrobial-resistant hypervirulent (hypermucoviscous) Klebsiella pneumoniae isolates in China Whole genome assembly and functional portrait of hypervirulent extensively drug-resistant NDM-1 and KPC-2 co-producing Klebsiella pneumoniae of capsular serotype K2 and ST86 OmpK35 and OmpK36 porin variants associated with specific sequence types of Klebsiella pneumoniae Allocation of Klebsiella pneumoniae bloodstream isolates into four distinct groups by ompK36 typing in a Taiwanese university hospital Multiplex PCR for detection of seven virulence factors and K1/K2 capsular serotypes of Klebsiella pneumoniae Comparison of prevalence of virulence factors for Klebsiella pneumoniae liver abscesses between isolates with capsular K1/K2 and non-K1/K2 serotypes Fecal carriage of extended-spectrum β-lactamases and AmpC-producing Escherichia coli in a Libyan community Circulation of imipenem-resistant Acinetobacter baumannii ST10, ST2 and ST3 in a university teaching hospital from Tehran Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations The authors would like to thank the personnel in the bacteriology department of the Pasteur Institute of Iran for their help. This research was supported by the Pasteur Institute of Iran. This project was done based on ethical guidelines as previously approved by the Pasteur Institute of Iran (IR. B-9427). The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.