An accessible, efficient and global approach for the large-scale sequencing of bacterial genomes

Analysis of global Aeromonas caviae genomes revealed that strains carrying T6SS are more common in human gastroenteritis than in environmental sources and are often phylogenetically related

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.

Characterizing subgenome recombination and chromosomal imbalances in banana varietal lineages

BACKGROUND

Bananas and plantains (Musa spp.) are among the most important crops worldwide. The cultivated varieties are vegetatively propagated, so their genetic diversity is essentially fixed over time. Musa acuminata, M. balbisiana and M. schizocarpa have provided the named A, B and S subgenomes that predominantly constitute these varieties. Here we aimed to characterize intergenetic recombination and chromosomal imbalances between these A/B/S subgenomes, which often result in copy-number variants (CNVs) leading to changes in gene dosage and phenotype, in a diverse panel of bananas and plantains. This will allow us to characterize varietal lineages better and identify sources of genetic variation.

METHODS

We delimited population structure and clonal lineages in a diverse panel of 188 banana and plantain accessions from the most common cultivars using admixture, principal component and phylogenetic analyses. We used new scalable alignment-based methods, Relative Averaged Alignment (RAA) and Relative Coverage, to infer subgenome composition (AA, AAB, etc.) and interspecific recombination.

RESULTS

In our panel, we identified ten varietal lineages composed of somatic clones, plus three groups of tetraploid accessions. We identified chromosomal exchanges resulting in gains/losses in chromosomal segments (CNVs), particularly in AAB and ABB varieties.

CONCLUSIONS

We demonstrated alignment-based RAA and Relative Coverage can identify subgenome composition and introgressions with similar results to more complex approaches based on single nucleotide polymorphism (SNP) databases. These ab initio species-agnostic methods can be used without sequencing a panel of wild ancestors to find private SNPs, or in recently diverged pools where private SNPs are uncommon. The extensive A/B/S exchanges and the variation in the length of some introgressions between lineages further support multiple foundational events of hybridization and residual backcrossing. Imbalances between A/B/S may have resulted in CNVs and gene dosage variation. Since most edible banana genomes are fixed on time, these CNVs are stable genetic variations probably associated with phenotypic variation for future genetic studies.

The genotype of barley cultivars influences multiple aspects of their associated microbiota via differential root exudate secretion

Plant-associated microbes play vital roles in promoting plant growth and health, with plants secreting root exudates into the rhizosphere to attract beneficial microbes. Exudate composition defines the nature of microbial recruitment, with different plant species attracting distinct microbiota to enable optimal adaptation to the soil environment. To more closely examine the relationship between plant genotype and microbial recruitment, we analysed the rhizosphere microbiomes of landrace (Chevallier) and modern (NFC Tipple) barley (Hordeum vulgare) cultivars. Distinct differences were observed between the plant-associated microbiomes of the 2 cultivars, with the plant-growth promoting rhizobacterial genus Pseudomonas substantially more abundant in the Tipple rhizosphere. Striking differences were also observed between the phenotypes of recruited Pseudomonas populations, alongside distinct genotypic clustering by cultivar. Cultivar-driven Pseudomonas selection was driven by root exudate composition, with the greater abundance of hexose sugars secreted from Tipple roots attracting microbes better adapted to growth on these metabolites and vice versa. Cultivar-driven selection also operates at the molecular level, with both gene expression and the abundance of ecologically relevant loci differing between Tipple and Chevallier Pseudomonas isolates. Finally, cultivar-driven selection is important for plant health, with both cultivars showing a distinct preference for microbes selected by their genetic siblings in rhizosphere transplantation assays.

Ancient and recent origins of shared polymorphisms in yeast

Shared genetic polymorphisms between populations and species can be ascribed to ancestral variation or to more recent gene flow. Here, we mapped shared polymorphisms in Saccharomyces cerevisiae and its sister species Saccharomyces paradoxus, which diverged 4-6 million years ago. We used a dense map of single-nucleotide diagnostic markers (mean distance 15.6 base pairs) in 1,673 sequenced S. cerevisiae isolates to catalogue 3,852 sequence blocks (≥5 consecutive markers) introgressed from S. paradoxus, with most being recent and clade-specific. The highly diverged wild Chinese S. cerevisiae lineages were depleted of introgressed blocks but retained an excess of individual ancestral polymorphisms derived from incomplete lineage sorting, perhaps due to less dramatic population bottlenecks. In the non-Chinese S. cerevisiae lineages, we inferred major hybridization events and detected cases of overlapping introgressed blocks across distinct clades due to either shared histories or convergent evolution. We experimentally engineered, in otherwise isogenic backgrounds, the introgressed PAD1-FDC1 gene pair that independently arose in two S. cerevisiae clades and revealed that it increases resistance against diverse antifungal drugs. Overall, our study retraces the histories of divergence and secondary contacts across S. cerevisiae and S. paradoxus populations and unveils a functional outcome.

Invasive non-typhoidal Salmonella (iNTS) aminoglycoside-resistant ST313 isolates feature unique pathogenic mechanisms to reach the bloodstream

Invasive non-typhoidal Salmonella (iNTS) from the clonal type ST313 (S. Typhimurium ST313) is the primary cause of invasive salmonellosis in Africa. Recently, in Brazil, iNTS ST313 strains have been isolated from different sources, but there is a lack of understanding of the mechanisms behind how these gut bacteria can break the gut barrier and reach the patient's bloodstream. Here, we compare 13 strains of S. Typhimurium ST313, previously unreported isolates, from human blood cultures, investigating aspects of virulence and mechanisms of resistance. Initially, RNAseq analyses between ST13-blood isolate and SL1344 (ST19) prototype revealed 15 upregulated genes directly related to cellular invasion and replication, such as sopD2, sifB, and pipB. Limited information is available about S. Typhimurium ST313 pathogenesis and epidemiology, especially related to the global distribution of strains. Herein, the correlation of strains isolated from different sources in Brazil was employed to compare clinical and non-clinical isolates, a total of 22 genomes were studied by single nucleotide polymorphism (SNPs). The epidemiological analysis of 22 genomes of S. Typhimurium ST313 strains grouped them into three distinct clusters (A, B, and C) by SNP analysis, where cluster A comprised five, group B six, and group C 11. The 13 clinical blood isolates were all resistant to streptomycin, 92.3% of strains were resistant to ampicillin and 15.39% were resistant to kanamycin. The resistance genes acrA, acrB, mdtK, emrB, emrR, mdsA, and mdsB related to the production of efflux pumps were detected in all (100%) strains studied, similar to pathogenic traits investigated. In conclusion, we evidenced that S. Typhimurium ST313 strains isolated in Brazil have unique epidemiology. The elevated frequencies of virulence genes such as sseJ, sopD2, and pipB are a major concern in these Brazilian isolates, showing a higher pathogenic potential.

The genomic epidemiology of shigellosis in South Africa

Shigellosis, a leading cause of diarrhoeal mortality and morbidity globally, predominantly affects children under five years of age living in low- and middle-income countries. While whole genome sequence analysis (WGSA) has been effectively used to further our understanding of shigellosis epidemiology, antimicrobial resistance, and transmission, it has been under-utilised in sub-Saharan Africa. In this study, we applied WGSA to large sub-sample of surveillance isolates from South Africa, collected from 2011 to 2015, focussing on Shigella flexneri 2a and Shigella sonnei. We find each serotype is epidemiologically distinct. The four identified S. flexneri 2a clusters having distinct geographical distributions, and antimicrobial resistance (AMR) and virulence profiles, while the four sub-Clades of S. sonnei varied in virulence plasmid retention. Our results support serotype specific lifestyles as a driver for epidemiological differences, show AMR is not required for epidemiological success in S. flexneri, and that the HIV epidemic may have promoted Shigella population expansion.

Whole-Genome Subtyping Reveals Population Structure and Host Adaptation of Salmonella Typhimurium from Wild Birds

Within-host evolution of bacterial pathogens can lead to host-associated variants of the same species or serovar. Identification and characterization of closely related variants from diverse host species are crucial to public health and host-pathogen adaptation research. However, the work remained largely underexplored at a strain level until the advent of whole-genome sequencing (WGS). Here, we performed WGS-based subtyping and analyses of Salmonella enterica serovar Typhimurium (n = 787) from different wild birds across 18 countries over a 75-year period. We revealed seven avian host-associated S. Typhimurium variants/lineages. These lineages emerged globally over short timescales and presented genetic features distinct from S. Typhimurium lineages circulating among humans and domestic animals. We further showed that, in terms of virulence, host adaptation of these variants was driven by genome degradation. Our results provide a snapshot of the population structure and genetic diversity of S. Typhimurium within avian hosts. We also demonstrate the value of WGS-based subtyping and analyses in unravelling closely related variants at the strain level.

A global genomic analysis of Salmonella Concord reveals lineages with high antimicrobial resistance in Ethiopia

Antimicrobial resistant Salmonella enterica serovar Concord (S. Concord) is known to cause severe gastrointestinal and bloodstream infections in patients from Ethiopia and Ethiopian adoptees, and occasional records exist of S. Concord linked to other countries. The evolution and geographical distribution of S. Concord remained unclear. Here, we provide a genomic overview of the population structure and antimicrobial resistance (AMR) of S. Concord by analysing genomes from 284 historical and contemporary isolates obtained between 1944 and 2022 across the globe. We demonstrate that S. Concord is a polyphyletic serovar distributed among three Salmonella super-lineages. Super-lineage A is composed of eight S. Concord lineages, of which four are associated with multiple countries and low levels of AMR. Other lineages are restricted to Ethiopia and horizontally acquired resistance to most antimicrobials used for treating invasive Salmonella infections in low- and middle-income countries. By reconstructing complete genomes for 10 representative strains, we demonstrate the presence of AMR markers integrated in structurally diverse IncHI2 and IncA/C2 plasmids, and/or the chromosome. Molecular surveillance of pathogens such as S. Concord supports the understanding of AMR and the multi-sector response to the global AMR threat. This study provides a comprehensive baseline data set essential for future molecular surveillance.

Novel Multilocus Sequence Typing and Global Sequence Clustering Schemes for Characterizing the Population Diversity of Streptococcus mitis

Streptococcus mitis is a common oral commensal and an opportunistic pathogen that causes bacteremia and infective endocarditis; however, the species has received little attention compared to other pathogenic streptococcal species. Effective and easy-to-use molecular typing tools are essential for understanding bacterial population diversity and biology, but schemes specific for S. mitis are not currently available. We therefore developed a multilocus sequence typing (MLST) scheme and defined sequence clusters or lineages of S. mitis using a comprehensive global data set of 322 genomes (148 publicly available and 174 newly sequenced). We used internal 450-bp sequence fragments of seven housekeeping genes (accA, gki, hom, oppC, patB, rlmN, and tsf) to define the MLST scheme and derived the global S. mitis sequence clusters using the PopPUNK clustering algorithm. We identified an initial set of 259 sequence types (STs) and 258 global sequence clusters. The schemes showed high concordance (100%), capturing extensive S. mitis diversity with strains assigned to multiple unique STs and global sequence clusters. The tools also identified extensive within- and between-host S. mitis genetic diversity among isolates sampled from a cohort of healthy individuals, together with potential transmission events, supported by both phylogeny and pairwise single nucleotide polymorphism (SNP) distances. Our novel molecular typing and strain clustering schemes for S. mitis allow for the integration of new strain data, are electronically portable at the PubMLST database (https://pubmlst.org/smitis), and offer a standardized approach to understanding the population structure of S. mitis. These robust tools will enable new insights into the epidemiology of S. mitis colonization, disease and transmission.

Multiplex PCR Assay for Clade Typing of Salmonella enterica Serovar Enteritidis

Salmonella enterica serovar Enteritidis is one of the most commonly reported serovars of nontyphoidal Salmonella causing human disease and is responsible for both gastroenteritis and invasive nontyphoidal Salmonella (iNTS) disease worldwide. Whole-genome sequence (WGS) comparison of Salmonella Enteritidis isolates from across the world has identified three distinct clades, global epidemic, Central/East African, and West African, all of which have been implicated in epidemics: the global epidemic clade was linked to poultry-associated gastroenteritis, while the two African clades were related to iNTS disease. However, the distribution and epidemiology of these clades across Africa are poorly understood because identification of these clades currently requires whole-genome sequencing capacity. Here, we report a sensitive, time- and cost-effective real-time PCR assay capable of differentiating between the Salmonella Enteritidis clades to facilitate surveillance and to inform public health responses. The assay described here is limited to previously confirmed S. Enteritidis isolates. IMPORTANCE Challenges in the diagnosis and treatment of invasive Salmonella Enteritidis bloodstream infections in sub-Saharan Africa are responsible for a case fatality rate of approximately 15%. It is important to identify distinct clades of S. Enteritidis in diagnostic laboratories in the African setting to determine the different health outcomes associated with particular outbreaks. Here, we describe the development of a high-quality molecular classification assay for clade typing of S. Enteritidis that is ideal for use in public health laboratories in resource-limited settings.

The evolutionary diversification of the Salmonella artAB toxin locus

Salmonella enterica is a diverse species of bacterial pathogens comprised of >2,500 serovars with variable host ranges and virulence properties. Accumulating evidence indicates that two AB₅-type toxins, typhoid toxin and ArtAB toxin, contribute to the more severe virulence properties of the Salmonella strains that encode them. It was recently discovered that there are two distinct types of artAB-like genetic elements in Salmonella: those that encode ArtAB toxins (artAB elements) and those in which the artA gene is degraded and the ArtB homolog, dubbed PltC, serves as an alternative delivery subunit for typhoid toxin (pltC elements). Here, we take a multifaceted approach to explore the evolutionary diversification of artAB-like genetic elements in Salmonella. We identify 7 subtypes of ArtAB toxins and 4 different PltC sequence groups that are distributed throughout the Salmonella genus. Both artAB and pltC are encoded within numerous diverse prophages, indicating a central role for phages in their evolutionary diversification. Genetic and structural analyses revealed features that distinguish pltC elements from artAB and identified evolutionary adaptations that enable PltC to efficiently engage typhoid toxin A subunits. For both pltC and artAB, we find that the sequences of the B subunits are especially variable, particularly amongst amino acid residues that fine tune the chemical environment of their glycan binding pockets. This study provides a framework to delineate the remarkably complex collection of Salmonella artAB/pltC-like genetic elements and provides a window into the mechanisms of evolution for AB₅-type toxins.

A low-cost genomics workflow enables isolate screening and strain-level analyses within microbiomes

Earth's environments harbor complex consortia of microbes that affect processes ranging from host health to biogeochemical cycles. Understanding their evolution and function is limited by an inability to isolate genomes in a high-throughput manner. Here, we present a workflow for bacterial whole-genome sequencing using open-source labware and the OpenTrons robotics platform, reducing costs to approximately $10 per genome. We assess genomic diversity within 45 gut bacterial species from wild-living chimpanzees and bonobos. We quantify intraspecific genomic diversity and reveal divergence of homologous plasmids between hosts. This enables population genetic analyses of bacterial strains not currently possible with metagenomic data alone.

WGS-Based Lineage and Antimicrobial Resistance Pattern of Salmonella Typhimurium Isolated during 2000-2017 in Peru

Salmonella Typhimurium is associated with foodborne diseases worldwide, including in Peru, and its emerging antibiotic resistance (AMR) is now a global public health problem. Therefore, country-specific monitoring of the AMR emergence is vital to control this pathogen, and in these aspects, whole genome sequence (WGS)—based approaches are better than gene-based analyses. Here, we performed the antimicrobial susceptibility test for ten widely used antibiotics and WGS-based various analyses of 90 S. Typhimurium isolates (human, animal, and environment) from 14 cities of Peru isolated from 2000 to 2017 to understand the lineage and antimicrobial resistance pattern of this pathogen in Peru. Our results suggest that the Peruvian isolates are of Typhimurium serovar and predominantly belong to sequence type ST19. Genomic diversity analyses indicate an open pan-genome, and at least ten lineages are circulating in Peru. A total of 48.8% and 31.0% of isolates are phenotypically and genotypically resistant to at least one antibiotic, while 12.0% are multi-drug resistant (MDR). Genotype−phenotype correlations for ten tested drugs show >80% accuracy, and >90% specificity. Sensitivity above 90% was only achieved for ciprofloxacin and ceftazidime. Two lineages exhibit the majority of the MDR isolates. A total of 63 different AMR genes are detected, of which 30 are found in 17 different plasmids. Transmissible plasmids such as lncI-gamma/k, IncI1-I(Alpha), Col(pHAD28), IncFIB, IncHI2, and lncI2 that carry AMR genes associated with third-generation antibiotics are also identified. Finally, three new non-synonymous single nucleotide variations (SNVs) for nalidixic acid and eight new SNVs for nitrofurantoin resistance are predicted using genome-wide association studies, comparative genomics, and functional annotation. Our analysis provides for the first time the WGS-based details of the circulating S. Typhimurium lineages and their antimicrobial resistance pattern in Peru.

Combination of Whole Genome Sequencing and Metagenomics for Microbiological Diagnostics

Whole genome sequencing (WGS) provides the highest resolution for genome-based species identification and can provide insight into the antimicrobial resistance and virulence potential of a single microbiological isolate during the diagnostic process. In contrast, metagenomic sequencing allows the analysis of DNA segments from multiple microorganisms within a community, either using an amplicon- or shotgun-based approach. However, WGS and shotgun metagenomic data are rarely combined, although such an approach may generate additive or synergistic information, critical for, e.g., patient management, infection control, and pathogen surveillance. To produce a combined workflow with actionable outputs, we need to understand the pre-to-post analytical process of both technologies. This will require specific databases storing interlinked sequencing and metadata, and also involves customized bioinformatic analytical pipelines. This review article will provide an overview of the critical steps and potential clinical application of combining WGS and metagenomics together for microbiological diagnosis.

The Notable Achievements and the Prospects of Bacterial Pathogen Genomics

Throughout the entirety of human history, bacterial pathogens have played an important role and even shaped the fate of civilizations. The application of genomics within the last 27 years has radically changed the way we understand the biology and evolution of these pathogens. In this review, we discuss how the short- (Illumina) and long-read (PacBio, Oxford Nanopore) sequencing technologies have shaped the discipline of bacterial pathogen genomics, in terms of fundamental research (i.e., evolution of pathogenicity), forensics, food safety, and routine clinical microbiology. We have mined and discuss some of the most prominent data/bioinformatics resources such as NCBI pathogens, PATRIC, and Pathogenwatch. Based on this mining, we present some of the most popular sequencing technologies, hybrid approaches, assemblers, and annotation pipelines. A small number of bacterial pathogens are of very high importance, and we also present the wealth of the genomic data for these species (i.e., which ones they are, the number of antimicrobial resistance genes per genome, the number of virulence factors). Finally, we discuss how this discipline will probably be transformed in the near future, especially by transitioning into metagenome-assembled genomes (MAGs), thanks to long-read sequencing.