Antimicrobial-resistant Neisseria gonorrhoeae in Europe in 2020 compared with in 2013 and 2018: a retrospective genomic surveillance study

BACKGROUND

Regular quality-assured whole-genome sequencing linked to antimicrobial resistance (AMR) and patient metadata is imperative to elucidate the shifting gonorrhoea epidemiology, both nationally and internationally. We aimed to examine the gonococcal population in the European Economic Area (EEA) in 2020, elucidate emerging and disappearing gonococcal lineages associated with AMR and patient metadata, compare with 2013 and 2018 whole-genome sequencing data, and explain changes in gonococcal AMR and gonorrhoea epidemiology.

METHODS

In this retrospective genomic surveillance study, we analysed consecutive gonococcal isolates that were collected in EEA countries through the European Gonococcal Antimicrobial Surveillance Programme (Euro-GASP) in 2020, and made comparisons with Euro-GASP data from 2013 and 2018. All isolates had linked AMR data (based on minimum inhibitory concentration determination) and patient metadata. We performed whole-genome sequencing and molecular typing and AMR determinants were derived from quality-checked whole-genome sequencing data. Links between genomic lineages, AMR, and patient metadata were examined.

FINDINGS

1932 gonococcal isolates collected in 2020 in 21 EEA countries were included. The majority (81·2%, 147 of 181 isolates) of azithromycin resistance (present in 9·4%, 181 of 1932) was explained by the continued expansion of the Neisseria gonorrhoeae sequence typing for antimicrobial resistance (NG-STAR) clonal complexes (CCs) 63, 168, and 213 (with mtrD/mtrR promoter mosaic 2) and the novel NG-STAR CC1031 (semi-mosaic mtrD variant 13), associated with men who have sex with men and anorectal or oropharyngeal infections. The declining cefixime resistance (0·5%, nine of 1932) and negligible ceftriaxone resistance (0·1%, one of 1932) was largely because of the progressive disappearance of NG-STAR CC90 (with mosaic penA allele), which was predominant in 2013. No known resistance determinants for novel antimicrobials (zoliflodacin, gepotidacin, and lefamulin) were found.

INTERPRETATION

Azithromycin-resistant clones, mainly with mtrD mosaic or semi-mosaic variants, appear to be stabilising at a relatively high level in the EEA. This mostly low-level azithromycin resistance might threaten the recommended ceftriaxone-azithromycin therapy, but the negligible ceftriaxone resistance is encouraging. The decreased genomic population diversity and increased clonality could be explained in part by the COVID-19 pandemic resulting in lower importation of novel strains into Europe.

FUNDING

European Centre for Disease Prevention and Control and Örebro University Hospital.

Capturing Data on Antimicrobial Resistance Patterns and Trends in Use in Regions of Asia (CAPTURA)

BACKGROUND

In 2015, the UK government established the Fleming Fund with the aim to address critical gaps in surveillance of antimicrobial resistance (AMR) in low- and middle-income countries in Asia and Africa. Among a large portfolio of grants, the Capturing Data on Antimicrobial Resistance Patterns and Trends in Use in Regions of Asia (CAPTURA) project was awarded with the specific objective of expanding the volume of historical data on AMR, consumption (AMC), and use (AMU) in the human healthcare sector across 12 countries in South and Southeast Asia.

METHODS

Starting in early 2019, the CAPTURA consortium began working with local governments and >100 relevant data-holding facilities across the region to identify, assess for quality, prioritize, and subsequently retrieve data on AMR, AMC, and AMU. Relevant and shared data were collated and analyzed to provide local overviews for national stakeholders as well as regional context, wherever possible.

RESULTS

From the vast information resource generated on current surveillance capacity and data availability, the project has highlighted gaps and areas for quality improvement and supported comprehensive capacity-building activities to optimize local data-collection and -management practices.

CONCLUSIONS

The project has paved the way for expansion of surveillance networks to include both the academic and private sector in several countries and has actively engaged in discussions to promote data sharing at the local, national, and regional levels. This paper describes the overarching approach to, and emerging lessons from, the CAPTURA project, and how it contributes to other ongoing efforts to strengthen national AMR surveillance in the region and globally.

Methodological Approach to Identify and Expand the Volume of Antimicrobial Resistance (AMR) Data in the Human Health Sector in Low- and Middle-Income Countries in Asia: Implications for Local and Regional AMR Surveillance Systems Strengthening

Antimicrobial resistance (AMR) is a multifaceted global health problem disproportionately affecting low- and middle-income countries (LMICs). The Capturing data on Antimicrobial resistance Patterns and Trends in Use in Regions of Asia (CAPTURA) project was tasked to expand the volume of AMR and antimicrobial use data in Asia. The CAPTURA project used 2 data-collection streams: facility data and project metadata. Project metadata constituted information collected to map out data sources and assess data quality, while facility data referred to the retrospective data collected from healthcare facilities. A down-selection process, labelled "the funnel approach" by the project, was adopted to use the project metadata in prioritizing and selecting laboratories for retrospective AMR data collection. Moreover, the metadata served as a guide for understanding the AMR data once they were collected. The findings from CAPTURA's metadata add to the current discourse on the limitation of AMR data in LMICs. There is generally a low volume of AMR data generated as there is a lack of microbiology laboratories with sufficient antimicrobial susceptibility testing capacity. Many laboratories in Asia are still capturing data on paper, resulting in scattered or unused data not readily accessible or shareable for analyses. There is also a lack of clinical and epidemiological data captured, impeding interpretation and in-depth understanding of the AMR data. CAPTURA's experience in Asia suggests that there is a wide spectrum of capacity and capability of microbiology laboratories within a country and region. As local AMR surveillance is a crucial instrument to inform context-specific measures to combat AMR, it is important to understand and assess current capacity-building needs while implementing activities to enhance surveillance systems.

Genomics for public health and international surveillance of antimicrobial resistance

Historically, epidemiological investigation and surveillance for bacterial antimicrobial resistance (AMR) has relied on low-resolution isolate-based phenotypic analyses undertaken at local and national reference laboratories. Genomic sequencing has the potential to provide a far more high-resolution picture of AMR evolution and transmission, and is already beginning to revolutionise how public health surveillance networks monitor and tackle bacterial AMR. However, the routine integration of genomics in surveillance pipelines still has considerable barriers to overcome. In 2022, a workshop series and online consultation brought together international experts in AMR and pathogen genomics to assess the status of genomic applications for AMR surveillance in a range of settings. Here we focus on discussions around the use of genomics for public health and international AMR surveillance, noting the potential advantages of, and barriers to, implementation, and proposing recommendations from the working group to help to drive the adoption of genomics in public health AMR surveillance. These recommendations include the need to build capacity for genome sequencing and analysis, harmonising and standardising surveillance systems, developing equitable data sharing and governance frameworks, and strengthening interactions and relationships among stakeholders at multiple levels.

Evidence review and recommendations for the implementation of genomics for antimicrobial resistance surveillance: reports from an international expert group

Nearly a century after the beginning of the antibiotic era, which has been associated with unparalleled improvements in human health and reductions in mortality associated with infection, the dwindling pipeline for new antibiotic classes coupled with the inevitable spread of antimicrobial resistance (AMR) poses a major global challenge. Historically, surveillance of bacteria with AMR typically relied on phenotypic analysis of isolates taken from infected individuals, which provides only a low-resolution view of the epidemiology behind an individual infection or wider outbreak. Recent years have seen increasing adoption of powerful new genomic technologies with the potential to revolutionise AMR surveillance by providing a high-resolution picture of the AMR profile of the bacteria causing infections and providing real-time actionable information for treating and preventing infection. However, many barriers remain to be overcome before genomic technologies can be adopted as a standard part of routine AMR surveillance around the world. Accordingly, the Surveillance and Epidemiology of Drug-resistant Infections Consortium convened an expert working group to assess the benefits and challenges of using genomics for AMR surveillance. In this Series, we detail these discussions and provide recommendations from the working group that can help to realise the massive potential benefits for genomics in surveillance of AMR.

Novel multidrug-resistant sublineages of Staphylococcus aureus clonal complex 22 discovered in India

Staphylococcus aureus is a major pathogen in India causing community and nosocomial infections, but little is known about its molecular epidemiology and mechanisms of resistance in hospital settings. Here, we use whole-genome sequencing (WGS) to characterize 478 S. aureus clinical isolates (393 methicillin-resistant Staphylococcus aureus (MRSA) and 85 methicilin-sensitive Staphylococcus aureus (MSSA) collected from 17 sentinel sites across India between 2014 and 2019. Sequencing results confirmed that sequence type 22 (ST22) (142 isolates, 29.7%), ST239 (74 isolates, 15.48%), and ST772 (67 isolates, 14%) were the most common clones. An in-depth analysis of 175 clonal complex (CC) 22 Indian isolates identified two novel ST22 MRSA lineages, both Panton-Valentine leukocidin+, both resistant to fluoroquinolones and aminoglycosides, and one harboring the the gene for toxic shock syndrome toxin 1 (tst). A temporal analysis of 1797 CC22 global isolates from 14 different studies showed that the two Indian ST22 lineages shared a common ancestor in 1984 (95% highest posterior density [HPD]: 1982-1986), as well as evidence of transmission to other parts of the world. Moreover, the study also gives a comprehensive view of ST2371, a sublineage of CC22, as a new emerging lineage in India and describes it in relationship with the other Indian ST22 isolates. In addition, the retrospective identification of a putative outbreak of multidrug-resistant (MDR) ST239 from a single hospital in Bangalore that persisted over a period of 3 years highlights the need for the implementation of routine surveillance and simple infection prevention and control measures to reduce these outbreaks. To our knowledge, this is the first WGS study that characterized CC22 in India and showed that the Indian clones are distinct from the EMRSA-15 clone. Thus, with the improved resolution afforded by WGS, this study substantially contributed to our understanding of the global population of MRSA. IMPORTANCE The study conducted in India between 2014 and 2019 presents novel insights into the prevalence of MRSA in the region. Previous studies have characterized two dominant clones of MRSA in India, ST772 and ST239, using whole-genome sequencing. However, this study is the first to describe the third dominant clone, ST22, using the same approach. The ST22 Indian isolates were analyzed in-depth, leading to the discovery of two new sublineages of hospital-acquired Staphylococcus aureus in India, both carrying antimicrobial resistance genes and mutations, which limit treatment options for patients. One of the newly characterized sublineages, second Indian cluster, carries the tsst-1 virulence gene, increasing the risk of severe infections. The geographic spread of the two novel lineages, both within India and internationally, could pose a global public health threat. The study also sheds light on ST2371 in India, a single-locus variant of ST22. The identification of a putative outbreak of MDR ST239 in a single hospital in Bangalore emphasizes the need for routine surveillance and simple infection prevention and control measures to reduce these outbreaks. Overall, this study significantly contributes to our understanding of the global population of MRSA, thanks to the improved resolution afforded by WGS.

Global diversity and antimicrobial resistance of typhoid fever pathogens: Insights from a meta-analysis of 13,000 Salmonella Typhi genomes

Background

The Global Typhoid Genomics Consortium was established to bring together the typhoid research community to aggregate and analyse Salmonella enterica serovar Typhi (Typhi) genomic data to inform public health action. This analysis, which marks 22 years since the publication of the first Typhi genome, represents the largest Typhi genome sequence collection to date (n=13,000).

Methods

This is a meta-analysis of global genotype and antimicrobial resistance (AMR) determinants extracted from previously sequenced genome data and analysed using consistent methods implemented in open analysis platforms GenoTyphi and Pathogenwatch.

Results

Compared with previous global snapshots, the data highlight that genotype 4.3.1 (H58) has not spread beyond Asia and Eastern/Southern Africa; in other regions, distinct genotypes dominate and have independently evolved AMR. Data gaps remain in many parts of the world, and we show the potential of travel-associated sequences to provide informal 'sentinel' surveillance for such locations. The data indicate that ciprofloxacin non-susceptibility (>1 resistance determinant) is widespread across geographies and genotypes, with high-level ciprofloxacin resistance (≥3 determinants) reaching 20% prevalence in South Asia. Extensively drug-resistant (XDR) typhoid has become dominant in Pakistan (70% in 2020) but has not yet become established elsewhere. Ceftriaxone resistance has emerged in eight non-XDR genotypes, including a ciprofloxacin-resistant lineage (4.3.1.2.1) in India. Azithromycin resistance mutations were detected at low prevalence in South Asia, including in two common ciprofloxacin-resistant genotypes.

Conclusions

The consortium's aim is to encourage continued data sharing and collaboration to monitor the emergence and global spread of AMR Typhi, and to inform decision-making around the introduction of typhoid conjugate vaccines (TCVs) and other prevention and control strategies.

Funding

No specific funding was awarded for this meta-analysis. Coordinators were supported by fellowships from the European Union (ZAD received funding from the European Union's Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No 845681), the Wellcome Trust (SB, Wellcome Trust Senior Fellowship), and the National Health and Medical Research Council (DJI is supported by an NHMRC Investigator Grant [GNT1195210]).

High Genetic Diversity of Carbapenem-Resistant Acinetobacter baumannii Isolates Recovered in Nigerian Hospitals in 2016 to 2020

Acinetobacter baumannii causes difficult-to-treat infections mostly among immunocompromised patients. Clinically relevant A. baumannii lineages and their carbapenem resistance mechanisms are sparsely described in Nigeria. This study aimed to characterize the diversity and genetic mechanisms of carbapenem resistance among A. baumannii strains isolated from hospitals in southwestern Nigeria. We sequenced the genomes of all A. baumannii isolates submitted to Nigeria's antimicrobial resistance surveillance reference laboratory between 2016 and 2020 on an Illumina platform and performed in silico genomic characterization. Selected strains were sequenced using the Oxford Nanopore technology to characterize the genetic context of carbapenem resistance genes. The 86 A. baumannii isolates were phylogenetically diverse and belonged to 35 distinct Oxford sequence types (oxfSTs), 16 of which were novel, and 28 Institut Pasteur STs (pasSTs). Thirty-eight (44.2%) isolates belonged to none of the known international clones (ICs). Over 50% of the isolates were phenotypically resistant to 10 of 12 tested antimicrobials. The majority (n = 54) of the isolates were carbapenem resistant, particularly the IC7 (pasST25; 100%) and IC9 (pasST85; >91.7%) strains. blaOXA-23 (34.9%) and blaNDM-1 (27.9%) were the most common carbapenem resistance genes detected. All blaOXA-23 genes were carried on Tn2006 or Tn2006-like transposons. Our findings suggest that a 10-kb Tn125 composite transposon is the primary means of blaNDM-1 dissemination. Our findings highlight an increase in blaNDM-1 prevalence and the widespread transposon-facilitated dissemination of carbapenemase genes in diverse A. baumannii lineages in southwestern Nigeria. We make the case for improving surveillance of these pathogens in Nigeria and other understudied settings. IMPORTANCE Acinetobacter baumannii bacteria are increasingly clinically relevant due to their propensity to harbor genes conferring resistance to multiple antimicrobials, as well as their ability to persist and disseminate in hospital environments and cause difficult-to-treat nosocomial infections. Little is known about the molecular epidemiology and antimicrobial resistance profiles of these organisms in Nigeria, largely due to limited capacity for their isolation, identification, and antimicrobial susceptibility testing. Our study characterized the diversity and antimicrobial resistance profiles of clinical A. baumannii in southwestern Nigeria using whole-genome sequencing. We also identified the key genetic elements facilitating the dissemination of carbapenem resistance genes within this species. This study provides key insights into the clinical burden and population dynamics of A. baumannii in hospitals in Nigeria and highlights the importance of routine whole-genome sequencing-based surveillance of this and other previously understudied pathogens in Nigeria and other similar settings.

Genomic epidemiology of Staphylococcus aureus isolated from bloodstream infections in South America during 2019 supports regional surveillance

Staphylococcus aureus remains one of the leading causes of infections worldwide and a common cause of bacteraemia. However, studies documenting the epidemiology of S. aureus in South America using genomics are scarce. We hereby report on the largest genomic epidemiology study to date of both methicillin-resistant S. aureus (MRSA) and methicillin-susceptible S. aureus (MSSA) in South America, conducted by the StaphNET-SA network. We characterised 404 genomes recovered from a prospective observational study of S. aureus bacteraemia in 58 hospitals from Argentina, Bolivia, Brazil, Paraguay and Uruguay between April and October 2019. We show that a minority of S. aureus isolates are phenotypically multi-drug resistant (5.2%), but more than a quarter are resistant to macrolide-lincosamide-streptogramin B (MLSb). MSSA were more genetically diverse than MRSA. Lower rates of associated antimicrobial resistance in community-associated(CA)-MRSA versus hospital-associated (HA)-MRSA were found in association with three S. aureus genotypes dominating the MRSA population: CC30-MRSA-IVc-t019-lukS/F-PV+, CC5-MRSA-IV-t002-lukS/F-PV- and CC8-MRSA-IVc-t008-lukS/F-PV+-COMER+. These are historically from a CA origin, carry on average fewer antimicrobial resistance determinants, and often lack key virulence genes. Surprisingly, CC398-MSSA-t1451-lukS/F-PV- related to the CC398 human-associated lineage is widely disseminated throughout the region, and is described here for the first time as the most prevalent MSSA lineage in South America. Moreover, CC398 strains carrying ermT (largely responsible for the MLSb resistance rates of MSSA strains: inducible iMLSb phenotype) and sh_fabI (related to triclosan resistance) were recovered from both CA and HA origin. The frequency of MRSA and MSSA lineages differed between countries but the most prevalent S. aureus genotypes are high-risk clones widely distributed in the South American region without a clear country-specific phylogeographical structure. Therefore, our findings underline the need for continuous genomic surveillance by regional networks such as StaphNET-SA. This article contains data hosted by Microreact.

Comparing genomic variant identification protocols for Candida auris

Genomic analyses are widely applied to epidemiological, population genetic and experimental studies of pathogenic fungi. A wide range of methods are employed to carry out these analyses, typically without including controls that gauge the accuracy of variant prediction. The importance of tracking outbreaks at a global scale has raised the urgency of establishing high-accuracy pipelines that generate consistent results between research groups. To evaluate currently employed methods for whole-genome variant detection and elaborate best practices for fungal pathogens, we compared how 14 independent variant calling pipelines performed across 35 Candida auris isolates from 4 distinct clades and evaluated the performance of variant calling, single-nucleotide polymorphism (SNP) counts and phylogenetic inference results. Although these pipelines used different variant callers and filtering criteria, we found high overall agreement of SNPs from each pipeline. This concordance correlated with site quality, as SNPs discovered by a few pipelines tended to show lower mapping quality scores and depth of coverage than those recovered by all pipelines. We observed that the major differences between pipelines were due to variation in read trimming strategies, SNP calling methods and parameters, and downstream filtration criteria. We calculated specificity and sensitivity for each pipeline by aligning three isolates with chromosomal level assemblies and found that the GATK-based pipelines were well balanced between these metrics. Selection of trimming methods had a greater impact on SAMtools-based pipelines than those using GATK. Phylogenetic trees inferred by each pipeline showed high consistency at the clade level, but there was more variability between isolates from a single outbreak, with pipelines that used more stringent cutoffs having lower resolution. This project generated two truth datasets useful for routine benchmarking of C. auris variant calling, a consensus VCF of genotypes discovered by 10 or more pipelines across these 35 diverse isolates and variants for 2 samples identified from whole-genome alignments. This study provides a foundation for evaluating SNP calling pipelines and developing best practices for future fungal genomic studies.

Global scenario of the RmtE pan-aminoglycoside-resistance mechanism: emergence of the rmtE4 gene in South America associated with a hospital-related IncL plasmid

Antimicrobial resistance (AMR) mechanisms, especially those conferring resistance to critically important antibiotics, are a great concern for public health. 16S rRNA methyltransferases (16S-RMTases) abolish the effectiveness of most clinically used aminoglycosides, but some of them are considered sporadic, such as RmtE. The main goals of this work were the genomic analysis of bacteria producing 16S-RMTases from a 'One Health' perspective in Venezuela, and the study of the epidemiological and evolutionary scenario of RmtE variants and their related mobile genetic elements (MGEs) worldwide. A total of 21 samples were collected in 2014 from different animal and environmental sources in the Cumaná region (Venezuela). Highly aminoglycoside-resistant Enterobacteriaceae isolates were selected, identified and screened for 16S-RMTase genes. Illumina and Nanopore whole-genome sequencing data were combined to obtain hybrid assemblies and analyse their sequence type, resistome, plasmidome and pan-genome. Genomic collections of rmtE variants and their associated MGEs were generated to perform epidemiological and phylogenetic analyses. A single 16S-RMTase, the novel RmtE4, was identified in five Klebsiella isolates from wastewater samples of Cumaná. This variant possessed three amino acid modifications with respect to RmtE1-3 (Asn152Asp, Val216Ile and Lys267Ile), representing the most genetic distant among all known and novel variants described in this work, and the second most prevalent. rmtE variants were globally spread, and their geographical distribution was determined by the associated MGEs and the carrying bacterial species. Thus, rmtE4 was found to be confined to Klebsiella isolates from South America, where it was closely related to ISVsa3 and an uncommon IncL plasmid related with hospital environments. This work uncovered the global scenario of RmtE and the existence of RmtE4, which could potentially emerge from South America. Surveillance and control measures should be developed based on these findings in order to prevent the dissemination of this AMR mechanism and preserve public health worldwide.

Nanopore-only assemblies for genomic surveillance of the global priority drug-resistant pathogen, Klebsiella pneumoniae

Oxford Nanopore Technologies (ONT) sequencing has rich potential for genomic epidemiology and public health investigations of bacterial pathogens, particularly in low-resource settings and at the point of care, due to its portability and affordability. However, low base-call accuracy has limited the reliability of ONT data for critical tasks such as antimicrobial resistance (AMR) and virulence gene detection and typing, serotype prediction, and cluster identification. Thus, Illumina sequencing remains the standard for genomic surveillance despite higher capital and running costs. We tested the accuracy of ONT-only assemblies for common applied bacterial genomics tasks (genotyping and cluster detection, implemented via Kleborate, Kaptive and Pathogenwatch), using data from 54 unique Klebsiella pneumoniae isolates. ONT reads generated via MinION with R9.4.1 flowcells were basecalled using three alternative models [Fast, High-accuracy (HAC) and Super-accuracy (SUP), available within ONT's Guppy software], assembled with Flye and polished using Medaka. Accuracy of typing using ONT-only assemblies was compared with that of Illumina-only and hybrid ONT+Illumina assemblies, constructed from the same isolates as reference standards. The most resource-intensive ONT-assembly approach (SUP basecalling, with or without Medaka polishing) performed best, yielding reliable capsule (K) type calls for all strains (100 % exact or best matching locus), reliable multi-locus sequence type (MLST) assignment (98.3 % exact match or single-locus variants), and good detection of acquired AMR genes and mutations (88-100 % correct identification across the various drug classes). Distance-based trees generated from SUP+Medaka assemblies accurately reflected overall genetic relationships between isolates. The definition of outbreak clusters from ONT-only assemblies was problematic due to inflation of SNP counts by high base-call errors. However, ONT data could be reliably used to 'rule out' isolates of distinct lineages from suspected transmission clusters. HAC basecalling + Medaka polishing performed similarly to SUP basecalling without polishing. Therefore, we recommend investing compute resources into basecalling (SUP model), wherever compute resources and time allow, and note that polishing is also worthwhile for improved performance. Overall, our results show that MLST, K type and AMR determinants can be reliably identified with ONT-only R9.4.1 flowcell data. However, cluster detection remains challenging with this technology.

Conservation of vaccine antigen sequences encoded by sequenced strains of Streptococcus equi subsp. equi

BACKGROUND

Streptococcus equi subspecies equi (S equi) is the cause of Strangles, one of the most prevalent diseases of horses worldwide. Variation within the immunodominant SeM protein has been documented, but a new eight-component fusion protein vaccine, Strangvac, does not contain live S equi or SeM and conservation of the antigens it contains have not been reported.

OBJECTIVE

To define the diversity of the eight Strangvac antigens across a diverse S equi population.

STUDY DESIGN

Genomic description.

METHODS

Antigen sequences from the genomes of 759 S equi isolates from 19 countries, recovered between 1955 and 2018, were analysed. Predicted amino acid sequences in the antigen fragments of SEQ0256(Eq5), SEQ0402(Eq8), SEQ0721(EAG), SEQ0855(SclF), SEQ0935(CNE), SEQ0999(IdeE), SEQ1817(SclI) and SEQ2101(SclC) in Strangvac and SeM were extracted from the 759 assembled genomes and compared.

RESULTS

The predicted amino acid sequences of SclC, SclI and IdeE were identical across all 759 genomes. CNE was truncated in the genome of five (0.7%) isolates. SclF was absent from one genome and another encoded a single amino acid substitution. EAG was truncated in two genomes. Eq5 was truncated in four genomes and 123 genomes encoded a single amino acid substitution. Eq8 was truncated in three genomes, one genome encoded four amino acid substitutions and 398 genomes encoded a single amino acid substitution at the final amino acid of the Eq8 antigen fragment. Therefore, at least 1579 (99.9%) of 1580 amino acids in Strangvac were identical in 743 (97.9%) genomes, and all genomes encoded identical amino acid sequences for at least six of the eight Strangvac antigens.

MAIN LIMITATIONS

Three hundred and seven (40.4%) isolates in this study were recovered from horses in the UK.

CONCLUSIONS

The predicted amino acid sequences of antigens in Strangvac were highly conserved across this collection of S equi.

An ST131 clade and a phylogroup A clade bearing an O101-like O-antigen cluster predominate among bloodstream Escherichia coli isolates from South-West Nigeria hospitals

Escherichia coli bloodstream infections are typically attributed to a limited number of lineages that carry virulence factors associated with invasiveness. In Nigeria, the identity of circulating clones is largely unknown and surveillance of their antimicrobial resistance has been limited. We verified and whole-genome sequenced 68 2016-2018 bloodstream E. coli isolates from three sentinel sites in South-Western Nigeria and susceptibility tested 67 of them. Resistance to antimicrobials commonly used in Nigeria was high, with 67 (100 %), 62 (92.5 %), 53 (79.1 %) and 37 (55.2 %) showing resistance to trimethoprim, ampicillin, ciprofloxacin and aminoglycosides, respectively. Thirty-five (51 %) isolates carried extended-spectrum β-lactamase genes and 32 (91 %) of these were multidrug resistant. All the isolates were susceptible to carbapenems and colistin. The strain set included globally disseminated high-risk clones from sequence type (ST)12 (2), ST131 (12) and ST648 (4). Twenty-three (33.8 %) of the isolates clustered within two clades. The first of these consisted of ST131 strains, comprising O16:H5 and O25:H4 sub-lineages. The second was an ST10-ST167 complex clade comprising strains carrying O-antigen and capsular genes of likely Klebsiella origin, identical to those of avian pathogenic E. coli Sanji, and serotyped in silico as O89, O101 or ONovel32, depending on the tool used. Four temporally associated ST90 strains from one sentinel were closely related enough to suggest that at least some of them represented a retrospectively detected outbreak cluster. Our data implicate a broad repertoire of E. coli isolates associated with bloodstream infections in South-West Nigeria. Continued genomic surveillance is valuable for tracking clones of importance and for outbreak identification.

Context-specific emergence and growth of the SARS-CoV-2 Delta variant

The SARS-CoV-2 Delta (Pango lineage B.1.617.2) variant of concern spread globally, causing resurgences of COVID-19 worldwide1,2. The emergence of the Delta variant in the UK occurred on the background of a heterogeneous landscape of immunity and relaxation of non-pharmaceutical interventions. Here we analyse 52,992 SARS-CoV-2 genomes from England together with 93,649 genomes from the rest of the world to reconstruct the emergence of Delta and quantify its introduction to and regional dissemination across England in the context of changing travel and social restrictions. Using analysis of human movement, contact tracing and virus genomic data, we find that the geographic focus of the expansion of Delta shifted from India to a more global pattern in early May 2021. In England, Delta lineages were introduced more than 1,000 times and spread nationally as non-pharmaceutical interventions were relaxed. We find that hotel quarantine for travellers reduced onward transmission from importations; however, the transmission chains that later dominated the Delta wave in England were seeded before travel restrictions were introduced. Increasing inter-regional travel within England drove the nationwide dissemination of Delta, with some cities receiving more than 2,000 observable lineage introductions from elsewhere. Subsequently, increased levels of local population mixing-and not the number of importations-were associated with the faster relative spread of Delta. The invasion dynamics of Delta depended on spatial heterogeneity in contact patterns, and our findings will inform optimal spatial interventions to reduce the transmission of current and future variants of concern, such as Omicron (Pango lineage B.1.1.529).

Phylogenetic and phylodynamic approaches to understanding and combating the early SARS-CoV-2 pandemic

Determining the transmissibility, prevalence and patterns of movement of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections is central to our understanding of the impact of the pandemic and to the design of effective control strategies. Phylogenies (evolutionary trees) have provided key insights into the international spread of SARS-CoV-2 and enabled investigation of individual outbreaks and transmission chains in specific settings. Phylodynamic approaches combine evolutionary, demographic and epidemiological concepts and have helped track virus genetic changes, identify emerging variants and inform public health strategy. Here, we review and synthesize studies that illustrate how phylogenetic and phylodynamic techniques were applied during the first year of the pandemic, and summarize their contributions to our understanding of SARS-CoV-2 transmission and control.

Genomic characterization of invasive typhoidal and non-typhoidal Salmonella in southwestern Nigeria

Background

Salmonellosis causes significant morbidity and mortality in Africa. Information on lineages of invasive Salmonella circulating in Nigeria is sparse.

Methods

Salmonella enterica isolated from blood (n = 60) and cerebrospinal fluid (CSF, n = 3) between 2016 and 2020 from five tertiary hospitals in southwest Nigeria were antimicrobial susceptibility-tested and Illumina-sequenced. Genomes were analysed using publicly-available bioinformatic tools.

Results

Isolates and sequence types (STs) from blood were S. Typhi [ST1, n = 1 and ST2, n = 43] and invasive non-typhoidal Salmonella (iNTS) (S. Enteritidis [ST11, n = 7], S. Durham [ST10, n = 2], S. Rissen [ST8756, n = 2], S. Chester [ST2063, n = 1], S. Dublin [ST10, n = 1], S. Infantis [ST603, n = 1], S. Telelkebir [ST8757, n = 1] and S. Typhimurium [ST313, n = 1]). S. Typhi ST2 (n = 2) and S. Adabraka ST8757 (n = 1) were recovered from CSF. Most S. Typhi belonged to genotype 3.1.1 (n = 44), carried an IncY plasmid, had several antibiotic resistance genes (ARGs) including bla_TEM-1 (n = 38), aph(6)-Id (n = 32), tet(A) (n = 33), sul2 (n = 32), dfrA14 (n = 30) as well as quinolone resistance-conferring gyrA_S83Y single-nucleotide polymorphisms (n = 37). All S. Enteritidis harboured aph(3”)-Ib, bla_TEM-1, catA1, dfrA7, sul1, sul2, tet(B) genes, and a single ARG, qnrB19, was detected in S. Telelkebir. Typhoidal toxins cdtB, pltA and pltB were detected in S. Typhi, Rissen, Chester, and Telelkebir.

Conclusion

Most invasive salmonelloses in southwest Nigeria are vaccine-preventable infections due to multidrug-resistant, West African dominant S. Typhi lineage 3.1.1. Invasive NTS serovars, including some harbouring typhoidal toxin or resistance genes, represented a third of the isolates emphasizing the need for better diagnosis and surveillance.

Whole genome sequencing of 63 invasive Salmonella from 5 tertiary hospitals in Nigeria revealed multiple serovars including a dominant antibiotic-resistance-gene harbouring S. Typhi 3.1.1 genotype comprising a gyrA_S83Y and IncY plasmid. We also report invasive non-typhoidal Salmonella harbouring typhoidal toxins.

Genomic surveillance of Salmonella spp. in the Philippines during 2013-2014

BACKGROUND

Increasing antimicrobial resistance (AMR) in Salmonella has been observed in the Philippines. We aimed to characterise the population and AMR mechanisms of Salmonella with whole genome sequencing (WGS) and compare it with laboratory surveillance methods.

METHODS

The serotype, multilocus sequence type, AMR genes and relatedness between isolates were determined from the genomes of 148 Salmonella Typhi (S. Typhi) and 65 non-typhoidal Salmonella (NTS) collected by the Antimicrobial Resistance Surveillance Program during 2013-2014. Genotypic serotypes and AMR prediction were compared with phenotypic data.

RESULTS

AMR rates in S. Typhi were low, with sparse acquisition of mutations associated with reduced susceptibility to fluoroquinolones or extended-spectrum beta-lactamases (ESBL) genes. By contrast, 75% of NTS isolates were insusceptible to at least one antimicrobial, with more than half carrying mutations and/or genes linked to fluoroquinolone resistance. ESBL genes were detected in five genomes, which also carried other AMR determinants. The population of S. Typhi was dominated by likely endemic genotype 3.0, which caused a putative local outbreak. The main NTS clades were global epidemic S. Enteritidis ST11 and S. Typhimurium monophasic variant (I,4,[5],12: i: -) ST34.

CONCLUSION

We provide the first genomic characterisation of Salmonella from the Philippines and evidence of WGS utility for ongoing surveillance.

Establishing a national reference laboratory for antimicrobial resistance using a whole-genome sequencing framework: Nigeria's experience

Whole-genome sequencing (WGS) is finding important applications in the surveillance of antimicrobial resistance (AMR), providing the most granular data and broadening the scope of niches and locations that can be surveilled. A common but often overlooked application of WGS is to replace or augment reference laboratory services for AMR surveillance. WGS has supplanted traditional strain subtyping in many comprehensive reference laboratories and is now the gold standard for rapidly ruling isolates into or out of suspected outbreak clusters. These and other properties give WGS the potential to serve in AMR reference functioning where a reference laboratory did not hitherto exist. In this perspective, we describe how we have employed a WGS approach, and an academic-public health system collaboration, to provide AMR reference laboratory services in Nigeria, as a model for leapfrogging to national AMR surveillance.

Widespread emergence of OmpK36 loop 3 insertions among multidrug-resistant clones of Klebsiella pneumoniae

Mutations in outer membrane porins act in synergy with carbapenemase enzymes to increase carbapenem resistance in the important nosocomial pathogen, Klebsiella pneumoniae (KP). A key example is a di-amino acid insertion, Glycine-Aspartate (GD), in the extracellular loop 3 (L3) region of OmpK36 which constricts the pore and restricts entry of carbapenems into the bacterial cell. Here we combined genomic and experimental approaches to characterise the diversity, spread and impact of different L3 insertion types in OmpK36. We identified L3 insertions in 3588 (24.1%) of 14,888 KP genomes with an intact ompK36 gene from a global collection. GD insertions were most common, with a high concentration in the ST258/512 clone that has spread widely in Europe and the Americas. Aspartate (D) and Threonine-Aspartate (TD) insertions were prevalent in genomes from Asia, due in part to acquisitions by KP sequence types ST16 and ST231 and subsequent clonal expansions. By solving the crystal structures of novel OmpK36 variants, we found that the TD insertion causes a pore constriction of 41%, significantly greater than that achieved by GD (10%) or D (8%), resulting in the highest levels of resistance to selected antibiotics. We show that in the absence of antibiotics KP mutants harbouring these L3 insertions exhibit both an in vitro and in vivo competitive disadvantage relative to the isogenic parental strain expressing wild type OmpK36. We propose that this explains the reversion of GD and TD insertions observed at low frequency among KP genomes. Finally, we demonstrate that strains expressing L3 insertions remain susceptible to drugs targeting carbapenemase-producing KP, including novel beta lactam-beta lactamase inhibitor combinations. This study provides a contemporary global view of OmpK36-mediated resistance mechanisms in KP, integrating surveillance and experimental data to guide treatment and drug development strategies.

Rapidly rising rates of antibiotic resistance among Klebsiella pneumoniae (KP) necessitate a comprehensive understanding of the diversity, spread and clinical impact of resistance mutations. In KP, mutations in outer membrane porins play an important role in mediating resistance to carbapenems, a key class of antibiotics. Here we show that resistance mutations in the extracellular loop 3 (L3) region of the OmpK36 porin are found at high prevalence among clinical genomes and we characterise their diversity and impact on resistance and virulence. They include amino acid insertions of Aspartate (D), Glycine-Aspartate (GD) and Threonine-Aspartate (TD), which act by decreasing the pore size and restricting entry of carbapenems into the bacterial cell. We show that these L3 insertions are associated with large clonal expansions of resistant lineages and impose a fitness cost evident during in vivo competition. Critically, strains harbouring L3 insertions remain susceptible to novel drugs, including beta lactam-beta lactamase inhibitor combinations. This study highlights the importance of monitoring the emergence and spread of strains with OmpK36 L3 insertions for the control of resistant KP infections and provides crucial data for drug development and treatment strategies.

Rectal Colonization and Nosocomial Transmission of Carbapenem-Resistant Acinetobacter baumannii in an Intensive Care Unit, Southwest Nigeria

Background

Acinetobacter baumannii are of major human health importance because they cause life-threatening nosocomial infections and often are highly resistant to antimicrobials. Specific multidrug-resistant A. baumannii lineages are implicated in hospital outbreaks globally. We retrospectively investigated a suspected outbreak of carbapenem-resistant A. baumannii (CRAB) colonizing patients in an intensive care unit (ICU) of a tertiary hospital in Southwest Nigeria where genomic surveillance of Acinetobacter has hitherto not been conducted.

Methods

A prospective observational study was conducted among all patients admitted to the ICU between August 2017 and June 2018. Acinetobacter species were isolated from rectal swabs and verified phenotypically with the Biomerieux Vitek 2 system. Whole genome sequencing (WGS) was performed on the Illumina platform to characterize isolates from a suspected outbreak during the study period. Phylogenetic analysis, multilocus sequence typing, and antimicrobial resistance gene prediction were carried out in silico.

Results

Acinetobacter isolates belonging to the A. baumannii complex were recovered from 20 (18.5%) ICU patients. Single nucleotide polymorphism (SNP) analysis and epidemiological information revealed a putative outbreak clone comprising seven CRAB strains belonging to the globally disseminated international clone (IC) 2. These isolates had ≤2 SNP differences, identical antimicrobial resistance and virulence genes, and were all ST1114/1841.

Conclusion

We report a carbapenem-resistant IC2 A. baumannii clone causing an outbreak in an ICU in Nigeria. The study findings underscore the need to strengthen the capacity to detect A. baumannii in human clinical samples in Nigeria and assess which interventions can effectively mitigate CRAB transmission in Nigerian hospital settings.

Future-proofing and maximizing the utility of metadata: The PHA4GE SARS-CoV-2 contextual data specification package

Background

The Public Health Alliance for Genomic Epidemiology (PHA4GE) (https://pha4ge.org) is a global coalition that is actively working to establish consensus standards, document and share best practices, improve the availability of critical bioinformatics tools and resources, and advocate for greater openness, interoperability, accessibility, and reproducibility in public health microbial bioinformatics. In the face of the current pandemic, PHA4GE has identified a need for a fit-for-purpose, open-source SARS-CoV-2 contextual data standard.

Results

As such, we have developed a SARS-CoV-2 contextual data specification package based on harmonizable, publicly available community standards. The specification can be implemented via a collection template, as well as an array of protocols and tools to support both the harmonization and submission of sequence data and contextual information to public biorepositories.

Conclusions

Well-structured, rich contextual data add value, promote reuse, and enable aggregation and integration of disparate datasets. Adoption of the proposed standard and practices will better enable interoperability between datasets and systems, improve the consistency and utility of generated data, and ultimately facilitate novel insights and discoveries in SARS-CoV-2 and COVID-19. The package is now supported by the NCBI’s BioSample database.

Context-specific emergence and growth of the SARS-CoV-2 Delta variant

The Delta variant of concern of SARS-CoV-2 has spread globally causing large outbreaks and resurgences of COVID-19 cases 1-3 . The emergence of Delta in the UK occurred on the background of a heterogeneous landscape of immunity and relaxation of non-pharmaceutical interventions 4,5 . Here we analyse 52,992 Delta genomes from England in combination with 93,649 global genomes to reconstruct the emergence of Delta, and quantify its introduction to and regional dissemination across England, in the context of changing travel and social restrictions. Through analysis of human movement, contact tracing, and virus genomic data, we find that the focus of geographic expansion of Delta shifted from India to a more global pattern in early May 2021. In England, Delta lineages were introduced >1,000 times and spread nationally as non-pharmaceutical interventions were relaxed. We find that hotel quarantine for travellers from India reduced onward transmission from importations; however the transmission chains that later dominated the Delta wave in England had been already seeded before restrictions were introduced. In England, increasing inter-regional travel drove Delta's nationwide dissemination, with some cities receiving >2,000 observable lineage introductions from other regions. Subsequently, increased levels of local population mixing, not the number of importations, was associated with faster relative growth of Delta. Among US states, we find that regions that previously experienced large waves also had faster Delta growth rates, and a model including interactions between immunity and human behaviour could accurately predict the rise of Delta there. Delta's invasion dynamics depended on fine scale spatial heterogeneity in immunity and contact patterns and our findings will inform optimal spatial interventions to reduce transmission of current and future VOCs such as Omicron.

Context-specific emergence and growth of the SARS-CoV-2 Delta variant

The Delta variant of concern of SARS-CoV-2 has spread globally causing large outbreaks and resurgences of COVID-19 cases. The emergence of Delta in the UK occurred on the background of a heterogeneous landscape of immunity and relaxation of non-pharmaceutical interventions. Here we analyse 52,992 Delta genomes from England in combination with 93,649 global genomes to reconstruct the emergence of Delta, and quantify its introduction to and regional dissemination across England, in the context of changing travel and social restrictions. Through analysis of human movement, contact tracing, and virus genomic data, we find that the focus of geographic expansion of Delta shifted from India to a more global pattern in early May 2021. In England, Delta lineages were introduced >1,000 times and spread nationally as non-pharmaceutical interventions were relaxed. We find that hotel quarantine for travellers from India reduced onward transmission from importations; however the transmission chains that later dominated the Delta wave in England had been already seeded before restrictions were introduced. In England, increasing inter-regional travel drove Delta's nationwide dissemination, with some cities receiving >2,000 observable lineage introductions from other regions. Subsequently, increased levels of local population mixing, not the number of importations, was associated with faster relative growth of Delta. Among US states, we find that regions that previously experienced large waves also had faster Delta growth rates, and a model including interactions between immunity and human behaviour could accurately predict the rise of Delta there. Delta’s invasion dynamics depended on fine scale spatial heterogeneity in immunity and contact patterns and our findings will inform optimal spatial interventions to reduce transmission of current and future VOCs such as Omicron.

Clones and Clusters of Antimicrobial-Resistant Klebsiella From Southwestern Nigeria

BACKGROUND

Klebsiella pneumoniae is a World Health Organization high-priority antibiotic-resistant pathogen. However, little is known about Klebsiella lineages circulating in Nigeria.

METHODS

We performed whole-genome sequencing (WGS) of 141 Klebsiella isolated between 2016 and 2018 from clinical specimens at 3 antimicrobial-resistance (AMR) sentinel surveillance tertiary hospitals in southwestern Nigeria. We conducted in silico multilocus sequence typing; AMR gene, virulence gene, plasmid, and K and O loci profiling; as well as phylogenetic analyses, using publicly available tools and Nextflow pipelines.

RESULTS

Phylogenetic analysis revealed that the majority of the 134 K. pneumoniae and 5 K. quasipneumoniae isolates from Nigeria characterized are closely related to globally disseminated multidrug-resistant clones. Of the 39 K. pneumoniae sequence types (STs) identified, the most common were ST307 (15%), ST5241 (12%), ST15 (~9%), and ST25 (~6%). ST5241, 1 of 10 novel STs detected, is a single locus variant of ST636 carrying dfrA14, tetD, qnrS, and oqxAB resistance genes. The extended-spectrum β-lactamase (ESBL) gene blaCTX_M-15 was seen in 72% of K. pneumoniae genomes, while 8% encoded a carbapenemase. No isolate carried a combination of carbapenemase-producing genes. Four likely outbreak clusters from 1 facility, within STs 17, 25, 307, and 5241, were ESBL but not carbapenemase-bearing clones.

CONCLUSIONS

This study uncovered known and novel K. pneumoniae lineages circulating in 3 hospitals in Southwest Nigeria that include multidrug-resistant ESBL producers. Carbapenemase-producing isolates remain uncommon. WGS retrospectively identified outbreak clusters, pointing to the value of genomic approaches in AMR surveillance for improving infection prevention and control in Nigerian hospitals.

Genome Sequencing Identifies Previously Unrecognized Klebsiella pneumoniae Outbreaks in Neonatal Intensive Care Units in the Philippines

BACKGROUND

Klebsiella pneumoniae is a critically important pathogen in the Philippines. Isolates are commonly resistant to at least 2 classes of antibiotics, yet mechanisms and spread of its resistance are not well studied.

METHODS

A retrospective sequencing survey was performed on carbapenem-, extended spectrum beta-lactam-, and cephalosporin-resistant Klebsiella pneumoniae isolated at 20 antimicrobial resistance (AMR) surveillance sentinel sites from 2015 through 2017. We characterized 259 isolates using biochemical methods, antimicrobial susceptibility testing, and whole-genome sequencing (WGS). Known AMR mechanisms were identified. Potential outbreaks were investigated by detecting clusters from epidemiologic, phenotypic, and genome-derived data.

RESULTS

Prevalent AMR mechanisms detected include blaCTX-M-15 (76.8%) and blaNDM-1 (37.5%). An epidemic IncFII(Yp) plasmid carrying blaNDM-1 was also detected in 46 isolates from 6 sentinel sites and 14 different sequence types (STs). This plasmid was also identified as the main vehicle of carbapenem resistance in 2 previously unrecognized local outbreaks of ST348 and ST283 at 2 different sentinel sites. A third local outbreak of ST397 was also identified but without the IncFII(Yp) plasmid. Isolates in each outbreak site showed identical STs and K- and O-loci, and similar resistance profiles and AMR genes. All outbreak isolates were collected from blood of children aged < 1 year.

CONCLUSION

WGS provided a better understanding of the epidemiology of multidrug resistant Klebsiella in the Philippines, which was not possible with only phenotypic and epidemiologic data. The identification of 3 previously unrecognized Klebsiella outbreaks highlights the utility of WGS in outbreak detection, as well as its importance in public health and in implementing infection control programs.

Rapid Genomic Characterization and Global Surveillance of Klebsiella Using Pathogenwatch

BACKGROUND

Klebsiella species, including the notable pathogen K. pneumoniae, are increasingly associated with antimicrobial resistance (AMR). Genome-based surveillance can inform interventions aimed at controlling AMR. However, its widespread implementation requires tools to streamline bioinformatic analyses and public health reporting.

METHODS

We developed the web application Pathogenwatch, which implements analytics tailored to Klebsiella species for integration and visualization of genomic and epidemiological data. We populated Pathogenwatch with 16 537 public Klebsiella genomes to enable contextualization of user genomes. We demonstrated its features with 1636 genomes from 4 low- and middle-income countries (LMICs) participating in the NIHR Global Health Research Unit (GHRU) on AMR.

RESULTS

Using Pathogenwatch, we found that GHRU genomes were dominated by a small number of epidemic drug-resistant clones of K. pneumoniae. However, differences in their distribution were observed (eg, ST258/512 dominated in Colombia, ST231 in India, ST307 in Nigeria, ST147 in the Philippines). Phylogenetic analyses including public genomes for contextualization enabled retrospective monitoring of their spread. In particular, we identified hospital outbreaks, detected introductions from abroad, and uncovered clonal expansions associated with resistance and virulence genes. Assessment of loci encoding O-antigens and capsule in K. pneumoniae, which represent possible vaccine candidates, showed that 3 O-types (O1-O3) represented 88.9% of all genomes, whereas capsule types were much more diverse.

CONCLUSIONS

Pathogenwatch provides a free, accessible platform for real-time analysis of Klebsiella genomes to aid surveillance at local, national, and global levels. We have improved representation of genomes from GHRU participant countries, further facilitating ongoing surveillance.

Train-the-Trainer as an Effective Approach to Building Global Networks of Experts in Genomic Surveillance of Antimicrobial Resistance (AMR)

Advanced genomics and sequencing technologies are increasingly becoming critical for global health applications such as pathogen and antimicrobial resistance (AMR) surveillance. Limited resources challenge capacity development in low- and middle-income countries (LMICs), with few countries having genomics facilities and adequately trained staff. Training research and public health experts who are directly involved in the establishment of such facilities offers an effective, but limited, solution to a growing need. Instead, training them to impart their knowledge and skills to others provides a sustainable model for scaling up the much needed capacity and capability for genomic sequencing and analysis locally with global impact. We designed and developed a Train-the-Trainer course integrating pedagogical aspects with genomic and bioinformatics activities. The course was delivered to 18 participants from 12 countries in Africa, Asia, and Latin America. A combination of teaching strategies culminating in a group project created a foundation for continued development at home institutions. Upon follow-up after 6 months, at least 40% of trainees had initiated training programs and collaborations to build capacity at local, national, and regional level. This work provides a framework for implementing a training and capacity building program for the application of genomics tools and resources in AMR surveillance.

Overcoming Data Bottlenecks in Genomic Pathogen Surveillance

Performing whole genome sequencing (WGS) for the surveillance of antimicrobial resistance offers the ability to determine not only the antimicrobials to which rates of resistance are increasing, but also the evolutionary mechanisms and transmission routes responsible for the increase at local, national, and global scales. To derive WGS-based outputs, a series of processes are required, beginning with sample and metadata collection, followed by nucleic acid extraction, library preparation, sequencing, and analysis. Throughout this pathway there are many data-related operations required (informatics) combined with more biologically focused procedures (bioinformatics). For a laboratory aiming to implement pathogen genomics, the informatics and bioinformatics activities can be a barrier to starting on the journey; for a laboratory that has already started, these activities may become overwhelming. Here we describe these data bottlenecks and how they have been addressed in laboratories in India, Colombia, Nigeria, and the Philippines, as part of the National Institute for Health Research Global Health Research Unit on Genomic Surveillance of Antimicrobial Resistance. The approaches taken include the use of reproducible data parsing pipelines and genome sequence analysis workflows, using technologies such as Data-flo, the Nextflow workflow manager, and containerization of software dependencies. By overcoming barriers to WGS implementation in countries where genome sampling for some species may be underrepresented, a body of evidence can be built to determine the concordance of antimicrobial sensitivity testing and genome-derived resistance, and novel high-risk clones and unknown mechanisms of resistance can be discovered.

Implementing Whole-Genome Sequencing for Ongoing Surveillance of Antimicrobial Resistance: Exemplifying Insights Into Klebsiella pneumoniae

In this Supplement, we detail outputs of the National Institute for Health Research Global Health Research Unit on Genomic Surveillance of Antimicrobial Resistance project, covering practical implementation of whole-genome sequencing across our consortium, which consists of laboratories in Colombia, India, Nigeria, and the Philippines.

Complexity of Genomic Epidemiology of Carbapenem-Resistant Klebsiella pneumoniae Isolates in Colombia Urges the Reinforcement of Whole Genome Sequencing-Based Surveillance Programs

BACKGROUND

Carbapenem-resistant Klebsiella pneumoniae (CRKP) is an emerging public health problem. This study explores the specifics of CRKP epidemiology in Colombia based on whole genome sequencing (WGS) of the National Reference Laboratory at Instituto Nacional de Salud (INS)'s 2013-2017 sample collection.

METHODS

A total of 425 CRKP isolates from 21 departments were analyzed by HiSeq-X10®Illumina high-throughput sequencing. Bioinformatic analysis was performed, primarily using the pipelines developed collaboratively by the National Institute for Health Research Global Health Research Unit (GHRU) on Genomic Surveillance of Antimicrobial Resistance (AMR), and AGROSAVIA.

RESULTS

Of the 425 CRKP isolates, 91.5% were carbapenemase-producing strains. The data support a recent expansion and the endemicity of CRKP in Colombia with the circulation of 7 high-risk clones, the most frequent being CG258 (48.39% of isolates). We identified genes encoding carbapenemases blaKPC-3, blaKPC-2, blaNDM-1, blaNDM-9, blaVIM-2, blaVIM-4, and blaVIM-24, and various mobile genetic elements (MGE). The virulence of CRKP isolates was low, but colibactin (clb3) was present in 25.2% of isolates, and a hypervirulent CRKP clone (CG380) was reported for the first time in Colombia. ST258, ST512, and ST4851 were characterized by low levels of diversity in the core genome (ANI > 99.9%).

CONCLUSIONS

The study outlines complex CRKP epidemiology in Colombia. CG258 expanded clonally and carries specific carbapenemases in specific MGEs, while the other high-risk clones (CG147, CG307, and CG152) present a more diverse complement of carbapenemases. The specifics of the Colombian situation stress the importance of WGS-based surveillance to monitor evolutionary trends of sequence types (STs), MGE, and resistance and virulence genes.

High-Resolution Genomic Profiling of Carbapenem-Resistant Klebsiella pneumoniae Isolates: A Multicentric Retrospective Indian Study

BACKGROUND

Carbapenem-resistant Klebsiella pneumoniae (CRKP) is a threat to public health in India because of its high dissemination, mortality, and limited treatment options. Its genomic variability is reflected in the diversity of sequence types, virulence factors, and antimicrobial resistance (AMR) mechanisms. This study aims to characterize the clonal relationships and genetic mechanisms of resistance and virulence in CRKP isolates in India.

MATERIALS AND METHODS

We characterized 344 retrospective K. pneumoniae clinical isolates collected from 8 centers across India collected in 2013-2019. Susceptibility to antibiotics was tested with VITEK 2. Capsular types, multilocus sequence type, virulence genes, AMR determinants, plasmid replicon types, and a single-nucleotide polymorphism phylogeny were inferred from their whole genome sequences.

RESULTS

Phylogenetic analysis of the 325 Klebsiella isolates that passed quality control revealed 3 groups: K. pneumoniae sensu stricto (n = 307), K. quasipneumoniae (n = 17), and K. variicola (n = 1). Sequencing and capsular diversity analysis of the 307 K. pneumoniae sensu stricto isolates revealed 28 sequence types, 26 K-locus types, and 11 O-locus types, with ST231, KL51, and O1V2 being predominant. blaOXA-48-like and blaNDM-1/5 were present in 73.2% and 24.4% of isolates, respectively. The major plasmid replicon types associated with carbapenase genes were IncF (51.0%) and Col group (35.0%).

CONCLUSION

Our study documents for the first time the genetic diversity of K and O antigens circulating in India. The results demonstrate the practical applicability of genomic surveillance and its utility in tracking the population dynamics of CRKP. It alerts us to the urgency for longitudinal surveillance of these transmissible lineages.

Integrating Scalable Genome Sequencing Into Microbiology Laboratories for Routine Antimicrobial Resistance Surveillance

Antimicrobial resistance (AMR) is considered a global threat, and novel drug discovery needs to be complemented with systematic and standardized epidemiological surveillance. Surveillance data are currently generated using phenotypic characterization. However, due to poor scalability, this approach does little for true epidemiological investigations. There is a strong case for whole-genome sequencing (WGS) to enhance the phenotypic data. To establish global AMR surveillance using WGS, we developed a laboratory implementation approach that we applied within the NIHR Global Health Research Unit (GHRU) on Genomic Surveillance of Antimicrobial Resistance. In this paper, we outline the laboratory implementation at 4 units: Colombia, India, Nigeria, and the Philippines. The journey to embedding WGS capacity was split into 4 phases: Assessment, Assembly, Optimization, and Reassessment. We show that on-boarding WGS capabilities can greatly enhance the real-time processing power within regional and national AMR surveillance initiatives, despite the high initial investment in laboratory infrastructure and maintenance. Countries looking to introduce WGS as a surveillance tool could begin by sequencing select Global Antimicrobial Resistance Surveillance System (GLASS) priority pathogens that can demonstrate the standardization and impact genome sequencing has in tackling AMR.

Good Financial Grant Practice: A Tool for Developing and Demonstrating Institutional Financial and Grant Management Capacity in Global Health

The administration and governance of grant funding across global health organizations presents enormous challenges. Meeting these challenges is crucial to ensuring that funds are used in the most effective way to improve health outcomes, in line with the United Nations' Sustainable Development Goal 3, "Ensure healthy lives and promote well-being for all at all ages." The Good Financial Grant Practice (GFGP) Standard (ARS 1651) is the world's first and, currently, only international standard for the financial governance and management of grant funding. Through consensus building and global harmonization between both low- and middle-income and high-income country players, the GFGP Standard has achieved a leveling impact: GFGP applies equally to, and can be implemented by, all types of organization, regardless of location, size, or whether they predominantly give or receive funding. GFGP can be used as a tool for addressing some of the challenges of the current funding model. Here, we describe our experiences and lessons learned from implementing GFGP across 4 diverse research institutions in India, Nigeria, Colombia, and the Philippines as part of our National Institute for Health Research Global Health Research Unit on Genomic Surveillance of Antimicrobial Resistance.

Circulation of third-generation cephalosporin resistant Salmonella Typhi in Mumbai, India

We report the persistent circulation of third-generation cephalosporin resistant Salmonella Typhi in Mumbai, linked to the acquisition and maintenance of a previously characterized IncX3 plasmid carrying the ESBL gene bla_SHV-12 and the fluoroquinolone resistance gene qnrB7 in the genetic context of a triple mutant also associated with fluoroquinolone resistance.

High-resolution sweep metagenomics using fast probabilistic inference

Determining the composition of bacterial communities beyond the level of a genus or species is challenging because of the considerable overlap between genomes representing close relatives. Here, we present the mSWEEP pipeline for identifying and estimating the relative sequence abundances of bacterial lineages from plate sweeps of enrichment cultures. mSWEEP leverages biologically grouped sequence assembly databases, applying probabilistic modelling, and provides controls for false positive results. Using sequencing data from major pathogens, we demonstrate significant improvements in lineage quantification and detection accuracy. Our pipeline facilitates investigating cultures comprising mixtures of bacteria, and opens up a new field of plate sweep metagenomics.

Genomic surveillance of Acinetobacter baumannii in the Philippines, 2013-2014

OBJECTIVE

Acinetobacter baumannii is an opportunistic nosocomial pathogen that has increasingly become resistant to carbapenems worldwide. In the Philippines, rates of carbapenem resistance and multidrug resistance are above 50%. We undertook a genomic study of carbapenem-resistant A. baumannii in the Philippines to characterize the population diversity and antimicrobial resistance mechanisms.

METHODS

We sequenced the whole genomes of 117 A. baumannii isolates recovered by 16 hospitals in the Philippines between 2013 and 2014. From the genome sequences, we determined the multilocus sequence type, presence of acquired determinants of antimicrobial resistance and relatedness between isolates. We also compared the phenotypic and genotypic resistance results.

RESULTS

Carbapenem resistance was mainly explained by acquisition of the class-D β-lactamase gene blaOXA-23. The concordance between phenotypic and genotypic resistance to imipenem was 98.15%, and it was 94.97% overall for the seven antibiotics analysed. Twenty-two different sequence types were identified, including 7 novel types. The population was dominated by the high-risk international clone 2 (i.e. clonal complex 92), in particular by ST195 and ST208 and their single locus variants. Using whole-genome sequencing, we identified local clusters representing potentially undetected nosocomial outbreaks, as well as multihospital clusters that indicated interhospital dissemination. Comparison with global genomes suggested that the establishment of carbapenem-resistant international clone 2 in the Philippines is likely the result of clonal expansion and geographical dissemination, and at least partly explained by inadequate hospital infection control and prevention.

DISCUSSION

This is the first extensive genomic study of carbapenem-resistant A. baumannii in the Philippines, and it underscores the importance of hospital infection control and prevention measures to contain high-risk clones.

Spatiotemporal invasion dynamics of SARS-CoV-2 lineage B.1.1.7 emergence

Understanding the causes and consequences of the emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern is crucial to pandemic control yet difficult to achieve because they arise in the context of variable human behavior and immunity. We investigated the spatial invasion dynamics of lineage B.1.1.7 by jointly analyzing UK human mobility, virus genomes, and community-based polymerase chain reaction data. We identified a multistage spatial invasion process in which early B.1.1.7 growth rates were associated with mobility and asymmetric lineage export from a dominant source location, enhancing the effects of B.1.1.7's increased intrinsic transmissibility. We further explored how B.1.1.7 spread was shaped by nonpharmaceutical interventions and spatial variation in previous attack rates. Our findings show that careful accounting of the behavioral and epidemiological context within which variants of concern emerge is necessary to interpret correctly their observed relative growth rates.

Assignment of epidemiological lineages in an emerging pandemic using the pangolin tool

The response of the global virus genomics community to the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has been unprecedented, with significant advances made towards the 'real-time' generation and sharing of SARS-CoV-2 genomic data. The rapid growth in virus genome data production has necessitated the development of new analytical methods that can deal with orders of magnitude of more genomes than previously available. Here, we present and describe Phylogenetic Assignment of Named Global Outbreak Lineages (pangolin), a computational tool that has been developed to assign the most likely lineage to a given SARS-CoV-2 genome sequence according to the Pango dynamic lineage nomenclature scheme. To date, nearly two million virus genomes have been submitted to the web-application implementation of pangolin, which has facilitated the SARS-CoV-2 genomic epidemiology and provided researchers with access to actionable information about the pandemic's transmission lineages.

CLIMB-COVID: continuous integration supporting decentralised sequencing for SARS-CoV-2 genomic surveillance

In response to the ongoing SARS-CoV-2 pandemic in the UK, the COVID-19 Genomics UK (COG-UK) consortium was formed to rapidly sequence SARS-CoV-2 genomes as part of a national-scale genomic surveillance strategy. The network consists of universities, academic institutes, regional sequencing centres and the four UK Public Health Agencies. We describe the development and deployment of CLIMB-COVID, an encompassing digital infrastructure to address the challenge of collecting and integrating both genomic sequencing data and sample-associated metadata produced across the COG-UK network.

The changing epidemiology of carbapenemase-producing Klebsiella pneumoniae in Italy: toward polyclonal evolution with emergence of high-risk lineages

BACKGROUND

Previous studies showed that the epidemic of carbapenem-resistant Klebsiella pneumoniae (CR-KP) observed in Italy since 2010 was sustained mostly by strains of clonal group (CG) 258 producing KPC-type carbapenemases. In the framework of the National Antibiotic-Resistance Surveillance (AR-ISS), a countrywide survey was conducted in 2016 to explore the evolution of the phenotypic and genotypic characteristics of CR-KP isolates.

METHODS

From March to July 2016, hospital laboratories participating in AR-ISS were requested to provide consecutive, non-duplicated CR-KP (meropenem and/or imipenem MIC >1 mg/L) from invasive infections. Antibiotic susceptibility was determined according to EUCAST recommendations. A WGS approach was adopted to characterize the isolates by investigating phylogeny, resistome and virulome.

RESULTS

Twenty-four laboratories provided 157 CR-KP isolates, of which 156 were confirmed as K. pneumoniae sensu stricto by WGS and found to carry at least one carbapenemase-encoding gene, corresponding in most cases (96.1%) to blaKPC. MLST- and SNP-based phylogeny revealed that 87.8% of the isolates clustered in four major lineages: CG258 (47.4%), with ST512 as the most common clone, CG307 (19.9%), ST101 (15.4%) and ST395 (5.1%). A close association was identified between lineages and antibiotic resistance phenotypes and genotypes, virulence traits and capsular types. Colistin resistance, mainly associated with mgrB mutations, was common in all major lineages except ST395.

CONCLUSIONS

This WGS-based survey showed that, although CG258 remained the most common CR-KP lineage in Italy, a polyclonal population has emerged with the spread of the new high-risk lineages CG307, ST101 and ST395, while KPC remained the most common carbapenemase.

A community-driven resource for genomic epidemiology and antimicrobial resistance prediction of Neisseria gonorrhoeae at Pathogenwatch

BACKGROUND

Antimicrobial-resistant (AMR) Neisseria gonorrhoeae is an urgent threat to public health, as strains resistant to at least one of the two last-line antibiotics used in empiric therapy of gonorrhoea, ceftriaxone and azithromycin, have spread internationally. Whole genome sequencing (WGS) data can be used to identify new AMR clones and transmission networks and inform the development of point-of-care tests for antimicrobial susceptibility, novel antimicrobials and vaccines. Community-driven tools that provide an easy access to and analysis of genomic and epidemiological data is the way forward for public health surveillance.

METHODS

Here we present a public health-focussed scheme for genomic epidemiology of N. gonorrhoeae at Pathogenwatch ( https://pathogen.watch/ngonorrhoeae ). An international advisory group of experts in epidemiology, public health, genetics and genomics of N. gonorrhoeae was convened to inform on the utility of current and future analytics in the platform. We implement backwards compatibility with MLST, NG-MAST and NG-STAR typing schemes as well as an exhaustive library of genetic AMR determinants linked to a genotypic prediction of resistance to eight antibiotics. A collection of over 12,000 N. gonorrhoeae genome sequences from public archives has been quality-checked, assembled and made public together with available metadata for contextualization.

RESULTS

AMR prediction from genome data revealed specificity values over 99% for azithromycin, ciprofloxacin and ceftriaxone and sensitivity values around 99% for benzylpenicillin and tetracycline. A case study using the Pathogenwatch collection of N. gonorrhoeae public genomes showed the global expansion of an azithromycin-resistant lineage carrying a mosaic mtr over at least the last 10 years, emphasising the power of Pathogenwatch to explore and evaluate genomic epidemiology questions of public health concern.

CONCLUSIONS

The N. gonorrhoeae scheme in Pathogenwatch provides customised bioinformatic pipelines guided by expert opinion that can be adapted to public health agencies and departments with little expertise in bioinformatics and lower-resourced settings with internet connection but limited computational infrastructure. The advisory group will assess and identify ongoing public health needs in the field of gonorrhoea, particularly regarding gonococcal AMR, in order to further enhance utility with modified or new analytic methods.

Population genomics and antimicrobial resistance dynamics of Escherichia coli in wastewater and river environments

Aquatic environments are key niches for the emergence, evolution and dissemination of antimicrobial resistance. However, the population diversity and the genetic elements that drive the dynamics of resistant bacteria in different aquatic environments are still largely unknown. The aim of this study was to understand the population genomics and evolutionary events of Escherichia coli resistant to clinically important antibiotics including aminoglycosides, in anthropogenic and natural water ecosystems. Here we show that less different E. coli sequence types (STs) are identified in wastewater than in rivers, albeit more resistant to antibiotics, and with significantly more plasmids/cell (6.36 vs 3.72). However, the genomic diversity within E. coli STs in both aquatic environments is similar. Wastewater environments favor the selection of conserved chromosomal structures associated with diverse flexible plasmids, unraveling promiscuous interplasmidic resistance genes flux. On the contrary, the key driver for river E. coli adaptation is a mutable chromosome along with few plasmid types shared between diverse STs harboring a limited resistance gene content.

Genomic Epidemiology of CC30 Methicillin-Resistant Staphylococcus aureus Strains from Argentina Reveals Four Major Clades with Distinctive Genetic Features

Staphylococcus aureus clonal complex 30 (CC30) has given rise to epidemics worldwide and is one of the most prevalent lineages in Argentina, represented by sequence type 30 methicillin-resistant S. aureus SCCmec type IV (ST30-MRSA-IV). ST30-MRSA-IV has displaced previous prevalent clones in the country and demonstrated increased virulence. Despite the burden of infections caused by ST30-MRSA-IV both in hospitals and in communities in Argentina, no detailed genome-based characterization of this clone is available to date. In this study, we used whole-genome sequencing (WGS) to evaluate the genetic diversity, population structure, and genomic characteristics of 190 CC30-MRSA strains circulating in Argentina between 2004 and 2015. Phylogenetic analysis revealed the existence of 4 major clades: ARG-1 (CC30-MRSA-IVc-spa t012), ARG-2 (ST30-MRSA-IVc-spa t021 related), ARG-3 (ST30-MRSA-IVh/j-spa t021 and related), and ARG-4 (CC30-MRSA-IVc-spa t019 and related). The clades were characterized by different distributions of antimicrobial resistance determinants, virulence genes, and mobile genetic elements (MGEs). While ARG-1 and ARG-4 were related to global epidemic MRSA-16 (EMRSA-16) and South West Pacific (SWP) clones, respectively, ARG-3 was phylogenetically distinct from previously defined CC30 epidemic clones. ARG-4, the most prevalent and geographically disseminated in the collection (N = 164), was characterized by specific MGEs and chromosomal mutations that might have contributed to its virulence and success. To our knowledge, this is the first genomic epidemiology study of CC30-MRSA in Argentina, which will serve as baseline genomic data going forward to inform public health measures for infection prevention and control.IMPORTANCE The rise in prevalence of community-associated methicillin-resistant Staphylococcus aureus (CA-MRSA) is of public health concern. In Argentina, several studies documented a shift in the epidemiology of CA-MRSA since 2009, with clonal complex 30 (CC30) and, in particular, sequence type 30 MRSA SCCmec type IV (ST30-MRSA-IV) replacing other clones both in communities and in hospitals and possibly displaying increased virulence. By sequencing the whole genomes of 190 CC30 MRSA isolates recovered from Argentina between 2005 and 2015, we showed that they represented a diverse population composed of 4 major clades. The predominant clade evolved from the South West Pacific clone but has acquired a distinct repertoire of mobile genetic elements, virulence genes, and chromosomal mutations that might play a role in its success. Our work is the first extensive genomic study of CC30 S. aureus in Argentina and will contribute not only to the development of genomic surveillance in the region but also to our understanding of the global epidemiology of this pathogen.

Protocol for an interdisciplinary cross-sectional study investigating the social, biological and community-level drivers of antimicrobial resistance (AMR): Holistic Approach to Unravel Antibacterial Resistance in East Africa (HATUA)

INTRODUCTION

Antimicrobial resistance (AMR) is a global health threat that requires urgent research using a multidisciplinary approach. The biological drivers of AMR are well understood, but factors related to treatment seeking and the social contexts of antibiotic (AB) use behaviours are less understood. Here we describe the Holistic Approach to Unravel Antibacterial Resistance in East Africa, a multicentre consortium that investigates the diverse drivers of drug resistance in urinary tract infections (UTIs) in East Africa.

METHODS AND ANALYSIS

This study will take place in Uganda, Kenya and Tanzania. We will conduct geospatial mapping of AB sellers, and conduct mystery client studies and in-depth interviews (IDIs) with drug sellers to investigate AB provision practices. In parallel, we will conduct IDIs with doctors, alongside community focus groups. Clinically diagnosed patients with UTI will be recruited from healthcare centres, provide urine samples and complete a questionnaire capturing retrospective treatment pathways, sociodemographic characteristics, attitudes and knowledge. Bacterial isolates from urine and stool samples will be subject to culture and antibiotic sensitivity testing. Genomic DNA from bacterial isolates will be extracted with a subset being sequenced. A follow-up household interview will be conducted with 1800 UTI-positive patients, where further environmental samples will be collected. A subsample of patients will be interviewed using qualitative tools. Questionnaire data, microbiological analysis and qualitative data will be linked at the individual level. Quantitative data will be analysed using statistical modelling, including Bayesian network analysis, and all forms of qualitative data analysed through iterative thematic content analysis.

ETHICS AND DISSEMINATION

Approvals have been obtained from all national and local ethical review bodies in East Africa and the UK. Results will be disseminated in communities, with local and global policy stakeholders, and in academic circles. They will have great potential to inform policy, improve clinical practice and build regional pathogen surveillance capacity.