- Volume 10, Issue 5, 2024

Wastewater-based surveillance (WBS) is an important epidemiological and public health tool for tracking pathogens across the scale of a building, neighbourhood, city, or region. WBS gained widespread adoption globally during the SARS-CoV-2 pandemic for estimating community infection levels by qPCR. Sequencing pathogen genes or genomes from wastewater adds information about pathogen genetic diversity, which can be used to identify viral lineages (including variants of concern) that are circulating in a local population. Capturing the genetic diversity by WBS sequencing is not trivial, as wastewater samples often contain a diverse mixture of viral lineages with real mutations and sequencing errors, which must be deconvoluted computationally from short sequencing reads. In this study we assess nine different computational tools that have recently been developed to address this challenge. We simulated 100 wastewater sequence samples consisting of SARS-CoV-2 BA.1, BA.2, and Delta lineages, in various mixtures, as well as a Delta–Omicron recombinant and a synthetic ‘novel’ lineage. Most tools performed well in identifying the true lineages present and estimating their relative abundances and were generally robust to variation in sequencing depth and read length. While many tools identified lineages present down to 1 % frequency, results were more reliable above a 5 % threshold. The presence of an unknown synthetic lineage, which represents an unclassified SARS-CoV-2 lineage, increases the error in relative abundance estimates of other lineages, but the magnitude of this effect was small for most tools. The tools also varied in how they labelled novel synthetic lineages and recombinants. While our simulated dataset represents just one of many possible use cases for these methods, we hope it helps users understand potential sources of error or bias in wastewater sequencing analysis and to appreciate the commonalities and differences across methods.

Staphylococcus aureus asymptomatically colonises 30 % of humans but can also cause a range of diseases, which can be fatal. In 2017 S. aureus was associated with 20 000 deaths in the USA alone. Dividing S. aureus isolates into smaller sub-groups can reveal the emergence of distinct sub-populations with varying potential to cause infections. Despite multiple molecular typing methods categorising such sub-groups, they do not take full advantage of S. aureus genome sequences when describing the fundamental population structure of the species. In this study, we developed Staphylococcus aureus Lineage Typing (SaLTy), which rapidly divides the species into 61 phylogenetically congruent lineages. Alleles of three core genes were identified that uniquely define the 61 lineages and were used for SaLTy typing. SaLTy was validated on 5000 genomes and 99.12 % (4956/5000) of isolates were assigned the correct lineage. We compared SaLTy lineages to previously calculated clonal complexes (CCs) from BIGSdb (n=21 173). SALTy improves on CCs by grouping isolates congruently with phylogenetic structure. SaLTy lineages were further used to describe the carriage of Staphylococcal chromosomal cassette containing mecA (SCCmec) which is carried by methicillin-resistant S. aureus (MRSA). Most lineages had isolates lacking SCCmec and the four largest lineages varied in SCCmec over time. Classifying isolates into SaLTy lineages, which were further SCCmec typed, allowed SaLTy to describe high-level MRSA epidemiology. We provide SaLTy as a simple typing method that defines phylogenetic lineages (https://github.com/LanLab/SaLTy). SaLTy is highly accurate and can quickly analyse large amounts of S. aureus genome data. SaLTy will aid the characterisation of S. aureus populations and ongoing surveillance of sub-groups that threaten human health.

Lyme disease (LD), caused by spirochete bacteria of the genus Borrelia burgdorferi sensu lato, remains the most common vector-borne disease in the northern hemisphere. Borrelia outer surface protein A (OspA) is an integral surface protein expressed during the tick cycle, and a validated vaccine target. There are at least 20 recognized Borrelia genospecies, that vary in OspA serotype. This study presents a new in silico sequence-based method for OspA typing using next-generation sequence data. Using a compiled database of over 400 Borrelia genomes encompassing the 4 most common disease-causing genospecies, we characterized OspA diversity in a manner that can accommodate existing and new OspA types and then defined boundaries for classification and assignment of OspA types based on the sequence similarity. To accommodate potential novel OspA types, we have developed a new nomenclature: OspA in silico type (IST). Beyond the ISTs that corresponded to existing OspA serotypes 1–8, we identified nine additional ISTs that cover new OspA variants in B. bavariensis (IST9–10), B. garinii (IST11–12), and other Borrelia genospecies (IST13–17). The IST typing scheme and associated OspA variants are available as part of the PubMLST Borrelia spp. database. Compared to traditional OspA serotyping methods, this new computational pipeline provides a more comprehensive and broadly applicable approach for characterization of OspA type and Borrelia genospecies to support vaccine development.

The Genome Taxonomy Database (GTDB) provides a species to domain classification of publicly available genomes based on average nucleotide identity (ANI) (for species) and a concatenated gene phylogeny normalized by evolutionary rates (for genus to phylum), which has been widely adopted by the scientific community. Here, we use the Genome UNClutterer (GUNC) software to identify putatively contaminated genomes in GTDB release 07-RS207. We found that GUNC reported 35,723 genomes as putatively contaminated, comprising 11.25 % of the 317,542 genomes in GTDB release 07-RS207. To assess the impact of this high level of inferred contamination on the delineation of taxa, we created ‘clean’ versions of the 34,846 putatively contaminated bacterial genomes by removing the most contaminated half. For each clean half, we re-calculated the ANI and concatenated gene phylogeny and found that only 77 (0.22 %) of the genomes were not consistent with their original classification. We conclude that the delineation of taxa in GTDB is robust to the putative contamination detected by GUNC.

The genus Catenibacillus (family Lachnospiraceae, phylum Bacillota) includes only one cultivated species so far, Catenibacillus scindens, isolated from human faeces and capable of deglycosylating dietary polyphenols and degrading flavonoid aglycones. Another human intestinal Catenibacillus strain not taxonomically resolved at that time was recently genome-sequenced. We analysed the genome of this novel isolate, designated Catenibacillus decagia, and showed its ability to deglycosylate C-coupled flavone and xanthone glucosides and O-coupled flavonoid glycosides. Most of the resulting aglycones were further degraded to the corresponding phenolic acids. Including the recently sequenced genome of C. scindens and ten faecal metagenome-assembled genomes assigned to the genus Catenibacillus, we performed a comparative genome analysis and searched for genes encoding potential C-glycosidases and other polyphenol-converting enzymes. According to genome data and physiological characterization, the core metabolism of Catenibacillus strains is based on a fermentative lifestyle with butyrate production and hydrogen evolution. Both C. scindens and C. decagia encode a flavonoid O-glycosidase, a flavone reductase, a flavanone/flavanonol-cleaving reductase and a phloretin hydrolase. Several gene clusters encode enzymes similar to those of the flavonoid C-deglycosylation system of Dorea strain PUE (DgpBC), while separately located genes encode putative polyphenol-glucoside oxidases (DgpA) required for C-deglycosylation. The diversity of dgpA and dgpBC gene clusters might explain the broad C-glycoside substrate spectrum of C. scindens and C. decagia. The other Catenibacillus genomes encode only a few potential flavonoid-converting enzymes. Our results indicate that several Catenibacillus species are well-equipped to deglycosylate and degrade dietary plant polyphenols and might inhabit a corresponding, specific niche in the gut.

Natural products from Actinomycetota have served as inspiration for many clinically relevant therapeutics. Despite early triumphs in natural product discovery, the rate of unearthing new compounds has decreased, necessitating inventive approaches. One promising strategy is to explore environments where survival is challenging. These harsh environments are hypothesized to lead to bacteria developing chemical adaptations (e.g. natural products) to enable their survival. This investigation focuses on ore-forming environments, particularly fluoride mines, which typically have extreme pH, salinity and nutrient scarcity. Herein, we have utilized metagenomics, metabolomics and evolutionary genome mining to dissect the biodiversity and metabolism in these harsh environments. This work has unveiled the promising biosynthetic potential of these bacteria and has demonstrated their ability to produce bioactive secondary metabolites. This research constitutes a pioneering endeavour in bioprospection within fluoride mining regions, providing insights into uncharted microbial ecosystems and their previously unexplored natural products.

Paramyxoviruses are a group of single-stranded, negative-sense RNA viruses, some of which are responsible for acute human disease, including parainfluenza virus, measles virus, Nipah virus and Hendra virus. In recent years, a large number of novel paramyxoviruses, particularly members of the genus Jeilongvirus, have been discovered in wild mammals, suggesting that the diversity of paramyxoviruses may be underestimated. Here we used hemi-nested reverse transcription PCR to obtain 190 paramyxovirus sequences from 969 small mammals in Hubei Province, Central China. These newly identified paramyxoviruses were classified into four clades: genera Jeilongvirus, Morbillivirus, Henipavirus and Narmovirus, with most of them belonging to the genus Jeilongvirus. Using Illumina sequencing and Sanger sequencing, we successfully recovered six near-full-length genomes with different genomic organizations, revealing the more complex genome content of paramyxoviruses. Co-divergence analysis of jeilongviruses and their known hosts indicates that host-switching occurred more frequently in the evolutionary histories of the genus Jeilongvirus. Together, our findings demonstrate the high prevalence of paramyxoviruses in small mammals, especially jeilongviruses, and highlight the diversity of paramyxoviruses and their genome content, as well as the evolution of jeilongviruses.

Clostridioides difficile infection (CDI) remains a significant public health threat globally. New interventions to treat CDI rely on an understanding of the evolution and epidemiology of circulating strains. Here we provide longitudinal genomic data on strain diversity, transmission dynamics and antimicrobial resistance (AMR) of C. difficile ribotypes (RTs) 014/020 (n=169), 002 (n=77) and 056 (n=36), the three most prominent C. difficile strains causing CDI in Australia. Genome scrutiny showed that AMR was uncommon in these lineages, with resistance-conferring alleles present in only 15/169 RT014/020 strains (8.9 %), 1/36 RT056 strains (2.78 %) and none of 77 RT002 strains. Notably, ~90 % of strains were resistant to MLSB agents in vitro, but only ~5.9 % harboured known resistance alleles, highlighting an incongruence between AMR genotype and phenotype. Core genome analyses revealed all three RTs contained genetically heterogeneous strain populations with limited evidence of clonal transmission between CDI cases. The average number of pairwise core genome SNP (cgSNP) differences within each RT group ranged from 23.3 (RT056, ST34, n=36) to 115.6 (RT002, ST8, n=77) and 315.9 (RT014/020, STs 2, 13, 14, 49, n=169). Just 19 clonal groups (encompassing 40 isolates), defined as isolates differing by ≤2 cgSNPs, were identified across all three RTs (RT014/020, n=14; RT002, n=3; RT056, n=2). Of these clonal groups, 63 % (12/19) comprised isolates from the same Australian State and 37 % (7/19) comprised isolates from different States. The low number of plausible transmission events found for these major RTs (and previously documented populations in animal and environmental sources/reservoirs) points to widespread and persistent community sources of diverse C. difficile strains as opposed to ongoing nationwide healthcare outbreaks dominated by a single clone. Together, these data provide new insights into the evolution of major lineages causing CDI in Australia and highlight the urgent need for enhanced surveillance, and for public health interventions to move beyond the healthcare setting and into a One Health paradigm to effectively combat this complex pathogen.

Invasive fungal pathogens pose a substantial threat to widely cultivated crop species, owing to their capacity to adapt to new hosts and new environmental conditions. Gaining insights into the demographic history of these pathogens and unravelling the mechanisms driving coevolutionary processes are crucial for developing durably effective disease management programmes. Pyrenophora teres is a significant fungal pathogen of barley, consisting of two lineages, Ptt and Ptm, with global distributions and demographic histories reflecting barley domestication and spread. However, the factors influencing the population structure of P. teres remain poorly understood, despite the varietal and environmental heterogeneity of barley agrosystems. Here, we report on the population genomic structure of P. teres in France and globally. We used genotyping-by-sequencing to show that Ptt and Ptm can coexist in the same area in France, with Ptt predominating. Furthermore, we showed that differences in the vernalization requirement of barley varieties were associated with population differentiation within Ptt in France and at a global scale, with one population cluster found on spring barley and another population cluster found on winter barley. Our results demonstrate how cultivation conditions, possibly associated with genetic differences between host populations, can be associated with the maintenance of divergent invasive pathogen populations coexisting over large geographic areas. This study not only advances our understanding of the coevolutionary dynamics of the Pt-barley pathosystem but also prompts further research on the relative contributions of adaptation to the host versus adaptation to abiotic conditions in shaping Ptt populations.

Evidence is accumulating in the literature that the horizontal spread of antimicrobial resistance (AMR) genes mediated by bacteriophages and bacteriophage-like plasmid (phage-plasmid) elements is much more common than previously envisioned. For instance, we recently identified and characterized a circular P1-like phage-plasmid harbouring a bla CTX-M-15 gene conferring extended-spectrum beta-lactamase (ESBL) resistance in Salmonella enterica serovar Typhi. As the prevalence and epidemiological relevance of such mechanisms has never been systematically assessed in Enterobacterales, in this study we carried out a follow-up retrospective analysis of UK Salmonella isolates previously sequenced as part of routine surveillance protocols between 2016 and 2021. Using a high-throughput bioinformatics pipeline we screened 47 784 isolates for the presence of the P1 lytic replication gene repL, identifying 226 positive isolates from 25 serovars and demonstrating that phage-plasmid elements are more frequent than previously thought. The affinity for phage-plasmids appears highly serovar-dependent, with several serovars being more likely hosts than others; most of the positive isolates (170/226) belonged to S. Typhimurium ST34 and ST19. The phage-plasmids ranged between 85.8 and 98.2 kb in size, with an average length of 92.1 kb; detailed analysis indicated a high amount of diversity in gene content and genomic architecture. In total, 132 phage-plasmids had the p0111 plasmid replication type, and 94 the IncY type; phylogenetic analysis indicated that both horizontal and vertical gene transmission mechanisms are likely to be involved in phage-plasmid propagation. Finally, phage-plasmids were present in isolates that were resistant and non-resistant to antimicrobials. In addition to providing a first comprehensive view of the presence of phage-plasmids in Salmonella, our work highlights the need for a better surveillance and understanding of phage-plasmids as AMR carriers, especially through their characterization with long-read sequencing.

Aeromonas spp. are commonly found in the aquatic environment and have been responsible for motile Aeromonas septicemia (MAS) in striped catfish, resulting in significant economic loss. These organisms also cause a range of opportunistic infections in humans with compromised immune systems. Here, we conducted a genomic investigation of 87 Aeromonas isolates derived from diseased catfish, healthy catfish and environmental water in catfish farms affected by MAS outbreaks in eight provinces in Mekong Delta (years: 2012–2022), together with 25 isolates from humans with bloodstream infections (years: 2010–2020). Genomics-based typing method precisely delineated Aeromonas species while traditional methods such as aerA PCR and MALDI-TOF were unable identify A. dhakensis. A. dhakensis was found to be more prevalent than A. hydrophila in both diseased catfish and human infections. A. dhakensis sequence type (ST) 656 followed by A. hydrophila ST251 were the predominant virulent species-lineages in diseased catfish (43.7 and 20.7 %, respectively), while diverse STs were found in humans with bloodstream infections. There was evidence of widespread transmission of ST656 and ST251 on striped catfish in the Mekong Delta region. ST656 and ST251 isolates carried a significantly higher number of acquired antimicrobial resistance (AMR) genes and virulence factors in comparison to other STs. They, however, exhibited several distinctions in key virulence factors (i.e. lack of type IV pili and enterotoxin ast in A. dhakensis), AMR genes (i.e. presence of imiH carbapenemase in A. dhakensis), and accessory gene content. To uncover potential conserved proteins of Aeromonas spp. for vaccine development, pangenome analysis has unveiled 2202 core genes between ST656 and ST251, of which 78 proteins were in either outer membrane or extracellular proteins. Our study represents one of the first genomic investigations of the species distribution, genetic landscape, and epidemiology of Aeromonas in diseased catfish and human infections in Vietnam. The emergence of antimicrobial resistant and virulent A. dhakensis strains underscores the needs of enhanced genomic surveillance and strengthening vaccine research and development in preventing Aeromonas diseases in catfish and humans, and the search for potential vaccine candidates could focus on Aeromonas core genes encoded for membrane and secreted proteins.

Aeromonas caviae is an emerging human enteric pathogen. However, the genomic features and virulence genes of A. caviae strains from human gastroenteritis and other sources have not been fully elucidated. Here, we conducted a genomic analysis of 565 global A. caviae strains isolated from different sources, including 261 strains isolated from faecal samples of gastroenteritis patients, of which 18 genomes were sequenced in this study. The presence of bacterial virulence genes and secretion systems in A. caviae strains from different sources was compared, and the phylogenetic relationship of A. caviae strains was assessed based on the core genome. The complete genome of A. caviae strain A20-9 isolated from a gastroenteritis patient was obtained in this study, from which 300 putative virulence factors and a T4SS-encoding plasmid, pAC, were identified. Genes encoding T4SS were also identified in a novel genomic island, ACI-1, from other T4SS-positive strains. The prevalence of T4SS was significantly lower in A. caviae strains from gastroenteritis patients than in environmental strains (3 %, P<0.0001 vs 14 %, P<0.01). Conversely, the prevalence of T6SS was significantly higher in A. caviae strains isolated from gastroenteritis patients than in environmental strains (25 %, P<0.05 vs 13  %, P<0.01). Four phylogenetic clusters were formed based on the core genome of 565 A. caviae strains, and strains carrying T6SS often showed close phylogenetic relationships. T3SS, aerolysin and thermostable cytotonic enterotoxin were absent in all 565 A. caviae strains. Our findings provide novel information on the genomic features of A. caviae and suggest that T6SS may play a role in A. caviae-induced human gastroenteritis.

Cutaneous ulcers are common in yaws-endemic areas. Although often attributed to 'Treponema pallidum subsp. pertenue' and Haemophilus ducreyi, quantitative PCR has highlighted a significant proportion of these ulcers are negative for both pathogens and are considered idiopathic. This is a retrospective analysis utilising existing 16S rRNA sequencing data from two independent yaws studies that took place in Ghana and the Solomon Islands. We characterized bacterial diversity in 38 samples to identify potential causative agents for idiopathic cutaneous ulcers. We identified a diverse bacterial profile, including Arcanobacterium haemolyticum, Campylobacter concisus, Corynebacterium diphtheriae, Staphylococcus spp. and Streptococcus pyogenes, consistent with findings from previous cutaneous ulcer microbiome studies. No single bacterial species was universally present across all samples. The most prevalent bacterium, Campylobacter ureolyticus, appeared in 42% of samples, suggesting a multifactorial aetiology for cutaneous ulcers in yaws-endemic areas. This study emphasizes the need for a nuanced understanding of potential causative agents. The findings prompt further exploration into the intricate microbial interactions contributing to idiopathic yaw-like ulcers, guiding future research toward comprehensive diagnostic and therapeutic strategies.