Unravelling patient pathways in the context of antibacterial resistance in East Africa

BACKGROUND

A key factor driving the development and maintenance of antibacterial resistance (ABR) is individuals' use of antibiotics (ABs) to treat illness. To better understand motivations and context for antibiotic use we use the concept of a patient treatment-seeking pathway: a treatment journey encompassing where patients go when they are unwell, what motivates their choices, and how they obtain antibiotics. This paper investigates patterns and determinants of patient treatment-seeking pathways, and how they intersect with AB use in East Africa, a region where ABR-attributable deaths are exceptionally high.

METHODS

The Holistic Approach to Unravelling Antibacterial Resistance (HATUA) Consortium collected quantitative data from 6,827 adult outpatients presenting with urinary tract infection (UTI) symptoms in Kenya, Tanzania, and Uganda between February 2019- September 2020, and conducted qualitative in-depth patient interviews with a subset (n = 116). We described patterns of treatment-seeking visually using Sankey plots and explored explanations and motivations using mixed-methods. Using Bayesian hierarchical regression modelling, we investigated the associations between socio-demographic, economic, healthcare, and attitudinal factors and three factors related to ABR: self-treatment as a first step, having a multi-step treatment pathway, and consuming ABs.

RESULTS

Although most patients (86%) sought help from medical facilities in the first instance, many (56%) described multi-step, repetitive treatment-seeking pathways, which further increased the likelihood of consuming ABs. Higher socio-economic status patients were more likely to consume ABs and have multi-step pathways. Reasons for choosing providers (e.g., cost, location, time) were conditioned by wider structural factors such as hybrid healthcare systems and AB availability.

CONCLUSION

There is likely to be a reinforcing cycle between complex, repetitive treatment pathways, AB consumption and ABR. A focus on individual antibiotic use as the key intervention point in this cycle ignores the contextual challenges patients face when treatment seeking, which include inadequate access to diagnostics, perceived inefficient public healthcare and ease of purchasing antibiotics without prescription. Pluralistic healthcare landscapes may promote more complex treatment seeking and therefore inappropriate AB use. We recommend further attention to healthcare system factors, focussing on medical facilities (e.g., accessible diagnostics, patient-doctor interactions, information flows), and community AB access points (e.g., drug sellers).

Pan-Resistome Characterization of Uropathogenic Escherichia coli and Klebsiella pneumoniae Strains Circulating in Uganda and Kenya, Isolated from 2017-2018

Urinary tract infection (UTI) develops after a pathogen adheres to the inner lining of the urinary tract. Cases of UTIs are predominantly caused by several Gram-negative bacteria and account for high morbidity in the clinical and community settings. Of greater concern are the strains carrying antimicrobial resistance (AMR)-conferring genes. The gravity of a UTI is also determined by a spectrum of other virulence factors. This study represents a pilot project to investigate the burden of AMR among uropathogens in East Africa. We examined bacterial samples isolated in 2017–2018 from in- and out-patients in Kenya (KY) and Uganda (UG) that presented with clinical symptoms of UTI. We reconstructed the evolutionary history of the strains, investigated their population structure, and performed comparative analysis their pangenome contents. We found 55 Escherichia coli and 19 Klebsiella pneumoniae strains confirmed uropathogenic following screening for the prevalence of UTI virulence genes including fimH, iutA, feoA/B/C, mrkD, and foc. We identified 18 different sequence types in E. coli population while all K. pneumoniae strains belong to ST11. The most prevalent E. coli sequence types were ST131 (26%), ST335/1193 (10%), and ST10 (6%). Diverse plasmid types were observed in both collections such as Incompatibility (IncF/IncH/IncQ1/IncX4) and Col groups. Pangenome analysis of each set revealed a total of 2862 and 3464 genes comprised the core genome of E. coli and K. pneumoniae population, respectively. Among these are acquired AMR determinants including fluoroquinolone resistance-conferring genes aac(3)-Ib-cr and other significant genes: aad, tet, sul1, sul2, and cat, which are associated with aminoglycoside, tetracycline, sulfonamide, and chloramphenicol resistance, respectively. Accessory genomes of both species collections were detected several β-lactamase genes, bla_CTX-M, bla_TEM and bla_OXA, or bla_NDM. Overall, 93% are multi-drug resistant in the E. coli collection while 100% of the K. pneumoniae strains contained genes that are associated with resistance to three or more antibiotic classes. Our findings illustrate the abundant acquired resistome and virulome repertoire in uropathogenic E. coli and K. pneumoniae, which are mainly disseminated via clonal and horizontal transfer, circulating in the East African region. We further demonstrate here that routine genomic surveillance is necessary for high-resolution bacterial epidemiology of these important AMR pathogens.

Plasmids shape the diverse accessory resistomes of Escherichia coli ST131

The human gut microbiome includes beneficial, commensal and pathogenic bacteria that possess antimicrobial resistance (AMR) genes and exchange these predominantly through conjugative plasmids. Escherichia coli is a significant component of the gastrointestinal microbiome and is typically non-pathogenic in this niche. In contrast, extra-intestinal pathogenic E. coli (ExPEC) including ST131 may occupy other environments like the urinary tract or bloodstream where they express genes enabling AMR and host cell adhesion like type 1 fimbriae. The extent to which commensal E. coli and uropathogenic ExPEC ST131 share AMR genes remains understudied at a genomic level, and we examined this here using a preterm infant resistome. We found that individual ST131 had small differences in AMR gene content relative to a larger shared resistome. Comparisons with a range of plasmids common in ST131 showed that AMR gene composition was driven by conjugation, recombination and mobile genetic elements. Plasmid pEK499 had extended regions in most ST131 Clade C isolates, and it had evidence of a co-evolutionary signal based on protein-level interactions with chromosomal gene products, as did pEK204 that had a type IV fimbrial pil operon. ST131 possessed extensive diversity of selective type 1, type IV, P and F17-like fimbriae genes that was highest in subclade C2. The structure and composition of AMR genes, plasmids and fimbriae vary widely in ST131 Clade C and this may mediate pathogenicity and infection outcomes.

Genomic surveillance of Escherichia coli ST131 identifies local expansion and serial replacement of subclones

Escherichia coli sequence type 131 (ST131) is a pandemic clone that is evolving rapidly with increasing levels of antimicrobial resistance. Here, we investigated an outbreak of E. coli ST131 producing extended spectrum β-lactamases (ESBLs) in a long-term care facility (LTCF) in Ireland by combining data from this LTCF (n=69) with other Irish (n=35) and global (n=690) ST131 genomes to reconstruct the evolutionary history and understand changes in population structure and genome architecture over time. This required a combination of short- and long-read genome sequencing, de novo assembly, read mapping, ESBL gene screening, plasmid alignment and temporal phylogenetics. We found that Clade C was the most prevalent (686 out of 794 isolates, 86 %) of the three major ST131 clades circulating worldwide (A with fimH41, B with fimH22, C with fimH30), and was associated with the presence of different ESBL alleles, diverse plasmids and transposable elements. Clade C was estimated to have emerged in c. 1985 and subsequently acquired different ESBL gene variants (blaCTX-M-14 vs blaCTX-M-15). An ISEcp1-mediated transposition of the blaCTX-M-15 gene further increased the diversity within Clade C. We discovered a local clonal expansion of a rare C2 lineage (C2_8) with a chromosomal insertion of blaCTX-M-15 at the mppA gene. This was acquired from an IncFIA plasmid. The C2_8 lineage clonally expanded in the Irish LTCF from 2006, displacing the existing C1 strain (C1_10), highlighting the potential for novel ESBL-producing ST131 with a distinct genetic profile to cause outbreaks strongly associated with specific healthcare environments.

An Escherichia coli ST131 pangenome atlas reveals population structure and evolution across 4,071 isolates

Escherichia coli ST131 is a major cause of infection with extensive antimicrobial resistance (AMR) facilitated by widespread beta-lactam antibiotic use. This drug pressure has driven extended-spectrum beta-lactamase (ESBL) gene acquisition and evolution in pathogens, so a clearer resolution of ST131's origin, adaptation and spread is essential. E. coli ST131's ESBL genes are typically embedded in mobile genetic elements (MGEs) that aid transfer to new plasmid or chromosomal locations, which are mobilised further by plasmid conjugation and recombination, resulting in a flexible ESBL, MGE and plasmid composition with a conserved core genome. We used population genomics to trace the evolution of AMR in ST131 more precisely by extracting all available high-quality Illumina HiSeq read libraries to investigate 4,071 globally-sourced genomes, the largest ST131 collection examined so far. We applied rigorous quality-control, genome de novo assembly and ESBL gene screening to resolve ST131's population structure across three genetically distinct Clades (A, B, C) and abundant subclades from the dominant Clade C. We reconstructed their evolutionary relationships across the core and accessory genomes using published reference genomes, long read assemblies and k-mer-based methods to contextualise pangenome diversity. The three main C subclades have co-circulated globally at relatively stable frequencies over time, suggesting attaining an equilibrium after their origin and initial rapid spread. This contrasted with their ESBL genes, which had stronger patterns across time, geography and subclade, and were located at distinct locations across the chromosomes and plasmids between isolates. Within the three C subclades, the core and accessory genome diversity levels were not correlated due to plasmid and MGE activity, unlike patterns between the three main clades, A, B and C. This population genomic study highlights the dynamic nature of the accessory genomes in ST131, suggesting that surveillance should anticipate genetically variable outbreaks with broader antibiotic resistance levels. Our findings emphasise the potential of evolutionary pangenomics to improve our understanding of AMR gene transfer, adaptation and transmission to discover accessory genome changes linked to novel subtypes.

Complete Assembly of Escherichia coli Sequence Type 131 Genomes Using Long Reads Demonstrates Antibiotic Resistance Gene Variation within Diverse Plasmid and Chromosomal Contexts

The incidence of infections caused by extraintestinal Escherichia coli (ExPEC) is rising globally, which is a major public health concern. ExPEC strains that are resistant to antimicrobials have been associated with excess mortality, prolonged hospital stays, and higher health care costs. E. coli sequence type 131 (ST131) is a major ExPEC clonal group worldwide, with variable plasmid composition, and has an array of genes enabling antimicrobial resistance (AMR). ST131 isolates frequently encode the AMR genes blaCTX-M-14, blaCTX-M-15, and blaCTX-M-27, which are often rearranged, amplified, and translocated by mobile genetic elements (MGEs). Short DNA reads do not fully resolve the architecture of repetitive elements on plasmids to allow MGE structures encoding blaCTX-M genes to be fully determined. Here, we performed long-read sequencing to decipher the genome structures of six E. coli ST131 isolates from six patients. Most long-read assemblies generated entire chromosomes and plasmids as single contigs, in contrast to more fragmented assemblies created with short reads alone. The long-read assemblies highlighted diverse accessory genomes with blaCTX-M-15, blaCTX-M-14, and blaCTX-M-27 genes identified in three, one, and one isolates, respectively. One sample had no blaCTX-M gene. Two samples had chromosomal blaCTX-M-14 and blaCTX-M-15 genes, and the latter was at three distinct locations, likely transposed by the adjacent MGEs: ISEcp1, IS903B, and Tn2 This study showed that AMR genes exist in multiple different chromosomal and plasmid contexts, even between closely related isolates within a clonal group such as E. coli ST131.IMPORTANCE Drug-resistant bacteria are a major cause of illness worldwide, and a specific subtype called Escherichia coli ST131 causes a significant number of these infections. ST131 bacteria become resistant to treatments by modifying their DNA and by transferring genes among one another via large packages of genes called plasmids, like a game of pass-the-parcel. Tackling infections more effectively requires a better understanding of what plasmids are being exchanged and their exact contents. To achieve this, we applied new high-resolution DNA sequencing technology to six ST131 samples from infected patients and compared the output to that of an existing approach. A combination of methods shows that drug resistance genes on plasmids are highly mobile because they can jump into ST131's chromosomes. We found that the plasmids are very elastic and undergo extensive rearrangements even in closely related samples. This application of DNA sequencing technologies illustrates at a new level the highly dynamic nature of ST131 genomes.