Hypervirulent Klebsiella pneumoniae in a South African tertiary hospital—Clinical profile, genetic determinants, and virulence in Caenorhabditis elegans

1 Introduction

The hypervirulent Klebsiella pneumoniae (hvKp), a causative agent of fulminant and invasive diseases relating to community-acquired pneumoniae (CAP), capable of bearing plasmids of hypervirulence or carbapenem resistance, is recognized as another circulating pathotype in addition to classical K. pneumoniae (cKp), which is a frequent pathogenic agent relating to hospital-acquired pneumoniae (HAP), considered to have limited virulence capability compared to hvKp (Karampatakis et al., 2023). This pathotype is associated with high pathogenicity and mortality and was initially identified in Taiwan in the mid-1980s (Russo and Marr, 2019; Wang et al., 2019). An increasing number of reports across the globe has indicated the geographical spread of hvKp (Shon et al., 2013). Liu et al. (2014) and Nahavandinejad and Asadpour (2017) reported an incidence of 31.4% (22/70) and 33.8% (22/65) of hypervirulent cases in China and Northern Iran, respectively. In South Africa (SA), limited studies on hvKp have been conducted, with only referred reports (collaborative studies conducted abroad that included samples collected in SA). Consequently, little is known about the true prevalence, clinical significance and presentation of hvKp infections.

Previously, hvKp was clinically associated with the ability to cause invasive pyogenic liver abscess (PLA) in a host (Pomakova et al., 2012; Russo and Marr, 2019). However, this pathotype has been clinically reported to cause community-acquired infections (CAI) which spreads to the eyes (endophthalmitis) and brain (meningitis) (Ullmann and Podschun, 1998; Pomakova et al., 2012; Russo et al., 2015; Paczosa and Mecsas, 2016). Unlike cKp infections, hvKp infections usually affect young and healthy individuals (Russo and Marr, 2019; Choby et al., 2020). However, current studies have reported the intrusion of hvKp in the healthcare settings (Li et al., 2018; Lan et al., 2019; Mukherjee et al., 2021). The mortality rate of hvKp infections ranges between 3 and 32% in healthy individuals in the community, suggesting that hvKp is a variant of concern and should be monitored (Han, 1995; Wang et al., 1998; Ko et al., 2002; Fang et al., 2007; Pomakova et al., 2012). In China, the incidence rate of hvKp infections has been reported to be ~74%, and a high mortality rate of 60%, with community transmission being further promoted via colonization of gastrointestinal system (Choby et al., 2020; Su et al., 2021). Several genetic determinants are associated with hypervirulence that distinguish hvKp from cKp (Wu et al., 2017; Lam et al., 2018a,b; Wyres et al., 2020).

The most common feature of hvKp is the increased production of a capsule, mainly capsule K1 and K2 (Ullmann and Podschun, 1998; Remya et al., 2018; Russo et al., 2018). A relationship between the presence of capsule type and serotypes (sequence type) has been reported, which has shed some additional light in the virulence of hvKp (Conlan et al., 2012; Hyun et al., 2019). The sequence type (ST) 23 is believed to be associated with K1 capsule type, while distinct STs (e.g., ST 65, ST 86, ST 375, and ST 380) are associated with K2 (Follador et al., 2016; Catalán-Nájera et al., 2017; Martin and Bachman, 2018). Iron chelators (siderophores), such as aerobactin and salmochelin, have also been reported to be associated with increased virulence in hvKp (Wu et al., 2017; Lam et al., 2018a,b; Wyres et al., 2020). These genetic determinants are located on large virulence plasmids along with the regulator of mucoid phenotype A (rmpA, and rmpA2) that enhances the chromosomal capsule polysaccharide (CPS) (Choby et al., 2020). The result of enhanced expression of CPS is a hyper-capsule, a phenomenon known as a hypermucoviscous phenotype detected by a string test (Russo and Marr, 2019; Choby et al., 2020).

The occurrence of hvKp varies and on average ranged between 12 and 45% in regions of China including Hong Kong, Beijing, Changsha, and Anhui (Zhang et al., 2016; Liu et al., 2018, 2019; Lan et al., 2019). However, in developing countries, this variant is often unidentified. Therefore, this study aimed to characterize K. pneumoniae isolates through clinical, molecular, and genomic studies to provide a broader knowledge about the newly emerging hvKp, and whether it is a public concern in developing countries similar to what is seen globally.

2 Materials and methods 2.1 Study design and location

This was a descriptive, cross-sectional study using isolates previously identified as K. pneumoniae routinely using the VITEK® 2 system (bioMérieux, France) by the National Health Laboratory Service (NHLS), Pelonomi Tertiary Hospital, in Bloemfontein, Free State Province, South Africa. Over a period of 12 months (September 2020 to October 2021), a total of 74 K. pneumoniae isolates were obtained from the laboratory. For patient confidentiality, all patient names were anonymized and a study number was assigned to each isolate. The selected K. pneumoniae isolates were from different wards in the hospital, including the medical ward, maternity, casualty, multidisciplinary intensive care unit (MICU), trauma unit, outpatient department (OPD) and neonatal unit. Furthermore, these isolates were obtained from different sample types including blood cultures, tracheal aspirate, urine, and sputum.

2.2 Patient demographics and antimicrobial susceptibility profiles

The laboratory database, NHLS TrakCare Lab WebView system, was used to extract antibiotic susceptibility profiles for each isolate. Patients demographic information and clinical data, including type of sample taken, age, gender, and the day the sample was taken, were also obtained from the NHLS TrakCare WebView system. The date of admission for each patient and clinical presentation were extracted from the Meditech Medical Record System. Hospital-acquired infections were defined as K. pneumoniae infection acquired by patients 48 h after hospitalization, while community-acquired infections were defined as infection occurring 48 h prior to hospitalization. Mild symptoms were treated on an outpatient basis. Moderate to severe disease was defined as patients presenting with organ impairment due to K. pneumoniae infection were patients who were admitted to the intensive care unit (ICU), while mild infections were patients who were admitted to normal wards or outpatients.

2.3 String test for detection of hypermucoviscous phenotypes

All the K. pneumoniae isolates obtained (n = 74) were sub-cultured onto 5% blood agar and incubated at 37°C for 24 h. A string test was performed on all isolates as described by Guo et al. (2017). A colony length of ≥5 mm was defined as positive for the hypermucoviscous (hmv) phenotype.

2.4 Molecular detection of Klebsiella pneumoniae strains

Crude DNA was extracted from all 74 K. pneumoniae isolates. Briefly, samples were heated at 95°C for 30 min and frozen at −80°C for 30 min. The samples were then centrifuged at 13 000 rpm for 10 min and the supernatant was used after centrifugation. A total of six separate conventional single-plex PCR reactions were performed on all isolates. The target sites iuc transporter (iutA), rmpA, capsule K1/K2 (magA/k2A, respectively), yersiniabactin (irp2) and iroB siderophores were used (Table 1).

Table 1. The list of primers used for detection of K1/K2 (magA/k2A) serotypes, rmpA, iutA, iroB, and irp2.

The PCR amplicons were separated using 1.5% SeaKem™ LE Agarose (Lonza™, USA) and the nucleic acid stain gelRed (Biotum GelRed™, Australia), as per the manufacturer's instructions, for 40 min at 100 volts (V). The gel was viewed using the Geldoc™ EZ system (Bio-Rad, USA) for the detection of specific band sizes. Two reference molecular markers were used, including 200 and 500 + 100 bp (Thermo Fisher Scientific, USA). The hvKp was defined as isolates bearing the iuc transport genetic determinant, iutA (Russo et al., 2018).

2.5 In vivo Caenorhabditis elegans killing assay

Seven Klebsiella pneumoniae (K. pneumoniae) representatives of different spectra of hypervirulent Klebsiella pneumoniae (hvKp) genetic determinants and classical Klebsiella pneumoniae (cKp) were selected for in vivo assay (shown in Figure 3A). Summarily, the selection criteria for this assay included two (2) cKp that possessed irp2 (Isolate 54) and K2 (Isolate 60) gene, respectively. Five hvKp that possessed all five genes with K1/K2 (Isolate 1, 5, and 17, respectively), isolate 24 (possess four virulent genes/K1), isolate 51 (had only irp2 gene). The Caenorhabditis elegans (C. elegans) glp-4; sek-1 hermaphrodites used in this study was obtained from the Caenorhabditis Genetic Center, College of Biological Sciences, University of Minnesota, USA. Escherichia coli OP50 was used as a food source for the nematodes (Brenner, 1974). The C. elegans were infected with the seven selected strains of K. pneumoniae according to a modified protocol described by Kamaladevi and Balamurugan (2016). Briefly, after moving Synchronized L4 C. elegans to brain heart infusion (BHI) plates (EMD Millipore, Germany) seeded with selected strains of K. pneumoniae, the plates were incubated at 25°C for 4 h, followed by washing of C. elegans with 9 ml of M9 buffer (Sigma-Aldrich, USA). The plates were incubated at 25°C and the nematodes were scored daily as alive or dead. The C. elegans were regarded as dead when no movement in response to prickling mechanical stimulation was observed and were removed from the media and discarded. The experiment was done in triplicate.

2.6 DNA library preparation, sequencing, editing, genome assembly and phylogeny analysis

The extracted DNA of five K. pneumoniae was sent for sequencing at the Agricultural Research Council (ARC), Pretoria, South Africa. The selected isolates were representatives of both hvKp (n = 4), and cKp (n = 1) and their selection was based on the C. elegans killing assay results. To construct a DNA library for the sequencing of the whole genome of K. pneumoniae, a Nextera® XT library preparation kit (Illumina®, San Diego, CA, USA) was used to construct a DNA library adhering to the manufacturer's procedures. This involved the fragmentation of DNA followed by the addition of dual barcodes to the DNA fragments. For the purification of the barcoded libraries, Agencourt AMPure magnetic beads (Beckman Coulter, Indianapolis, Indiana, USA) were utilized and the selection of an average insert of 300 bp (range 200–400 bp) was carried out simultaneously. In addition, validation of the library and quantification were conducted prior to sequencing using a 2100 Bioanalyzer platform (Agilent Technologies, Santa Clara, CA, USA) and the Qubit™ 3.0 fluorometer (Invitrogen, Carlsbad, CA, USA), respectively. This was followed by the pooling of the validated and quantified libraries, and the whole genome sequencing was executed, using a 5% PhiX DNA control spike-in, on an Illumina® MiSeq platform with using the Illumina MiSeq platform using a v3 standard sequencing kit (Illumina®, San Diego, CA, USA) for 500 cycles.

The raw sequence files (fastq format) were retrieved from Illumina Miseq instrument.

Quality control of sequences was performed using FastQC (v.0.12.1; https://www.bioinformatics.babraham.ac.uk/projects/fastqc/), whereafter trimming and removal of barcodes was achieved with BBDuk which is a component of the BBTools package (https://sourceforge.net/projects/bbmap/). Sequence assembly was performed with Spades (v.3.15.5; Prjibelski et al., 2020). The assembled genomes were assessed for completeness using QUAST (https://github.com/ablab/quast). An additional analysis of genome completeness was performed using BUSCO (v.5.4.7; https://busco.ezlab.org/). The identification of putative virulence genes was achieved through a BLAST comparison of assembled contigs against a database of representative sequences for all virulence genes as identified from the Institut Pasteur Klebsiella pneumoniae virulence genes scheme (https://bigsdb.pasteur.fr/cgi-bin/bigsdb/bigsdb.pl?db=pubmlst_klebsiella_seqdef&page=downloadAlleles). Phylogeny analysis was conducted using a 1,000 replicate bootstrapping and the distance-based neighbor-joining algorithm method as implemented in Mega 11 software was adopted to compute the evolutionary distances of the aligned sequences (Tamura et al., 2021). Annotation and plasmid visualization of the plasmid were conducted using SnapGene® software version 7.1.1 (available at www.snapgene.com).

2.7 Statistical analysis

Genetic determinants of these isolates were analyzed using logistic regression in MS Excel (version 2016) to identify variables associated with hvKp. The presence of iutA was used as an independent variable for defining hvKp. Statistical significance between the virulence levels of distinct strains (hvKp, cKp and E. coli OP50) was analyzed using the online application for survival analysis version 2 (OASIS 2) as described by Han et al. (2016). All experiments performed in the C. elegans killing assay were done in triplicate and were expressed as a mean with standard deviation. A p-value of < 0.05 was considered statistically significant using logistic regression in MS Excel.

2.8 Ethical considerations

Ethics approval to conduct the research was obtained from the Health Sciences Research Ethics Committee (HSREC) of the University of the Free State in Bloemfontein, South Africa (UFS-HSD2020/1579/2302) and Environmental and Biosafety Research Ethics Committee (UFS-ESD2020/0148). Further approval from the Free State Department of Health was obtained.

3 Results 3.1 Patient demographics and genotypic biomarker

Hypervirulence was found in nine (12.2%) of the 74 K. pneumoniae isolates (Supplementary Figure 1). The average patient age in hypervirulent Klebsiella pneumoniae (hvKp) strains was 35 years (19–57 years) while cKp was 34 years (Newborn–85).

The Clinical presentations and antimicrobial profiles of hvKp are summarized in Tables 2, 3 respectively. The antimicrobial susceptibility profile of cKp isolates is shown in the supplementary Table 2. The overall antimicrobial resistance of all isolates (n = 74) indicated extended-spectrum beta-lactamases (ESBLs) of (19%) whereby none of the hvKp were ESBL. Regarding Carbapenem-resistant Enterobacteriaceae (CRE), the overall incidence was 4%, none of which were hvKp.

Table 2. Clinical presentation of patients with hypervirulent Klebsiella pneumoniae (hvKp) infections.

Table 3. The antimicrobial susceptibility profile of hypervirulent Klebsiella pneumoniae isolates.

Hypermucoviscous phenotype (hmv) was found in 7/9 (78%) hvKp strains. An additional (1/64) classical Klebsiella pneumoniae (cKp) strain also possessed this phenotype. However, there was no statistical significance associated with the severity of the disease and Klebsiella pathotypes (Figure 1).

Figure 1. The graphical representation of the spectrum of Klebsiella pneumoniae disease amongst the studied isolates.

3.2 Virulence-associated factors of hvKp by PCR

The nine patients with hvKp in Pelonomi Tertiary Hospital varied regarding admission, as two patients came from the main medical ward and another two from the maternity ward. Other hvKp cases were from casualty, the multidisciplinary intensive care unit (MICU), the antenatal clinic and the out-patient department (OPD; n = 1 each). The presence of five genotypic traits [iroB, rmpA, K1 (magA), and K2 (k2A)] and a single phenotypic characteristic (hmv phenotype) were statistically associated p < 0.05) with hvKp (Table 4). Moreover, the virulence genes rmpA, K1 (magA), K2 (k2A) and iroB were detected in 88.9% (n = 8), 44.4% (n = 49), 55.6% (n = 5) and 88.9% (n = 8) hvKp strains, respectively. The presence of the irp2 gene was not significantly associated with hvKp (p = 0.17) as 19 (29.2%) of the 65 cKp isolates possessed this gene. There was no statistically significant association between hvKp/cKp (p < 0.05) and the severity of disease. However, in 44.4% (n = 4/9) of infections caused by hvKp isolates, the patients had a moderate to severe spectrum of disease, whereas most patients (n = 40/65; 61.5%) with cKp infections had mild symptoms.

Table 4. Comparison of phenotypic and virulence-associated genes between hypervirulent Klebsiella pneumoniae (hvKp) and classical K. pneumoniae (cKp).

3.3 Klebsiella pneumoniae virulence in Caenorhabditis elegans

Live worms moved freely and responded to the prickling mechanical stimulation (movement when touched; Figure 2A), while dead C. elegans were rigid, straight, and floated on the media (Figure 2B). All seven selected K. pneumoniae strains (Figure 3B) showed a statistically significant (p < 0.05) increase in virulence compared to C. elegans exposed to E. coli OP50 only. The Two hvKp strains expressing all five genetic virulence markers showed enhanced virulence (p < 0.05), as none of the worms survived between day 9 and day 10, as shown in Figure 3B. In contrast, hvKp isolate 17 which had all five genetic determinants and belonged to K2 capsular type, had a reduced virulence compared to isolate 1 which had enhanced virulence (presence of five genetic determinants) and belonged to K1. There was no significant difference (p = 0.0795) in the virulence observed in a single iutA-positive hvKp (isolate 51) that also possessed the irp2 gene.

Figure 2. Light visualization of Caenorhabditis elegans under stereomicroscope (×1.28 magnification). (A) Live C. elegans that responded to prickling mechanical stimulation. (B) Dead, straight and rigid C. elegans that did not respond to prickling mechanical stimulation.

Figure 3. In vivo Caenorhabditis elegans infection model comparing hvKp and cKp. (A) Virulence genes present (green) or absent (light blue) in the strains selected for the infection model. Green: virulence determinant present; blue: virulence determinant absent. (B) The survival rate of C. elegans fed with different hvKp or cKp strains.

3.4 Klebsiella pneumoniae sequence types and K-serotype

The population structure of the isolates was established by determining multi-locus sequence typing (MLST) and the K-serotype grouping (Larsen et al., 2012; Wyres et al., 2016). The sequences were identified for isolate 5, 7, 24, 51, and 54 as ST20, ST23, ST65, ST985, and ST3430, which is associated with the K loci groups KL1, KL2, KL28, KL39, and KL52, respectively, having a very high match confidence of ≥99% confidence to the reference locus, with no expected genes missing and no unexpected locus genes.

3.5 Genome sequence analysis

The analysis confirmed all isolates to be K. pneumoniae. The genome sequence revealed distinct virulence genes including siderophores (such as aerobactin, salmochelin, yersiniabactin, and enterobactin), iron uptake, allantoin metabolism, and colibactin with an identity range of 95%−100% (as shown in Supplementary Table 1). The prototype (Figure 4) shows K. pneumoniae plasmid (125,269 bp). The plasmid encodes rmpA, rmpA2, iucABCD, iutA, and iroBCDN amongst others.

Figure 4. The complete sequence of K.pneumoniae plasmid (125,269 bp) showing the position of the hvKp-associated genes (iucABCDiutA, iroBCDN, rmpA, and rmpA2) adopted for the study (left). The region of the iutA adopted as a biomarker in this study is amplified (right), and corresponding Open Reading Frames (ORFs) are indicated on the prototype isolate. This visualization was generated using SnapGene software.

The phylogenetic analysis of the curated sequences of K. pneumoniae further authenticated the virulence-associated factor of hvKp with >80% bootstrap values in the genes analyzed and their clustering patterns as shown in Figure 5, Sample 5 (hvKp) was found to cluster in the same clade with reference strain (Accession number: CP1392441) from the United States, having ~99% similarity, as well as reference strain (Accession number: CP071163) from West China. Also, sample 17 (hvKp) clustered unambiguously with other reference sequences; accession numbers: CP139682 and CP140295 from China and the United Kingdom respectively, while sample 24 was found to cluster with reference sequences CP137400 and CP137360 with 99% identity. Lastly, sample 54 (cKp) clustered unambiguously with other reference sequences (Accession number: CP124750, CP074539, OW969925, and CP132634) from NCBI with 98% similarity in one clade (Figure 5). The curated sequences from this study have been submitted to the GenBank and assigned with the following accession numbers: PP296963, PP296964, PP296965, and PP296966.

Figure 5. The evolutionary history was inferred using the Neighbor-Joining method (Saitou and Nei, 1987). The percentage of replicate trees in which the associated K.pneumoniae clustered together in the bootstrap test (1,000 replicates) are shown next to the branches (Felsenstein, 1985). The evolutionary distances were computed using the Maximum Composite Likelihood method and are in the units of the number of base substitutions per site. Evolutionary analyses were conducted in MEGA11 (Tamura et al., 2021).

4 Discussion

This study aimed to investigate the occurrence of hvKp in Pelonomi Tertiary Hospital, one of the largest academic hospitals in the Free State Province. In addition, the aim was to assess the virulence of hvKp vs. cKp isolates and determine if any association between hypervirulent genetic determinates and virulence could be established. In total, 12.2% (n = 9/74) of the isolates collected over the 1-year period were hvKp. Based on available data, the majority of infections caused by these hvKp strains were community-acquired. These findings correlate with the literature which has documented a relationship between hypervirulence and community-acquired infections (CAIs) (Fang et al., 2007; Shankar et al., 2018; Wyres et al., 2020).

Several studies have described the importance of primary pyogenic liver abscess (PLA) or secondary infections (endophthalmitis and meningitis) as invasive infections in community-acquired hvKp patients (Fang et al., 2007; Guo et al., 2017; Li et al., 2021). This includes a study by Liu et al. (2014), who reported that hvKp strains found in the community caused severe infections such as liver abscess and pneumonia. Other studies reported that no statistically significant difference (p > 0.05) had been observed between hvKp and cKp regarding the severity of disease (Yan et al., 2016). Our findings support the latter, as no statistically significant association between hvKp/cKp and severity of disease was found. However, four of the nine hvKp patients had moderate to severe infections associated with organ impairment. In contrast, most patients with cKp infections had mild symptoms and were treated as outpatients. All hvKp isolates from this study were sensitive to antimicrobial drugs, except intrinsic resistance to penicillin and none (0/9) were extended-spectrum beta-lactamase (ESBL) producers.

The phenotypic, and genotypic differences between cKp and hvKp regarding virulence-associated genes are summarized in Table 4. The string test based on hypermucoviscous (hmv) phenotype (string ≥5 mm) has widely been used as a marker for hvKp, with ~90% predicted accuracy for clinical hvKp strains (Russo et al., 2018). However, this semi-qualitative method is easily influenced by colony conditions and the user's technique (Tan et al., 2014). The hmv phenotype from this study was found in 7/9 (78%) hvKp strains. An additional cKp (1/64) strain also possessed this phenotype. Similar results were reported by Zhang et al. (2016) who reported an average 75% (65/87) of hmv phenotype. Moreover, six hmv phenotypes were further detected in cKp strains and suggested that the hmv phenotype is indeed not a suitable biomarker for hypervirulence as some cKp strains also possesses this phenotype (Zhang et al., 2016).

The use of iutA as a biomarker in this current study confirmed that iutA is a reliable indication of hypervirulence (Hsieh et al., 2008; Hsu et al., 2011; Russo et al., 2014; Lam et al., 2018a,b). This is supported by current findings that the iutA gene was only detected in K. pneumoniae isolates that also possessed the virulent determinants rmpA, iroB and rmpA2 (Russo et al., 2014; Russo and Marr, 2019). These three genetic determinants along with iutA are located on large virulence plasmids that cKp do not possess (Russo et al., 2014; Nahavandinejad and Asadpour, 2017; Russo and Marr, 2019; Choby et al., 2020). The results from the current study further showed that these genetic determinants are located on large virulence plasmids, and are associated with hypervirulence as none of these genetic determinants were found in isolate 54 (cKp) in this study, this is supported by the findings in a recent study conducted on hvKp from hospital and community settings that showed the presence of iutA, iucB, iucC, and iutA in all hvKp isolates (Raj et al., 2022). Recent studies have reported that cKp strains can acquire hvKp virulence-like plasmids that encode virulence-associated determinants (Gu et al., 2018; Huang et al., 2018; Russo and Marr, 2019). Therefore, more studies are needed to understand in depth the acquired virulence by cKp strains, e.g., through plasmids, and whether they result in enhanced virulence in the host.

Using the C. elegans nematode model, hvKp isolates were more virulent than cKp, in particularly isolate 1 (that contained all five hypervirulent genetic determinants including K1) and isolate 5 (which had five genes as well as K1 and K2). These findings are in keeping with Pomakova et al. (2012), who noted the presence of genetic determinants (K1/K2, iroB and rmpA) in hvKp strains (defined by positive hmv phenotype). Pomakova et al. (2012) further reported increased virulence in a rat abscess model when all hvKp genetic determinants were present. Other studies found that despite an isolate being hmv-negative, it still exhibited high virulence when possessing rmpA and iucABCDiutA genes (Li et al., 2019).

The presence of the hmv phenotype and irp2 gene did not have any significant effect in virulence of K. pneumoniae strains. Isolate 24 which lacked the hmv phenotype was more virulent than isolate 17, which had the hmv phenotype (Figure 3B), while isolate 54 that possessed only the irp2 was not significantly associated with a higher rate of killing C. elegans, as >60% of the worms survived by day 10. These results are in line with the literature, as various authors reported similar findings that these two factors could be found in cKp strains and therefore are not suitable markers for hypervirulence and not statistically significantly virulent in in vivo models (Wu et al., 2017; Lam et al., 2018a,b; Wyres et al., 2020).

A significant difference between isolate 17 (belonging to capsule type K2) and isolate 1 (belonging to K1) was noted, as isolate 1 was more virulent than isolate 17. Isolate 5 (belonged to K1) was also more virulent than isolate 17 (belonged to K2). Studies have shed some additional light on the virulence of these two capsule types. It is believed that the association of hvKp with K1/K2 serotypes depends on the presence of distinct sequence types (STs) that are conserved or diverse between the two serotypes (Russo et al., 2014; Lam et al., 2018a,b). The ST 23 is associated with the K1 capsule type, while distinct STs are associated with K2 (e.g., ST 65, ST 86, ST 375 and ST 380) (Follador et al., 2016; Catalán-Nájera et al., 2017; Martin and Bachman, 2018). Findings from multi-locus sequence typing (MLST) in this study support the above as it revealed that capsule K1 indeed belonged to ST23, while K2 belonged to ST 65. The capsule K1 serotype is more virulent compared to K2 (Russo and Marr, 2019). The C. elegans killing assay from the current study showed that isolate 5 (ST 23, K1) was more virulent than isolate 17 (ST 65, K2). Furthermore, a gene encoding an alginate lyase isozyme was exclusively found in strains of the K1 serotype, which was associated with enhanced virulence and co-infection in the host (Nahavandinejad and Asadpour, 2017). More genomic sequencing is still required to understand the genetic determinants and the association between the capsule type and invasive infections.

The phylogenetic analysis from this study showed a high percentage of relatedness (89%−99%) with other reference sequences from global endemic countries such as West China, United Kingdom, and the United States (Figure 5). Also, the clustering pattern of the cKp and hvKp from this study showed a 100% similarity of relatedness based on the evolutionary tree (Figure 5), owing to the indisputable ubiquitous nature and ability of contagious pathogens to cross populations.

This study had limitations, including a small sample size. Secondly, the clinical history of some patients was incomplete we were unable to establish the underlying conditions of patients.

In conclusion, to the best of our knowledge, this is the first study to confirm the presence of hvKp in one of the largest hospitals in an African developing country. The results from this study suggest that hvKp strains are more virulent relative to cKp in in vivo models which we considered as the strength of this study amongst others. Enhanced virulence is associated with the presence of virulence determinants, particularly K1/K2, iroB and rmpA. Currently, the occurrence of hvKp in the Free State Province is not a major concern, although further studies in different settings are recommended to deduce whether hypervirulent Klebsiella pneumoniae is a public concern as it is globally.

Data availability statement

The data presented in the study are deposited in the GenBank repository, accession number: PP296963, PP296964, PP296965, and PP296966.

Ethics statement

The studies involving humans were approved by Health Sciences Research Ethics Committee, University of the Free State. The studies were conducted in accordance with the local legislation and institutional requirements. The human samples used in this study were acquired from the National Health Laboratory Service (NHLS), Pelonomi Tertiary Hospital, in Bloemfontein, Free State Province, South Africa. Written informed consent for participation was not required from the participants or the participants' legal guardians/next of kin in accordance with the national legislation and institutional requirements. The manuscript presents research on animals that do not require ethical approval for their study.

Author contributions

LD: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Visualization, Writing—original draft, Project administration, Writing—review & editing. OA: Conceptualization, Formal analysis, Visualization, Writing—review & editing. EA-C: Conceptualization, Formal analysis, Investigation, Writing—review & editing. CP: Formal analysis, Investigation, Writing—review & editing. NM: Formal analysis, Investigation, Writing—review & editing. MD: Conceptualization, Formal analysis, Writing—review & editing. JM: Conceptualization, Investigation, Supervision, Writing—review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was funded by National Research Foundation (NRF) reference number: MND190620449605 and National Health Laboratory Services (NHLS) grant number 94012.

Acknowledgments

We would like to acknowledge the Caenorhabditis Genetics Center (CGC), which was funded by National Institutes of Health (NIH)—Office of Research Infrastructure Programs (P40OD010440) and located in Minnesota, USA, for providing us with Caenorhabditis elegans specimens. We further extend our gratitude to the Free State Province Department of Health for the approval of this study, and National Health Laboratory Services (NHLS).

Conflict of interest

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

Publisher's note

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

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2024.1385724/full#supplementary-material

References

Catalán-Nájera, J. C., Garza-Ramos, U., and Barrios-Camacho, H. (2017). Hypervirulence and hypermucoviscosity: two different but complementary Klebsiella spp. phenotypes? Virulence 7, 1111–1123. doi: 10.1080/21505594.2017.1317412

PubMed Abstract | Crossref Full Text | Google Scholar

Choby, J. E., Howard-Anderson, J., and Weiss, D. S. (2020). Hypervirulent Klebsiella pneumoniae – clinical and molecular perspectives. J. Intern. Med. 287, 283–300. doi: 10.1111/joim.13007

PubMed Abstract | Crossref Full Text | Google Scholar

Compain, F., Babosan, A., Brisse, S., Genel, N., Audo, J., Ailloud, F., et al. (2014). Multiplex PCR for detection of seven virulence factors and K1/K2 capsular serotypes of Klebsiella pneumoniae. J. Clin. Microbiol. 52, 4377–4380. doi: 10.1128/JCM.02316-14

PubMed Abstract | Crossref Full Text | Google Scholar

Conlan, S., Kong, H. H., and Segre, J. A. (2012). Species-level analysis of DNA sequence data from the NIH Human Microbiome Project. PLoS ONE 7:e47075. doi: 10.1371/journal.pone.0047075

PubMed Abstract | Crossref Full Text | Google Scholar

Fang, C. T., Lai, S. Y., Yi, W. C., Hsueh, P. R., Liu, K. L., Chang, S. C., et al. (2007). Klebsiella pneumoniae genotype K1: an emerging pathogen that causes septic ocular or central nervous system complications from pyogenic liver abscess. Clin. Infect. Dis. 45, 284–293. doi: 10.1086/519262