Extensive Expression of the Virulome Related to Antibiotic Genotyping in Nosocomial Strains of Klebsiella pneumoniae

Abstract

The emergence of hyper-virulent and multidrug-resistant (MDR) strains of Klebsiella pneumoniae isolated from patients with hospital- and community-acquired infections is a serious health problem that increases mortality. The molecular analysis of virulome expression related to antimicrobial-resistant genotype and infection type in K. pneumoniae strains isolated from patients with hospital- and community-acquired infections has been poorly studied. In this study, we analyzed the overall expression of the virulence genotype associated with the antimicrobial resistance genotype and pulse field gel electrophoresis (PFGE) type (PFtype) in K. pneumoniae. We studied 25 strains of K. pneumoniae isolated from patients who developed bacteremia and pneumonia during their hospital stay and 125 strains from outpatients who acquired community-acquired infections. Susceptibility to 12 antimicrobials was determined by Kirby–Bauer. The identification of K. pneumoniae and antibiotic-resistance genes was performed using polymerase chain reaction (PCR). To promote the expression of the virulence genes of K. pneumoniae, an in vitro infection model was used in human epithelial cell lines A549 and A431. Bacterial RNA was extracted with the QIAcube robotic workstation, and reverse transcription to cDNA was performed with the Reverse Transcription QuantiTect kit (Qiagen). The determination of the expression of the virulence genes was performed by real-time PCR. In addition, 57.3% (n = 86) of the strains isolated from patients with hospital- and community-acquired infections were multidrug-resistant (MDR), mainly to beta-lactam antibiotics (CB, AM, CFX, and CF), aminoglycosides (GE), quinolones (CPF and NOF), nitrofurantoin (NF), and sulfamethoxazole/trimethoprim (SXT). The most frequently expressed genes among strains isolated from hospital- and community-acquired infections were adhesion-type, ycfm (80%), mrkD (51.3%), and fimH (30.7%); iron uptake, irp2 (84%), fyuA (68.7%), entB (64.7%), and irp1 (56.7%); and protectins, rpmA (26%), which were related to antibiotic-resistance genes, bla_TEM (96%), bla_SHV (64%), bla_CITM (52.6%), bla_CTXM-1 (44.7%), tetA (74%), sul1 (57.3%), aac(3)-IV (40.7%), and aadA1 (36%). The results showed the existence of different patterns of expression of virulome related to the genotype of resistance to antimicrobials and to the PFtypes in the strains of K. pneumoniae that cause hospital- and community-acquired infections. These findings are important and may contribute to improving medical treatment strategies against infections caused by K. pneumoniae.

1. Introduction

The emergence of multidrug-resistant Klebsiella pneumoniae strains associated with hospital- and community-acquired infections, such as bacteremia, meningitis, pneumonia, and urinary tract infections (UTIs) [1,2,3], is a major concern for choosing the appropriate medical treatment [4]. The multidrug resistance of K. pneumoniae occurs mainly through the horizontal transfer of resistance genes through mobile genetic elements, such as plasmids, transposons, and integrons [5,6]. The frequency of genes encoding extended-spectrum beta–lactamases (ESBL) in K. pneumoniae strains has increased [7], as have genes for resistance to aminoglycosides, quinolones, tetracycline, and carbapenemases [8,9,10].

Recently, an increase in multidrug-resistant (MDR) and hyper-virulent (hvKpn) variants of K. pneumoniae [11,12], which cause highly invasive infections such as liver abscesses, endophthalmitis, meningitis, and septic arthritis [13], have been described in both healthy and immuno-compromised individuals [3]. The hypervirulence of K. pneumoniae has been associated with the hypermucoviscous phenotype, the K1 and K2 capsular serotype, the presence of the magA (mucoviscosisty-associated gene A) and rmpA (regulator of mucoid phenotype A) genes, and the salmochelin (iro) and aerobactin (iuc) siderophore systems [14].

The pathogenic capacity of K. pneumoniae is associated with the expression of a large number of virulence markers, which include adhesins (fimbriae Type 1 and Type 3), iron-acquisition systems, protectins, and toxins [15]. The genes that encode the fimbriae (fimH-1, kpn, ycfM, and mrkD) are considered the main virulence markers of K. pneumoniae that mediate adhesion to host cells. Therefore, their expression favors bacterial colonization [16], while the expression of siderophore genes, such as enterobactin (entB), yersiniabactin (irp1, irp2, fyuA, and ybtS), salmochelin (iroN) and aerobactin (iutA), allowing for iron uptake and facilitating survival and bacterial multiplication during infections [17], while the expression of the cnf-1 and hlyA toxin genes promote the degradation of the host tissues, facilitating bacterial dissemination [18].

The increase in mortality rates in patients with hospital- and community-acquired infections caused by hyper-virulent and multidrug-resistant strains of K. pneumoniae is a critical health situation [19], so it is currently important to analyze the properties molecular virulence and antimicrobial resistance of the new variants of K. pneumoniae with the purpose of improving medical treatment alternatives.

Molecular characterization of hyper-virulent hospital strains of MDR K. pneumoniae in Mexico [20,21], and in other parts of the world, have been studied [22,23,24]. However, reports on the expression of virulence markers, and their correlation with multidrug-resistance and clinical origin, are scarce. Therefore, in this study, we aimed to establish an in vitro model of infection using human cell lines to globally and comprehensively determine the virulome expression related to the antimicrobial resistance genotype and pulse field gel electrophoresis (PFGE) type in a group of K. pneumoniae strains from different infections.

2. Results

2.1. Multidrug-Resistance Phenotype

Utilizing the Kirby–Bauer disk diffusion method, our analysis revealed notable levels of antibiotic resistance among the tested strains, with the highest percentages of resistance observed for carbenicillin, ampicillin, and cefotaxime, as well as for nitrofurantoin (Table 1). The 57.3% (n = 86) of the K. pneumoniae strains isolated from patients with hospital- and community-acquired infections were multidrug-resistant in the range of 4–12 antimicrobials, within which more than 40% (n = 60) of the strains presented the phenotype of multi-resistance to nine antibiotics from different groups; beta-lactams (CB, AM, CFX, and CF), aminoglycosides (NET), quinolones (CPF and NOF), nitrofurans (NF), and sulfonamide/trimethoprim (SXT).

2.2. Origin of the Strains and In Vitro Infection of Human Epithelial Cell Lines to Determine Viruloma Expression

The majority of hospital-acquired K. pneumoniae strains were isolated from patients with bacteremia (21/150; Table 2) and community-acquired infections from patients with UTIs (61/150) and respiratory infection (53/150). Using the real-time PCR method, we found that the most frequently expressed genes among the strains after the in vitro infection model of the human epithelial cell lines A549 and A431 were the genes encoding adhesins (ycfm, mrkD, and fimH), the genes encoding iron-acquisition systems (irp2, fyuA, entB, and irp-1) and protectins (rpmA; Table 2). The frequency of expression of most virulence genes in the strains according to clinical origin was very similar, except for adhesin (mrkD) and enterobactin (entB), whose percentages were higher in pneumonia and bacteremia (hospital-acquired infections), respectively, while yersiniabactin (ybtS) and protectin (rpmA) percentages were higher in others (community-acquired infections). The genes expressed with the least frequency after the in vitro model of infection of the cell lines with K. pneumoniae strains were hlyA (n = 6) and magA (n = 1), both from strains isolated from patients with infections acquired in the community. The gene that codes for the CNF-1 toxin was not expressed in any of the strains studied.

2.3. Detection of Antibiotic Resistance Genes by Conventional PCR

The most frequently detected antibiotic-resistance genes by the conventional PCR method among K. pneumoniae strains were beta-lactams (bla_TEM, bla_SHV, bla_CITM, and bla_CTXM-1), tetracycline (tetA), sulfonamide (sul1), gentamicin (aac(3)-IV), streptomycin (aadA1), and trimethoprim (dfrA1). No statistically significant differences were found between the frequency of most antibiotic-resistance genes in strains according to clinical origin, except for bla_CTX-1 and aac(3)-IV, where frequencies were higher in strains isolated from bacteremia and pneumonia of patients with hospital acquired infections than in others (Table 3), respectively. The antibiotic resistance genes bla_CTXM-2 (n = 4) and qnr (n = 2) in strains isolated from patients with hospital-acquired infections (UTI and respiratory infection) had the lowest frequency.

2.4. Virulence Genes Expression According to Antibiotic Resistance Genes

Overall, in the K. pneumoniae strains isolated from patients with hospital-acquired infections, our real-time PCR assays showed an elevated expression of adhesion genes (mrkD and ycfM) and iron-acquisition systems (entB, irp1, irp2, and fyuA) associated with the presence of antibiotic-resistance genes to beta-lactam antibiotics (bla_SHV, bla_CITM, bla_TEM, and bla_CTXM-1), sulfonamides (sul1), and tetracycline (tetA) (Table 4). However, statistically, the individual frequency of expression of each virulence gene related to the detection of each of the antibiotic-resistance genes in the strains was different (Table 4).

Similarly, among the K. pneumoniae strains isolated from patients with community-acquired infections (Table 5), we detected a high expression by real-time PCR of adhesion genes (mrkD and ycfM), and iron acquisition systems (entB, irp1, irp2, and fyuA) related to detection by conventional PCR of resistance genes to beta–lactam antibiotics (bla_SHV and bla_TEM) and tetracycline (tetA). However, statistically, the individual frequency of expression by real-time PCR of each virulence gene related to the detection by conventional PCR of each of the antibiotic-resistance genes in the strains was different (Table 5).

2.5. PFGE and Expression Patterns of Virulence Genes

The PFGE analysis using XbaI showed two main clades (A and B) (Figure 1). Clade A, comprising 92 strains, was divided into two compact subgroups, A1 (n = 41) and A2 (n = 51). Subgroup A1 showed a Dice similarity index > 61.1% among the strains, and Subgroup A2 showed a similarity index of 65.1%. Subgroup A1 showed 10 different PFtypes (pulse field types), where the PFtype comprises four strains (41, 48, 84, and 154) isolated from patients with community infections (UTI (n = 2) and respiratory infection (n = 1)), and hospital infections (pneumonia (n = 1)) showed different expression profiles of virulence genes associated with the antibiotic resistance genotype, with strain 84 showing the broadest profile (mrkD/ycfM/rpmA/entB/irp-1/irp-2/fyuA/iutA/hlyA/aac(3)-IV/bla_SHV/tetA/bla_CTXM-1/bla_TEM). Similarly, in Subgroup A2, a PFtype comprising nine strains (12, 61, 73, 76, 94, 101, 112, 124, and 133) of different origins was identified, all of which had different virulence gene-expression profiles associated with the antibiotic-resistance genotype.

Clade B, comprising 57 strains, was distributed into two subgroups: B1 (n = 6) and B2 (n = 51) (Figure 1). Subgroup B1 showed a Dice similarity index > 65% among strains, and Subgroup B2 showed a similarity index of 67.9%. Subgroup B1 had only one PFtype, comprising 4 strains (80, 146, 102, and 81) of different origins and with different expression profiles of virulence genes associated with the antibiotic-resistance genotype. In Subgroup B2, 12 different PFtypes were identified, some of which comprised paired strains (119–116; 1–4; and 54–150), and one PFtype that comprised eight strains (50, 72, 122, 140, 142, 149, 151, and 153), all from community infections (UTI = 5; respiratory infection = 2, and other = 1), and with different virulence gene expression profile associated to the antibiotic resistance genotype.

2.6. Unsupervised Hierarchical Clustering

Unsupervised hierarchical clustering analysis showed three main groups based on similarities between K. pneumoniae strains (Figure 2). Group 1 (strain range 40–108) had the highest number of expressed virulence genes and was characterized by the simultaneous expression of one-to-seven virulence genes. Group 2 (strains between 70–147) had an intermediate number of expressed virulence genes (four to nine genes), and Group 3 (strains between 138–48) had the lowest number and was defined by the expression of 3–10 virulence genes.

Based on the expression of the virulome (VE) by real-time PCR, two different clades were identified: (A) Highly expressed among the strains (ymfM, irp-2, fyuA, entB, irp-1, and mrkD), and (B) Lowly expressed between strains (fimH, iutA, rpmA, kpn, ybtS, iron, hlyA, magA, and cnf). The antibiotic-resistance genotype (ARG) was grouped into two main clades according to their frequency of detection by conventional PCR. Clade C had a large number of genes detected (bla_TEM, tet(A), sul-1, bla_SHV, bla_CTXM-1, bla_CITM, aadA1, and aac(3)-IV), and Clade D was composed by a lower number of gene (dfrA-1, cmlA, cat1, tet(B), bla_CTXM-9, bla_CTXM-2, and qnr).

The antibiotic-resistance phenotype (ARP) was characterized by a high frequency of resistance to beta–lactams (AM, CB, CF, and CFX), as well as NF, SXT, and GE. The percentages of the resistance phenotype to beta–lactam antibiotics, aminoglycosides, SXT, and CL in the strains coincided with the resistance genotype for beta-lactams (bla_SHV, bla_TEM, bla_CITM, and bla_CTXM-1), aminoglycosides (aac(3)-IV)), sulfonamide (sul-1), trimethoprim (dfrA-1), and chloramphenicol (cat1 and cmlA) (Figure 2).

Overall, the unsupervised hierarchical analysis showed in the three groups (1, 2, and 3) strains with the same pattern of virolome expression. For example, in Group 2, it is observed in the cladogram that Strains 101, 98, 90, 25, 9, and 21 presented the same expression profile of the virulence genes (ycfM, irp-1, fyuA, entB, irp-1, and mrkD) associated with the frequency of the genotype and phenotype of resistance to antibiotics bla_TEM-TET(A)-AM-CB, and with the clinical origin. Interestingly, this shared pattern of expression was independent of the strain origin because Strains 101, 98, 90, and 9 were isolated from bacteremia of patients with hospital-acquired infections, as well as Strains 25 (infected ulcer) and 21 (respiratory infection) from patients with community-acquired infections. The Strains 113-50-79; 45-7-43 from Group 1, and 102-134, 42-104, and 2-12 from Group 3, also had the same virulome expression pattern.

3. Discussion

The emergence of hypervirulent and MDR hospital- and community-acquired strains of K. pneumoniae is a serious concern that increases mortality, mainly in immune-compromised patients [14,25]. Numerous virulence factors contribute to the pathogenicity of K. pneumoniae, thus increasing the severity of hospital and community infections [26]. Our working group previously determined the frequency of virulence genes by conventional PCR in these K. pneumoniae strains (n = 150), as well as its relationship with the hypermucoviscosity phenotype, capsular serotypes, biofilm formation, and the resistance phenotype to antibiotics [27]. In this work, we continue our research by characterizing in these K. pneumoniae strains the expression of the virulome by real-time PCR related to the antimicrobial resistance genotype, to the pulse field gel electrophoresis (PFGE) type, and their clinical origin.

In this study, the overall percentages for the expression of different virulence genes in K. pneumoniae strains isolated from patients with hospital- and community-acquired infections following the in vitro infection of human epithelial cell lines were very similar, except for the frequency of expression of mrkD, entB, ybtS, and rmpA genes, where there were statistically significant differences. We found a higher expression of the adhesion marker mrkD in strains from hospital-acquired infections (pneumonia) than in strains from community-acquired infections. The adhesion gene ycfM was expressed in all strains isolated from bacteremia (21/21), pneumonia (4/4), and stool culture and tumor biopsy (3/3; other). The co-expression of the adhesion genes fimH, mrkD, and ycfM detected in some strains of K. pneumoniae, following the infection of A549 and A431 cell lines, suggests that this expression combination could increase bacterial colonization and persistence in infected tissues. Type 1 fimbriae (fimH) was previously found to be frequently expressed during the infection of urinary tissue and participate in the invasion of bladder cells [16], while Type 3 fimbriae (mrkD) was found to promote biofilm formation and, consequently, the evasion of the host immune response [28].

Iron is an essential element for electron transport, DNA synthesis, and peroxide reduction, favoring bacterial survival and multiplication during infections [29]. Among the expression of iron uptake genes, entB expression was higher in bacteremia and pneumonia (hospital-acquired infections), and in respiratory infection (community-acquired infection) as well. High expression percentages of irp1 (bacteremia, respiratory infection, and others), irp2, and fyuA (bacteremia, UTI, respiratory infection, and others) were also found. The simultaneous high expression of the iron uptake genes entB, irp1, irp2, ybtS, and fyuA by some of the strains after infection of A549 and A431 cell lines highlights the virulence of the strains to cause acute infections because entB (enterobactin) and irp-1 (yersiniabactin) expression increases biofilm formation in Klebsiella pneumoniae, causing liver abscess [30]; ybtS (yersiniabactin) promotes respiratory tract infection [31], while irp and fyuA (yersiniabactin) reduce the bactericidal capacity of innate immune cells [32].

The expression of the rmpA gene (protectin) was higher in strains from community-acquired infection [respiratory infection (23/53) and others (2/3)] than in that from hospital-acquired infections. The toxin gene hlyA (hemolysin A) was only expressed in strains isolated from community-acquired infection [UTI (4/61) and respiratory infection (2/53)]. The expression of the mucoid phenotype A (rmpA) gene was detected mainly in the strains from respiratory infection and UTI. This proves the ability of the strains to cause more acute infections because rmpA has been previously associated with hvKpn strains, causing liver abscesses and bacteremia [33], while hlyA, which was more frequently expressed in strains originating from UTI (community-acquired infection), favors bladder inflammation during UTIs [34].

More than half of the K. pneumoniae strains (n = 86) isolated from patients with hospital- and community-acquired infections were multidrug-resistant (MDR), mainly to beta-lactam antibiotics (CB, AM, CFX and CF), aminoglycosides (GE), quinolones (CPF and NOF), nitrofurantoin (NF), and sulfamethoxazole/trimethoprim (SXT), which coincides with the multidrug-resistance described in strains of K. pneumoniae isolated from patients with hospital-acquired infections in Mexico [20].

The frequency of the detection of antibiotic-resistance genes in K. pneumoniae strains isolated from patients with hospital- and community-acquired infections was very similar, except for bla_CTXM-1 (beta-lactam) and aac(3)-IV (gentamicin), where the frequency was higher in strains isolated from bacteremia and pneumonia, respectively, both from hospital-acquired infections. High percentages for beta–lactam (bla_SHV and bla_TEM), tetracycline (tet(A)), and sulfonamide (sul1)-resistance genes were found in strains isolated from patients with hospital-acquired infections (bacteremia and pneumonia) and community-acquired infections (UTI, respiratory infection, infected ulcer, and others). High percentages for the streptomycin- (aadA1) and gentamicin- (aac(3)-IV) resistance gene were also detected in strains from patients with hospital-acquired infections (pneumonia) and community-acquired infections (UTI, respiratory infection, and infected ulcer). The overall number of bla_TEM (144/150) and bla_CXTX-9 (10/150) coincide with those described in K. pneumoniae strains isolated from Mexican patients in a health facility [35], while the numbers of bla_SHV (96/150), bla_CTXM-1 (67/150), and bla_TEM (144/150) are higher than those described in K. pneumoniae strains isolated from hospital-acquired infections in other parts of the world [36,37]. The wide distribution of the resistance genotype to other non-beta–lactam antibiotics, such as tetracycline (tet(A)), sulfonamide (sul1), gentamicin (aac(3)-IV), and streptomycin (aadA1) detected in K. pneumoniae strains from patients with hospital- and community-acquired infections coincides with the high MDR phenotype found in Mexico in hospital strains of K. pneumoniae from bacteremia, UTIs, respiratory infection, and other infections [35,38].

In this study, the correlation of the expression percentages of adhesion genes (mrkD and ycfM) and iron-acquisition systems (entB, irp-1, irp2, and fyuA) with the genotype of resistance to beta–lactam antibiotics (bla_SHV, bla_TEM, and bla_CTXM-1), sulfonamide, and tetracycline was higher in K. pneumoniae strains isolated from patients with hospital-acquired infections compared with strains from community-acquired infections. These results are substantial and suggest that the expression of the virulence genotype related to the antimicrobial-resistance genotype in strains from bacteremia and pneumonia (hospital-acquired infections) may cause more acute infections, increasing patient mortality. Mortality in patients with hypervirulent K. pneumoniae pneumonia has been found to be 23.1% [39], 54.3% in bacteremia, 48.9% in intensive care, 43.1% in organ transplantation, and 13.5% in urinary tract infection [40]. The prevalence of the cnf-1 gene (n = 22) and the magA gene (n = 25), which had been previously identified by our research group within the genome of these K. pneumoniae strains [27], does not correspond with the expression results in this study using an in vitro infection model in A549 and A431 cell lines with the same strains. Interestingly, cnf-1 was not expressed in any of the strains, whereas magA was expressed in only one strain. This observed inconsistency in detection and expression is probably attributable to the cell lines used in our study since prior reports have demonstrated differential levels of expression according to the infected cell type. One example is that cnf-1 tends to have a higher frequency of expression during urinary tract infections [41], whereas magA expression has been linked to hypermucoviscosity in strains associated with liver abscesses [42]. Furthermore, the temporal dynamics of the infectious process have a profound effect on gene expression. Some genes are known to be expressed early in the infection, while others have later expression. For instance, the early-expressed type 1 fimbrial genes are regulated by an invertible element located in the promoter region [43]. This element permits gene expression to be repressed or activated based on the prevailing conditions of the infection process.

PFGE is a frequently used method for clonal typing of isolates in outbreaks and allows for source identification [44]. In this study, 13 different PFtypes distributed in Subtypes B1 and B2 of Clade B were identified, some of them integrated by two, three, four, and up to eight strains. In Clade A, 19 PFtypes distributed in Subtypes A1 and A2 were identified. These comprehensive findings showed different virulome expression patterns associated with the antibiotic-resistance genotype and clinical origin in K. pneumoniae strains from the same PFtype, where the largest PFtype (subtype A2) was found to comprise nine strains (12, 61, 73, 76, 94, 101, 112, 124, and 133) isolated from hospital-acquired (bacteremia [n = 3]) and community-acquired infections (respiratory infection [n = 3], and infected ulcers [n = 2] and UTI [n = 1]), all with different virulence gene expression profiles associated with the antibiotic resistance genotype.

The results showed that K. pneumoniae strains from the same PFtype, isolated from different patients with hospital or community infections, are possibly found in the hospital environment. Therefore, some patients with community-acquired strains could have acquired these strains during their hospital stay.

The unsupervised hierarchical clustering analysis, similar to the PFGE method, revealed a diverse spectrum of virulome expression patterns associated with the genotype-phenotype characteristics of antibiotic resistance and the origin of K. pneumoniae strains isolated from patients with either community- or hospital-acquired infections. High throughput technologies, such as whole genome sequencing, could provide a more comprehensive understanding of the global composition of the virulence and antimicrobial-resistance genes. These analyses will be performed in future works. To the best of our knowledge, this is the first study performed in Mexico where the wide distribution of different virulome expression patterns associated with the antibiotic resistance genotype in K. pneumoniae strains isolated from patients with different hospital- and community-acquired infections is shown, suggesting the high pathogenic potential of these strains during infections.

4. Materials and Methods

4.1. Origin of the Strains

A total of 150 strains of K. pneumoniae were analyzed, which were collected from September 2019 to March 2020 in the Microbiology laboratory of the Hospital General Regional No. 72 (Instituto Mexicano del Seguro Social, State of Mexico, Mexico) located in the municipality of Tlalnepantla de Baz, Edo. from Mexico, Mexico. The strains were isolated from samples of hospitalized patients with ongoing infections, such as bacteremia (n = 21) and pneumonia (n = 4), which they acquired during their hospital stay after admission for managing complications arising from other comorbidities such as diabetes, high blood pressure, chronic obstructive pulmonary disease, and obesity. Additionally, strains were collected from the non-hospitalized external community, encompassing UTIs (urinary tract infection; n = 61), respiratory (n = 53), ulcers (n = 8) and other infections (stool culture (n = 2), and tumor biopsy (n = 1)). All patients provided informed consent for inclusion in this study. The study was approved by the institution’s Ethics Committee (identification code: R-2022-1406-033).

4.2. Bacterial DNA Extraction

The DNA of the K. pneumoniae strains was extracted by the boiling method described previously [45]. Each of the strains was seeded on eosin methylene blue agar (EMB; DIBICO, Edo. de México, Mexico) at 37 °C for 12 h under constant agitation. Four isolated 2-mm (in diameter) colonies were collected, suspended in Eppendorf tubes with 200 µL of sterile water, incubated at 100 °C for 10 min, and centrifuged at 10,000× g for 5 min. DNA obtained in the supernatant was stored at −20 °C. In order to use 100 ng/µL DNA for the PCR reaction, DNA concentration and purity were measured using a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA).

4.3. Identification of K. pneumoniae through Polymerase Chain Reaction (PCR)

K. pneumoniae strains were identified using PCR amplification of 16S–23S rDNA internal transcribed spacer [46]. The final volume per reaction mixture for each uniplex PCR assay was 20 µL; 1 µL of forward primer and 1 µL of reverse primer (10 pmol, Integrated DNA Technologies, San Diego, CA, USA), 12 µL of Taq DNA Polymerase 2× Master Mix RED (AMPLIQON, Copenhagen, Denmark), 3 µL of nuclease-free water, and 3 µL of DNA template (100 ng). The PCR amplification conditions were performed using a T100™ Thermal Cycler (Bio-Rad, Feldkirchen, Germany) instrument as follows; initial denaturation at 94 °C for 5 min, 30 cycles of denaturation at 94 °C for 30 s, annealing at 55 °C for 30 s, and extension at 72 °C for 50 s; with a final extension of 72 °C for 5 min. The K. pneumoniae ATCC 700721 strain was used as a positive control. The amplified DNA fragments (130 pb) were visualized using 2% agarose gels stained with Midori Green Direct (Nippon Genetics, Düren, Germany).

4.4. Multidrug-Resistance Phenotype

The Kirby–Bauer disk diffusion method (Diagnostic Research, CDMX, Mexico) was used to determine resistance to the following 12 antimicrobials; carbenicillin (CB; 100 µg), ampicillin (AM; 10 µg), cefotaxime (CFX; 30 µg), cephalothin (CF; 30 µg), gentamicin (GE; 10 µg), amikacin (AK; 30 µg), netilmicin (NET; 30 µg), ciprofloxacin (CPF; 5 µg), norfloxacin (NOF; 10 µg), nitrofurantoin (NF; 300 µg), chloramphenicol (CL; 30 µg), and trimethoprim with sulfamethoxazole (SXT; 25 µg). For the reproducibility of the method, the Escherichia coli ATCC 25922 strain was used as a control. The criteria established by the Clinical and Laboratory Standards Institute were used for the interpretation of the results [47]. K. pneumoniae strains were classified as multidrug resistant (MDR) when they presented resistance to more than three antibiotics from different groups.

4.5. Detection of Antibiotic-Resistance Genes Using Conventional PCR

The conventional polymerase chain reaction (PCR) method was used to determine on the chromosome of the K. pneumoniae strains; the presence or absence of the genes encode resistance to beta–lactams (bla_SHV, bla_CITM, bla_TEM, bla_CTX-1, bla_CTXM-2, and bla_CTXM-9), streptomycin (aadA1), gentamicin (aac(3)-IV), sulfonamide (sul1), chloramphenicol (cat1 and cmlA), tetracycline (tet(A) and tet(B)), trimethoprim (dfrA1), and quinolones (qnr). The primers and PCR conditions used for detection of antibiotic resistance genes were as described previously [48,49]. The final volume-per-reaction mixture for each uniplex PCR assay was 20 µL; 1 µL of forward primer and 1 µL of reverse primer (10 pmol, Integrated DNA Technologies, San Diego, CA, USA), 12 µL of Taq DNA Polymerase 2× Master Mix RED (AMPLIQON, Copenhagen, Denmark), 3 µL of nuclease- free water, and 3 µL of DNA template (100 ng). The K. pneumoniae strains from clinical isolates Kp7, Kp17, Kp37, Kp39, and Kp74, which harbor the antibiotic-resistance genes studied, were used as positive controls.

4.6. Bacterial Dilution of Strains for Infection of Human Epithelial Cell Lines

One colony (2 mm in diameter) of every K. pneumoniae strain was seeded separately in brain heart infusion broth (MCD Lab, Edo. de México, Mexico) and incubated for 12 h at 37 °C under constant shaking. Following the instructions in the RNA protect Bacteria Reagent Manual (Qiagen, Hilden, Germany), a 1:4 dilution of each K. pneumoniae culture was performed using phosphate-buffered saline solution to obtain an optical density at 600 nm = 1.0 using a Beckman DU-7400 spectrophotometer (Laguna Hills, CA, USA), which corresponded to a concentration of 1 × 10⁹ cells/mL. Subsequently, dilutions were performed in order to obtain a concentration of 2 × 10⁶ cells/mL.

4.7. Preparation of Human Epithelial Cell Lines and In Vitro Infection Model

The in vitro infection model using human epithelial cell lines to determine virulence gene expression was performed as previously described [50,51]. With the aim of emulating a human infection process, we prepared cultures of human lung cancer A549 (ATCC CCL-185) and human epidermoid A431 (CRL-1555, Manassas, VA, USA) cell lines. The cell lines were plated in DMEM medium (Dulbecco’s Modified Eagle Medium, Corning NY, USA) and supplemented with 10% FBS (Fetal Bovine Serum, Corning NY, USA) in a 24-well cell culture plate, until a confluency of 1.8 × 10⁵ cells/well was reached for the in vitro infection experiment. All cultures were incubated at 37 °C under 5% of CO₂ atmosphere.

To stimulate the bacterial infection of the human cell lines, each dilution of K. pneumoniae culture isolated from patients with bacteremia, pneumonia, and respiratory, ulcer, and other infections were inoculated (50 µL; 2 × 10⁶ cells/mL) onto the surface of a monolayer of the cultured 1.8 × 10⁵ cells of the epithelial cell line A549 (ATCC CCL-185) derived from lung cancer, while dilutions (50 µL; 2 × 10⁶ cells/mL) of K. pneumoniae from patients with UTIs were inoculated on the monolayer of the human epidermoid cell line A431 (CRL-1555, Manassas, VA, USA). The 24-well plates were incubated with 1 mL of F12K plus 10% fetal bovine serum at 37 °C for 48 h under a 5% CO₂ atmosphere with saturated humidity. The maintenance medium (F12K plus 10% fetal bovine serum) was changed every 24 h. The K. pneumoniae strain, ATCC 700721, was used as a control.

4.8. K. pneumoniae RNA Extraction and Reverse-Transcription to cDNA

To establish the expression of virulence genes, RNA was extracted from K. pneumoniae strains after in vitro infection of human epithelial cell lines, for which K. pneumoniae strains were harvested from the surface of A549 and A431 cell line cultures and suspended in 1000 µL of RNA Protect Bacteria reagent (Qiagen). Samples were centrifuged at 8000× g for 10 min to obtain bacterial cell sediment. RNA was extracted from the bacterial sediment using a QIAcube automated extraction equipment (Qiagen, Hilden, Germany) with the commercial RNeasy Mini Kit (Qiagen), which involved bacterial lysis with TE buffer (10 mM Tris-HCl, 1 mM EDTA, pH 8) containing 1 mg/mL lysozyme. A NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA, EE. UU) was used to determine the concentration and purity of total RNA. The Reverse Transcription QuantiTect kit (Qiagen) was used to perform first-strand cDNA synthesis.

4.9. Determination of K. pneumoniae Virulome Expression Using Real-Time PCR

To determine virulome expression in K. pneumoniae strains by real-time PCR, the Rotor-Gene Q 5plex HRM device (Qiagen, Hilden, Germany) was used. The primers used to assess the expression of the virulence genes of K. pneumoniae virulence genes encoding for adhesins, fimH-1 (type 1 fimbriae), mrkD (type 3 fimbriae), kpn (fimH-type adhesin), and ycfM (outer membrane lipoprotein); iron-acquisition systems, entB (enterobactin biosynthesis), irp1, irp2, ybtS (yersiniabactin biosynthesis), fyuA (yersiniabactin receptor), iutA (aerobactin receptor), and iroN (catecholate siderophore receptor); protectins, magA (mucoviscosity-associated gene A) and rmpA (mucoid regulator phenotype A), and toxins; hlyA (hemolysin) and cnf-1 (necrotizing cytotoxic factor 1) were as previously described [15]. The final volume-per-reaction mixture for each real-time PCR assay was 20 µL: 12.5 µL of the Rotor–Gene SYBR Green PCR kit master mix (Qiagen), 1 µL of each forward and reverse primer (1 µM), 8.5 µL of RNAase-free water, and 2 µL of the cDNA (100 ng). Amplification was performed at 95 °C for 5 min, followed by 40 cycles of 95 °C for 5 s, and an extension at 60 °C for 10 s. For each quantitative real-time PCR (qPCR) assay, we performed a standard curve prepared from three dilutions of cDNA (100, 200, and 300 ng/µL) from a strain harboring the gene of interest (positive control). From these three cDNA dilutions, the Rotor–Gene Q 5plex HRM System software version 2.3.1.49 (Qiagen) calculated the cutoff for the definition of the threshold cycle (TC) of each strain. In this way, with the CT values of each strain obtained with respect to the CT of the standard curve of the control strain, the arithmetic mean was obtained, and the expression percentage of each gene was determined. Each real-time PCR assay included a melting curve, a housekeeping gene (rpoB), and a non-templated control (NTC). The K. pneumoniae strains from clinical isolates Kp1, Kp6, and Kp12, harboring the 15 virulence genes studied (fimH-1, mrkD, kpn, ycfM, entB, irp1, irp2, ybtS, fyuA, iutA, iroN, magA rmpA, hlyA, and cnf-1) were used as positive controls for the preparation of the standard curve.

4.10. Pulse Field Gel Electrophoresis (PFGE)

The preparation of agarose blocks with K. pneumoniae genomic DNA was performed using the method previously described by the standardized PFGE protocol PulseNet, USA [52]. To separate XbaI fragments, electrophoresis was performed on 1% of agarose gels and 0.5X TBE buffer at 14 °C with an increased pulse time from 2.2 to 54.2 s for 20 h and 6.0 V/cm using the CHEF–MAPPER (Bio-Rad, Hercules, CA, USA) device. XbaI fragment sizes were estimated using XbaI fragments of Salmonella enterica serotype Braenderup global standard H9812. The Gel Logic 112 image photo documentation system (Kodak, Rochester, NY, USA) was used to digitize the images. The BioNumerics v.7.1 software package (Applied Maths, Sint Martens-Latem, Belgium) was used to analyze the fingerprint profile on the PFGE gel. After background subtraction and gel normalization, the fingerprints in the profiles were typed based on band similarity and dissimilarity, using Dice’s similarity coefficient and the unweighted pairwise clustering method with arithmetic mean according to average linkage clustering methods. PFGE patterns were interpreted according to the criteria previously described for the typing of bacterial strains [53]. (1) Identical: Strains are designated as genetically indistinguishable if their restriction patterns have the same number of bands. (2) Closely related: The PFGE patterns of the strains present two or three bands of difference. (3) Possibly related: The PFGE pattern of the strains differs between four and six bands. (4) Unrelated: The PFGE pattern differs in seven or more bands.

4.11. Unsupervised Hierarchical Clustering

K. pneumoniae strains were systematically grouped according to virulome expression, antibiotic-resistance genotype, antibiotic-resistance phenotype, diagnosis, and clinical origin from patients with hospital- and community-acquired infections using unsupervised hierarchical clustering based on Euclidean distances for categorical variables. A categorical data matrix that included virulence gene expression, antibiotic-resistance genotype, antibiotic-resistance phenotype (resistance and sensible), and the clinical origin of patients with hospital-acquired infections and community-acquired infections, as well as diagnosis, was created in R (v.3.6.2). The distance of each strain was calculated based on the overall similarity coefficient, which estimates the maximum possible absolute discrepancy between each combined pair of strains. Strains were visualized in a distribution diagram with a dendrogram constructed using ComplexHeatmap (v3.6.2, R core).

4.12. Statistical Analysis

The χ² test was applied using the SPSS statistical software (version 20.0; SPSS Inc., Chicago, IL, USA) (p < 0.05) to establish differences between the frequency of expression of virulence genes related to the antibiotic-resistance genotype and clinical origin in K. pneumoniae strains isolated from patients with hospital- and community-acquired infections.

5. Conclusions

Our results highlighted the existence of different PFtypes of K. pneumoniae causing hospital- and community-acquired infections, with wide virulome expression profiles related to the antimicrobial-resistance genotype and infection type. Therefore, it is important to implement new monitoring and control programs on the emergence of new hypervirulent and MDR K. pneumoniae strains to improve medical treatment and reduce mortality rates.

Author Contributions

Conceptualization, G.L.P.-C. and E.M.-P.; Methodology, A.B.-C., G.D., C.E.D.-V., A.H.d.l.C.-M., L.R.G.-C. and M.P.S.-Y.; Software, G.D., A.B.-C. and R.M.-E.; Validation, G.L.P.-C. and E.M.-P.; Formal Analysis, G.L.P.-C. and E.M.-P.; Investigation, G.L.P.-C., E.M.-P., R.M.-E., A.B.-C. and F.V.-P.; Resources, G.L.P.-C. and E.M.-P.; Data Curation, G.L.P.-C., E.M.-P., A.B.-C. and M.P.S.-Y.; Writing—Original Draft Preparation, G.L.P.-C. and E.M.-P.; Writing—Review and Editing, G.L.P.-C., E.M.-P. and F.V.-P.; Visualization, G.L.P.-C. and E.M.-P.; Supervision, G.L.P.-C. and E.M.-P.; Project Administration, G.L.P.-C. and E.M.-P.; Funding Acquisition, G.L.P.-C. and E.M.-P. All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement

The study was approved by the Ethics and Research Committee of the Regional General Hospital No. 72 (Instituto Mexicano del Seguro Social) (protocol code R-2022-1406-033).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Not applicable.

Acknowledgments

We want to thank the UNAM-DGAPA-PAPIIT IN225020 Program.

Conflicts of Interest

The authors declare no conflict of interest.

Footnotes

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Figures and Tables

Figure 1. Pulse-field gel electrophoresis (PFGE) profile dendrogram and expression patterns of virulence genes and antibiotic-resistance genotyping in K. pneumoniae strains from hospital- and community-acquired infections. The dendrogram was produced through the Dice similarity coefficient and unweighted pairwise clustering method with arithmetic mean clustering methods using PFGE imaging of XbaI-digested genomic DNA. The scale bar shows the correlation coefficient (%). Capital letters (A,B) represent the two main clades produced. The letters B1, B2, A1, and A2 represent subgroups of strains. The strain number, type of infection, origin, and virulence gene expression pattern associated with the antibiotic resistance genotype are shown on the right.

Figure 2. Hierarchical grouping of K. pneumoniae strains according to virulome expression related to antibiotic-resistance genotype–phenotype, diagnosis, and clinical origins. Positivity and negativity for a particular genotype is represented by a red and grey rectangle, respectively. Upper panel: Virulome expression (VE). Middle panel: Antibiotic-resistance genotype (ARG). Lower panel: Antibiotic-resistance phenotype (ARP). Left axis: Expression of virulence factors (Clades A and B) and detected antibiotic-resistance genes (Clades D and C). Right axis: Absolute frequency of expression by gene and antibiotic resistance. Top: Cladogram total of VE, ARG, and ARP of the strains. Bottom: Diagnosis and origin of the strains. ARP section: AM = ampicillin, CB = carbenicillin, NF = nitrofurantoin, CF = cephalothin, SXT = trimethoprim-sulfamethoxazole, CFX = cefotaxime, GE = gentamicin, CPF = ciprofloxacin, NOF = norfloxacin, NET = netilmicin, AK = amikacin and CL = chloramphenicol.

Table 1

Antibiotic resistance phenotype in strains of Klebsiella pneumoniae isolated from patients with hospital- and community-acquired infections.

Antibiotics	Antibiotic Group	Number of Resistant Strains (%)
Carbenicillin (CB)	beta-lactams	148 (98.6)
Ampicillin (AM)		148 (98.6)
Cefotaxime (CFX)		66 (44)
Cephalothin (CF)		81 (54)
Gentamicin (GE)	Aminoglycosides	63 (42)
Amikacin (AK)		28 (18.6)
Netilmycin (NET)		42 (28)
Ciprofloxacin (CPF)	Quinolones	64 (42.6)
Norfloxacin (NOF)		61 (40.6)
Nitrofurantoin (NF)	Nitrofurans	86 (57.3)
Chloramphenicol (CL)	Phenicols	31 (20.6)
Sulfamethoxazole/trimethoprim (SXT)	Sulfonamide/Trimethoprim	67 (44.6)

Table 2

Frequency of expression of virulence genes in strains of Klebsiella pneumoniae isolated from patients with hospital- and community-acquired infections.

Origin of the Strains		Number of Strains (%)
		Adhesins				Iron Acquisition Systems							Protectins		Toxins
		fimH	mrkD	kpn	ycfM	entB	irp1	irp2	ybtS	fyuA	iutA	iroN	magA	rmpA	hlyA	cnf-1
Hospital-acquired(n = 25)	Bacteremia(n = 21)	5 (23.8)	15 (71.4)	5(23.8)	21 (100)	19 (90.5)	14 (66.7)	19 (90.5)	6(28.6)	17 (81)	2 (9.5)	0	0	0	0	0
	Pneumonia(n = 4)	1 (25)	4 (100)	0	4 (100)	3 (75)	2 (50)	2 (50)	1 (25)	3 (75)	3 (75)	0	0	1 (25)	0	0
Community-acquired (n = 125)	UTI(n = 61)	22 (36.1)	13 (21.3)	6 (9.8)	45 (73.8)	27 (44.3)	27 (44.3)	52 (85.2)	2 (3.3)	38 (62.3)	17 (27.9)	7 (11.5)	0	11 (18)	4 (6.6)	0
	Respiratory infection(n = 53)	16 (30.2)	37 (69.8)	6(11.3)	40 (75.5)	41 (77.4)	36 (67.9)	45 (84.9)	5 (9.4)	38 (71.7)	20 (37.7)	5 (9.4)	1 (1.9)	23 (43.4)	2 (3.8)	0
	Infected ulcer(n = 8)	2 (25)	6 (75)	1 (12.5)	7(87.5)	5 (62.5)	4 (50)	5 (62.5)	0	4 (50)	1(12.5)	1 (12.5)	0	2 (25)	0	0
	* Other(n = 3)	0	2 (66.7)	1 (33.3)	3 (100)	2 (66.7)	2 (66.7)	3 (100)	1(33.3)	3 (100)	1(33.3)	1(33.3)	0	2 (66.7)	0	0
p-value		0.732	1.8 × 10−8	0.264	0.054	0.0001	0.103	0.422	0.009	0.265	0.058	0.210	0.593	0.0001	0.763	-
Total (n = 150)		46 (30.7)	77 (51.3)	19 (12.7)	120(80)	97 (64.7)	85 (56.7)	126 (84)	15 (10)	103 (68.7)	44 (29.3)	14 (9.3)	1 (0.7)	39 (26)	6 (4)	0

Significant p-values (<0.05) are shown in bold. * Other: stool culture (n = 2) and tumor biopsy (n = 1).

Table 3

Frequencies of antibiotic-resistance genes in strains of Klebsiella pneumoniae isolated from patients with hospital- and community-acquired infections.

Origin of the Strains		Number of Strains (%)
		Beta–Lactams						Streptomycin	Gentamicin	Sulfonamide	Chloramphenicol		Tetracycline		Trimethoprim	Quinolones
		blaSHV	blaCITM	blaTEM	blaCTXM-1	blaCTXM-2	blaCTXM-9	aadA1	aac(3)-IV	sul1	cat1	cmlA	tet(A)	tet(B)	dfrA1	qnr
Hospital-acquired(n = 25)	Bacteremia(n = 21)	16 (76.2)	10 (47.6)	20 (95.2)	17 (80.9)	0	1 (4.7)	5 (23.8)	2 (9.5)	12 (57.1)	5 (23.8)	2 (9.52)	20 (95.2)	1 (4.7)	1 (4.7)	0
	Pneumonia(n = 4)	2 (50)	3 (75)	4 (100)	2 (50)	0	0	2 (50)	3 (75)	3 (75)	1 (25)	1 (25)	3 (75)	1 (25)	0	0
Community-acquired(n = 125)	UTI(n = 61)	39 (63.9)	36 (59)	56 (91.8)	30 (49)	3 (4.9)	4 (6.5)	20 (32.7)	32 (52.4)	32 (52.4)	3 (4.9)	10 (16.3)	46 (75.4)	8 (13.1)	14 (22.9)	2 (3.2)
	Respiratory infection(n = 53)	32 (60.3)	26 (49)	53 (100)	13 (24.5)	1 (1.8)	5 (9.4)	24 (45.2)	19 (35.8)	34 (64.1)	7 (13.2)	7 (13.2)	35 (66)	4 (7.5)	10 (18.8)	0
	Infected ulcer(n = 8)	5 (62.5)	4 (50)	8 (100)	3 (37.5)	0	0	3 (37.5)	5 (62.5)	4 (50)	1 (12.5)	1 (12.5)	5 (62.5)	1 (12.5)	2 (25)	0
	Other(n = 3)	2 (66.6)	0	3 (100)	2 (66.6)	0	0	0	0	1 (33.3)	0	0	2 (66.6)	0	1 (33.3)	0
p-value		0.840	0.389	0.363	0.0001	0.765	0.941	0.363	0.001	0.709	0.125	0.906	0.094	0.565	0.334	0.710
Total (n = 150)		96 (64)	79 (52.6)	144 (96)	67 (44.7)	4 (2.7)	10(6.7)	54 (36)	61 (40.7)	86 (57.3)	17 (11.3)	21 (14)	111 (74)	15 (10)	28 (18.7)	2(1.3)

Significant p-values (<0.05) are shown in bold.

Table 4

Distribution of virulence genes expression according to antibiotic resistance genes in strains of Klebsiella pneumoniae isolated from patients with hospital-acquired infections.

Hospital-Acquired Infections (n = 25)
Antibiotic Resistance Gene		Number of Strains (%)
		Adhesins				Iron Acquisition Systems							Protectins		Toxins
		fimH	mrkD	kpn	ycfM	entB	irp1	irp2	ybtS	fyuA	iutA	iroN	magA	rmpA	hlyA	cnf-1
Beta–lactams	blaSHV	6 (24)	18 (72)	4 (16)	24 (96)	18 (72)	15 (60)	20 (80)	7 (28)	18 (72)	5 (20)	0	0	1 (4)	0	0
	blaCITM	4 (16)	12 (48)	2 (8)	13 (52)	10 (40)	8 (32)	12 (48)	5 (20)	10 (20)	3 (12)	0	0	1 (4)	0	0
	blaTEM	5 (20)	22 (88)	5 (20)	24 (96)	22 (88)	16 (64)	20 (80)	7 (28)	20 (80)	5 (20)	0	0	1 (4)	0	0
	blaCTXM-1	6 (24)	13 (52)	5 (20)	19 (76)	16 (64)	11 (44)	17 (68)	5 (20)	15 (60)	3 (12)	0	0	0	0	0
	blaCTXM-2	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
	blaCTXM-9	0	1 (4)	1 (4)	0	1 (4)	0	1 (4)	0	1 (4)	0	0	0	0	0	0
Streptomycin	aadA1	1 (4)	6 (24)	1 (4)	7 (28)	6 (24)	5 (20)	6 (24)	3 (12)	6 (24)	1 (4)	0	0	1 (4)	0	0
Gentamicin	aac(3)-IV	1 (4)	4 (16)	0	5 (20)	4 (16)	2 (8)	3 (12)	1 (4)	3 (12)	2 (8)	0	0	1 (4)	0	0
Sulfonamide	sul1	4 (16)	13 (52)	2 (8)	15 (60)	13 (52)	11 (44)	12 (48)	3 (12)	11 (44)	3 (12)	0	0	1 (4)	0	0
Chloramphenicol	cat1	1 (4)	4 (16)	1 (4)	5 (20)	5 (20)	5 (20)	5 (20)	2 (8)	5 (20)	2 (8)	0	0	1 (4)	0	0
	cmlA	1 (4)	2 (8)	1 (4)	3(12)	3 (12)	1 (4)	3 (12)	1 (4)	2 (8)	1 (4)	0	0	1 (4)	0	0
Tetracycline	tet(A)	6 (24)	14 (56)	5 (20)	21 (84)	15 (60)	14 (56)	20 (80)	5 (20)	18 (72)	4 (16)	0	0	3 (12)	0	0
	tet(B)	0	2 (8)	0	2 (8)	2 (8)	1(4)	1 (4)	1 (4)	1 (4)	0	0	0	1 (4)	0	0
Trimethoprim	dfrA1	0	0	1 (4)	1 (4)	1 (4)	0	1 (4)	0	1 (4)	0	0	0	0	0	0
Quinolones	qnr	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
p-value		7.1 × 10−11	0.0004	1.2 × 10−7	0.0004	0.0004	0.0004	0.0004	0.0004	0.0004	0.0004	-	-	1	-	-

Significant p-values (<0.05) are shown in bold.

Table 5

Distribution of virulence genes expression according to antibiotic resistance genes in strains of Klebsiella pneumoniae isolated from patients with community-acquired infections.

Community-Acquired Infections (n = 125)
Antibiotic Resistance Gene		Number of Strains (%)
		Adhesins				Iron-Acquisition Systems							Protectins		Toxins
		fimH	mrkD	kpn	ycfM	entB	irp1	irp2	ybtS	fyuA	iutA	iroN	magA	rmpA	hlyA	cnf-1
Beta-lactams	blaSHV	31 (24.8)	40 (32)	8 (6.4)	72 (57.6)	56 (44.8)	51 (40.8)	76 (60.8)	5 (4)	58 (46.4)	28 (22.4)	12 (9.6)	0	26 (20.8)	3 (2.4)	0
	blaCITM	18 (14.4)	30 (24)	6 (4.8)	56 (44.8)	41 (32.8)	38 (30.4)	59 (47.2)	3 (2.4)	45 (36)	23 (18.4)	8 (6.4)	0	18 (14.4)	3 (2.4)	0
	blaTEM	38 (30.4)	57 (45.6)	14 (11.2)	91 (72.8)	73 (58.4)	65(52)	98 (78.4)	8 (6.4)	75 (60)	39 (31.2)	12 (9.6)	1 (0.8)	35 (28)	5 (4)	0
	blaCTXM-1	14 (11.2)	20 (16)	6 (4.8)	36 (28.8)	26 (20.8)	27 (21.6)	42 (33.6)	3 (2.4)	30 (24)	12 (9.6)	6 (4.8)	0	9 (7.2)	2 (1.6)	0
	blaCTXM-2	1 (0.8)	1 (0.8)	0	1 (0.8)	2 (1.6)	2 (1.6)	2 (1.6)	0	2 (1.6)	1 (0.8)	0	0	1 (0.8)	0	0
	blaCTXM-9	3 (2.4)	3 (2.4)	1 (0.8)	10 (8)	5 (4)	4 (3.2)	5 (4)	0	5 (4)	4 (3.2)	2 (1.6)	0	1 (0.8)	0	0
Streptomycin	aadA1	15 (12)	23 (18.4)	4 (3.2)	39 (31.2)	31 (24.8)	26 (20.8)	40 (32)	2 (1.6)	32 (25.6)	17 (13.6)	5 (4)	0	10 (8)	1 (0.8)	0
Gentamicin	aac(3)-IV	17 (13.6)	23 (18.4)	6 (4.8)	41 (32.8)	35 (28)	29 (23.2)	48 (38.4)	1 (0.8)	36 (28.8)	18 (14.4)	8 (6.4)	1 (0.8)	8 (6.4)	3 (2.4)	0
Sulfonamide	sul1	24 (19.2)	32 (25.6)	9 (7.2)	52 (41.6)	42 (33.6)	41(32.8)	58 (46.4)	2 (1.6)	46 (36.8)	18 (14.4)	3 (2.4)	0	15 (12)	3 (2.4)	0
Chloramphenicol	cat1	4 (3.2)	4 (3.2)	1 (0.8)	6 (4.8)	6 (4.8)	6 (4.8)	8 (6.4)	1 (0.8)	6 (4.8)	1 (0.8)	1 (0.8)	0	1 (0.8)	0	0
	cmlA	6 (4.8)	7 (5.6)	3 (2.4)	12 (9.6)	8 (6.4)	14 (11.2)	16 (12.8)	1 (0.8)	17 (13.6)	5 (4)	4 (3.2)	0	4 (3.2)	0	0
Tetracycline	tet(A)	26 (20.8)	35 (28)	12 (9.6)	64 (51.2)	48 (38.4)	48 (38.4)	71 (56.8)	5 (4)	52 (41.6)	24 (19.2)	9 (7.2)	1 (0.8)	21 (16.8)	5	0
	tet(B)	4 (3.2)	2 (1.6)	2 (1.6)	7 (5.6)	4 (3.2)	7 (5.6)	10 (8)	1 (0.8)	5 (4)	3 (2.4)	0	0	0	0	0
Trimethoprim	dfrA1	7 (5.6)	10 (8)	1 (0.8)	20 (16)	18(14.4)	15 (12)	19 (15.2)	2 (1.6)	14 (11.2)	9 (7.2)	3 (2.4)	0	5 (4)	1 (0.8)	0
Quinolones	qnr	0	0	0	1 (0.8)	1 (0.8)	0	1 (0.8)	0	1 (0.8)	0	0	0	0	0	0
p-value		0.0004	0.0004	0.0004	0.0004	0.0004	0.0004	0.0004	0.0004	0.0004	0.0004	0.0004	1	0.0004	0.0004	-

Significant p-values (<0.05) are shown in bold.

References

1. Patil, S.; Dong, S.; Francisco, N.M.; Liu, S.; Wen, F. Emergence of multidrug-resistant Klebsiella pneumoniae in hospitalised young children. Lancet Microbe; 2023; 4, 390. [DOI: https://dx.doi.org/10.1016/S2666-5247(23)00094-0]

2. Davoudabadi, S.; Goudarzi, M.; Hashemi, A. Detection of Virulence Factors and Antibiotic Resistance among Klebsiella pneumoniae Isolates from Iran. Biomed. Res. Int.; 2023; 3624497. [DOI: https://dx.doi.org/10.1155/2023/3624497] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/36825037]

3. Kochan, T.J.; Nozick, S.H.; Medernach, R.L.; Cheung, B.H.; Gatesy, S.W.M.; Lebrun-Corbin, M.; Mitra, S.D.; Khalatyan, N.; Krapp, F.; Qi, C. et al. Genomic surveillance for multidrug-resistant or hypervirulent Klebsiella pneumoniae among United States bloodstream isolates. BMC Infect. Dis.; 2022; 22, 603. [DOI: https://dx.doi.org/10.1186/s12879-022-07558-1]

4. Gu, D.; Dong, N.; Zheng, Z.; Lin, D.; Huang, M.; Wang, L.; Chan, E.W.; Shu, L.; Yu, J.; Zhang, R. et al. A fatal outbreak of ST11 carbapenem-resistant hypervirulent Klebsiella pneumoniae in a Chinese hospital: A molecular epidemiological study. Lancet Infect. Dis.; 2018; 18, pp. 37-46. [DOI: https://dx.doi.org/10.1016/S1473-3099(17)30489-9] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28864030]

5. Toledano-Tableros, J.E.; Gayosso-Vázquez, C.; Jarillo-Quijada, M.D.; Fernández-Vázquez, J.L.; Morfin-Otero, R.; Rodríguez-Noriega, E.; Giono-Cerezo, S.; Gutkind, G.; Di Conza, J.; Santos-Preciado, J.I. et al. Dissemination of bla _{NDM- 1} Gene Among Several Klebsiella pneumoniae Sequence Types in Mexico Associated With Horizontal Transfer Mediated by IncF-Like Plasmids. Front. Microbiol.; 2021; 12, 611274. [DOI: https://dx.doi.org/10.3389/fmicb.2021.611274] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33841344]

6. Cortés-Ortíz, I.A.; Juárez-Gómez, J.C.; Cu-Quijano, C.; Flores-Paz, R.; Durán-Manuel, E.M.; Cruz-Cruz, C.; Gutiérrez-Muñoz, V.H.; Sosa-Hernández, O.; Escobar-Escamilla, N.; Bravata-Alcántara, J.C. et al. Klebsiella pneumoniae blaNDM-1 carrying a class 1 integron causing a hospital outbreak in a Mexican attention center. J. Infect. Dev. Ctries; 2021; 15, pp. 657-664. [DOI: https://dx.doi.org/10.3855/jidc.12996]

7. Dehshiri, M.; Khoramrooz, S.S.; Zoladl, M.; Khosravani, S.A.; Parhizgari, N.; Motazedian, M.H.; Jahedi, S.; Sharifi, A. The frequency of Klebsiella pneumonia encoding genes for CTX-M, TEM-1 and SHV-1 extended-spectrum beta lactamases enzymes isolated from urinary tract infection. Ann. Clin. Microbiol. Antimicrob.; 2018; 17, 4. [DOI: https://dx.doi.org/10.1186/s12941-018-0256-y]

8. Troyano-Hernáez, P.; Gutiérrez-Arroyo, A.; Gómez-Gil, R.; Mingorance, J.; Lázaro-Perona, F. Emergence of Klebsiella pneumoniae Harboring the aac(6ʹ)-Ian Amikacin Resistance Gene. Antimicrob. Agents Chemother.; 2018; 62, e01952-18. [DOI: https://dx.doi.org/10.1128/AAC.01952-18]

9. Ferreira, R.L.; da Silva, B.C.M.; Rezende, G.S.; Nakamura-Silva, R.; Pitondo-Silva, A.; Campanini, E.B.; Brito, M.C.A.; da Silva, E.M.L.; Freire, C.C.M.; da Cunha, A.F. et al. High Prevalence of Multidrug-Resistant Klebsiella pneumoniae Harboring Several Virulence and β-Lactamase Encoding Genes in a Brazilian Intensive Care Unit. Front. Microbiol.; 2019; 9, 3198. [DOI: https://dx.doi.org/10.3389/fmicb.2018.03198]

10. Yang, X.; Dong, N.; Chan, E.W.; Zhang, R.; Chen, S. Carbapenem Resistance-Encoding and Virulence-Encoding Conjugative Plasmids in Klebsiella pneumoniae. Trends Microbiol.; 2021; 29, pp. 65-83. [DOI: https://dx.doi.org/10.1016/j.tim.2020.04.012]

11. Hua, Y.; Wang, J.; Huang, M.; Huang, Y.; Zhang, R.; Bu, F.; Yang, B.; Chen, J.; Lin, X.; Hu, X. et al. Outer membrane vesicles-transmitted virulence genes mediate the emergence of new antimicrobial-resistant hypervirulent Klebsiella pneumoniae. Emerg. Microbes Infect.; 2022; 11, pp. 1281-1292. [DOI: https://dx.doi.org/10.1080/22221751.2022.2065935] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35437096]

12. Yang, X.; Sun, Q.; Li, J.; Jiang, Y.; Li, Y.; Lin, J.; Chen, K.; Chan, E.W.; Zhang, R.; Chen, S. Molecular epidemiology of carbapenem-resistant hypervirulent Klebsiella pneumoniae in China. Emerg. Microbes Infect.; 2022; 11, pp. 841-849. [DOI: https://dx.doi.org/10.1080/22221751.2022.2049458]

13. Zhu, J.; Wang, T.; Chen, L.; Du, H. Virulence Factors in Hypervirulent Klebsiella pneumoniae. Front. Microbiol.; 2021; 12, 642484. [DOI: https://dx.doi.org/10.3389/fmicb.2021.642484] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33897652]

14. Namikawa, H.; Oinuma, K.I.; Yamada, K.; Kaneko, Y.; Kakeya, H.; Shuto, T. Predictors of hypervirulent Klebsiella pneumoniae infections: A systematic review and meta-analysis. J. Hosp. Infect.; 2023; 134, pp. 153-160. [DOI: https://dx.doi.org/10.1016/j.jhin.2023.02.005] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/36813165]

15. El Fertas, A.R.; Messai, Y.; Alouache, S.; Bakour, R. Virulence profiles and antibiotic susceptibility patterns of Klebsiella pneumoniae strains isolated from different clinical specimens. Pathol. Biol. Paris; 2013; 61, pp. 209-216. [DOI: https://dx.doi.org/10.1016/j.patbio.2012.10.004]

16. Struve, C.; Bojer, M.; Krogfelt, K.A. Characterization of Klebsiella pneumoniae type 1 fimbriae by detection of phase variation during colonization and infection and impact on virulence. Infect. Immun.; 2008; 76, pp. 4055-4065. [DOI: https://dx.doi.org/10.1128/IAI.00494-08]

17. Paczosa, M.K.; Mecsas, J. Klebsiella pneumoniae: Going on the Offense whit a Strong Defense. Microbiol. Mol. Biol. Rev.; 2016; 80, pp. 629-661. [DOI: https://dx.doi.org/10.1128/MMBR.00078-15]

18. Subashchandrabose, S.; Mobley, H.L.T. Virulence and Fitness Determinants of Uropathogenic Escherichia coli. Microbiol. Spectr.; 2015; 3, [DOI: https://dx.doi.org/10.1128/microbiolspec.UTI-0015-2012]

19. Han, Y.L.; Wen, X.H.; Zhao, W.; Cao, X.S.; Wen, J.X.; Wang, J.R.; Hu, Z.D.; Zheng, W.Q. Epidemiological characteristics and molecular evolution mechanisms of carbapenem-resistant hypervirulent Klebsiella pneumoniae. Front. Microbiol.; 2022; 13, 1003783. [DOI: https://dx.doi.org/10.3389/fmicb.2022.1003783]

20. Ostria-Hernandez, M.L.; Juárez-de la Rosa, K.C.; Arzate-Barbosa, P.; Lara-Hernández, A.; Sakai, F.; Ibarra, J.A.; Castro-Escarpulli, G.; Vidal, J.E. Nosocomial, Multidrug-Resistant Klebsiella pneumoniae Strains Isolated from Mexico City Produce Robust Biofilms on Abiotic Surfaces but Not on Human Lung Cells. Microb. Drug Resist.; 2018; 24, pp. 422-433. [DOI: https://dx.doi.org/10.1089/mdr.2017.0073]

21. Catalán-Nájera, J.C.; Barrios-Camacho, H.; Duran-Bedolla, J.; Sagal-Prado, A.; Hernández-Castro, R.; García-Méndez, J.; Morfín-OteroMorfín-Otero, R.; Velázquez-Larios, M.D.R.; Ortíz-Navarrete, V.; Gutierrez-Xicotencatl, L. et al. Molecular characterization and pathogenicity determination of hypervirulent Klebsiella pneumoniae clinical isolates serotype K2 in Mexico. Diagn. Microbiol. Infect. Dis.; 2019; 94, pp. 316-319. [DOI: https://dx.doi.org/10.1016/j.diagmicrobio.2019.01.013] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30857917]

22. Remya, P.A.; Shanthi, M.; Sekar, U. Characterisation of virulence genes associated with pathogenicity in Klebsiella pneumoniae. Indian J. Med. Microbiol.; 2019; 37, pp. 210-218. [DOI: https://dx.doi.org/10.4103/ijmm.IJMM_19_157] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31745021]

23. Vargas, J.M.; Moreno Mochi, M.P.; Nuñez, J.M.; Cáceres, M.; Mochi, S.; Del Campo Moreno, R.; Jure, M.A. Virulence factors and clinical patterns of multiple-clone hypermucoviscous KPC-2 producing K. pneumoniae. Heliyon; 2019; 5, e01829. [DOI: https://dx.doi.org/10.1016/j.heliyon.2019.e01829] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31286076]

24. Gołębiewska, J.E.; Krawczyk, B.; Wysocka, M.; Ewiak, A.; Komarnicka, J.; Bronk, M.; Rutkowski, B.; Dębska-Ślizień, A. Host and pathogen factors in Klebsiella pneumoniae upper urinary tract infections in renal transplant patients. J. Med. Microbiol.; 2019; 68, pp. 382-394. [DOI: https://dx.doi.org/10.1099/jmm.0.000942] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30747620]

25. Hallal Ferreira Raro, O.; Nordmann, P.; Dominguez Pino, M.; Findlay, J.; Poirel, L. Emergence of Carbapenemase-Producing Hypervirulent Klebsiella pneumoniae in Switzerland. Antimicrob. Agents Chemother.; 2023; 67, e0142422. [DOI: https://dx.doi.org/10.1128/aac.01424-22] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/36853006]

26. Wang, G.; Zhao, G.; Chao, X.; Xie, L.; Wang, H. The Characteristic of Virulence, Biofilm and Antibiotic Resistance of Klebsiella pneumoniae. Int. J. Environ. Res. Public Health; 2020; 17, 6278. [DOI: https://dx.doi.org/10.3390/ijerph17176278]

27. Bautista-Cerón, A.; Monroy-Pérez, E.; García-Cortés, L.R.; Rojas-Jiménez, E.A.; Vaca-Paniagua, F.; Paniagua-Contreras, G.L. Hypervirulence and Multiresistance to Antibiotics in Klebsiella pneumoniae Strains Isolated from Patients with Hospital- and Community-Acquired Infections in a Mexican Medical Center. Microorganisms; 2022; 10, 2043. [DOI: https://dx.doi.org/10.3390/microorganisms10102043]

28. Bakhtiari, R.; Javadi, A.; Aminzadeh, M.; Molaee-Aghaee, E.; Shaffaghat, Z. Association between Presence of RmpA, MrkA and MrkD Genes and Antibiotic Resistance in Clinical Klebsiella pneumoniae Isolates from Hospitals in Tehran, Iran. Iran. J. Public Health; 2021; 50, pp. 1009-1016. [DOI: https://dx.doi.org/10.18502/ijph.v50i5.6118]

29. Braun, V. Iron uptake mechanisms and their regulation in pathogenic bacteria. Int. J. Med. Microbiol.; 2001; 291, pp. 67-79. [DOI: https://dx.doi.org/10.1078/1438-4221-00103]

30. Chen, T.; Dong, G.; Zhang, S.; Zhang, X.; Zhao, Y.; Cao, J.; Zhou, T.; Wu, Q. Effects of iron on the growth, biofilm formation and virulence of Klebsiella pneumoniae causing liver abscess. BMC Microbiol.; 2020; 20, 36. [DOI: https://dx.doi.org/10.1186/s12866-020-01727-5]

31. Bachman, M.A.; Oyler, J.E.; Burns, S.H.; Caza, M.; Lepine, F.; Dozois, C.M.; Weiser, J.N. Klebsiella pneumoniae yersiniabactin promotes respiratory tract infection through evasion of lipocalin 2. Infect. Immun.; 2011; 79, pp. 3309-3316. [DOI: https://dx.doi.org/10.1128/IAI.05114-11] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21576334]

32. Schubert, S.; Picard, B.; Gouriou, S.; Heesemann, J.; Denamur, E. Yersinia high-pathogenicity island contributes to virulence in Escherichia coli causing extraintestinal infections. Infect. Immun.; 2002; 70, pp. 5335-5337. [DOI: https://dx.doi.org/10.1128/IAI.70.9.5335-5337.2002] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/12183596]

33. Kikuchi, S.; Kosai, K.; Ota, K.; Mitsumoto-Kaseida, F.; Sakamoto, K.; Hasegawa, H.; Izumikawa, K.; Mukae, H.; Yanagihara, K. Clinical and microbiological characteristics of bloodstream infection caused by Klebsiella pneumoniae harboring rmpA in Japanese adults. Sci. Rep.; 2023; 13, 6571. [DOI: https://dx.doi.org/10.1038/s41598-023-33265-1]

34. Schulz, E.; Schumann, M.; Schneemann, M.; Dony, V.; Fromm, A.; Nagel, O.; Schulzke, J.D.; Bücker, R. Escherichia coli Alpha-Hemolysin HlyA Induces Host Cell Polarity Changes, Epithelial Barrier Dysfunction and Cell Detachment in Human Colon Carcinoma Caco-2 Cell Model via PTEN-Dependent Dysregulation of Cell Junctions. Toxins; 2021; 13, 520. [DOI: https://dx.doi.org/10.3390/toxins13080520] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34437391]

35. Gallegos-Miranda, V.; Garza-Ramos, U.; Bolado-Martínez, E.; Navarro-Navarro, M.; Félix-Murray, K.R.; Candia-Plata, M.D.C.; Sanchez-Martinez, G.; Dúran-Bedolla, J.; Silva-Sánchez, J. ESBL-producing Escherichia coli and Klebsiella pneumoniae from health-care institutions in Mexico. J. Chemother.; 2021; 33, pp. 122-127. [DOI: https://dx.doi.org/10.1080/1120009X.2020.1861815]

36. Al-Sheboul, S.A.; Al-Madi, G.S.; Brown, B.; Hayajneh, W.A. Prevalence of Extended-Spectrum β-Lactamases in Multidrug-Resistant Klebsiella pneumoniae Isolates in Jordanian Hospitals. J. Epidemiol. Glob. Health; 2023; 13, pp. 180-190. [DOI: https://dx.doi.org/10.1007/s44197-023-00096-2] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/37095370]

37. Dan, B.; Dai, H.; Zhou, D.; Tong, H.; Zhu, M. Relationship Between Drug Resistance Characteristics and Biofilm Formation in Klebsiella pneumoniae Strains. Infect. Drug Resist.; 2023; 16, pp. 985-998. [DOI: https://dx.doi.org/10.2147/IDR.S396609]

38. Córdova-Espinoza, M.G.; Giono-Cerezo, S.; Sierra-Atanacio, E.G.; Escamilla-Gutiérrez, A.; Carrillo-Tapia, E.; Carrillo-Vázquez, L.I.; Mendoza-Pérez, F.; Leyte-Lugo, M.; González-Vázquez, R.; Mayorga-Reyes, L. et al. Isolation and Identification of Multidrug-Resistant Klebsiella pneumoniae Clones from the Hospital Environment. Pathogens; 2023; 12, 634. [DOI: https://dx.doi.org/10.3390/pathogens12050634]

39. Hwang, J.H.; Handigund, M.; Hwang, J.H.; Cho, Y.G.; Kim, D.S.; Lee, J. Clinical Features and Risk Factors Associated With 30-Day Mortality in Patients With Pneumonia Caused by Hypervirulent Klebsiella pneumoniae (hvKP). Ann. Lab. Med.; 2020; 40, pp. 481-487. [DOI: https://dx.doi.org/10.3343/alm.2020.40.6.481]

40. Xu, L.; Sun, X.; Ma, X. Systematic review and meta-analysis of mortality of patients infected with carbapenem-resistant Klebsiella pneumoniae. Ann. Clin. Microbiol. Antimicrob.; 2017; 16, 18. [DOI: https://dx.doi.org/10.1186/s12941-017-0191-3]

41. Rippere-Lampe, K.E.; O’Brien, A.D.; Conran, R.; Lockman, H.A. Mutation of the gene encoding cytotoxic necrotizing factor type 1 (cnf-1) attenuates the virulence of uropathogenic Escherichia coli. Infect. Immun.; 2001; 69, pp. 3954-3964. [DOI: https://dx.doi.org/10.1128/IAI.69.6.3954-3964.2001]

42. Yeh, K.M.; Lin, J.C.; Yin, F.Y.; Fung, C.P.; Hung, H.C.; Siu, L.K.; Chang, F.Y. Revisiting the importance of virulence determinant magA and its surrounding genes in Klebsiella pneumoniae causing pyogenic liver abscesses: Exact role in serotype K1 capsule formation. J. Infect. Dis.; 2010; 201, pp. 1259-1267. [DOI: https://dx.doi.org/10.1086/606010]

43. Müller, C.M.; Aberg, A.; Straseviçiene, J.; Emody, L.; Uhlin, B.E.; Balsalobre, C. Type 1 fimbriae, a colonization factor of uropathogenic Escherichia coli, are controlled by the metabolic sensor CRP- cAMP. PLoS Pathog.; 2009; 5, e1000303. [DOI: https://dx.doi.org/10.1371/journal.ppat.1000303] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19229313]

44. Lopez-Canovas, L.; Martinez Benitez, M.B.; Herrera Isidron, J.A.; Flores Soto, E. Pulsed Field Gel Electrophoresis: Past, present, and future. Anal. Biochem.; 2019; 573, pp. 17-29. [DOI: https://dx.doi.org/10.1016/j.ab.2019.02.020] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30826351]

45. Ooka, T.; Terajima, J.; Kusumoto, M.; Iguchi, A.; Kurokawa, K.; Ogura, Y.; Asadulghani, M.; Nakayama, K.; Murase, K.; Ohnishi, M. et al. Development of a multiplex PCR-based rapid typing method for enterohemorrhagic Escherichia coli O157 strains. J. Clin. Microbiol.; 2009; 47, pp. 2888-2894. [DOI: https://dx.doi.org/10.1128/JCM.00792-09] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19641072]

46. Liu, Y.; Liu, C.; Zheng, W.; Zhang, X.; Yu, J.; Gao, Q.; Hou, Y.; Huang, X. PCR detection of Klebsiella pneumoniae in infant formula based on 16S-23S internal transcribed spacer. Int. J. Food Microbiol.; 2008; 125, pp. 230-235. [DOI: https://dx.doi.org/10.1016/j.ijfoodmicro.2008.03.005] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18579248]

47. CLSI. Performance standards for antimicrobial susceptibility testing. CLSI, Twenty-Third Informational Supplement M100-S23; Carpenter, D.E. Clinical and Laboratory Standards Institute: Wayne, PA, USA, 2013; 74.

48. Dallenne, C.; Da Costa, A.; Decré, D.; Favier, C.; Arlet, G. Development of a set of multiplex PCR assays for the detection of genes encoding important beta-lactamases in Enterobacteriaceae. J. Antimicrob. Chemother.; 2010; 65, pp. 490-495. [DOI: https://dx.doi.org/10.1093/jac/dkp498] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20071363]

49. Momtaz, H.; Karimian, A.; Madani, M.; Safarpoor Dehkordi, F.; Ranjbar, R.; Sarshar, M.; Souod, N. Uropathogenic Escherichia coli in Iran: Serogroup distributions, virulence factors and antimicrobial resistance properties. Ann. Clin. Microbiol. Antimicrob.; 2013; 12, 8. [DOI: https://dx.doi.org/10.1186/1476-0711-12-8]

50. Cheng, G.; Wozniak, K.; Wallig, M.A.; Fidel, P.L., Jr.; Trupin, S.R.; Hoyer, L.L. Comparison between Candida albicans Agglutinin-Like Sequence Gene Expression Patterns in Human Clinical Specimens and Models of Vaginal Candidiasis. Infect. Immun.; 2005; 73, pp. 1656-1663. [DOI: https://dx.doi.org/10.1128/IAI.73.3.1656-1663.2005]

51. Naglik, J.R.; Moyes, D.; Makwana, J.; Kanzaria, P.; Tsichlaki, E.; Weindl, G.; Tappuni, A.R.; Rodgers, C.A.; Woodman, A.J.; Challacombe, S.J. et al. Quantitative expression of the Candida albicans secreted aspartyl proteinase gene family in human oral and vaginal candidiasis. Microbiology; 2008; 154, pp. 3266-3280. [DOI: https://dx.doi.org/10.1099/mic.0.2008/022293-0]

52. Ribot, E.M.; Fair, M.A.; Gautom, R.; Cameron, D.N.; Hunter, S.B.; Swaminathan, B.; Barrett, T.J. Standardization of pulsed- field gel electrophoresis protocols for the subtyping of Escherichia coli O157:H7, Salmonella and Shigella for PulseNet. Foodborne Pathog. Dis.; 2006; 3, pp. 59-67. [DOI: https://dx.doi.org/10.1089/fpd.2006.3.59] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16602980]

53. Tenover, F.C.; Arbeit, R.D.; Goering, R.V.; Mickelsen, P.A.; Murray, B.E.; Persing, D.H.; Swaminathan, B. Interpreting chromosomal DNA restriction patterns produced by pulsed-field gel electrophoresis: Criteria for bacterial strain typing. J. Clin. Microbiol.; 1995; 33, pp. 2233-2239. [DOI: https://dx.doi.org/10.1128/jcm.33.9.2233-2239.1995] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/7494007]