Salmonella Typhi Haplotype 58 (H58) Biofilm Formation and Genetic Variation in Typhoid Fever Patients with Gallstones in an Endemic Setting in Kenya

Genotype Identification and Clustering Tree

Since the S. Typhi population is highly structured, with dozens of subclades being associated with specific geographical regions and antimicrobial resistance patterns, the genotypes causing asymptomatic carriage in the current study settings were identified. All 22 bacteria isolates were genotype 4.3.1 (S. Typhi Haplotype 58 [H58]). Within this genotype, 11 of the isolates (from 2/4 of the households) fell under the lineage 4.3.1.1 (from households A and B) while the remaining 11 isolates (from households C and D) belong to lineage 4.3.1.2 (Figure 1). Lineage 4.3.1.2 strains were further grouped into two sub-lineages, 4.3.1.2EA2 (4 isolates from the index case in household C) and 4.3.1.2EA3 (7 isolates from an index case in household D). Each patient shed typhoid bacilli belonging to only one lineage/sub-lineage. From household A, isolates (i)-(v) were from and index case while isolate (vi) was from a household contact. Isolates (i) and (ii), household B, were from index case while (iii), (iv) and (v) were from an asymptomatic household contact living with the index case. All the S. Typhi strains were isolated from stool samples apart from two, isolate (i) from household B, and isolate (i) from household D, which were isolated from blood samples. The time of isolation/shedding of each isolate is shown in Table 1. Assembled genomes can be accessed in NCBI database, BioProject ID PRJNA1101423, GenBank accession numbers are shown in the data availability section.

• Download figure
• Open in new tab

Figure 1.

Clustering tree for S. Typhi strains isolated from the study participants living in the different households as generated using the pathogenwatch database and visualized using microreact. Roman numbers indicate specific isolates from the different households (A, B, C and D). Isolates from each household are shown by specific leaf nodes color. Branch tips are color-coded according to the household of isolation. The bar represents branch length scale bar. The tree can be visualized here https://microreact.org/project/s3biz8PE9LDnFayjco22MZ-s-typhi-kenya-2023.

Antimicrobial Resistance Genes

Different antimicrobial resistance patterns were observed in the isolated S. Typhi strains. The seven isolates belonging to sub-lineage 4.3.1.2EA3 (from Household D) were multidrug resistant, all expressed the following acquired antimicrobial resistance genes; sul1, dfrA7, catA1, aph(6)-Id, aph(3’’)-Ib, sul2 and bla_TEM-1, and a point mutation in the Quinolone Resistance Determining Region (QRDR) of gyrA (gyrA S83Y). Phenotypic susceptibility data showed that these seven isolates were resistant to ampicillin, chloramphenicol, trimethoprim-sulfamethoxazole and nalidixic acid but non-susceptible to ciprofloxacin. The four sub-lineage 4.3.1.2EA2 isolates (S. Typhi strains from Household C) had a gyrB S464F mutation in the QRDR, and all were non-susceptible to nalidixic acid and ciprofloxacin according to phenotypic susceptibility results. The third group, lineage 4.3.1.1, had six strains (Household A isolates) with a gyrB S464F mutation, also demonstrating non-susceptibility to nalidixic acid and ciprofloxacin. The other five lineage 4.3.1.1 strains (household B isolates) had gyrA S83F mutations in the QRDR, and showed resistance to nalidixic acid and non-susceptibility to ciprofloxacin (Table 2).

Variant Calling

When comparing genome sequences among the isolates, missense mutations were observed in a gene that codes for Trehalose Phosphotransferase System (PTS), trehalose transporter subunit IIBC, treB gene (treB A383T), in 3/7 sub-lineage 4.3.1.2EA3 (Household D) isolates and in tviE gene (tviE P263S), which codes for the gene coding for biosynthesis of Vi polysaccharide biosynthesis protein TviE, in 4/7 isolates in the same group, while 2/7 had both mutations (treB A383T and tviE P263S), (Table 3). A nonsense mutation was observed in 1/7 of the sub-lineage 4.3.1.2EA3 isolates in the BCKKFA_12280 locus (E237*) that codes for integrase/transposase family protein. The first follow-up isolate (S. Typhi sub-lineage 4.3.1.2EA2) from the index case living in household C (C ii) had a total of 16 mutations, of which 7 were silent, 8 missense mutations and a deletion in mutL gene that codes for DNA mismatch repair endonuclease MutL (Table 3 and Supplementary files 2 and 3). The second and third follow-up isolates from the same patient (C iii and C iv) had two silent mutations, one in the yccC gene that codes for putative membrane protein YccC and tehA gene that codes for dicarboxylate transporter/tellurite-resistance protein TehA (tehA A124A and yccC R199R mutations respectively) and missense mutations in the amiA gene that codes for N-acetylmuramoyl-L-alanine amidase AmiA (amiA V145A mutation). An additional mutation was observed in the LEJNAJ_18700 locus (K124E) that encodes for 4-hydroxyphenylacetate permease. The second follow-up isolate from the patient in household C (C iii) had a third missense mutation in the LEJNAJ_10515 locus (D78G) that encodes phage baseplate assembly protein V. The follow-up lineage 4.3.1.1 strains, which were isolated from households A and B, also had a few mutations. While no mutation was detected in the first follow-up sample from the index case in household A, the second follow-up sample (A iii) had a single nucleotide polymorphism in the crl gene (crl L38P), crl codes for sigma factor-binding protein Crl. The third and fourth follow-up S. Typhi isolates from the same patient (A iv and A v) had a silent mutation in the tnpA gene (tnpA Y41Y) that codes for IS200/IS605 family transposase (Supplementary file 1). From the same household, an S. Typhi strain isolated from a household contact’s stool sample (A vi) had a missense mutation in the treB gene (treB A383T). S. Typhi bacteria isolated from the stool of the index case living in household B before treatment (isolate ii), had a missense mutation in the waaK gene (waaK P167L) that codes for lipopolysaccharide N-acetylglucosaminyltransferase, compared to the strain from blood sample collected on the same day (Supplementary file 2).

Plasmid Identification

S. Typhi genotype 4.3.1 lineages/sub-lineages identified in this study were found to contain plasmids. A total of 22 different plasmids were identified, 12/22 of these had sequences corresponding to known Salmonella species plasmids, 6/22 had sequences corresponding to Escherichia coli plasmids, 2/22 corresponding to Klebsiella pneumoniae plasmids, 1/22 corresponding to Enterobacter kobei and 1/22 an Escherichia albertii plasmid (Table 4). Out of the 22 plasmids, 8/22 were detected in all 22 S. Typhi strains. The eight plasmids include Salmonella Infantis strain CFSAN003307 plasmid pCFSAN003307, Salmonella Senftenberg strain NCTC10384 plasmid 4, E. coli strain 13TMH22 plasmid p13TMH22-1, Salmonella enterica subsp. enterica serovar 1,4,[5],12:i:-strain PNCS014880 plasmid p16-6773.2, S. Typhi strain ERL11909 plasmid 3, Salmonella Typhimurium strain AUSMDU00027951 isolate AUSMDU00027951 plasmid P02, Salmonella Weltevreden plasmid pSH17G0407 and S. Senftenberg strain NCTC10384 plasmid 3. Plasmids with sequences corresponding to E. coli strain 88COLEC plasmid p88COLEC-1 and E. coli strain 13P477T plasmid p13P477T-2 were only detected in lineage 4.3.1.2EA2 (Household C) isolates, in 4/4 and 2/4 of the strains respectively. The E. albertii strain 2014C-4356 plasmid unnamed6 was only found in 1/11 of lineage 4.3.1.1 isolates while Salmonella Uganda strain CFSAN006208 plasmid pCFSAN006208_1 was detected in 1/11 of lineage 4.3.1.1 isolates. Sequences corresponding to five known plasmids, Salmonella Enteritidis strain 81-1705 plasmid pSE81-1705-3, S. Typhi strain 311189_252186 plasmid pHCM1, E. coli NES6 plasmid pN-ES-6-1, K. pneumoniae strain RIVM_C018860 plasmid pRIVM_C018860_2 and K. pneumoniae strain 03108465-40B plasmid pKpQIL, were detected in 7/7 of sub-group 4.3.1.2EA3 strains but not in the other sub-groups.

Genes conferring multidrug resistance in S. Typhi were detected in two different plasmids identified in sub-group 4.3.1.2EA3 (Household D) strains. The antimicrobial resistance genes sul1, dfrA, catA1, aph(6)-Id, aph(3’’)-Ib and sul2 were detected in a plasmid with sequences corresponding to S. Typhi strain 311189_252186 plasmid pHCM1, while bla_TEM-1 was detected in a plasmid with sequences corresponding to E. coli NES6 plasmid pN-ES-6-1. None of the plasmids in isolates from households A, B and C had identifiable antimicrobial resistance genes.

Biofilm Formation

We wanted to determine if biofilm-forming ability is correlated with the stage of typhoid fever/carriage, antimicrobial resistance or with the presence of plasmids. Biofilms were examined in gallbladder-simulating conditions. There was varying ability to form biofilms under in vitro conditions in the S. Typhi strains tested. All isolates formed weak biofilms in absence of both cholesterol and bile (OD₅₇₀ below 0.3) and significantly strong biofilms in presence of cholesterol and 2.5% human bile (Figure 2A and 2B). Differences in biofilm forming ability was observed across the identified genotype 4.3.1 lineages/sub lineages. The sub-lineage 4.3.1.2EA2 formed the strongest biofilms (OD₅₇₀ slightly above 2.0) while sub-lineage 4.3.1.2EA3 formed relatively weak biofilms (OD₅₇₀ below 0.5) even in the presence of cholesterol and human bile (Figure 2A and 2B). The weak biofilm forming isolates, sub-lineage 4.3.1.2EA3, had AMR genes encoded in 2/19 of the plasmids identified in this sub-group. Sub-lineages 4.3.1.2EA2 and 4.3.1.1EA1 had fewer identified plasmids, (14 and 15 respectively) (Table 4), None of these plasmids had identifiable AMR genes, but isolates in these subgroups formed significantly stronger biofilms in the presence of both cholesterol and human bile (Figure 2C). However, there was no statistical significance in biofilm forming ability in strains isolated during the symptomatic vs. asymptomatic stage in each of the three sub-groups of S. Typhi (Figure 2D).

• Download figure
• Open in new tab

• Download figure
• Open in new tab

Figure 2.

Salmonella Typhi biofilms. (A and B) Quantity of biofilms after growth in presence/absence of cholesterol and/or bile. (C) Biofilm formation by MDR S. Typhi strains vs. non-MDR strains. (D) Comparison of biofilms formed by S. Typhi strains isolated before treatment vs. after treatment. Error bars represent SEM, ****, P<0.001; **, P<0.05; ^nsP>0.5.

Source of Bacteria strains

The whole genome sequences of 22 S. Typhi strains isolated from blood and stool samples of six patients living in four different households in Mukuru informal settlement, a typhoid endemic area in Nairobi, Kenya, were analyzed. Nineteen of these strains were from four patients who were also diagnosed with cholelithiasis, of whom, three had typhoid symptoms and continued to shed S. Typhi after treatment. The presence of gallstones was confirmed through an ultrasound scan, the primary imaging modality performed by a radiologist used to evaluate patients suspected of having gallbladder disease. One of the patients with gallstones was asymptomatic but shedding S. Typhi and living in the same household (household B) with an acute typhoid fever case (not diagnosed with gallstones). In a different household (household A), a study participant (household contact) not diagnosed with cholelithiasis shed S. Typhi once. The index case from this household was diagnosed with cholelithiasis and continued to shed S. Typhi after treatment with antibiotics. Laboratory methods on isolation and identification of the isolates are as described in our previous publication (38).

Whole Genome Sequencing

DNA extracted from S. Typhi strains using GenElute™ Bacterial Genomic DNA Kit (Missouri, United States) was prepared for whole genome sequencing by SeqCoast Genomics using an Illumina DNA Prep tagmentation kit and unique dual indexes. Sequencing was performed on the Illumina NextSeq2000 platform using a 300-cycle flow cell kit to produce 2×150bp paired reads as previously described (39). PhiX control, 1-2%, was spiked into the run to support optimal base calling. Read demultiplexing, read trimming, and run analytics were performed using DRAGEN v3.10.12, an on-board analysis software on the NextSeq2000. Raw sequencing data generated during this study are available in the National Center for Biotechnology Information (NCBI) public data, BioProject ID PRJNA1101423.

Genome Assembly and Annotation

Quality-trimming of the reads was done using Trimmomatic (40). Read error was corrected using SPAdes (41). Short (and long, if available) reads were assembled into contigs using SPAdes (wrapped in Unicycler). Read mapping was done using Bowtie2 and SAMtools (wrapped in Unicycler) (42–44) (https://github.com/rrwick/Porechop). Pilon (wrapped in Unicycler) was used in polishing of each assembly (45). Gene prediction and functional annotation was performed using BAKTA (46). The annotation pipeline was as follows; prediction of protein-coding genes using Prodigal, tRNA identification using tRNAscan-SE (47), tRNA and tmRNA identification using Aragorn (48), prediction of rRNA sequences using Infernal and the Rfam database (49), CRISPR prediction using PILER-CR (50), antimicrobial resistance gene identification using AMRFinderPlus (51), prediction of signal peptides using DeepSig (52), prediction of transposases using ISFinder (53) and computation of codon usage biases for each amino acid.

Genotype Identification and Bacteria Clustering

Identification and clustering of S. Typhi genotypes was performed using Pathogenwatch (https://pathogen.watch/), a web-based platform with several different components (54,55). The platform provides compatibility with S. Typhi typing information for MLST (56), in silico serotyping (SISTR) (57), and a SNP genotyping scheme (GenoTyphi) (55). S. Typhi Assemblies (in fasta format) were uploaded to the platform https://pathogen.watch/upload/fasta for the analysis.

Antimicrobial Susceptibility Testing

Antimicrobial susceptibility testing was performed using the disk diffusion technique for all antimicrobials commonly used in Kenya for typhoid fever treatment including ampicillin (10 µg), tetracycline (30 µg), co-trimoxazole (25 µg), chloramphenicol (30 µg), amoxicillin-clavulanate (20/10 μg), cefpodoxime 30 µg, ceftazidime (30 µg), ceftriaxone (30 µg), cefotaxime (30 µg), azithromycin (15 μg), ciprofloxacin (5 µg), nalidixic acid (10 µg), kanamycin (30 μg) and gentamicin (10 μg). The diameter of the zone of inhibition was measured after 18-24 hours and results were interpreted according the Clinical and Laboratory Standards Institute (CLSI) guidelines for Salmonella (CLSI 2023).

Screening of Antimicrobial Resistance Genes

NCBI Antimicrobial Resistance Gene Finder Plus (AMRFinderPlus) (https://github.com/ncbi/amr/wiki), was used to identify acquired antimicrobial resistance genes and known resistance-associated point mutations in S. Typhi assembled nucleotide sequences (51). The results were compared with phenotypic susceptibility data.

Variant Calling

To identify genetic variations between the longitudinal clinical isolates, S. Typhi strains isolated before each typhoid index case was treated with antibiotics were used as the reference genome, and were compared with those isolated after treatment (follow-up isolates), and/or those isolated from household contacts using breseq (58). The variant calling pipeline was as follows; quality-filtering of raw reads using Trimmomatic (40), mapping of reads against a reference genome (first strain isolated from the index case), analysis of possible mutations based on mapping data, identification of mutations and graphical and tabular summaries of mutation profile across samples (42).

Plasmid Identification

To identify plasmids in the isolated S. Typhi strains, a plasmid detection tool PLASMe was used (https://github.com/HubertTang/PLASMe). The tool uses the alignment component in PLASMe to identify closely related plasmids while diverged plasmids are predicted using order-specific Transformer models (59).

In Vitro Biofilm Formation Assays

Because of the importance of biofilm formation on gallstones in chronic carriage (17,19) the biofilm forming ability of all the 22 human S. Typhi isolates were tested under gallbladder simulating conditions. S. Typhi biofilms were grown on non-treated polystyrene 96-well plates (Corning, Kennebunkport, ME). To simulate growth conditions on gallstones, wells in two plates were pre-coated with cholesterol by adding a solution of 5 mg/mL in 1:1 isopropanol:ethanol and air-dried overnight. A pure colony of S. Typhi on an XLD agar plate was cultured in Tryptone Soy Broth (TSB). Overnight (O/N) cultures in broth were normalized to OD₆₀₀=0.8, diluted 1:2500 in TSB or TSB containing 2.5% human bile, and 100 µL/well were dispensed into the plates. The plates were incubated at 25°C in a Fisherbrand™ nutating mixer (Thermo Fisher Scientific; Hampton, NH) at 24 rpm for 96 hours. Media (TSB or TSB containing bile) was changed after every 24 hours for consistent S. Typhi biofilm growth. Plates were emptied on the fourth day and washed twice before heat fixing at 60°C for 1 hour. The biofilms were stained using a crystal violet solution and acetic acid (33%) used to elute crystal violet before reading the OD₅₇₀. GraphPad prism 9.5 was used to analyze the biofilm formation results. One way Analysis of variance (ANOVA) was used to test level of significance in biofilm formation between the different S. Typhi sub-lineages and in different conditions, i.e., biofilms in absence of cholesterol and bile, in cholesterol coated plates in absence of bile, and in presence of cholesterol and bile. Student’s t-test was used to test level of significance in biofilm formation in strains isolated before treatment vs last strains shed by the patient, P-values less than 0.05 (P<0.05) were considered significant.