Advancing pathogen genomics in resource-limited settings.

Genomic sequencing has emerged as a powerful tool to enhance early pathogen detection and characterization with implications for public health and clinical decision making. Although widely available in developed countries, the application of pathogen genomics among low-resource, high-disease burden settings remains at an early stage. In these contexts, tailored approaches for integrating pathogen genomics within infectious disease control programs will be essential to optimize cost efficiency and public health impact. We propose a framework for embedding pathogen genomics within national surveillance plans across a spectrum of surveillance and laboratory capacities. We adopt a public health approach to genomics and examine its application to high-priority diseases relevant in resource-limited settings. For each grouping, we assess the value proposition for genomics to inform public health and clinical decision-making, alongside its contribution toward research and development of novel diagnostics, therapeutics, and vaccines.

Exploring programmatic indicators of tuberculosis control that incorporate routine Mycobacterium tuberculosis sequencing in low incidence settings: a comprehensive (2017-2021) patient cohort analysis.

Background: Routine whole genome sequencing of Mycobacterium tuberculosis has been implemented with increasing frequency. However, its value for tuberculosis (TB) control programs beyond individual case management and enhanced drug resistance detection has not yet been explored. Methods: We analysed routine sequencing data of culture-confirmed TB cases notified between 1st January 2017 and 31st December 2021 in New South Wales (NSW), Australia. Genomic surveillance included evidence of local TB transmission, defined by single nucleotide polymorphism (SNP) clustering over a variable (0-25) SNP threshold, and drug resistance conferring mutations. Findings: M. tuberculosis sequences from 1831 patients were examined, representing 64.8% of all notified TB cases and 96.2% of culture-confirmed cases. Applying a traditional 5-SNP cluster threshold identified 62 transmission clusters with 183 clustered cases; 101/183 (55.2%) had 0 SNP differences. Cluster assessment over a 5-year period, using a 5-SNP threshold, provided a comprehensive overview of likely recent transmission within NSW, Australia, as an indicator of local TB control. Genotypic drug susceptibility testing (DST) was highly concordant with phenotypic DST and provided a 6.8% increase in antimycobacterial resistance detection. Importantly, it detected mutations missed by routine molecular tests. Lineage 2 strains were more likely to be drug resistant (p < 0.0001) and locally transmitted if drug resistant (p < 0.0001). Interpretation: Performing routine prospective WGS in a low incidence country like Australia, provides genomically informed programmatic indicators of local TB control. A rolling 5-year cluster assessment reflects epidemic containment and progress towards 'zero TB transmission'. Genomic DST also provides valuable information for clinical care and drug resistance surveillance. Funding: NHMRC Centre for Research Excellence in Tuberculosis (www.tbcre.org.au) and NSW Health Prevention Research Support Program.

Genomic markers of drug resistance in Mycobacterium tuberculosis populations with minority variants.

Minority variants of Mycobacterium tuberculosis harboring mutations conferring resistance can become dominant populations during tuberculosis (TB) treatment, leading to treatment failure. Our understanding of drug-resistant within-host subpopulations and the frequency of resistance-conferring mutations in minority variants remains limited. M. tuberculosis sequences recovered from liquid cultures of culture-confirmed TB cases notified between January 2017 and December 2021 in New South Wales, Australia were examined. Potential drug resistance-conferring minority variants were identified using LoFreq, and mixed populations of different M. tuberculosis strains (≥100 SNPs apart) were examined using QuantTB. A total of 1831 routinely sequenced M. tuberculosis strains were included in the analysis. Drug resistance-conferring minority variants were detected in 3.5% (65/1831) of sequenced cultures; 84.6% (55/65) had majority strains that were drug susceptible and 15.4% (10/65) had majority strains that were drug resistant. Minority variants with high-confidence drug resistance-conferring mutations were 1.5 times more common when the majority strains were drug resistant. Mixed M. tuberculosis strain populations were documented in 10.0% (183/1831) of specimens. Minority variants with high-confidence drug resistance-conferring mutations were more frequently detected in mixed M. tuberculosis strain populations (2.7%, 5/183) than in single strain populations (0.6%, 10/1648; P = 0.01). Drug-resistant minority variants require monitoring in settings that implement routine M. tuberculosis sequencing. The frequency with which drug-resistant minority variants are detected is likely influenced by pre-culture requirement. Culture-independent sequencing methods should provide a more accurate reflection of drug-resistant subpopulations.

Persistent Salmonella enterica serovar Typhi sub-populations within host interrogated by whole genome sequencing and metagenomics.

Salmonella enterica serovar Typhi (S. Typhi) causes typhoid fever and, in some cases, chronic carriage after resolution of acute disease. This study examined sequential isolates of S. Typhi from a single host with persistent asymptomatic infection. These isolates, along with another S. Typhi isolate recovered from a household contact with typhoid fever, were subjected to whole genome sequencing and analysis. In addition, direct sequencing of the bile fluid from the host with persistent infection was also performed. Comparative analysis of isolates revealed three sub-populations of S. Typhi with distinct genetic patterns. Metagenomic sequencing recognised only two of the three sub-populations within the bile fluid. The detection and investigation of insertion sequences IS10R and associated deletions complemented analysis of single nucleotide polymorphisms. These findings improve our understanding of within-host dynamics of S. Typhi in cases of persistent infection and inform epidemiological investigations of transmission events associated with chronic carriers.

SABRes: in silico detection of drug resistance conferring mutations in subpopulations of SARS-CoV-2 genomes.

The emergence of resistance to antiviral drugs increasingly used to treat SARS-CoV-2 infections has been recognised as a significant threat to COVID-19 control. In addition, some SARS-CoV-2 variants of concern appear to be intrinsically resistant to several classes of these antiviral agents. Therefore, there is a critical need for rapid recognition of clinically relevant polymorphisms in SARS-CoV-2 genomes associated with significant reduction of drug activity in virus neutralisation experiments. Here we present SABRes, a bioinformatic tool, which leverages on expanding public datasets of SARS-CoV-2 genomes and allows detection of drug resistance mutations in consensus genomes as well as in viral subpopulations. We have applied SABRes to detect resistance-conferring mutations in 25,197 genomes generated over the course of the SARS-CoV-2 pandemic in Australia and identified 299 genomes containing resistance conferring mutations to the five antiviral therapeutics that retain effectiveness against currently circulating strains of SARS-CoV-2 - Sotrovimab, Bebtelovimab, Remdesivir, Nirmatrelvir and Molnupiravir. These genomes accounted for a 1.18% prevalence of resistant isolates discovered by SABRes, including 80 genomes with resistance conferring mutations found in viral subpopulations. Timely recognition of these mutations within subpopulations is critical as these mutations can provide an advantage under selective pressure and presents an important step forward in our ability to monitor SARS-CoV-2 drug resistance.

The evolution and international spread of extensively drug resistant Shigella sonnei.

Shigella sonnei causes shigellosis, a severe gastrointestinal illness that is sexually transmissible among men who have sex with men (MSM). Multidrug resistance in S. sonnei is common including against World Health Organisation recommended treatment options, azithromycin, and ciprofloxacin. Recently, an MSM-associated outbreak of extended-spectrum ß-lactamase producing, extensively drug resistant S. sonnei was reported in the United Kingdom. Here, we aimed to identify the genetic basis, evolutionary history, and international dissemination of the outbreak strain. Our genomic epidemiological analyses of 3,304 isolates from the United Kingdom, Australia, Belgium, France, and the United States of America revealed an internationally connected outbreak with a most recent common ancestor in 2018 carrying a low-fitness cost resistance plasmid, previously observed in travel associated sublineages of S. flexneri. Our results highlight the persistent threat of horizontally transmitted antimicrobial resistance and the value of continuing to work towards early and open international sharing of genomic surveillance data.

Pangenome Analysis of a Salmonella Enteritidis Population Links a Major Outbreak to a Gifsy-1-Like Prophage Containing Anti-Inflammatory Gene gogB.

A major outbreak of the globally significant Salmonella Enteritidis foodborne pathogen was identified within a large clinical data set by a program of routine WGS of clinical presentations of salmonellosis in New South Wales, Australia. Pangenome analysis helped to quantify and isolate prophage content within the accessory partition of the pangenome. A prophage similar to Gifsy-1 (henceforth GF-1L) was found to occur in all isolates of the outbreak core SNP cluster, and in three other isolates. Further analysis revealed that the GF-1L prophage carried the gogB virulence factor. These observations suggest that GF-1L may be an important marker of virulence for S. Enteritidis population screening and, that anti-inflammatory, gogB-mediated virulence currently associated with Salmonella Typhimurium may also be displayed by S. Enteritidis. IMPORTANCE We examined 5 years of genomic and epidemiological data for the significant global foodborne pathogen, Salmonella enterica. Although Salmonella enterica subspecies enterica serovar Enteritidis (S. Enteritidis) is the leading cause of salmonellosis in the USA and Europe, prior to 2018 it was not endemic in the southern states of Australia. However, in 2018 a large outbreak led to the endemicity of S. Enteritidis in New South Wales, Australia, and a unique opportunity to study this phenomenon. Using pangenome analysis we uncovered that this clone contained a Gifsy-1-like prophage harboring the known virulence factor gogB. The prophage reported has not previously been described in S. Enteritidis isolates.

Genomes of Vibrio cholerae O1 Serotype Ogawa Associated with Current Cholera Activity in Pakistan.

Here, this report presents two genomes of Vibrio cholerae O1 serotype Ogawa, recovered from cholera cases in Australia linked to travel to Pakistan in 2022. Their multidrug-resistant genotype represents the current activity of cholera within the seventh pandemic. One of the genome sequences was assembled using both short- and long-read sequences.

Case report: a genomics-guided reclassification of a blood culture isolate misassigned by MALDI-TOF as Yersinia pestis.

In this report, we describe a case where Gram-negative rods were isolated from a blood culture which subsequently presented a discordant Yersinia species result by MALDI-TOF. Rapid sequencing provided high-resolution identification of the isolate as Yersinia pseudotuberculosis , which was subsequently confirmed by biochemical tests.

Genomic and transcriptomic variation in Bordetella spp. following induction of erythromycin resistance.

BACKGROUND: The emergence of macrolide resistance in Bordetella pertussis, the causative agent of pertussis, due to mutations in the 23S rRNA gene has been recently recognized. However, resistance mechanisms to macrolides in Bordetella parapertussis and Bordetella holmesii remain unknown. OBJECTIVES: This study investigated genomic changes induced by in vitro exposure to erythromycin in these three main pathogens responsible for pertussis-like disease. METHODS: A set of 10 clinical and reference strains of B. pertussis, B. parapertussis and B. holmesii was exposed to erythromycin for 15 weeks or 30 subculture passages. Antibiotic pressure was achieved by growth on the selective media with erythromycin Etest strips or impregnated discs. Genome polymorphisms and transcriptomic profiles were examined by short- and long-read sequencing of passaged isolates. RESULTS: B. parapertussis and B. holmesii isolates developed significant in vitro resistance to erythromycin (MIC >256 mg/L) within 2 to 7 weeks and at 5 to 12 weeks, respectively. B. pertussis remained phenotypically susceptible to the antibiotic following 15 weeks of exposure, with the MIC between 0.032 to 0.38 mg/L. Genomic analysis revealed that B. holmesii developed resistance due to mutations in the 23S rRNA gene. The resistance mechanism in B. parapertussis was hypothesized as being due to upregulation of an efflux pump mechanism. CONCLUSIONS: These findings indicate that both B. holmesii and B. parapertussis can be more prone to induced resistance following exposure to treatment with erythromycin than B. pertussis. The surveillance of macrolide resistance in Bordetella isolates recovered from patients with pertussis, especially persistent disease, is warranted.

Co-infection with SARS-CoV-2 Omicron and Delta variants revealed by genomic surveillance.

Co-infections with different variants of SARS-CoV-2 are a key precursor to recombination events that are likely to drive SARS-CoV-2 evolution. Rapid identification of such co-infections is required to determine their frequency in the community, particularly in populations at-risk of severe COVID-19, which have already been identified as incubators for punctuated evolutionary events. However, limited data and tools are currently available to detect and characterise the SARS-CoV-2 co-infections associated with recognised variants of concern. Here we describe co-infection with the SARS-CoV-2 variants of concern Omicron and Delta in two epidemiologically unrelated adult patients with chronic kidney disease requiring maintenance haemodialysis. Both variants were co-circulating in the community at the time of detection. Genomic surveillance based on amplicon- and probe-based sequencing using short- and long-read technologies identified and quantified subpopulations of Delta and Omicron viruses in respiratory samples. These findings highlight the importance of integrated genomic surveillance in vulnerable populations and provide diagnostic pathways to recognise SARS-CoV-2 co-infection using genomic data.

Microbial Genomics as a Catalyst for Targeted Antivirulence Therapeutics.

Virulence arresting drugs (VAD) are an expanding class of antimicrobial treatment that act to "disarm" rather than kill bacteria. Despite an increasing number of VAD being registered for clinical use, uptake is hampered by the lack of methods that can identify patients who are most likely to benefit from these new agents. The application of pathogen genomics can facilitate the rational utilization of advanced therapeutics for infectious diseases. The development of genomic assessment of VAD targets is essential to support the early stages of VAD diffusion into infectious disease management. Genomic identification and characterization of VAD targets in clinical isolates can augment antimicrobial stewardship and pharmacovigilance. Personalized genomics guided use of VAD will provide crucial policy guidance to regulating agencies, assist hospitals to optimize the use of these expensive medicines and create market opportunities for biotech companies and diagnostic laboratories.

Added Value of Genomic Surveillance of Virulence Factors in Shiga Toxin-Producing Escherichia coli in New South Wales, Australia.

The disease caused by Shiga toxin-producing Escherichia coli (STEC) remains a significant public health challenge globally, but the incidence of human STEC infections in Australia remains relatively low. This study examined the virulence characteristics and diversity of STEC isolates in the state of New South Wales between December 2017 and May 2020. Utilisation of both whole and core genome multi-locus sequence typing (MLST) allowed for the inference of genomic diversity and detection of isolates that were likely to be epidemiologically linked. The most common STEC serotype and stx subtype detected in this study were O157:H7 and stx 1a, respectively. A genomic scan of other virulence factors present in STEC suggested interplay between iron uptake system and virulence factors that mediate either iron release or countermeasures against host defence that could result in a reduction of stx 1a expression. This reduced expression of the dominant stx genotype could contribute to the reduced incidence of STEC-related illness in Australia. Genomic surveillance of STEC becomes an important part of public health response and ongoing interrogation of virulence factors in STEC offers additional insights for the public health risk assessment.

Revealing COVID-19 transmission in Australia by SARS-CoV-2 genome sequencing and agent-based modeling.

In January 2020, a novel betacoronavirus (family Coronaviridae), named severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), was identified as the etiological agent of a cluster of pneumonia cases occurring in Wuhan City, Hubei Province, China1,2. The disease arising from SARS-CoV-2 infection, coronavirus disease 2019 (COVID-19), subsequently spread rapidly causing a worldwide pandemic. Here we examine the added value of near real-time genome sequencing of SARS-CoV-2 in a subpopulation of infected patients during the first 10 weeks of COVID-19 containment in Australia and compare findings from genomic surveillance with predictions of a computational agent-based model (ABM). Using the Australian census data, the ABM generates over 24 million software agents representing the population of Australia, each with demographic attributes of an anonymous individual. It then simulates transmission of the disease over time, spreading from specific infection sources, using contact rates of individuals within different social contexts. We report that the prospective sequencing of SARS-CoV-2 clarified the probable source of infection in cases where epidemiological links could not be determined, significantly decreased the proportion of COVID-19 cases with contentious links, documented genomically similar cases associated with concurrent transmission in several institutions and identified previously unsuspected links. Only a quarter of sequenced cases appeared to be locally acquired and were concordant with predictions from the ABM. These high-resolution genomic data are crucial to track cases with locally acquired COVID-19 and for timely recognition of independent importations once border restrictions are lifted and trade and travel resume.

Visualization of Phage Genomic Data: Comparative Genomics and Publication-Quality Diagrams.

The presentation of bacteriophage genomes as diagrams allows the location and organization of features to be communicated in a clear and effective manner. A wide range of software applications are available for the clear and accurate visualization of genomic data. Several of these applications incorporate comparative analysis tools, allowing for insertions, deletions, rearrangements and variations in syntenic regions to be visualized. In this chapter, freely available software and resources for the generation of high-quality graphical maps of bacteriophage genomes are listed and discussed.

Genomic Epidemiology of NDM-1-Encoding Plasmids in Latin American Clinical Isolates Reveals Insights into the Evolution of Multidrug Resistance.

Bacteria that produce the broad-spectrum Carbapenem antibiotic New Delhi Metallo-ß-lactamase (NDM) place a burden on health care systems worldwide, due to the limited treatment options for infections caused by them and the rapid global spread of this antibiotic resistance mechanism. Although it is believed that the associated resistance gene blaNDM-1 originated in Acinetobacter spp., the role of Enterobacteriaceae in its dissemination remains unclear. In this study, we used whole genome sequencing to investigate the dissemination dynamics of blaNDM-1-positive plasmids in a set of 21 clinical NDM-1-positive isolates from Colombia and Mexico (Providencia rettgeri, Klebsiella pneumoniae, and Acinetobacter baumannii) as well as six representative NDM-1-positive Escherichia coli transconjugants. Additionally, the plasmids from three representative P. rettgeri isolates were sequenced by PacBio sequencing and finished. Our results demonstrate the presence of previously reported plasmids from K. pneumoniae and A. baumannii in different genetic backgrounds and geographically distant locations in Colombia. Three new previously unclassified plasmids were also identified in P. rettgeri from Colombia and Mexico, plus an interesting genetic link between NDM-1-positive P. rettgeri from distant geographic locations (Canada, Mexico, Colombia, and Israel) without any reported epidemiological links was discovered. Finally, we detected a relationship between plasmids present in P. rettgeri and plasmids from A. baumannii and K. pneumoniae. Overall, our findings suggest a Russian doll model for the dissemination of blaNDM-1 in Latin America, with P. rettgeri playing a central role in this process, and reveal new insights into the evolution and dissemination of plasmids carrying such antibiotic resistance genes.

First Complete Providencia rettgeri Genome Sequence, the NDM-1-Producing Clinical Strain RB151.

Providencia rettgeri is an opportunistic bacterial pathogen of clinical significance due to its association with urinary tract infections and multidrug resistance. Here, we report the first complete genome sequence of P. rettgeri The genome of strain RB151 consists of a 4.8-Mbp chromosome and a 108-kbp blaNDM-1-positive plasmid.

Phylogenetically related Argentinean and Australian Escherichia coli O157 isolates are distinguished by virulence clades and alternative Shiga toxin 1 and 2 prophages.

Shiga toxigenic Escherichia coli O157 is the leading cause of hemolytic uremic syndrome (HUS) worldwide. The frequencies of stx genotypes and the incidences of O157-related illness and HUS vary significantly between Argentina and Australia. Locus-specific polymorphism analysis revealed that lineage I/II (LI/II) E. coli O157 isolates were most prevalent in Argentina (90%) and Australia (88%). Argentinean LI/II isolates were shown to belong to clades 4 (28%) and 8 (72%), while Australian LI/II isolates were identified as clades 6 (15%), 7 (83%), and 8 (2%). Clade 8 was significantly associated with Shiga toxin bacteriophage insertion (SBI) type stx(2) (locus of insertion, argW) in Argentinean isolates (P < 0.0001). In Argentinean LI/II strains, stx(2) is carried by a prophage inserted at argW, whereas in Australian LI/II strains the argW locus is occupied by the novel stx(1) prophage. In both Argentinean and Australian LI/II strains, stx(2c) is almost exclusively carried by a prophage inserted at sbcB. However, alternative q(933)- or q(21)-related alleles were identified in the Australian stx(2c) prophage. Argentinean LI/II isolates were also distinguished from Australian isolates by the presence of the putative virulence determinant ECSP_3286 and the predominance of motile O157:H7 strains. Characteristics common to both Argentinean and Australian LI/II O157 strains included the presence of putative virulence determinants (ECSP_3620, ECSP_0242, ECSP_2687, ECSP_2870, and ECSP_2872) and the predominance of the tir255T allele. These data support further understanding of O157 phylogeny and may foster greater insight into the differential virulence of O157 lineages.