How low can you go? Short-read polishing of Oxford Nanopore bacterial genome assemblies.

It is now possible to assemble near-perfect bacterial genomes using Oxford Nanopore Technologies (ONT) long reads, but short-read polishing is usually required for perfection. However, the effect of short-read depth on polishing performance is not well understood. Here, we introduce Pypolca (with default and careful parameters) and Polypolish v0.6.0 (with a new careful parameter). We then show that: (1) all polishers other than Pypolca-careful, Polypolish-default and Polypolish-careful commonly introduce false-positive errors at low read depth; (2) most of the benefit of short-read polishing occurs by 25× depth; (3) Polypolish-careful almost never introduces false-positive errors at any depth; and (4) Pypolca-careful is the single most effective polisher. Overall, we recommend the following polishing strategies: Polypolish-careful alone when depth is very low (<5×), Polypolish-careful and Pypolca-careful when depth is low (5-25×), and Polypolish-default and Pypolca-careful when depth is sufficient (>25×).

GROND: a quality-checked and publicly available database of full-length 16S-ITS-23S rRNA operon sequences.

Sequence comparison of 16S rRNA PCR amplicons is an established approach to taxonomically identify bacterial isolates and profile complex microbial communities. One potential application of recent advances in long-read sequencing technologies is to sequence entire rRNA operons and capture significantly more phylogenetic information compared to sequencing of the 16S rRNA (or regions thereof) alone, with the potential to increase the proportion of amplicons that can be reliably classified to lower taxonomic ranks. Here we describe GROND (Genome-derived Ribosomal Operon Database), a publicly available database of quality-checked 16S-ITS-23S rRNA operons, accompanied by multiple taxonomic classifications. GROND will aid researchers in analysis of their data and act as a standardised database to allow comparison of results between studies.

Deazaflavin metabolite produced by endosymbiotic bacteria controls fungal host reproduction.

The endosymbiosis between the pathogenic fungus Rhizopus microsporus and the toxin-producing bacterium Mycetohabitans rhizoxinica represents a unique example of host control by an endosymbiont. Fungal sporulation strictly depends on the presence of endosymbionts as well as bacterially produced secondary metabolites. However, an influence of primary metabolites on host control remained unexplored. Recently, we discovered that M. rhizoxinica produces FO and 3PG-F420, a derivative of the specialized redox cofactor F420. Whether FO/3PG-F420 plays a role in the symbiosis has yet to be investigated. Here, we report that FO, the precursor of 3PG-F420, is essential to the establishment of a stable symbiosis. Bioinformatic analysis revealed that the genetic inventory to produce cofactor 3PG-F420 is conserved in the genomes of eight endofungal Mycetohabitans strains. By developing a CRISPR/Cas-assisted base editing strategy for M. rhizoxinica, we generated mutant strains deficient in 3PG-F420 (M. rhizoxinica ΔcofC) and in both FO and 3PG-F420 (M. rhizoxinica ΔfbiC). Co-culture experiments demonstrated that the sporulating phenotype of apo-symbiotic R. microsporus is maintained upon reinfection with wild-type M. rhizoxinica or M. rhizoxinica ΔcofC. In contrast, R. microsporus is unable to sporulate when co-cultivated with M. rhizoxinica ΔfbiC, even though the fungus was observed by super-resolution fluorescence microscopy to be successfully colonized. Genetic and chemical complementation of the FO deficiency of M. rhizoxinica ΔfbiC led to restoration of fungal sporulation, signifying that FO is indispensable for establishing a functional symbiosis. Even though FO is known for its light-harvesting properties, our data illustrate an important role of FO in inter-kingdom communication.

Mycobacterium ulcerans challenge strain selection for a Buruli ulcer controlled human infection model.

Critical scientific questions remain regarding infection with Mycobacterium ulcerans, the organism responsible for the neglected tropical disease, Buruli ulcer (BU). A controlled human infection model has the potential to accelerate our knowledge of the immunological correlates of disease, to test prophylactic interventions and novel therapeutics. Here we present microbiological evidence supporting M. ulcerans JKD8049 as a suitable human challenge strain. This non-genetically modified Australian isolate is susceptible to clinically relevant antibiotics, can be cultured in animal-free and surfactant-free media, can be enumerated for precise dosing, and has stable viability following cryopreservation. Infectious challenge of humans with JKD8049 is anticipated to imitate natural infection, as M. ulcerans JKD8049 is genetically stable following in vitro passage and produces the key virulence factor, mycolactone. Also reported are considerations for the manufacture, storage, and administration of M. ulcerans JKD8049 for controlled human infection.

The 3' UTR of vigR is required for virulence in Staphylococcus aureus and has expanded through STAR sequence repeat insertions.

Infections caused by methicillin-resistant Staphylococcus aureus (MRSA) are alarmingly common, and treatment is confined to last-line antibiotics. Vancomycin is the treatment of choice for MRSA bacteremia, and treatment failure is often associated with vancomycin-intermediate S. aureus isolates. The regulatory 3' UTR of the vigR mRNA contributes to vancomycin tolerance and upregulates the autolysin IsaA. Using MS2-affinity purification coupled with RNA sequencing, we find that the vigR 3' UTR also regulates dapE, a succinyl-diaminopimelate desuccinylase required for lysine and peptidoglycan synthesis, suggesting a broader role in controlling cell wall metabolism and vancomycin tolerance. Deletion of the 3' UTR increased virulence, while the isaA mutant is completely attenuated in a wax moth larvae model. Sequence and structural analyses of vigR indicated that the 3' UTR has expanded through the acquisition of Staphylococcus aureus repeat insertions that contribute sequence for the isaA interaction seed and may functionalize the 3' UTR.

High performance Legionella pneumophila source attribution using genomics-based machine learning classification.

Fundamental to effective Legionnaires' disease outbreak control is the ability to rapidly identify the environmental source(s) of the causative agent, Legionella pneumophila. Genomics has revolutionized pathogen surveillance, but L. pneumophila has a complex ecology and population structure that can limit source inference based on standard core genome phylogenetics. Here, we present a powerful machine learning approach that assigns the geographical source of Legionnaires' disease outbreaks more accurately than current core genome comparisons. Models were developed upon 534 L. pneumophila genome sequences, including 149 genomes linked to 20 previously reported Legionnaires' disease outbreaks through detailed case investigations. Our classification models were developed in a cross-validation framework using only environmental L. pneumophila genomes. Assignments of clinical isolate geographic origins demonstrated high predictive sensitivity and specificity of the models, with no false positives or false negatives for 13 out of 20 outbreak groups, despite the presence of within-outbreak polyclonal population structure. Analysis of the same 534-genome panel with a conventional phylogenomic tree and a core genome multi-locus sequence type allelic distance-based classification approach revealed that our machine learning method had the highest overall classification performance-agreement with epidemiological information. Our multivariate statistical learning approach maximizes the use of genomic variation data and is thus well-suited for supporting Legionnaires' disease outbreak investigations.IMPORTANCEIdentifying the sources of Legionnaires' disease outbreaks is crucial for effective control. Current genomic methods, while useful, often fall short due to the complex ecology and population structure of Legionella pneumophila, the causative agent. Our study introduces a high-performing machine learning approach for more accurate geographical source attribution of Legionnaires' disease outbreaks. Developed using cross-validation on environmental L. pneumophila genomes, our models demonstrate excellent predictive sensitivity and specificity. Importantly, this new approach outperforms traditional methods like phylogenomic trees and core genome multi-locus sequence typing, proving more efficient at leveraging genomic variation data to infer outbreak sources. Our machine learning algorithms, harnessing both core and accessory genomic variation, offer significant promise in public health settings. By enabling rapid and precise source identification in Legionnaires' disease outbreaks, such approaches have the potential to expedite intervention efforts and curtail disease transmission.

Mosquitoes provide a transmission route between possums and humans for Buruli ulcer in southeastern Australia.

Buruli ulcer, a chronic subcutaneous infection caused by Mycobacterium ulcerans, is increasing in prevalence in southeastern Australia. Possums are a local wildlife reservoir for M. ulcerans and, although mosquitoes have been implicated in transmission, it remains unclear how humans acquire infection. We conducted extensive field survey analyses of M. ulcerans prevalence among mosquitoes in the Mornington Peninsula region of southeastern Australia. PCR screening of trapped mosquitoes revealed a significant association between M. ulcerans and Aedes notoscriptus. Spatial scanning statistics revealed overlap between clusters of M. ulcerans-positive Ae. notoscriptus, M. ulcerans-positive possum excreta and Buruli ulcer cases, and metabarcoding analyses showed individual mosquitoes had fed on humans and possums. Bacterial genomic analysis confirmed shared single-nucleotide-polymorphism profiles for M. ulcerans detected in mosquitoes, possum excreta and humans. These findings indicate Ae. notoscriptus probably transmit M. ulcerans in southeastern Australia and highlight mosquito control as a Buruli ulcer prevention measure.

High-performance enrichment-based genome sequencing to support the investigation of hepatitis A virus outbreaks.

IMPORTANCE: This proof-of-concept study introduces a hybrid capture oligo panel for whole-genome sequencing of all six human pathogenic hepatitis A virus (HAV) subgenotypes, exhibiting a higher sensitivity than some conventional genotyping assays. The ability of hybrid capture to enrich multiple targets allows for a single, streamlined workflow, thus facilitating the potential harmonization of molecular surveillance of HAV with other enteric viruses. Even challenging sample matrices can be accommodated, making them suitable for broad implementation in clinical and public health laboratories. This innovative approach has significant implications for enhancing multijurisdictional outbreak investigations as well as our understanding of the global diversity and transmission dynamics of HAV.

Transcription activator-like effectors from endosymbiotic bacteria control the reproduction of their fungal host.

IMPORTANCE: Interactions between fungi and bacteria are critically important in ecology, medicine, and biotechnology. In this study, we shed light on factors that promote the persistence of a toxin-producing, phytopathogenic Rhizopus-Mycetohabitans symbiosis that causes severe crop losses in Asia. We present an unprecedented case where bacterially produced transcription activator-like (TAL) effectors are key to maintaining a stable endosymbiosis. In their absence, fungal sporulation is abrogated, leading to collapse of the phytopathogenic alliance. The Mycetohabitans TAL (MTAL)-mediated mechanism of host control illustrates a unique role of bacterial effector molecules that has broader implications, potentially serving as a model to understand how prokaryotic symbionts interact with their eukaryotic hosts.

Bioinformatic investigation of discordant sequence data for SARS-CoV-2: insights for robust genomic analysis during pandemic surveillance.

The COVID-19 pandemic has necessitated the rapid development and implementation of whole-genome sequencing (WGS) and bioinformatic methods for managing the pandemic. However, variability in methods and capabilities between laboratories has posed challenges in ensuring data accuracy. A national working group comprising 18 laboratory scientists and bioinformaticians from Australia and New Zealand was formed to improve data concordance across public health laboratories (PHLs). One effort, presented in this study, sought to understand the impact of the methodology on consensus genome concordance and interpretation. SARS-CoV-2 WGS proficiency testing programme (PTP) data were retrospectively obtained from the 2021 Royal College of Pathologists of Australasia Quality Assurance Programmes (RCPAQAP), which included 11 participating Australian laboratories. The submitted consensus genomes and reads from eight contrived specimens were investigated, focusing on discordant sequence data and findings were presented to the working group to inform best practices. Despite using a variety of laboratory and bioinformatic methods for SARS-CoV-2 WGS, participants largely produced concordant genomes. Two participants returned five discordant sites in a high-Cτ replicate, which could be resolved with reasonable bioinformatic quality thresholds. We noted ten discrepancies in genome assessment that arose from nucleotide heterogeneity at three different sites in three cell-culture-derived control specimens. While these sites were ultimately accurate after considering the participants' bioinformatic parameters, it presented an interesting challenge for developing standards to account for intrahost single nucleotide variation (iSNV). Observed differences had little to no impact on key surveillance metrics, lineage assignment and phylogenetic clustering, while genome coverage <90 % affected both. We recommend PHLs bioinformatically generate two consensus genomes with and without ambiguity thresholds for quality control and downstream analysis, respectively, and adhere to a minimum 90 % genome coverage threshold for inclusion in surveillance interpretations. We also suggest additional PTP assessment criteria, including primer efficiency, detection of iSNVs and minimum genome coverage of 90 %. This study underscores the importance of multidisciplinary national working groups in informing guidelines in real time for bioinformatic quality acceptance criteria. It demonstrates the potential for enhancing public health responses through improved data concordance and quality control in SARS-CoV-2 genomic analysis during pandemic surveillance.

Clinical Staphylococcus aureus inhibits human T-cell activity through interaction with the PD-1 receptor.

IMPORTANCE: Therapies that target and aid the host immune defense to repel cancer cells or invading pathogens are rapidly emerging. Antibiotic resistance is among the largest threats to human health globally. Staphylococcus aureus (S. aureus) is the most common bacterial infection, and it poses a challenge to the healthcare system due to its significant ability to develop resistance toward current available therapies. In long-term infections, S. aureus further adapt to avoid clearance by the host immune defense. In this study, we discover a new interaction that allows S. aureus to avoid elimination by the immune system, which likely supports its persistence in the host. Moreover, we find that blocking the specific receptor (PD-1) using antibodies significantly relieves the S. aureus-imposed inhibition. Our findings suggest that therapeutically targeting PD-1 is a possible future strategy for treating certain antibiotic-resistant staphylococcal infections.

Comparison of contemporary invasive and non-invasive Streptococcus pneumoniae isolates reveals new insights into circulating anti-microbial resistance determinants.

Streptococcus pneumoniae is a major human pathogen with a high burden of disease. Non-invasive isolates (those found in non-sterile sites) are thought to be a key source of invasive isolates (those found in sterile sites) and a reservoir of anti-microbial resistance (AMR) determinants. Despite this, pneumococcal surveillance has almost exclusively focused on invasive isolates. We aimed to compare contemporaneous invasive and non-invasive isolate populations to understand how they interact and identify differences in AMR gene distribution. We used a combination of whole-genome sequencing and phenotypic anti-microbial susceptibility testing and a data set of invasive (n = 1,288) and non-invasive (n = 186) pneumococcal isolates, collected in Victoria, Australia, between 2018 and 2022. The non-invasive population had increased levels of antibiotic resistance to multiple classes of antibiotics including beta-lactam antibiotics penicillin and ceftriaxone. We identified genomic intersections between the invasive and non-invasive populations and no distinct phylogenetic clustering of the two populations. However, this analysis revealed sub-populations overrepresented in each population. The sub-populations that had high levels of AMR were overrepresented in the non-invasive population. We determined that WamR-Pneumo was the most accurate in silico tool for predicting resistance to the antibiotics tested. This tool was then used to assess the allelic diversity of the penicillin-binding protein genes, which acquire mutations leading to beta-lactam antibiotic resistance, and found that they were highly conserved (≥80% shared) between the two populations. These findings show the potential of non-invasive isolates to serve as reservoirs of AMR determinants.

Mycobacterium ulcerans not detected by PCR on human skin in Buruli ulcer endemic areas of south eastern Australia.

INTRODUCTION: Mycobacterium ulcerans (MU) causes Buruli ulcer (Buruli), a geographically restricted infection that can result in skin loss, contracture and permanent scarring. Lesion-location maps compiled from more than 640 cases in south eastern Australia suggest biting insects are likely involved in transmission, but it is unclear whether MU is brought by insects to humans or if MU is already on the skin and inoculation is an opportunistic event that need not be insect dependent. METHODS: We validated a PCR swab detection assay and defined its dynamic range using laboratory cultured M. ulcerans and fresh pigskin. We invited volunteers in Buruli-endemic and non-endemic areas to sample their skin surfaces with self-collected skin swabs tested by IS2404 quantitative PCR. RESULTS: Pigskin validation experiments established a limit-of-detection of 0.06 CFU/cm2 at a qPCR cycle threshold (Ct) of 35. Fifty-seven volunteers returned their self-collected kits of 4 swabs (bilateral ankles, calves, wrists, forearms), 10 from control areas and 47 from endemic areas. Collection was timed to coincide with the known peak-transmission period of Buruli. All swabs from human volunteers tested negative (Ct ≥35). CONCLUSIONS: M. ulcerans was not detected on the skin of humans from highly Buruli endemic areas.

The two-component system WalKR provides an essential link between cell wall homeostasis and DNA replication in Staphylococcus aureus.

Among the 16 two-component systems in the opportunistic human pathogen Staphylococcus aureus, only WalKR is essential. Like the orthologous systems in other Bacillota, S. aureus WalKR controls autolysins involved in peptidoglycan remodeling and is therefore intimately involved in cell division. However, despite the importance of WalKR in S. aureus, the basis for its essentiality is not understood and the regulon is poorly defined. Here, we defined a consensus WalR DNA-binding motif and the direct WalKR regulon by using functional genomics, including chromatin immunoprecipitation sequencing, with a panel of isogenic walKR mutants that had a spectrum of altered activities. Consistent with prior findings, the direct regulon includes multiple autolysin genes. However, this work also revealed that WalR directly regulates at least five essential genes involved in lipoteichoic acid synthesis (ltaS): translation (rplK), DNA compaction (hup), initiation of DNA replication (dnaA, hup) and purine nucleotide metabolism (prs). Thus, WalKR in S. aureus serves as a polyfunctional regulator that contributes to fundamental control over critical cell processes by coordinately linking cell wall homeostasis with purine biosynthesis, protein biosynthesis, and DNA replication. Our findings further address the essentiality of this locus and highlight the importance of WalKR as a bona fide target for novel anti-staphylococcal therapeutics. IMPORTANCE The opportunistic human pathogen Staphylococcus aureus uses an array of protein sensing systems called two-component systems (TCS) to sense environmental signals and adapt its physiology in response by regulating different genes. This sensory network is key to S. aureus versatility and success as a pathogen. Here, we reveal for the first time the full extent of the regulatory network of WalKR, the only staphylococcal TCS that is indispensable for survival under laboratory conditions. We found that WalKR is a master regulator of cell growth, coordinating the expression of genes from multiple, fundamental S. aureus cellular processes, including those involved in maintaining cell wall metabolism, protein biosynthesis, nucleotide metabolism, and the initiation of DNA replication.

A statistical genomics framework to trace bacterial genomic predictors of clinical outcomes in Staphylococcus aureus bacteremia.

Outcomes of severe bacterial infections are determined by the interplay between host, pathogen, and treatments. While human genomics has provided insights into host factors impacting Staphylococcus aureus infections, comparatively little is known about S. aureus genotypes and disease severity. Building on the hypothesis that bacterial pathoadaptation is a key outcome driver, we developed a genome-wide association study (GWAS) framework to identify adaptive mutations associated with treatment failure and mortality in S. aureus bacteremia (1,358 episodes). Our research highlights the potential of vancomycin-selected mutations and vancomycin minimum inhibitory concentration (MIC) as key explanatory variables to predict infection severity. The contribution of bacterial variation was much lower for clinical outcomes (heritability <5%); however, GWASs allowed us to identify additional, MIC-independent candidate pathogenesis loci. Using supervised machine learning, we were able to quantify the predictive potential of these adaptive signatures. Our statistical genomics framework provides a powerful means to capture adaptive mutations impacting severe bacterial infections.

Analysis of Rhizonin Biosynthesis Reveals Origin of Pharmacophoric Furylalanine Moieties in Diverse Cyclopeptides.

Rhizonin A and B are hepatotoxic cyclopeptides produced by bacterial endosymbionts (Mycetohabitans endofungorum) of the fungus Rhizopus microsporus. Their toxicity critically depends on the presence of 3-furylalanine (Fua) residues, which also occur in pharmaceutically relevant cyclopeptides of the endolide and bingchamide families. The biosynthesis and incorporation of Fua by non-ribosomal peptide synthetases (NRPS), however, has remained elusive. By genome sequencing and gene inactivation we elucidated the gene cluster responsible for rhizonin biosynthesis. A suite of isotope labeling experiments identified tyrosine and l-DOPA as Fua precursors and provided the first mechanistic insight. Bioinformatics, mutational analysis and heterologous reconstitution identified dioxygenase RhzB as necessary and sufficient for Fua formation. RhzB is a novel type of heme-dependent aromatic oxygenases (HDAO) that enabled the discovery of the bingchamide biosynthesis gene cluster through genome mining.

Population structure, serotype distribution and antibiotic resistance of Streptococcus pneumoniae causing invasive disease in Victoria, Australia.

Streptococcus pneumoniae is a major human pathogen and can cause a range of conditions from asymptomatic colonization to invasive pneumococcal disease (IPD). The epidemiology and distribution of IPD-causing serotypes in Australia has undergone large changes following the introduction of the 7-valent pneumococcal conjugate vaccine (PCV) in 2005 and the 13-valent PCV in 2011. In this study, to provide a contemporary understanding of the IPD causing population in Victoria, Australia, we aimed to examine the population structure and prevalence of antimicrobial resistance using whole-genome sequencing and comprehensive antimicrobial susceptibility data of 1288 isolates collected between 2018 and 2022. We observed high diversity among the isolates with 52 serotypes, 203 sequence types (STs) and 70 Global Pneumococcal Sequencing Project Clusters (GPSCs) identified. Serotypes contained in the 13v-PCV represented 35.3 % (n=405) of isolates. Antimicrobial resistance (AMR) to at least one antibiotic was identified in 23.8 % (n=358) of isolates with penicillin resistance the most prevalent (20.3 %, n=261 using meningitis breakpoints and 5.1 % n=65 using oral breakpoints). Of the AMR isolates, 28 % (n=101) were multidrug resistant (MDR) (resistant to three or more drug classes). Vaccination status of cases was determined for a subset of isolates with 34 cases classified as vaccine failure events (fully vaccinated IPD cases of vaccine serotype). However, no phylogenetic association with failure events was observed. Within the highly diverse IPD population, we identified six high-risk sub-populations of public health concern characterized by high prevalence, high rates of AMR and MDR, or serotype inclusion in vaccines. High-risk serotypes included serotypes 3, 19F, 19A, 14, 11A, 15A and serofamily 23. In addition, we present our data validating seroBA for in silico serotyping to facilitate ISO-accreditation of this test in routine use in a public health reference laboratory and have made this data set available. This study provides insights into the population dynamics, highlights non-vaccine serotypes of concern that are highly resistant, and provides a genomic framework for the ongoing surveillance of IPD in Australia which can inform next-generation IPD prevention strategies.

Discovery and Biosynthesis of the Cytotoxic Polyene Terpenomycin in Human Pathogenic Nocardia.

Nocardia are opportunistic human pathogens that can cause a range of debilitating and difficult to treat infections of the lungs, brain, skin, and soft tissues. Despite their close relationship to the well-known secondary metabolite-producing genus, Streptomyces, comparatively few natural products are known from the Nocardia, and even less is known about their involvement in the pathogenesis. Here, we combine chemistry, genomics, and molecular microbiology to reveal the production of terpenomycin, a new cytotoxic and antifungal polyene from a human pathogenic Nocardia terpenica isolate. We unveil the polyketide synthase (PKS) responsible for terpenomycin biosynthesis and show that it combines several unusual features, including "split", skipped, and iteratively used modules, and the use of the unusual extender unit methoxymalonate as a starter unit. To link genes to molecules, we constructed a transposon mutant library in N. terpenica, identifying a terpenomycin-null mutant with an inactivated terpenomycin PKS. Our findings show that the neglected actinomycetes have an unappreciated capacity for the production of bioactive molecules with unique biosynthetic pathways waiting to be uncovered and highlights these organisms as producers of diverse natural products.

A human model of Buruli ulcer: The case for controlled human infection and considerations for selecting a Mycobacterium ulcerans challenge strain.

Critical knowledge gaps regarding infection with Mycobacterium ulcerans, the cause of Buruli ulcer (BU), have impeded development of new therapeutic approaches and vaccines for prevention of this neglected tropical disease. Here, we review the current understanding of host-pathogen interactions and correlates of immune protection to explore the case for establishing a controlled human infection model of M. ulcerans infection. We also summarise the overarching safety considerations and present a rationale for selecting a suitable challenge strain.