Different evolutionary trends form the twilight zone of the bacterial pan-genome

Recent emergence of cephalosporin-resistant Salmonella Typhi in India due to the endemic clone acquiring IncFIB(K) plasmid encoding blaCTX-M-15 gene

The emergence and spread of Salmonella Typhi (S. Typhi) resistant to third-generation cephalosporins is a serious global health concern. In this study, we genomically characterized 142 cephalosporin-resistant S. Typhi strains isolated from India. Comparative genome analysis revealed the emergence of a new clone of ceftriaxone-resistant S. Typhi harboring three plasmids of the incompatibility groups IncFIB(K), IncX1, and IncFIB(pHCM2). Among these, the IncFIB(K) plasmid confers resistance to third-generation cephalosporins through the blaCTX-M-15 gene, along with other resistance determinants such as aph(3"), aph(6'), sul2, dfrA14, qnrS, and tet(A). Phylogenetic analysis showed that the isolates from Gujarat (n = 140/142) belong to a distinct subclade (genotype 4.3.1.2.2) within genotype 4.3.1.2 (H58 lineage II). Single nucleotide polymorphism-based phylogenetic analysis of the core genes in IncFIB(K) suggested a close relatedness of the plasmid backbone to that of IncFIB(K) from other Enterobacteriales, indicating that H58 lineage II possesses the capability to acquire MDR plasmids from these organisms. This could indicate the potential onset of a new wave of ceftriaxone-resistant S. Typhi in India. The implementation of control measures-such as vaccination and improved water, sanitation, and hygiene systems-is crucial in areas where MDR or extensively drug-resistant S. Typhi strains are prevalent to curb the spread and impact of these resistant strains.

IMPORTANCE

Typhoid fever remains a global health concern, especially in areas lacking sanitation and clean water. The rise of drug-resistant strains complicates treatment, increasing illness, death, and healthcare expenses. Travel facilitates the spread of these strains worldwide. Multidrug-resistant and extensively drug-resistant (XDR) strains, including those resistant to first-line antibiotics and fluoroquinolones, pose significant challenges. Azithromycin and third-generation cephalosporins are now preferred treatments. Recently, XDR typhoid emerged in Pakistan, resistant even to third-generation cephalosporins. India also faces challenges, with sporadic cases initially declining but now re-emerging. New strains in India show resistance to third-generation cephalosporins due to plasmid acquisition from other bacteria, particularly blaCTX-M-carrying IncFIB(K). Due to the ongoing nature of this outbreak, the data from this study deserve further consideration in order to control its spread in India.

Lineage-specific variation in frequency and hotspots of recombination in invasive Escherichia coli

BACKGROUND

The opportunistic bacterium Escherichia coli can invade normally sterile sites in the human body, potentially leading to life-threatening organ dysfunction and even death. However, our understanding of the evolutionary processes that shape its genetic diversity in this sterile environment remains limited. Here, we aim to quantify the frequency and characteristics of homologous recombination in E. coli from bloodstream infections.

RESULTS

Analysis of 557 short-read genome sequences revealed that the propensity to exchange DNA by homologous recombination varies within a distinct population (bloodstream) at narrow geographic (Dartmouth Hitchcock Medical Center, New Hampshire, USA) and temporal (years 2016 - 2022) scope. We identified the four largest monophyletic sequence clusters in the core genome phylogeny that are represented by prominent sequence types (ST): BAPS1 (mainly ST95), BAPS4 (mainly ST73), BAPS10 (mainly ST131), BAPS14 (mainly ST58). We show that the four dominant clusters vary in different characteristics of recombination: number of single nucleotide polymorphisms due to recombination, number of recombination blocks, cumulative bases in recombination blocks, ratio of probabilities that a given site was altered through recombination and mutation (r/m), and ratio of rates at which recombination and mutation occurred (ρ/θ). Each sequence cluster contains a unique set of antimicrobial resistance (AMR) and virulence genes that have experienced recombination. Common among the four sequence clusters were the recombined virulence genes with functions associated with the Curli secretion channel (csgG) and ferric enterobactin transport (entEF, fepEG). We did not identify any one recombined AMR gene that was present in all four sequence clusters. However, AMR genes mdtABC, baeSR, emrKY and tolC had experienced recombination in sequence clusters BAPS4, BAPS10, and BAPS14. These differences lie in part on the contributions of vertically inherited ancestral recombination and contemporary branch-specific recombination, with some genomes having relatively higher proportions of recombined DNA.

CONCLUSIONS

Our results highlight the variation in the propensity to exchange DNA via homologous recombination within a distinct population at narrow geographic and temporal ranges. Understanding the sources of the genetic variation in invasive E. coli will help inform the implementation of effective strategies to reduce the burden of disease and AMR.

A Syntenic Pangenome for Gardnerella Reveals Taxonomic Boundaries and Stratification of Metabolic and Virulence Potential across Species

Bacterial vaginosis (BV) is a prevalent condition associated with an imbalance in the vaginal microbiota, often involving species of Gardnerella . The taxonomic complexity and inconsistent nomenclature of Gardnerella have impeded progress in understanding the role of specific species in health and disease. In this study, we conducted a comprehensive genomic and pangenomic analysis to resolve taxonomic ambiguities and elucidate metabolic and virulence potential across Gardnerella species. We obtained complete, closed genomes for 42 Gardnerella isolates from women with BV and curated publicly available genome sequences (n = 291). Average nucleotide identity (ANI) analysis, digital DNA-DNA hybridization (dDDH), and the cpn60 gene sequences identified nine species and eleven subspecies within Gardnerella , for which we refined species and subspecies boundaries and proposed updated nomenclature. Pangenome analysis revealed species-specific gene clusters linked to metabolic pathways, virulence factors, and niche adaptations, distinguishing species specialized for mucin degradation in the vaginal environment from those potentially adapted to urinary tract colonization. Notably, we identified lineage-specific evolutionary divergence in gene clusters associated with biofilm formation, carbohydrate metabolism, and antimicrobial resistance. We further discovered the first cryptic plasmids naturally present within the Gardnerella genus. Our findings provide a unified framework for Gardnerella taxonomy and nomenclature, and enhance our understanding of species-specific functional capabilities, with implications for Gardnerella research, diagnostics, and targeted therapeutics in BV.

Competitive behaviors in Serratia marcescens are coordinately regulated by a lifestyle switch frequently inactivated in the clinical environment

Opportunistic bacterial pathogens must compete with other bacteria and switch between host- and environment-adapted states. Type VI secretion systems (T6SSs) occur widely in gram-negative bacteria and can efficiently kill neighboring competitors. We determined the distribution of T6SSs across the genus Serratia and observed that a highly conserved antibacterial T6SS is differentially active between closely related clinical isolates of Serratia marcescens. By combining genomic and experimental approaches, we identified a genus-core two-component system, BetR-Reg1-Reg2, that controls T6SS activity and exhibits frequent inactivating mutations, exclusively in S. marcescens isolates of clinical origin. This regulatory system controls a number of lifestyle-related traits at transcriptional and post-translational levels, including T6SS activity, antibiotic production, motility, and adhesion, with loss of BetR increasing virulence in an in vivo infection model. Our data support a model whereby this system represents a conserved, modular switch from sessile to pioneering and aggressive behavior, which is subject to selection pressure in clinical environments.

Niche-specific evolution and gene exchange of Salmonella in retail pork and chicken

Salmonella exhibits extensive genetic diversity, facilitated by horizontal gene transfer occurring within and between species, playing a pivotal role in this diversification. Nevertheless, most studies focus on clinical and farm animal isolates, and research on the pangenome dynamics of Salmonella isolates from retail stage of the animal food supply chain is limited. Here, we investigated the genomes of 950 Salmonella isolates recovered from retail chicken and pork meats in seven provinces and one municipality of China in 2018. We observed a strong correlation between Salmonella sublineage diversity and the accessory genome with meat type, revealing reduced diversity associated with increased resistance. Importantly, genes associated with antibiotic, biocide, and heavy metal resistance were unevenly distributed in Salmonella from retail chicken and pork. Pork Salmonella isolates showed a higher prevalence of copper and silver resistance genes, while chicken Salmonella isolates displayed a significant predominance of genetic determinants associated with cephalosporin and ciprofloxacin resistance. Moreover, co-occurrence patterns of resistance determinants and their interaction with mobile genetic elements also correlated with meat type. In summary, our findings shed light on how Salmonella achieves their ecological niche success driven by evolution and gene changes in the retail stage of the animal food supply chain.

MosAIC: An annotated collection of mosquito-associated bacteria with high-quality genome assemblies

Mosquitoes transmit medically important human pathogens, including viruses like dengue virus and parasites such as Plasmodium spp., the causative agent of malaria. Mosquito microbiomes are critically important for the ability of mosquitoes to transmit disease-causing agents. However, while large collections of bacterial isolates and genomic data exist for vertebrate microbiomes, the vast majority of work in mosquitoes to date is based on 16S rRNA gene amplicon data that provides limited taxonomic resolution and no functional information. To address this gap and facilitate future studies using experimental microbiome manipulations, we generated a bacterial Mosquito-Associated Isolate Collection (MosAIC) consisting of 392 bacterial isolates with extensive metadata and high-quality draft genome assemblies that are publicly available, both isolates and sequence data, for use by the scientific community. MosAIC encompasses 142 species spanning 29 bacterial families, with members of the Enterobacteriaceae comprising 40% of the collection. Phylogenomic analysis of 3 genera, Enterobacter, Serratia, and Elizabethkingia, reveal lineages of mosquito-associated bacteria isolated from different mosquito species in multiple laboratories. Investigation into species' pangenomes further reveals clusters of genes specific to these lineages, which are of interest for future work to test for functions connected to mosquito host association. Altogether, we describe the generation of a physical collection of mosquito-associated bacterial isolates, their genomic data, and analyses of selected groups in context of genome data from closely related isolates, providing a unique, highly valuable resource for research on bacterial colonisation and adaptation within mosquito hosts. Future efforts will expand the collection to include broader geographic and host species representation, especially from individuals collected from field populations, as well as other mosquito-associated microbes, including fungi, archaea, and protozoa.

Evolutionary dynamics of the accessory genomes of Staphylococcus aureus

Staphylococcus aureus is a ubiquitous commensal and opportunistic bacterial pathogen that can cause a wide gamut of infections, which are exacerbated by the presence of multidrug-resistant and methicillin-resistant S. aureus. S. aureus is genetically heterogeneous and consists of numerous distinct lineages. Using 558 complete genomes of S. aureus, we aim to determine how the accessory genome content among phylogenetic lineages of S. aureus is structured and has evolved. Bayesian hierarchical clustering identified 10 sequence clusters, of which seven contained major sequence types (ST 1, 5, 8, 30, 59, 239, and 398). The seven sequence clusters differed in their accessory gene content, including genes associated with antimicrobial resistance and virulence. Focusing on the two largest clusters, BAPS8 and BAPS10, and each consisting mostly of ST5 and ST8, respectively, we found that the structure and connected components in the co-occurrence networks of accessory genomes varied between them. These differences are explained, in part, by the variation in the rates at which the two sequence clusters gained and lost accessory genes, with the highest rate of gene accumulation occurring recently in their evolutionary histories. We also identified a divergent group within BAPS10 that has experienced high gene gain and loss early in its history. Together, our results show highly variable and dynamic accessory genomes in S. aureus that are structured by the history of the specific lineages that carry them.IMPORTANCEStaphylococcus aureus is an opportunistic, multi-host pathogen that can cause a variety of benign and life-threatening infections. Our results revealed considerable differences in the structure and evolution of the accessory genomes of major lineages within S. aureus. Such genomic variation within a species can have important implications on disease epidemiology, pathogenesis of infection, and interactions with the vertebrate host. Our findings provide important insights into the underlying genetic basis for the success of S. aureus as a highly adaptable and resistant pathogen, which will inform current efforts to control and treat staphylococcal diseases.

A distinct Fusobacterium nucleatum clade dominates the colorectal cancer niche

Fusobacterium nucleatum (Fn), a bacterium present in the human oral cavity and rarely found in the lower gastrointestinal tract of healthy individuals1, is enriched in human colorectal cancer (CRC) tumours2-5. High intratumoural Fn loads are associated with recurrence, metastases and poorer patient prognosis5-8. Here, to delineate Fn genetic factors facilitating tumour colonization, we generated closed genomes for 135 Fn strains; 80 oral strains from individuals without cancer and 55 unique cancer strains cultured from tumours from 51 patients with CRC. Pangenomic analyses identified 483 CRC-enriched genetic factors. Tumour-isolated strains predominantly belong to Fn subspecies animalis (Fna). However, genomic analyses reveal that Fna, considered a single subspecies, is instead composed of two distinct clades (Fna C1 and Fna C2). Of these, only Fna C2 dominates the CRC tumour niche. Inter-Fna analyses identified 195 Fna C2-associated genetic factors consistent with increased metabolic potential and colonization of the gastrointestinal tract. In support of this, Fna C2-treated mice had an increased number of intestinal adenomas and altered metabolites. Microbiome analysis of human tumour tissue from 116 patients with CRC demonstrated Fna C2 enrichment. Comparison of 62 paired specimens showed that only Fna C2 is tumour enriched compared to normal adjacent tissue. This was further supported by metagenomic analysis of stool samples from 627 patients with CRC and 619 healthy individuals. Collectively, our results identify the Fna clade bifurcation, show that specifically Fna C2 drives the reported Fn enrichment in human CRC and reveal the genetic underpinnings of pathoadaptation of Fna C2 to the CRC niche.

Genomic revisitation and reclassification of the genus Providencia

Members of Providencia, although typically opportunistic, can cause severe infections in immunocompromised hosts. Recent advances in genome sequencing provide an opportunity for more precise study of this genus. In this study, we first identified and characterized a novel species named Providencia zhijiangensis sp. nov. It has ≤88.23% average nucleotide identity (ANI) and ≤31.8% in silico DNA-DNA hybridization (dDDH) values with all known Providencia species, which fall significantly below the species-defining thresholds. Interestingly, we found that Providencia stuartii and Providencia thailandensis actually fall under the same species, evidenced by an ANI of 98.59% and a dDDH value of 90.4%. By fusing ANI with phylogeny, we have reclassified 545 genomes within this genus into 20 species, including seven unnamed taxa (provisionally titled Taxon 1-7), which can be further subdivided into 23 lineages. Pangenomic analysis identified 1,550 genus-core genes in Providencia, with coenzymes being the predominant category at 10.56%, suggesting significant intermediate metabolism activity. Resistance analysis revealed that most lineages of the genus (82.61%, 19/23) carry a high number of antibiotic-resistance genes (ARGs) and display diverse resistance profiles. Notably, the majority of ARGs are located on plasmids, underscoring the significant role of plasmids in the resistance evolution within this genus. Three species or lineages (P. stuartii, Taxon 3, and Providencia hangzhouensis L12) that possess the highest number of carbapenem-resistance genes suggest their potential influence on clinical treatment. These findings underscore the need for continued surveillance and study of this genus, particularly due to their role in harboring antibiotic-resistance genes.

IMPORTANCE

The Providencia genus, known to harbor opportunistic pathogens, has been a subject of interest due to its potential to cause severe infections, particularly in vulnerable individuals. Our research offers groundbreaking insights into this genus, unveiling a novel species, Providencia zhijiangensis sp. nov., and highlighting the need for a re-evaluation of existing classifications. Our comprehensive genomic assessment offers a detailed classification of 545 genomes into distinct species and lineages, revealing the rich biodiversity and intricate species diversity within the genus. The substantial presence of antibiotic-resistance genes in the Providencia genus underscores potential challenges for public health and clinical treatments. Our study highlights the pressing need for increased surveillance and research, enriching our understanding of antibiotic resistance in this realm.

Average Nucleotide Identity based Staphylococcus aureus strain grouping allows identification of strain-specific genes in the pangenome

Staphylococcus aureus causes both hospital and community acquired infections in humans worldwide. Due to the high incidence of infection S. aureus is also one of the most sampled and sequenced pathogens today, providing an outstanding resource to understand variation at the bacterial subspecies level. We processed and downsampled 83,383 public S. aureus Illumina whole genome shotgun sequences and 1,263 complete genomes to produce 7,954 representative substrains. Pairwise comparison of core gene Average Nucleotide Identity (ANI) revealed a natural boundary of 99.5% that could be used to define 145 distinct strains within the species. We found that intermediate frequency genes in the pangenome (present in 10-95% of genomes) could be divided into those closely linked to strain background ("strain-concentrated") and those highly variable within strains ("strain-diffuse"). Non-core genes had different patterns of chromosome location; notably, strain-diffuse associated with prophages, strain-concentrated with the vSaβ genome island and rare genes (<10% frequency) concentrated near the origin of replication. Antibiotic genes were enriched in the strain-diffuse class, while virulence genes were distributed between strain-diffuse, strain-concentrated, core and rare classes. This study shows how different patterns of gene movement help create strains as distinct subspecies entities and provide insight into the diverse histories of important S. aureus functions.

Antimicrobial resistance heterogeneity among multidrug-resistant Gram-negative pathogens: Phenotypic, genotypic, and proteomic analysis

Microbes evolve rapidly by modifying their genomes through mutations or through the horizontal acquisition of mobile genetic elements (MGEs) linked with fitness traits such as antimicrobial resistance (AMR), virulence, and metabolic functions. We conducted a multicentric study in India and collected different clinical samples for decoding the genome sequences of bacterial pathogens associated with sepsis, urinary tract infections, and respiratory infections to understand the functional potency associated with AMR and its dynamics. Genomic analysis identified several acquired AMR genes (ARGs) that have a pathogen-specific signature. We observed that blaCTX-M-15, blaCMY-42, blaNDM-5, and aadA(2) were prevalent in Escherichia coli, and blaTEM-1B, blaOXA-232, blaNDM-1, rmtB, and rmtC were dominant in Klebsiella pneumoniae. In contrast, Pseudomonas aeruginosa and Acinetobacter baumannii harbored blaVEB, blaVIM-2, aph(3'), strA/B, blaOXA-23, aph(3') variants, and amrA, respectively. Regardless of the type of ARG, the MGEs linked with ARGs were also pathogen-specific. The sequence type of these pathogens was identified as high-risk international clones, with only a few lineages being predominant and region-specific. Whole-cell proteome analysis of extensively drug-resistant K. pneumoniae, A. baumannii, E. coli, and P. aeruginosa strains revealed differential abundances of resistance-associated proteins in the presence and absence of different classes of antibiotics. The pathogen-specific resistance signatures and differential abundance of AMR-associated proteins identified in this study should add value to AMR diagnostics and the choice of appropriate drug combinations for successful antimicrobial therapy.

Pan-Genome Analysis of Staphylococcus aureus Reveals Key Factors Influencing Genomic Plasticity

The massive quantities of bacterial genomic data being generated have facilitated in-depth analyses of bacteria for pan-genomic studies. However, the pan-genome compositions of one species differed significantly between different studies, so we used Staphylococcus aureus as a model organism to explore the influences driving bacterial pan-genome composition. We selected a series of diverse strains for pan-genomic analysis to explore the pan-genomic composition of S. aureus at the species level and the actual contribution of influencing factors (sequence type [ST], source of isolation, country of isolation, and date of collection) to pan-genome composition. We found that the distribution of core genes in bacterial populations restrained under different conditions differed significantly and showed "local core gene regions" in the same ST. Therefore, we propose that ST may be a key factor driving the dynamic distribution of bacterial genomes and that phylogenetic analyses using whole-genome alignment are no longer appropriate in populations containing multiple ST strains. Pan-genomic analysis showed that some of the housekeeping genes of multilocus sequence typing (MLST) are carried at less than 60% in S. aureus strains. Consequently, we propose a new set of marker genes for the classification of S. aureus, which provides a reference for finding a new set of housekeeping genes to apply to MLST. In this study, we explored the role of driving factors influencing pan-genome composition, providing new insights into the study of bacterial pan-genomes. IMPORTANCE We sought to explore the impact of driving factors influencing pan-genome composition using Staphylococcus aureus as a model organism to provide new insights for the study of bacterial pan-genomes. We believe that the sequence type (ST) of the strains under consideration plays a significant role in the dynamic distribution of bacterial genes. Our findings indicate that there are a certain number of essential genes in Staphylococcus aureus; however, the number of core genes is not as high as previously thought. The new classification method proposed herein suggests that a new set of housekeeping genes more suitable for Staphylococcus aureus must be identified to improve the current classification status of this species.

Multiclonal human origin and global expansion of an endemic bacterial pathogen of livestock

Most new pathogens of humans and animals arise via switching events from distinct host species. However, our understanding of the evolutionary and ecological drivers of successful host adaptation, expansion, and dissemination are limited. Staphylococcus aureus is a major bacterial pathogen of humans and a leading cause of mastitis in dairy cows worldwide. Here we trace the evolutionary history of bovine S. aureus using a global dataset of 10,254 S. aureus genomes including 1,896 bovine isolates from 32 countries in 6 continents. We identified 7 major contemporary endemic clones of S. aureus causing bovine mastitis around the world and traced them back to 4 independent host-jump events from humans that occurred up to 2,500 y ago. Individual clones emerged and underwent clonal expansion from the mid-19th to late 20th century coinciding with the commercialization and industrialization of dairy farming, and older lineages have become globally distributed via established cattle trade links. Importantly, we identified lineage-dependent differences in the frequency of host transmission events between humans and cows in both directions revealing high risk clones threatening veterinary and human health. Finally, pangenome network analysis revealed that some bovine S. aureus lineages contained distinct sets of bovine-associated genes, consistent with multiple trajectories to host adaptation via gene acquisition. Taken together, we have dissected the evolutionary history of a major endemic pathogen of livestock providing a comprehensive temporal, geographic, and gene-level perspective of its remarkable success.

Distinct evolutionary trajectories in the Escherichia coli pangenome occur within sequence types

The Escherichia coli species contains a diverse set of sequence types and there remain important questions regarding differences in genetic content within this population that need to be addressed. Pangenomes are useful vehicles for studying gene content within sequence types. Here, we analyse 21 E. coli sequence type pangenomes using comparative pangenomics to identify variance in both pangenome structure and content. We present functional breakdowns of sequence type core genomes and identify sequence types that are enriched in metabolism, transcription and cell membrane biogenesis genes. We also uncover metabolism genes that have variable core classification, depending on which allele is present. Our comparative pangenomics approach allows for detailed exploration of sequence type pangenomes within the context of the species. We show that ongoing gene gain and loss in the E. coli pangenome is sequence type-specific, which may be a consequence of distinct sequence type-specific evolutionary drivers.

Vibrio cholerae O139 genomes provide a clue to why it may have failed to usher in the eighth cholera pandemic

Cholera is a life-threatening infectious disease that remains an important public health issue in several low and middle-income countries. In 1992, a newly identified O139 Vibrio cholerae temporarily displaced the O1 serogroup. No study has been able to answer why the potential eighth cholera pandemic (8CP) causing V. cholerae O139 emerged so successfully and then died out. We conducted a genomic study, including 330 O139 isolates, covering emergence of the serogroup in 1992 through to 2015. We noted two key genomic evolutionary changes that may have been responsible for the disappearance of genetically distinct but temporally overlapping waves (A-C) of O139. Firstly, as the waves progressed, a switch from a homogenous toxin genotype in wave-A to heterogeneous genotypes. Secondly, a gradual loss of antimicrobial resistance (AMR) with the progression of waves. We hypothesize that these two changes contributed to the eventual epidemiological decline of O139.

Prokaryote pangenomes are dynamic entities

Prokaryote pangenomes are influenced heavily by environmental factors and the opportunity for gene gain and loss events. As the field of pangenome analysis has expanded, so has the need to fully understand the complexity of how eco-evolutionary dynamics shape pangenomes. Here, we describe current models of pangenome evolution and discuss their suitability and accuracy. We suggest that pangenomes are dynamic entities under constant flux, highlighting the influence of two-way interactions between pangenome and environment. New classifications of core and accessory genes are also considered, underscoring the need for continuous evaluation of nomenclature in a fast-moving field. We conclude that future models of pangenome evolution should incorporate eco-evolutionary dynamics to fully encompass their dynamic, changeable nature.