- Volume 10, Issue 11, 2024

El Niño events, the warm phase of the El Niño Southern Oscillation, facilitate the movement of warm surface waters eastwards across the Pacific Ocean. Marine organisms transported by these waters can act as biological corridors for water-borne bacteria with attachment abilities. El Niño events have been hypothesized as driving the recent emergence of Vibrio parahaemolyticus (Vp) variants, marine bacterium causing gastroenteritis, in South America, but the lack of a robust methodological framework limited any further exploration. Here, we introduce two new analysis approaches to explore Vp dynamics in South America, which will be central to uncovering Vp dynamics in the future. Distributed non-linear lag models found that strong El Niño events increase the relative probability of Vp detection in Peru, with a 3–4-month lag time. Machine learning found that the presence of a specific gene (vopZ) involved in attachment to plankton in a pandemic Vp clone in South America was temporally associated with strong El Niño events, offering a possible strategy for survival over long-range dispersal, such as that offered by El Niño events. Robust surveillance of marine pathogens and methodological development are necessary to produce resolute conclusions on the effect of El Niño events on water-borne diseases.

Genome-based analysis allows for large-scale classification of diverse bacteria and has been widely adopted for delineating species. Unfortunately, for higher taxonomic ranks such as genus, establishing a generally accepted approach based on genome analysis is challenging. While core-genome phylogenies depict the evolutionary relationships among species, determining the correspondence between clades and genera may not be straightforward. For genotypic divergence, the percentage of conserved proteins and genome-wide average amino acid identity are commonly used, but often do not provide a clear threshold for classification. In this work, we investigated the utility of global comparisons and data visualization in identifying clusters of species based on their overall gene content and rationalized that such patterns can be integrated with phylogeny and other information such as phenotypes for improving taxonomy. As a proof of concept, we selected 177 representative genome sequences from the Mycoplasmatales–Entomoplasmatales clade within the class Mollicutes for a case study. We found that the clustering patterns corresponded to the current understanding of these organisms, namely the split into three above-genus groups: Hominis, Pneumoniae and Spiroplasma–Entomoplasmataceae–Mycoides. However, at the genus level, several important issues were found. For example, recent taxonomic revisions that split the Hominis group into three genera and Entomoplasmataceae into five genera are problematic, as those newly described or emended genera lack clear differentiations in gene content from one another. Moreover, several cases of misclassification were identified. These findings demonstrated the utility of this approach and its potential application to other bacteria.

Understanding the genomic diversity and functional implications of Ganoderma species is crucial for elucidating their evolutionary history and biotechnological potential. Here, we present the first pan-genomic analysis of Ganoderma spp., combining five newly sequenced genomes with ten publicly available genomes. Our comprehensive comparative study unveiled a rich genomic landscape, identifying core genes shared among all Ganoderma strains and species-specific gene sets. Additionally, we identified structural variants impacting the expression of key genes, including insights into the MSH4 gene involved in DNA repair and recombination processes, which exhibits a 440 bp insertion in the promoter region and a leucine-to-serine mutation in the gene body, potentially increasing spore production in the S3 strain. Overall, our study provides valuable insights into the genomic architecture and functional diversity of Ganoderma, paving the way for further research on its evolutionary dynamics, biotechnological applications and pharmaceutical potential.

Despite being a major human pathogen, limited studies have reported RNA modifications in Acinetobacter baumannii. These post-transcriptional modifications play crucial regulatory roles in bacteria and have also been shown to modulate bacterial virulence. Using nanopore sequencing, we characterized RNA modifications in a virulent A. baumannii strain (Ab-C98) under free-living (mid-exponential phase in vitro culture) and during an early stage of infection (3 h post-infection) in Galleria mellonella larvae. Analysis revealed that m5C methylations are essential for ribosome synthesis, while m6A and Ψ are involved in metabolic pathways and translation processes. Iron-chelating genes exbD (m5C and m6A) and feoB (m6A and Ψ) and RNA polymerase subunit rpoC (m6A and Ψ) were selectively modified during infection. This first transcriptome-wide study highlights the potential regulatory roles of m5C, m6A and Ψ modifications in A. baumannii during infection.

Psychrotrophic and thermoduric bacteria are the main reasons for the spoilage of dairy products. This study aims to address the composition and function of psychrotrophic and thermoduric bacteria in eight groups of raw milk samples obtained from Heilongjiang Province and Inner Mongolia (China). Microbial enumeration showed an average total bacterial count of 4.63 log c.f.u. ml−1 and psychrotrophic bacterial counts of 4.82 log c.f.u. ml−1. The mean counts of mesophilic and thermophilic thermoduric bacteria were 3.68 log and 1.81 log c.f.u. ml−1, respectively. Isolated psychrotrophic bacteria (26 genera and 50 species) and mesophilic thermoduric bacteria (20 genera and 32 species) showed high microbial diversity. Through metagenomic and proteomic analyses, significant disparities in the concentration and community structure of psychrotrophic and thermoduric bacteria were observed among different locations. A large number of peptidases were annotated by metagenomics, which may result in milk spoilage. They mainly come from some typical psychrotrophic and thermoduric bacteria, such as Chryseobacterium, Epilithonimonas, Pseudomonas, Psychrobacter, Acinetobacter, Lactococcus, Escherichia and Bacillus. However, the main proteins detected in fresh raw milk were associated with bacterial growth, reproduction and adaptation to cold environments. This investigation provides valuable insights into the microbial communities and protein profiles of raw milk, shedding light on the microbial factors contributing to milk deterioration.

Dehalobacter is a genus of organohalide-respiring bacteria that is recognized for its fastidious growth using reductive dehalogenases (RDases). In the SC05 culture, however, a Dehalobacter population also mineralizes dichloromethane (DCM) produced by chloroform dechlorination using the mec cassette, just downstream of its active RDase. A closed genome of this DCM-mineralizing lineage has previously evaded assembly. Here, we present the genomes of two novel Dehalobacter strains, each of which was assembled from the metagenome of a distinct subculture from SC05. A pangenomic analysis of the Dehalobacter genus, including RDase synteny and phylogenomics, reveals at least five species of Dehalobacter based on average nucleotide identity, RDase and core gene synteny, as well as differential functional genes. An integration hotspot is also pinpointed in the Dehalobacter genome, in which many recombinase islands have accumulated. This nested recombinase island encodes the active RDase and mec cassette in both SC05 Dehalobacter genomes, indicating the transfer of key functional genes between species of Dehalobacter. Horizontal gene transfer between these two novel Dehalobacter strains has implications for the evolutionary history within the SC05 subcultures and of the Dehalobacter genus as a whole, especially regarding adaptation to anthropogenic chemicals.

Vulnerable patients in an intensive care unit (ICU) setting are at high risk of infection from bacteria including gut-colonising Escherichia coli and Klebsiella species. Complex ICU procedures often depend on successful antimicrobial treatment, underscoring the importance of understanding the extent of patient colonisation by multi-drug-resistant organisms (MDROs) in large UK ICUs. Previous work on ICUs globally uncovered high rates of colonisation by transmission of MDROs, but the situation in UK ICUs is less understood. Here, we investigated the diversity and antibiotic resistance gene (ARG) carriage of bacteria present in one of the largest UK ICUs at the Queen Elizabeth Hospital Birmingham (QEHB), focusing primarily on E. coli as both a widespread commensal and a globally disseminated multi-drug-resistant pathogen. Samples were taken during highly restrictive coronavirus disease 2019 (COVID-19) control measures from May to December 2021. Whole-genome and metagenomic sequencing were used to detect and report strain-level colonisation of patients, focusing on E. coli sequence types (STs), their colonisation dynamics and antimicrobial resistance gene carriage. We found a lack of multi-drug resistance (MDR) in the QEHB. Only one carbapenemase-producing organism was isolated, a Citrobacter carrying bla KPC-2. There was no evidence supporting the spread of this strain, and there was little evidence overall of nosocomial acquisition or circulation of colonising E. coli. Whilst 22 different E. coli STs were identified, only 1 strain of the pandemic ST131 lineage was isolated. This ST131 strain was non-MDR and was found to be a clade A strain, associated with low levels of antibiotic resistance. Overall, the QEHB ICU had very low levels of pandemic or MDR strains, a result that may be influenced in part by the strict COVID-19 control measures in place at the time. Employing some of these infection prevention and control measures where reasonable in all ICUs might therefore assist in maintaining low levels of nosocomial MDR.

We used genomic and epidemiological data to assess and compare the population structure and origins of Campylobacter, a major foodborne pathogen, in two neighbouring countries with strong trade and cultural links, similar poultry production systems and frequent movement of people and food products. The most common sequence types (STs) differed between Australia and New Zealand, with many unique to each country. Over half of all STs were represented by a single isolate. Multidrug-resistant (MDR) genotypes were detected in 0.8% of all samples, with no MDR isolates detected in poultry. Quinolone and tetracycline resistant ST6964 was prevalent in New Zealand (10.6% of C. jejuni). Closely related isolates suggested some similar food sources or contacts. We have shown that there is little genetic overlap in human and poultry STs of Campylobacter between the countries, which highlights that this common foodborne pathogen has domestic origins in Australia and New Zealand.

Clustered regularly interspaced short palindromic repeat (CRISPR)-Cas systems are bacterial defences that target bacteriophages and mobile genetic elements. How these defences evolve in novel host environments remains largely unknown. We studied the evolution of the CRISPR-Cas system in Mycoplasma gallisepticum (also named Mycoplasmoides gallisepticum), a bacterial pathogen of poultry that jumped into a passerine host ~30 years ago. Over the decade following the host shift, all isolates displaying a functional CRISPR-Cas system were found not only to harbour completely new sets of spacers, but the DNA protospacer adjacent motif recognized by the main effector M. gallisepticum Cas9 (MgCas9) was also different. These changes in CRISPR-Cas diversity and specificity are consistent with a change in the community of phages and mobile elements infecting M. gallisepticum as it colonized the novel host. In the years following the host shift, we also detected a gradual rise in isolates displaying non-functional MgCas9. After 12 years, all circulating isolates harboured inactive forms only. This loss of CRISPR-Cas function comes at a time when the passerine host is known to have evolved widespread resistance, which in turn drove the evolution of increasing M. gallisepticum virulence through antagonistic coevolution. Such striking concordance in the rise of inactivated forms of CRISPR-Cas and the evolution of host resistance suggests that the inactivation of the CRISPR-Cas system was necessary for enabling adaptive bacterial responses to host-driven selection. We highlight the need to consider both host and pathogen selection pressures on bacteria for understanding the evolution of CRISPR-Cas systems and the key factors driving the emergence of a pathogenic bacterium in a novel host.

Shigellosis is a serious public health issue in many developing countries. The emergence of multidrug-resistant (MDR) Shigella isolates has deepened the treatment difficulty and health burden of shigellosis. China is the largest developing country in the world, but so far, the genome of MDR Shigella isolates has not been well characterized. In this study, 60 clinical isolates of Shigella spp. in Fujian Province, southeast China, from 2005 to 2019 were characterized for drug resistance phenotype, whole-genome sequencing and bioinformatics analysis. The results showed that the MDR rate of Shigella isolates was 100%, among which the resistance rates of cefotaxime, ciprofloxacin and azithromycin were 36.67, 21.67 and 10.00 %, respectively. The positive rate of extended-spectrum beta-lactamase (ESBL)-producing strains was 23.33%. The resistance profiles of Shigella flexneri and Shigella sonnei to some antimicrobials differed. The MDR isolates carried multiple antimicrobial resistance genes, among which blaCTX-M-14 and blaCTX-M-15 mediated ESBL resistance. ‘ISEcp1 -blaCTX−M -IS903’ (type I) and ‘ISEcp1 -blaCTX−M ’ (type II) were the most common genetic environments around the blaCTX-M genes, and plasmids containing these structures included IncFII, IncI1, IncI2 and IncN. The double gene mutation pattern of gyrA and parC resulted in a significant decrease in the sensitivity of S. flexneri to ciprofloxacin. The overall resistance phenotype and genotype concordance rate was 88.50%, and the sensitivity and specificity of the genotype antimicrobial susceptibility test (AST) were 93.35 and 82.53 %, respectively. However, inconsistency occurred between phenotypic and genotype profiles for a few antibiotics. Phylogenomic investigation with core genome multi-locus sequence typing (cgMLST) and SNPs were used to characterize the endemic transmission of these infections in Fujian and their evolutionary origin within the global context. For S. flexneri, Fujian isolates were all limited to PG3 and could be divided into three phylogenetic clusters. The ciprofloxacin-resistant strains were mainly F2a and FXv and assigned to the three clusters with different quinolone resistance-determining region mutation patterns. For S. sonnei, most Fujian strains belonged to Lineage III with genotype 3.7.6, except three isolates of Lineage I with genotype 1.3. The strains carrying the blaCTX-M genes were dispersed, indicating different origins of gene acquisition. Most of the circulating isolates in Fujian Province were not related to major international outbreak lineages and were only endemic to the country. In conclusion, multi-drug resistance of Shigella isolates in Fujian Province was serious, and genome-based laboratory surveillance will be crucial to the clinical treatment and public health measures for shigellosis.

Ascochyta blight caused by the ascomycete Ascochyta rabiei poses a major biotic threat to chickpea (Cicer arietinum) industries worldwide and incurs substantial costs to the Australian multimillion-dollar chickpea industry in both disease control and yield loss. The fungus was introduced to Australia in the 1970s from an unknown source population and, within a few decades, successfully established in all Australian agroecological chickpea-growing regions. Although genetically highly clonal, a broad range of phenotypic variation in terms of aggressiveness exists among the Australian A. rabiei isolates. More recently, highly aggressive isolates capable of causing severe disease symptoms on moderate to highly resistant chickpea cultivars have increased in frequency. To identify genetic loci potentially associated with A. rabiei aggressiveness on Australian chickpea cultivars, we performed deep genome sequencing of 230 isolates collected from a range of agroecological chickpea-growing regions between 2013 and 2020. Population genetic analyses using genome-wide SNP data identified three main clusters of genetically closely related isolates in Australia. Phylogenetic analyses showed that highly aggressive phenotypes developed multiple times independently throughout the phylogeny. The results point to a minor contribution of multiple genetic regions and most likely epigenomic variations to aggressiveness of A. rabiei isolates on Australian chickpea cultivars.

Natural products – small molecules generated by organisms to facilitate ecological interactions – are of great importance to society and are used as antibacterial, antiviral, antifungal and anticancer drugs. However, the role and evolution of these molecules and the fitness benefits they provide to their hosts in their natural habitat remain an outstanding question. In bacteria, the genes that encode the biosynthetic proteins that generate these molecules are organised into discrete loci termed biosynthetic gene clusters (BGCs). In this work, we asked the following question: How are biosynthetic gene clusters organised at the chromosomal level? We sought to answer this using publicly available high-quality assemblies of Micromonospora, an actinomycete genus with members responsible for biosynthesizing notable natural products, such as gentamicin and calicheamicin. By orienting the Micromonospora chromosome around the origin of replication, we demonstrated that Micromonospora has a conserved origin-proximal region, which becomes progressively more disordered towards the antipodes of the origin. We then demonstrated through genome mining of these organisms that the conserved origin-proximal region and the origin-distal region of Micromonospora have distinct populations of BGCs and, in this regard, parallel the organization of Streptomyces, which possesses linear chromosomes. Specifically, the origin-proximal region contains highly syntenous, conserved BGCs predicted to biosynthesize terpenes and a type III polyketide synthase. In contrast, the ori-distal region contains a highly diverse population of BGCs, with many BGCs belonging to unique gene cluster families. These data highlight that genomic plasticity in Micromonospora is locus-specific, and highlight the importance of using high-quality genome assemblies for natural product discovery and guide future natural product discovery by highlighting that biosynthetic novelty may be enriched in specific chromosomal neighbourhoods.