Comprehensive genomic and plasmid characterization of multidrug-resistant bacterial strains by R10.4.1 nanopore sequencing

The escalating prevalence of multidrug-resistant (MDR) bacteria pose a significant public health threat. Understanding the genomic features and deciphering the antibiotic resistance profiles of these pathogens is crucial for development of effective surveillance and treatment strategies. In this study, we employed the R10.4.1 nanopore sequencing technology, specifically through the use of the MinION platform, to analyze eight MDR bacterial strains originating from clinical, ecological and food sources. A single 72-hour sequencing run could yield approximately 12 million reads which covered a total of 34 gigabases (Gbp). The nanopore R10.4.1 data was processed using the Flye assembler, successfully assembling the genomes of eight bacterial strains and their 18 plasmids. Notably, the assemblies generated solely from R10.4.1 nanopore data closely matched those from next-generation sequencing data. Diverse antibiotic resistance patterns and specific resistance genes in the test strains were identified. Hospital strains that exhibited multidrug resistance were found to harbor various resistance genes that encode efflux pumps and extended-spectrum β-lactamases. Environmental and food sources were found to display resistance profiles in a species-specific manner. The composition of structurally complex plasmids in the test strains could also be revealed by analysis of nanopore long reads, which also suggested evidence of horizontal transfer of plasmids between different bacterial species. These findings provide valuable insights into the genetic characteristics of MDR bacteria and demonstrating the practicality of nanopore sequencing technology for detecting of resistance elements in bacterial pathogens.

Epidemiological and Genomic analysis of Vibrio parahaemolyticus isolated from imported travelers at the port of Shanghai, China (2017-2019)

BACKGROUND

Vibrio parahaemolyticus is the predominant etiological agent of seafood-associated foodborne illnesses on a global scale. It is essential to elucidate the mechanisms by which this pathogen disseminates. Given the existing research predominantly concentrates on localized outbreaks, there is a pressing necessity for a comprehensive investigation to capture strains of V. parahaemolyticus cross borders.

RESULTS

This study examined the frequency and genetic attributes of imported V. parahaemolyticus strains among travelers entering Shanghai Port, China, between 2017 and 2019.Through the collection of 21 strains from diverse countries and regions, Southeast Asia was pinpointed as a significant source for the emergence of V. parahaemolyticus. Phylogenetic analysis revealed clear delineation between strains originating from human and environmental sources, emphasizing that underlying genome data of foodborne pathogens is essential for environmental monitoring, food safety and early diagnosis of diseases. Furthermore, our study identified the presence of virulence genes (tdh and tlh) and approximately 120 antibiotic resistance-related genes in the majority of isolates, highlighting their crucial involvement in the pathogenesis of V. parahaemolyticus.

CONCLUSIONS

This research enhanced our comprehension of the worldwide transmission of V. parahaemolyticus and its antimicrobial resistance patterns. The findings have important implications for public health interventions and antimicrobial stewardship strategies, underscoring the necessity for epidemiological surveillance of pathogen at international travel hubs.

Expanding the flexibility of base editing for high-throughput genetic screens in bacteria

Genome-wide screens have become powerful tools for elucidating genotype-to-phenotype relationships in bacteria. Of the varying techniques to achieve knockout and knockdown, CRISPR base editors are emerging as promising options. However, the limited number of available, efficient target sites hampers their use for high-throughput screening. Here, we make multiple advances to enable flexible base editing as part of high-throughput genetic screening in bacteria. We first co-opt the Streptococcus canis Cas9 that exhibits more flexible protospacer-adjacent motif recognition than the traditional Streptococcus pyogenes Cas9. We then expand beyond introducing premature stop codons by mutating start codons. Next, we derive guide design rules by applying machine learning to an essentiality screen conducted in Escherichia coli. Finally, we rescue poorly edited sites by combining base editing with Cas9-induced cleavage of unedited cells, thereby enriching for intended edits. The efficiency of this dual system was validated through a conditional essentiality screen based on growth in minimal media. Overall, expanding the scope of genome-wide knockout screens with base editors could further facilitate the investigation of new gene functions and interactions in bacteria.

Overlapping Streptococcus pyogenes and Streptococcus dysgalactiae subspecies equisimilis household transmission and mobile genetic element exchange

Streptococcus dysgalactiae subspecies equisimilis (SDSE) and Streptococcus pyogenes share skin and throat niches with extensive genomic homology and horizontal gene transfer (HGT) possibly underlying shared disease phenotypes. It is unknown if cross-species transmission interaction occurs. Here, we conduct a genomic analysis of a longitudinal household survey in remote Australian First Nations communities for patterns of cross-species transmission interaction and HGT. Collected from 4547 person-consultations, we analyse 294 SDSE and 315 S. pyogenes genomes. We find SDSE and S. pyogenes transmission intersects extensively among households and show that patterns of co-occurrence and transmission links are consistent with independent transmission without inter-species interference. We identify at least one of three near-identical cross-species mobile genetic elements (MGEs) carrying antimicrobial resistance or streptodornase virulence genes in 55 (19%) SDSE and 23 (7%) S. pyogenes isolates. These findings demonstrate co-circulation of both pathogens and HGT in communities with a high burden of streptococcal disease, supporting a need to integrate SDSE and S. pyogenes surveillance and control efforts.

Uncovering the boundaries of Campylobacter species through large-scale phylogenetic and nucleotide identity analyses

Campylobacter species are typically helical shaped, Gram-negative, and non-spore-forming bacteria. Species in this genus include established foodborne and animal pathogens as well as emerging pathogens. The accumulation of genomic data from the Campylobacter genus has increased exponentially in recent years, accompanied by the discovery of putative new species. At present, the lack of a standardized species boundary complicates distinguishing established and novel species. We defined the Campylobacter genus core genome (500 loci) using publicly available Campylobacter complete genomes (n = 498) and constructed a core genome phylogeny using 2,193 publicly available Campylobacter genomes to examine inter-species diversity and species boundaries. Utilizing 8,440 Campylobacter genomes representing 33 species and 8 subspecies, we found species delineation based on an average nucleotide identity (ANI) cutoff of 94.2% is consistent with the core genome phylogeny. We identified 60 ANI genomic species that delineated Campylobacter species in concordance with previous comparative genetic studies. All pairwise ANI genomic species pairs had in silico DNA-DNA hybridization scores of less than 70%, supporting their delineation as separate species. We provide the tool Campylobacter Genomic Species typer (CampyGStyper) that assigns ANI genomic species to query genomes based on ANI similarities to medoid genomes from each ANI genomic species with an accuracy of 99.96%. The ANI genomic species definitions proposed here allow consistent species definition in the Campylobacter genus and will facilitate the detection of novel species in the future.IMPORTANCEIn recent years, Campylobacter has gained recognition as the leading cause of bacterial gastroenteritis worldwide, leading to a substantial rise in the collection of genomic data of the Campylobacter genus in public databases. Currently, a standardized Campylobacter species boundary at the genomic level is absent, leading to challenges in detecting emerging pathogens and defining putative novel species within this genus. We used a comprehensive representation of genomes of the Campylobacter genus to construct a core genome phylogenetic tree. Furthermore, we found an average nucleotide identity (ANI) of 94.2% as the optimal cutoff to define the Campylobacter species. Using this cutoff, we identified 60 ANI genomic species which provided a standardized species definition and nomenclature. Importantly, we have developed Campylobacter Genomic Species typer (CampyGStyper), which can robustly and accurately assign these ANI genomic species to Campylobacter genomes, thereby aiding pathogen surveillance and facilitating evolutionary and epidemiological studies of existing and emerging pathogens in the genus Campylobacter.

Novel mechanisms of diversity generation in Acinetobacter baumannii resistance islands driven by Tn7-like elements

Mobile genetic elements play an important role in the acquisition of antibiotic and biocide resistance, especially through the formation of resistance islands in bacterial chromosomes. We analyzed the contribution of Tn7-like transposons to island formation and diversification in the nosocomial pathogen Acinetobacter baumannii and identified four separate families that recognize different integration sites. One integration site is within the comM gene and coincides with the previously described Tn6022 elements suggested to account for the AbaR resistance island. We established Tn6022 in a heterologous E. coli host and confirmed basic features of transposition into the comM attachment site and the use of a novel transposition protein. By analyzing population features within Tn6022 elements we identified two potential novel transposon-encoded diversification mechanisms with this dynamic genetic island. The activities of these diversification features were confirmed in E. coli. One was a novel natural gain-of-activity allele that could function to broaden transposition targeting. The second was a transposon-encoded hybrid dif-like site that parasitizes the host dimer chromosome resolution system to function with its own tyrosine recombinase. This work establishes a highly active Tn7-like transposon that harnesses novel features allowing the spread and diversification of genetic islands in pathogenic bacteria.

Genetically encoded transcriptional plasticity underlies stress adaptation in Mycobacterium tuberculosis

Transcriptional regulation is a critical adaptive mechanism that allows bacteria to respond to changing environments, yet the concept of transcriptional plasticity (TP) - the variability of gene expression in response to environmental changes - remains largely unexplored. In this study, we investigate the genome-wide TP profiles of Mycobacterium tuberculosis (Mtb) genes by analyzing 894 RNA sequencing samples derived from 73 different environmental conditions. Our data reveal that Mtb genes exhibit significant TP variation that correlates with gene function and gene essentiality. We also find that critical genetic features, such as gene length, GC content, and operon size independently impose constraints on TP, beyond trans-regulation. By extending our analysis to include two other Mycobacterium species -- M. smegmatis and M. abscessus -- we demonstrate a striking conservation of the TP landscape. This study provides a comprehensive understanding of the TP exhibited by mycobacteria genes, shedding light on this significant, yet understudied, genetic feature encoded in bacterial genomes.

Clinical manifestations, antimicrobial resistance and genomic feature analysis of multidrug-resistant Elizabethkingia strains

BACKGROUND

Elizabethkingia is emerging as an opportunistic pathogen in humans. The aim of this study was to investigate the clinical epidemiology, antimicrobial susceptibility, virulence factors, and genome features of Elizabethkingia spp.

METHODS

Clinical data from 71 patients who were diagnosed with Elizabethkingia-induced pneumonia and bacteremia between August 2019 and September 2021 were analyzed. Whole-genome sequencing was performed on seven isolates, and the results were compared with a dataset of 83 available Elizabethkingia genomes. Genomic features, Kyoto Encyclopedia of Genes and Genomes (KEGG) results and clusters of orthologous groups (COGs) were analyzed.

RESULTS

The mean age of the patients was 56.9 ± 20.7 years, and the in-hospital mortality rate was 29.6% (21/71). Elizabethkingia strains were obtained mainly from intensive care units (36.6%, 26/71) and emergency departments (32.4%, 23/71). The majority of the strains were isolated from respiratory tract specimens (85.9%, 61/71). All patients had a history of broad-spectrum antimicrobial exposure. Hospitalization for invasive mechanical ventilation or catheter insertion was found to be a risk factor for infection. The isolates displayed a high rate of resistance to cephalosporins and carbapenems, but all were susceptible to minocycline and colistin. Genomic analysis identified five β-lactamase genes (blaGOB, blaBlaB, blaCME, blaOXA, and blaTEM) responsible for β-lactam resistance and virulence genes involved in stress adaptation (ureB/G, katA/B, and clpP), adherence (groEL, tufA, and htpB) and immune modulation (gmd, tviB, cps4J, wbtIL, cap8E/D/G, and rfbC). Functional analysis of the COGs revealed that "metabolism" constituted the largest category within the core genome, while "information storage and processing" was predominant in both the accessory and unique genomes. The unique genes in our 7 strains were mostly enriched in KEGG pathways related to microRNAs in cancer, drug resistance (β-lactam and vancomycin), ABC transporters, biological metabolism and biosynthesis, and nucleotide excision repair mechanisms.

CONCLUSION

The Elizabethkingia genus exhibits multidrug resistance and carries carbapenemase genes. This study presents a comparative genomic analysis of Elizabethkingia, providing knowledge that facilitates a better understanding of this microorganism.

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.

How do bacterial endosymbionts work with so few genes?

The move from a free-living environment to a long-term residence inside a host eukaryotic cell has profound effects on bacterial function. While endosymbioses are found in many eukaryotes, from protists to plants to animals, the bacteria that form these host-beneficial relationships are even more diverse. Endosymbiont genomes can become radically smaller than their free-living relatives, and their few remaining genes show extreme compositional biases. The details of how these reduced and divergent gene sets work, and how they interact with their host cell, remain mysterious. This Unsolved Mystery reviews how genome reduction alters endosymbiont biology and highlights a "tipping point" where the loss of the ability to build a cell envelope coincides with a marked erosion of translation-related genes.

Modeling the mosaic structure of bacterial genomes to infer their evolutionary history

The chronology and phylogeny of bacterial evolution are difficult to reconstruct due to a scarce fossil record. The analysis of bacterial genomes remains challenging because of large sequence divergence, the plasticity of bacterial genomes due to frequent gene loss, horizontal gene transfer, and differences in selective pressure from one locus to another. Therefore, taking advantage of the rich and rapidly accumulating genomic data requires accurate modeling of genome evolution. An important technical consideration is that loci with high effective mutation rates may diverge beyond the detection limit of the alignment algorithms used, biasing the genome-wide divergence estimates toward smaller divergences. In this article, we propose a novel method to gain insight into bacterial evolution based on statistical properties of genome comparisons. We find that the length distribution of sequence matches is shaped by the effective mutation rates of different loci, by the horizontal transfers, and by the aligner sensitivity. Based on these inputs, we build a model and show that it accounts for the empirically observed distributions, taking the Enterobacteriaceae family as an example. Our method allows to distinguish segments of vertical and horizontal origins and to estimate the time divergence and exchange rate between any pair of taxa from genome-wide alignments. Based on the estimated time divergences, we construct a time-calibrated phylogenetic tree to demonstrate the accuracy of the method.

bacLIFE: a user-friendly computational workflow for genome analysis and prediction of lifestyle-associated genes in bacteria

Bacteria have an extensive adaptive ability to live in close association with eukaryotic hosts, exhibiting detrimental, neutral or beneficial effects on host growth and health. However, the genes involved in niche adaptation are mostly unknown and their functions poorly characterized. Here, we present bacLIFE ( https://github.com/Carrion-lab/bacLIFE ) a streamlined computational workflow for genome annotation, large-scale comparative genomics, and prediction of lifestyle-associated genes (LAGs). As a proof of concept, we analyzed 16,846 genomes from the Burkholderia/Paraburkholderia and Pseudomonas genera, which led to the identification of hundreds of genes potentially associated with a plant pathogenic lifestyle. Site-directed mutagenesis of 14 of these predicted LAGs of unknown function, followed by plant bioassays, showed that 6 predicted LAGs are indeed involved in the phytopathogenic lifestyle of Burkholderia plantarii and Pseudomonas syringae pv. phaseolicola. These 6 LAGs encompassed a glycosyltransferase, extracellular binding proteins, homoserine dehydrogenases and hypothetical proteins. Collectively, our results highlight bacLIFE as an effective computational tool for prediction of LAGs and the generation of hypotheses for a better understanding of bacteria-host interactions.

Redefining the treponemal history through pre-Columbian genomes from Brazil

The origins of treponemal diseases have long remained unknown, especially considering the sudden onset of the first syphilis epidemic in the late 15th century in Europe and its hypothesized arrival from the Americas with Columbus' expeditions1,2. Recently, ancient DNA evidence has revealed various treponemal infections circulating in early modern Europe and colonial-era Mexico3-6. However, there has been to our knowledge no genomic evidence of treponematosis recovered from either the Americas or the Old World that can be reliably dated to the time before the first trans-Atlantic contacts. Here, we present treponemal genomes from nearly 2,000-year-old human remains from Brazil. We reconstruct four ancient genomes of a prehistoric treponemal pathogen, most closely related to the bejel-causing agent Treponema pallidum endemicum. Contradicting the modern day geographical niche of bejel in the arid regions of the world, the results call into question the previous palaeopathological characterization of treponeme subspecies and showcase their adaptive potential. A high-coverage genome is used to improve molecular clock date estimations, placing the divergence of modern T. pallidum subspecies firmly in pre-Columbian times. Overall, our study demonstrates the opportunities within archaeogenetics to uncover key events in pathogen evolution and emergence, paving the way to new hypotheses on the origin and spread of treponematoses.

Bacterial Mutation During Seasonal Epidemics

Rapidly evolving bacterial pathogens pose a unique challenge for long-term plant disease management. In this study, we investigated the types and rate of mutations in bacterial populations during seasonal disease epidemics. Two phylogenetically distinct strains of the bacterial spot pathogen, Xanthomonas perforans, were marked, released in tomato fields, and recaptured at several time points during the growing season. Genomic variations in recaptured isolates were identified by comparative analysis of their whole-genome sequences. In total, 180 unique variations (116 substitutions, 57 insertions/deletions, and 7 structural variations) were identified from 300 genomes, resulting in the overall host-associated mutation rate of ∼0.3 to 0.9/genome/week. This result serves as a benchmark for bacterial mutation during epidemics in similar pathosystems. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Transcription-replication interactions reveal bacterial genome regulation

Organisms determine the transcription rates of thousands of genes through a few modes of regulation that recur across the genome1. In bacteria, the relationship between the regulatory architecture of a gene and its expression is well understood for individual model gene circuits2,3. However, a broader perspective of these dynamics at the genome scale is lacking, in part because bacterial transcriptomics has hitherto captured only a static snapshot of expression averaged across millions of cells4. As a result, the full diversity of gene expression dynamics and their relation to regulatory architecture remains unknown. Here we present a novel genome-wide classification of regulatory modes based on the transcriptional response of each gene to its own replication, which we term the transcription-replication interaction profile (TRIP). Analysing single-bacterium RNA-sequencing data, we found that the response to the universal perturbation of chromosomal replication integrates biological regulatory factors with biophysical molecular events on the chromosome to reveal the local regulatory context of a gene. Whereas the TRIPs of many genes conform to a gene dosage-dependent pattern, others diverge in distinct ways, and this is shaped by factors such as intra-operon position and repression state. By revealing the underlying mechanistic drivers of gene expression heterogeneity, this work provides a quantitative, biophysical framework for modelling replication-dependent expression dynamics.

Resequencing and characterization of the first Corynebacterium pseudotuberculosis genome isolated from camel

Background

Corynebacterium pseudotuberculosis is a zoonotic Gram-positive bacterial pathogen known to cause different diseases in many mammals, including lymph node abscesses in camels. Strains from biovars equi and ovis of C. pseudotuberculosis can infect camels. Comparative genomics could help to identify features related to host adaptation, and currently strain Cp162 from biovar equi is the only one from camel with a sequenced genome.

Methods

In this work, we compared the quality of three genome assemblies of strain Cp162 that used data from the DNA sequencing platforms SOLiD v3 Plus, IonTorrent PGM, and Illumina HiSeq 2500 with an optical map and investigate the unique features of this strain. For this purpose, we applied comparative genomic analysis on the different Cp162 genome assembly versions and included other 129 genomes from the same species.

Results

Since the first version of the genome, there was an increase of 88 Kbp and 121 protein-coding sequences, a decrease of pseudogenes from 139 to 53, and two inversions and one rearrangement corrected. We identified 30 virulence genes, none associated to the camel host, and the genes rpob2 and rbpA predicted to confer resistance to rifampin. In comparison to 129 genomes of the same species, strain Cp162 has four genes exclusively present, two of them code transposases and two truncated proteins, and the three exclusively absent genes lysG, NUDIX domain protein, and Hypothetical protein. All 130 genomes had the rifampin resistance genes rpob2 and rbpA. Our results found no unique gene that could be associated with tropism to camel host, and further studies should include more genomes and genome-wide association studies testing for genes and SNPs.

Towards estimating the number of strains that make up a natural bacterial population

What a strain is and how many strains make up a natural bacterial population remain elusive concepts despite their apparent importance for assessing the role of intra-population diversity in disease emergence or response to environmental perturbations. To advance these concepts, we sequenced 138 randomly selected Salinibacter ruber isolates from two solar salterns and assessed these genomes against companion short-read metagenomes from the same samples. The distribution of genome-aggregate average nucleotide identity (ANI) values among these isolates revealed a bimodal distribution, with four-fold lower occurrence of values between 99.2% and 99.8% relative to ANI >99.8% or <99.2%, revealing a natural "gap" in the sequence space within species. Accordingly, we used this ANI gap to define genomovars and a higher ANI value of >99.99% and shared gene-content >99.0% to define strains. Using these thresholds and extrapolating from how many metagenomic reads each genomovar uniquely recruited, we estimated that -although our 138 isolates represented about 80% of the Sal. ruber population- the total population in one saltern pond is composed of 5,500 to 11,000 genomovars, the great majority of which appear to be rare in-situ. These data also revealed that the most frequently recovered isolate in lab media was often not the most abundant genomovar in-situ, suggesting that cultivation biases are significant, even in cases that cultivation procedures are thought to be robust. The methodology and ANI thresholds outlined here should represent a useful guide for future microdiversity surveys of additional microbial species.

gcPathogen: a comprehensive genomic resource of human pathogens for public health

Here, we present the manually curated Global Catalogue of Pathogens (gcPathogen), an extensive genomic resource designed to facilitate rapid and accurate pathogen analysis, epidemiological exploration and monitoring of antibiotic resistance features and virulence factors. The catalogue seamlessly integrates and analyzes genomic data and associated metadata for human pathogens isolated from infected patients, animal hosts, food and the environment. The pathogen list is supported by evidence from medical or government pathogenic lists and publications. The current version of gcPathogen boasts an impressive collection of 1 164 974 assemblies comprising 986 044 strains from 497 bacterial taxa, 4794 assemblies encompassing 4319 strains from 265 fungal taxa, 89 965 assemblies featuring 13 687 strains from 222 viral taxa, and 646 assemblies including 387 strains from 159 parasitic taxa. Through this database, researchers gain access to a comprehensive 'one-stop shop' that facilitates global, long-term public health surveillance while enabling in-depth analysis of genomes, sequence types, antibiotic resistance genes, virulence factors and mobile genetic elements across different countries, diseases and hosts. To access and explore the data and statistics, an interactive web interface has been developed, which can be accessed at https://nmdc.cn/gcpathogen/. This user-friendly platform allows seamless querying and exploration of the extensive information housed within the gcPathogen database.

RefSeq and the prokaryotic genome annotation pipeline in the age of metagenomes

The Reference Sequence (RefSeq) project at the National Center for Biotechnology Information (NCBI) contains over 315 000 bacterial and archaeal genomes and 236 million proteins with up-to-date and consistent annotation. In the past 3 years, we have expanded the diversity of the RefSeq collection by including the best quality metagenome-assembled genomes (MAGs) submitted to INSDC (DDBJ, ENA and GenBank), while maintaining its quality by adding validation checks. Assemblies are now more stringently evaluated for contamination and for completeness of annotation prior to acceptance into RefSeq. MAGs now account for over 17000 assemblies in RefSeq, split over 165 orders and 362 families. Changes in the Prokaryotic Genome Annotation Pipeline (PGAP), which is used to annotate nearly all RefSeq assemblies include better detection of protein-coding genes. Nearly 83% of RefSeq proteins are now named by a curated Protein Family Model, a 4.7% increase in the past three years ago. In addition to literature citations, Enzyme Commission numbers, and gene symbols, Gene Ontology terms are now assigned to 48% of RefSeq proteins, allowing for easier multi-genome comparison. RefSeq is found at https://www.ncbi.nlm.nih.gov/refseq/. PGAP is available as a stand-alone tool able to produce GenBank-ready files at https://github.com/ncbi/pgap.

Contingency, repeatability, and predictability in the evolution of a prokaryotic pangenome

Pangenomes exhibit remarkable variability in many prokaryotic species, much of which is maintained through the processes of horizontal gene transfer and gene loss. Repeated acquisitions of near-identical homologs can easily be observed across pangenomes, leading to the question of whether these parallel events potentiate similar evolutionary trajectories, or whether the remarkably different genetic backgrounds of the recipients mean that postacquisition evolutionary trajectories end up being quite different. In this study, we present a machine learning method that predicts the presence or absence of genes in the Escherichia coli pangenome based on complex patterns of the presence or absence of other accessory genes within a genome. Our analysis leverages the repeated transfer of genes through the E. coli pangenome to observe patterns of repeated evolution following similar events. We find that the presence or absence of a substantial set of genes is highly predictable from other genes alone, indicating that selection potentiates and maintains gene-gene co-occurrence and avoidance relationships deterministically over long-term bacterial evolution and is robust to differences in host evolutionary history. We propose that at least part of the pangenome can be understood as a set of genes with relationships that govern their likely cohabitants, analogous to an ecosystem's set of interacting organisms. Our findings indicate that intragenomic gene fitness effects may be key drivers of prokaryotic evolution, influencing the repeated emergence of complex gene-gene relationships across the pangenome.

Bioactive glycans in a microbiome-directed food for children with malnutrition

Evidence is accumulating that perturbed postnatal development of the gut microbiome contributes to childhood malnutrition1-4. Here we analyse biospecimens from a randomized, controlled trial of a microbiome-directed complementary food (MDCF-2) that produced superior rates of weight gain compared with a calorically more dense conventional ready-to-use supplementary food in 12-18-month-old Bangladeshi children with moderate acute malnutrition4. We reconstructed 1,000 bacterial genomes (metagenome-assembled genomes (MAGs)) from the faecal microbiomes of trial participants, identified 75 MAGs of which the abundances were positively associated with ponderal growth (change in weight-for-length Z score (WLZ)), characterized changes in MAG gene expression as a function of treatment type and WLZ response, and quantified carbohydrate structures in MDCF-2 and faeces. The results reveal that two Prevotella copri MAGs that are positively associated with WLZ are the principal contributors to MDCF-2-induced expression of metabolic pathways involved in utilizing the component glycans of MDCF-2. The predicted specificities of carbohydrate-active enzymes expressed by their polysaccharide-utilization loci are correlated with (1) the in vitro growth of Bangladeshi P. copri strains, possessing varying degrees of polysaccharide-utilization loci and genomic conservation with these MAGs, in defined medium containing different purified glycans representative of those in MDCF-2, and (2) the levels of faecal carbohydrate structures in the trial participants. These associations suggest that identifying bioactive glycan structures in MDCFs metabolized by growth-associated bacterial taxa will help to guide recommendations about their use in children with acute malnutrition and enable the development of additional formulations.

Phosphate-related genomic islands as drivers of environmental adaptation in the streamlined marine alphaproteobacterial HIMB59

IMPORTANCE

These results shed light on the evolutionary strategies of microbes with streamlined genomes to adapt and survive in the oligotrophic conditions that dominate the surface waters of the global ocean. At the individual level, these microbes have been subjected to evolutionary constraints that have led to a more efficient use of nutrients, removing non-essential genes named as "streamlining theory." However, at the population level, they conserve a highly diverse gene pool in flexible genomic islands resulting in polyclonal populations on the same genomic background as an evolutionary response to environmental pressures. Localization of these islands at equivalent positions in the genome facilitates horizontal transfer between clonal lineages. This high level of environmental genomic heterogeneity could explain their cosmopolitan distribution. In the case of the order HIMB59 within the class Alphaproteobacteria, two factors exert evolutionary pressure and determine this intraspecific diversity: phages and the concentration of P in the environment.

Genome-wide transformation reveals extensive exchange across closely related Bacillus species

Bacterial transformation is an important mode of horizontal gene transfer that helps spread genetic material across species boundaries. Yet, the factors that pose barriers to genome-wide cross-species gene transfer are poorly characterized. Here, we develop a replacement accumulation assay to study the effects of genomic distance on transfer dynamics. Using Bacillus subtilis as recipient and various species of the genus Bacillus as donors, we find that the rate of orthologous replacement decreases exponentially with the divergence of their core genomes. We reveal that at least 96% of the B. subtilis core genes are accessible to replacement by alleles from Bacillus spizizenii. For the more distantly related Bacillus atrophaeus, gene replacement events cluster at genomic locations with high sequence identity and preferentially replace ribosomal genes. Orthologous replacement also creates mosaic patterns between donor and recipient genomes, rearranges the genome architecture, and governs gain and loss of accessory genes. We conclude that cross-species gene transfer is dominated by orthologous replacement of core genes which occurs nearly unrestricted between closely related species. At a lower rate, the exchange of accessory genes gives rise to more complex genome dynamics.

The Helicobacter pylori Genome Project: insights into H. pylori population structure from analysis of a worldwide collection of complete genomes

Helicobacter pylori, a dominant member of the gastric microbiota, shares co-evolutionary history with humans. This has led to the development of genetically distinct H. pylori subpopulations associated with the geographic origin of the host and with differential gastric disease risk. Here, we provide insights into H. pylori population structure as a part of the Helicobacter pylori Genome Project (HpGP), a multi-disciplinary initiative aimed at elucidating H. pylori pathogenesis and identifying new therapeutic targets. We collected 1011 well-characterized clinical strains from 50 countries and generated high-quality genome sequences. We analysed core genome diversity and population structure of the HpGP dataset and 255 worldwide reference genomes to outline the ancestral contribution to Eurasian, African, and American populations. We found evidence of substantial contribution of population hpNorthAsia and subpopulation hspUral in Northern European H. pylori. The genomes of H. pylori isolated from northern and southern Indigenous Americans differed in that bacteria isolated in northern Indigenous communities were more similar to North Asian H. pylori while the southern had higher relatedness to hpEastAsia. Notably, we also found a highly clonal yet geographically dispersed North American subpopulation, which is negative for the cag pathogenicity island, and present in 7% of sequenced US genomes. We expect the HpGP dataset and the corresponding strains to become a major asset for H. pylori genomics.

Complex interactions between diverse mobile genetic elements drive the evolution of metal-resistant bacterial genomes

In this study, we compared the genomes of three metal-resistant bacteria isolated from mercury-contaminated soil. We identified diverse and novel MGEs with evidence of multiple LGT events shaping their genomic structure and heavy metal resistance. Among the three metal-resistant strains, Sphingobium sp SA2 and Sphingopyxis sp SE2 were resistant to multiple metals including mercury, cadmium, copper, zinc and lead. Pseudoxanthomonas sp SE1 showed resistance to mercury only. Whole genome sequencing by Illumina and Oxford Nanopore technologies was undertaken to obtain comprehensive genomic data. The Sphingobium and Sphingopyxis strains contained multiple chromosomes and plasmids, whereas the Pseudoxanthomonas strain contained one circular chromosome. Consistent with their metal resistance profiles, the strains of Sphingobium and Sphingopyxis contained a higher quantity of diverse metal resistance genes across their chromosomes and plasmids compared to the single-metal resistant Pseudoxanthomonas SE1. In all three strains, metal resistance genes were principally associated with various novel MGEs including genomic islands (GIs), integrative conjugative elements (ICEs), transposons, insertion sequences (IS), recombinase in trio (RIT) elements and group II introns, indicating their importance in facilitating metal resistance adaptation in a contaminated environment. In the Pseudoxanthomonas strain, metal resistance regions were largely situated on a GI. The chromosomes of the strains of Sphingobium and Sphingopyxis contained multiple metal resistance regions, which were likely acquired by several GIs, ICEs, numerous IS elements, several Tn3 family transposons and RIT elements. Two of the plasmids of Sphingobium were impacted by Tn3 family transposons and ISs likely integrating metal resistance genes. The two plasmids of Sphingopyxis harboured transposons, IS elements, an RIT element and a group II intron. This study provides a comprehensive annotation of complex genomic regions of metal resistance associated with novel MGEs. It highlights the critical importance of LGT in the evolution of metal resistance of bacteria in contaminated environments.

Statistical prediction of microbial metabolic traits from genomes

The metabolic activity of microbial communities is central to their role in biogeochemical cycles, human health, and biotechnology. Despite the abundance of sequencing data characterizing these consortia, it remains a serious challenge to predict microbial metabolic traits from sequencing data alone. Here we culture 96 bacterial isolates individually and assay their ability to grow on 10 distinct compounds as a sole carbon source. Using these data as well as two existing datasets, we show that statistical approaches can accurately predict bacterial carbon utilization traits from genomes. First, we show that classifiers trained on gene content can accurately predict bacterial carbon utilization phenotypes by encoding phylogenetic information. These models substantially outperform predictions made by constraint-based metabolic models automatically constructed from genomes. This result solidifies our current knowledge about the strong connection between phylogeny and metabolic traits. However, phylogeny-based predictions fail to predict traits for taxa that are phylogenetically distant from any strains in the training set. To overcome this we train improved models on gene presence/absence to predict carbon utilization traits from gene content. We show that models that predict carbon utilization traits from gene presence/absence can generalize to taxa that are phylogenetically distant from the training set either by exploiting biochemical information for feature selection or by having sufficiently large datasets. In the latter case, we provide evidence that a statistical approach can identify putatively mechanistic genes involved in metabolic traits. Our study demonstrates the potential power for predicting microbial phenotypes from genotypes using statistical approaches.

Distributed genotyping and clustering of Neisseria strains reveal continual emergence of epidemic meningococcus over a century

Core genome multilocus sequence typing (cgMLST) is commonly used to classify bacterial strains into different types, for taxonomical and epidemiological applications. However, cgMLST schemes require central databases for the nomenclature of new alleles and sequence types, which must be synchronized worldwide and involve increasingly intensive calculation and storage demands. Here, we describe a distributed cgMLST (dcgMLST) scheme that does not require a central database of allelic sequences and apply it to study evolutionary patterns of epidemic and endemic strains of the genus Neisseria. We classify 69,994 worldwide Neisseria strains into multi-level clusters that assign species, lineages, and local disease outbreaks. We divide Neisseria meningitidis into 168 endemic lineages and three epidemic lineages responsible for at least 9 epidemics in the past century. According to our analyses, the epidemic and endemic lineages experienced very different population dynamics in the past 100 years. Epidemic lineages repetitively emerged from endemic lineages, disseminated worldwide, and apparently disappeared rapidly afterward. We propose a stepwise model for the evolutionary trajectory of epidemic lineages in Neisseria, and expect that the development of similar dcgMLST schemes will facilitate epidemiological studies of other bacterial pathogens.

A core-genome multilocus sequence typing scheme for the detection of genetically related Streptococcus pyogenes clusters

The recently observed increase in invasive Streptococcus pyogenes infections causes concern in Europe. However, conventional molecular typing methods lack discriminatory power to aid investigations of outbreaks caused by S. pyogenes. Therefore, there is an urgent need for high-resolution molecular typing methods to assess genetic relatedness between S. pyogenes isolates. In the current study, we aimed to develop a novel high-resolution core-genome multilocus sequence typing (cgMLST) scheme for S. pyogenes and compared its discriminatory power to conventional molecular typing methods. The cgMLST scheme was designed with the commercial Ridom SeqSphere+ software package. To define a cluster threshold, the scheme was evaluated using publicly available data from nine defined S. pyogenes outbreaks in the United Kingdom. The cgMLST scheme was then applied to 23 isolates from a suspected S. pyogenes outbreak and 117 S. pyogenes surveillance isolates both from the Netherlands. MLST and emm-typing results were used for comparison to cgMLST results. The allelic differences between isolates from defined outbreaks ranged between 6 and 31 for isolates with the same emm-type, resulting in a proposed cluster threshold of <5 allelic differences out of 1,095 target loci. Seven out of twenty-three (30%) isolates from the suspected outbreak had an allelic difference of <2, thereby identifying a potential cluster that could not be linked to other isolates. The proposed cgMLST scheme shows a higher discriminatory ability when compared to conventional typing methods. The rapid and simple analysis workflow allows for extended detection of clusters of potential outbreak isolates and surveillance and may facilitate the sharing of sequencing results between (inter)national laboratories.

Combining Fusion of Cells with CRISPR-Cas9 Editing for the Cloning of Large DNA Fragments or Complete Bacterial Genomes in Yeast

The genetic engineering of genome fragments larger than 100 kbp is challenging and requires both specific methods and cloning hosts. The yeast Saccharomyces cerevisiae is considered as a host of choice for cloning and engineering whole or partial genomes from viruses, bacteria, and algae. Several methods are now available to perform these manipulations, each with its own limitations. In order to extend the range of yeast cloning strategies, a new approach combining two already described methods, Fusion cloning and CReasPy-Cloning, was developed. The CReasPy-Fusion method allows the simultaneous cloning and engineering of megabase-sized genomes in yeast by the fusion of bacterial cells with yeast spheroplasts carrying the CRISPR-Cas9 system. With this new approach, we demonstrate the feasibility of cloning and editing whole genomes from several Mycoplasma species belonging to different phylogenetic groups. We also show that CReasPy-Fusion allows the capture of large genome fragments with high efficacy, resulting in the successful cloning of selected loci in yeast. We finally identify bacterial nuclease encoding genes as barriers for CReasPy-Fusion by showing that their removal from the donor genome improves the cloning efficacy.

Bacterial genome size and gene functional diversity negatively correlate with taxonomic diversity along a pH gradient

Bacterial gene repertoires reflect adaptive strategies, contribute to ecosystem functioning and are limited by genome size. However, gene functional diversity does not necessarily correlate with taxonomic diversity because average genome size may vary by community. Here, we analyse gene functional diversity (by shotgun metagenomics) and taxonomic diversity (by 16S rRNA gene amplicon sequencing) to investigate soil bacterial communities along a natural pH gradient in 12 tropical, subtropical, and temperate forests. We find that bacterial average genome size and gene functional diversity decrease, whereas taxonomic diversity increases, as soil pH rises from acid to neutral; as a result, bacterial taxonomic and functional diversity are negatively correlated. The gene repertoire of acid-adapted oligotrophs is enriched in functions of signal transduction, cell motility, secretion system, and degradation of complex compounds, while that of neutral pH-adapted copiotrophs is enriched in functions of energy metabolism and membrane transport. Our results indicate that a mismatch between taxonomic and functional diversity can arise when environmental factors (such as pH) select for adaptive strategies that affect genome size distributions.

Toward Minimal Transcription-Translation Machinery

With the advantages of simple genetic composition, low metabolic background, low energy waste, and high genetic stability, genome-reduced strains, as promising functional chassis, have become an intensive direction for constructing potent biosynthesis factories. Herein, an innovative Genome-Reduced strain-based Active Cell-free Easy-to-make-protein (GRACE) system is built as minimal transcription-translation machinery. In this study, two Escherichia coli genome-reduced strains, ΔW3110 and ΔMG1655, with genome reduction of 11.53% and 37.85%, are fused with the cell-free transcription-translation (CFTT) system. The GRACE systems perform better than the corresponding CFTT systems derived from their parental strains in representative valuable applications, such as the expression and solubilization of membrane proteins or protein polymers, biosensing of inorganic or organic molecules based on different principles, and unnatural amino acid embedding. Obviously, the GRACE system has provided a brand-new enabling platform for cell-free transcription-translation basic and applied studies and also would inspire the potential of genome-reduced strains for versatile applications.

VBCG: 20 validated bacterial core genes for phylogenomic analysis with high fidelity and resolution

BACKGROUND

Phylogenomic analysis has become an inseparable part of studies of bacterial diversity and evolution, and many different bacterial core genes have been collated and used for phylogenomic tree reconstruction. However, these genes have been selected based on their presence and single-copy ratio in all bacterial genomes, leaving out the gene's 'phylogenetic fidelity' unexamined.

RESULTS

From 30,522 complete genomes covering 11,262 species, we examined 148 bacterial core genes that have been previously used for phylogenomic analysis. In addition to the gene presence and single-copy rations, we evaluated the gene's phylogenetic fidelity by comparing each gene's phylogeny with its corresponding 16S rRNA gene tree. Out of the 148 bacterial genes, 20 validated bacterial core genes (VBCG) were selected as the core gene set with the highest bacterial phylogenetic fidelity. Compared to the larger gene set, the 20-gene core set resulted in more species having all genes present and fewer species with missing data, thereby enhancing the accuracy of phylogenomic analysis. Using Escherichia coli strains as examples of prominent bacterial foodborne pathogens, we demonstrated that the 20 VBCG produced phylogenies with higher fidelity and resolution at species and strain levels while 16S rRNA gene tree alone could not.

CONCLUSION

The 20 validated core gene set improves the fidelity and speed of phylogenomic analysis. Among other uses, this tool improves our ability to explore the evolution, typing and tracking of bacterial strains, such as human pathogens. We have developed a Python pipeline and a desktop graphic app (available on GitHub) for users to perform phylogenomic analysis with high fidelity and resolution. Video Abstract.

A large-scale comparative genomics study reveals niche-driven and within-sample intra-species functional diversification in Lacticaseibacillus rhamnosus

Lacticaseibacillus rhamnosus (L. rhamnosus) is widely recognized as a probiotic species, and it exists in a variety of environments including host gut and dairy products. This work aimed at conducting a large-scale comparative genomics analysis of 384 L. rhamnosus genomes (257 whole-sequence or metagenomic-assembled genomes from gut-associated isolates [122 and 135 retrieved from the UHGG and NCBI databases, respectively] and 127 genomes from dairy isolates [34 from the NCBI database; 93 isolated from a cheese sample and sequenced here]). Our results showed that L. rhamnosus had a large and open pan-genome (15,253 pan-genes identified from all 384 genomes; 15,028 pan-genes if the 93 cheese-originated isolates were excluded). The core-gene phylogenetic tree constructed from the 384 L. rhamnosus genomes comprised five phylogenetic branches, with a random distribution of dairy and gut-associated isolates/genomes across the tree. No significant difference was identified in the overall profile of metabolism-related genes between dairy and gut-associated genomes; however, notably, the gut-associated strains/isolates contained more genes coding for specific metabolic pathways and carbohydrate-active enzymes, e.g., lacto-N-biosidase (EC 3.2.1.140; GT20) and lacto-N-biose phosphorylase/galacto-N-biose phosphorylase (EC 2.4.1.211; GH112). Further, we found that there was obvious intra-species diversification of the 93 cheese-originated L. rhamnosus isolates, forming three clades (Clades A, B, and C) in the reconstructed core-gene phylogenetic tree. There were numerous single nucleotide variations (over 10,000) across the three clades. Moreover, significant differences were observed in the content of metabolism-related genes across clades (p < 0.05, Adonis test), characterized by the enrichment in glycoside hydrolases in Clade C and the possession of unique metabolic pathways in each clade. These results implicated genomics/functional diversification of L. rhamnosus in a single food matrix and niche-driven adaptive evolution of isolates from dairy and host gut-associated origins. Our study shed insights into the selection of candidate strains for food industry applications.

Two species with a peculiar evolution within the genus Pectobacterium suggest adaptation to a new environmental niche

Historically, research on Soft Rot Pectobacteriacea (SRP) has focused on economically important crops and ornamentals and knowledge of these bacteria outside the plant context remains poorly investigated. Recently, two closely related species Pectobacterium aquaticum and Pectobacterium quasiaquaticum were isolated from water and have not been isolated from any plant yet. To identify the distinctive characteristics of these two species, we performed a comparative genomic analysis of 80 genomes representing 19 Pectobacterium species and performed an evolutionary reconstruction. Both water species underwent a reduction in genome size associated with a high pseudogene content. A high gene loss was predicted at the emergence of both species. Among the 199 gene families missing from both P. aquaticum and P. quasiaquaticum genomes but present in at least 80% of other Pectobacterium genomes, COG analysis identified many genes involved in nutrient transport systems. In addition, many type II secreted proteins were also missing in both species. Phenotypic analysis revealed that both species had reduced pectinolytic activity, a biofilm formation defect, were highly motile and had reduced virulence on several plants. These genomic and phenotypic data suggest that the ecological niche of P. aquaticum and P. quasiaquaticum may differ from that of other Pectobacterium species.

The heterogenous and diverse population of prophages in Mycobacterium genomes

IMPORTANCE

Mycobacterium species include several human pathogens and mycobacteriophages show potential for therapeutic use to control Mycobacterium infections. However, phage infection profiles vary greatly among Mycobacterium abscessus clinical isolates and phage therapies must be personalized for individual patients. Mycobacterium phage susceptibility is likely determined primarily by accessory parts of bacterial genomes, and we have identified the prophage and phage-related genomic regions across sequenced Mycobacterium strains. The prophages are numerous and diverse, especially in M. abscessus genomes, and provide a potentially rich reservoir of new viruses that can be propagated lytically and used to expand the repertoire of therapeutically useful phages.

Evaluation of Fourier transform-infrared spectroscopy (FT-IR) as a control measure for nosocomial outbreak investigations

Whole-genome sequencing (WGS) provides greater resolution than other molecular epidemiology strategies and is emerging as a new gold standard approach for microbial strain typing. The Bruker IR Biotyper is designed as a screening tool to identify bacterial isolates that require WGS to establish accurate relationships, but its performance and utility in nosocomial outbreak investigations have not been thoroughly investigated. Here, we evaluated the IR Biotyper by retrospectively examining isolates tested by WGS during investigations of potential nosocomial transmission events or outbreaks. Ninety-eight clinical isolates from 14 different outbreak investigations were examined: three collections of Acinetobacter baumannii (n = 2, n = 9, n = 5 isolates in each collection), one of Escherichia coli (n = 16), two of Pseudomonas aeruginosa (n = 2 and n = 5), two of Serratia marcescens (n = 9 and n = 7), five of Staphylococcus aureus (n = 8, n = 4, n = 3, n = 3, n = 17), and one of Stenotrophomonas maltophilia (n = 8). Linear regression demonstrated a weak, positive correlation between the number of pairwise genome-wide single-nucleotide polymorphisms (SNPs) and IR Biotyper spectral distance values for Gram-positive (r = 0.43, P ≤ 0.0001), Gram-negative (r = 0.1554, P = 0.0639), and all organisms combined (r = 0.342, P ≤ 0.0001). Overall, the IR Biotyper had a positive predictive value (PPV) of 55.81% for identifying strains that were closely related by genomic identity, but a negative predictive value (NPV) of 86.79% for identifying unrelated isolates. When experimentally adjusted cut-offs were applied to A. baumannii, P. aeruginosa, and E. coli, the PPV was 62% for identifying strains that were closely related and the NPV was 100% for identifying unrelated isolates. Implementation of the IR Biotyper as a screening tool in this cohort would have reduced the number of Gram-negative isolates requiring further WGS analysis by 50% and would reduce the number of S. aureus isolates needing WGS resolution by 48%.

Genomic footprinting uncovers global transcription factor responses to amino acids in Escherichia coli

Our knowledge of transcriptional responses to changes in nutrient availability comes primarily from few well-studied transcription factors (TFs), often lacking an unbiased genome-wide perspective. Leveraging recent advances allowing bacterial genomic footprinting, we comprehensively mapped the genome-wide regulatory responses of Escherichia coli to exogenous leucine, methionine, alanine, and lysine. The global TF Lrp was found to individually sense three amino acids and mount three different target gene responses. Overall, 531 genes had altered RNA polymerase occupancy, and 32 TFs responded directly or indirectly to the presence of amino acids, including regulators of membrane and osmotic pressure homeostasis. About 70% of the detected TF-DNA interactions had not been reported before. We thus identified 682 previously unknown TF-binding locations, for a subset of which the involved TFs were identified by affinity purification. This comprehensive map of amino acid regulation illustrates the incompleteness of the known transcriptional regulation network, even in E. coli.

Classification and prediction of Klebsiella pneumoniae strains with different MLST allelic profiles via SERS spectral analysis

The Gram-negative non-motile Klebsiella pneuomoniae is currently a major cause of hospital-acquired (HA) and community-acquired (CA) infections, leading to great public health concern globally, while rapid identification and accurate tracing of the pathogenic bacterium is essential in facilitating monitoring and controlling of K. pneumoniae outbreak and dissemination. Multi-locus sequence typing (MLST) is a commonly used typing approach with low cost that is able to distinguish bacterial isolates based on the allelic profiles of several housekeeping genes, despite low resolution and labor intensity of the method. Core-genome MLST scheme (cgMLST) is recently proposed to sub-type and monitor outbreaks of bacterial strains with high resolution and reliability, which uses hundreds or thousands of genes conserved in all or most members of the species. However, the method is complex and requires whole genome sequencing of bacterial strains with high costs. Therefore, it is urgently needed to develop novel methods with high resolution and low cost for bacterial typing. Surface enhanced Raman spectroscopy (SERS) is a rapid, sensitive and cheap method for bacterial identification. Previous studies confirmed that classification and prediction of bacterial strains via SERS spectral analysis correlated well with MLST typing results. However, there is currently no similar comparative analysis in K. pneumoniae strains. In this pilot study, 16 K. pneumoniae strains with different sequencing typings (STs) were selected and a phylogenetic tree was constructed based on core genome analysis. SERS spectra (N = 45/each strain) were generated for all the K. pneumoniae strains, which were then comparatively classified and predicted via six representative machine learning (ML) algorithms. According to the results, SERS technique coupled with the ML algorithm support vector machine (SVM) could achieve the highest accuracy (5-Fold Cross Validation = 100%) in terms of differentiating and predicting all the K. pneumoniae strains that were consistent to corresponding MLSTs. In sum, we show in this pilot study that the SERS-SVM based method is able to accurately predict K. pneumoniae MLST types, which has the application potential in clinical settings for tracing dissemination and controlling outbreak of K. pneumoniae in hospitals and communities with low costs and high rapidity.

A systematic approach to classify and characterize genomic islands driven by conjugative mobility using protein signatures

Genomic islands (GIs) play a crucial role in the spread of antibiotic resistance, virulence factors and antiviral defense systems in a broad range of bacterial species. However, the characterization and classification of GIs are challenging due to their relatively small size and considerable genetic diversity. Predicting their intercellular mobility is of utmost importance in the context of the emerging crisis of multidrug resistance. Here, we propose a large-scale classification method to categorize GIs according to their mobility profile and, subsequently, analyze their gene cargo. We based our classification decision scheme on a collection of mobility protein motif definitions available in publicly accessible databases. Our results show that the size distribution of GI classes correlates with their respective structure and complexity. Self-transmissible GIs are usually the largest, except in Bacillota and Actinomycetota, accumulate antibiotic and phage resistance genes, and favour the use of a tyrosine recombinase to insert into a host's replicon. Non-mobilizable GIs tend to use a DDE transposase instead. Finally, although tRNA genes are more frequently targeted as insertion sites by GIs encoding a tyrosine recombinase, most GIs insert in a protein-encoding gene. This study is a stepping stone toward a better characterization of mobile GIs in bacterial genomes and their mechanism of mobility.

Genomic typing and virulence gene profile analysis of Salmonella Derby from different sources

Salmonella enterica serovar Derby (S. Derby) is one of the most common Salmonella serovars which can infect poultry, swine, and humans. With the reduction of the sequencing cost and the improvement of sequencing technology, whole genome sequencing (WGS) has become an important method for bacterial determination, molecular investigation, and pathogenic tracing analysis. In this study, we investigated S. Derby isolates from different sources in China using in-silico multilocus sequence typing (MLST), core genome MLST (cgMLST) and whole genome MLST (wgMLST) analysis based on WGS. The results showed that 21 S. Derby strains were divided into 3 STs using MLST analysis, including ST40 (n = 19, accounting for 90.48%), ST71 (n = 1, accounting for 4.76%) and ST8016 (n = 1, accounting for 4.76%). cgMLST and wgMLST analysis categorized the tested strains into 13 cgSTs and 21 wgSTs, respectively. The minimum spanning trees of cgMLST and wgMLST both divided these strains into 3 clusters and 4 singletons. In addition, virulence gene profiles of S. Derby isolates were also analyzed, and a total of 174 virulence genes belonged to 8 categories were identified. In summary, we studied genomic typing, phylogenetic relationship and virulence gene profiles of S. Derby strains from different sources in China. These findings were beneficial for the epidemiology and pathogenesis of Salmonella.

Predicting variable gene content in Escherichia coli using conserved genes

Having the ability to predict the protein-encoding gene content of an incomplete genome or metagenome-assembled genome is important for a variety of bioinformatic tasks. In this study, as a proof of concept, we built machine learning classifiers for predicting variable gene content in Escherichia coli genomes using only the nucleotide k-mers from a set of 100 conserved genes as features. Protein families were used to define orthologs, and a single classifier was built for predicting the presence or absence of each protein family occurring in 10%-90% of all E. coli genomes. The resulting set of 3,259 extreme gradient boosting classifiers had a per-genome average macro F1 score of 0.944 [0.943-0.945, 95% CI]. We show that the F1 scores are stable across multi-locus sequence types and that the trend can be recapitulated by sampling a smaller number of core genes or diverse input genomes. Surprisingly, the presence or absence of poorly annotated proteins, including "hypothetical proteins" was accurately predicted (F1 = 0.902 [0.898-0.906, 95% CI]). Models for proteins with horizontal gene transfer-related functions had slightly lower F1 scores but were still accurate (F1s = 0.895, 0.872, 0.824, and 0.841 for transposon, phage, plasmid, and antimicrobial resistance-related functions, respectively). Finally, using a holdout set of 419 diverse E. coli genomes that were isolated from freshwater environmental sources, we observed an average per-genome F1 score of 0.880 [0.876-0.883, 95% CI], demonstrating the extensibility of the models. Overall, this study provides a framework for predicting variable gene content using a limited amount of input sequence data. IMPORTANCE Having the ability to predict the protein-encoding gene content of a genome is important for assessing genome quality, binning genomes from shotgun metagenomic assemblies, and assessing risk due to the presence of antimicrobial resistance and other virulence genes. In this study, we built a set of binary classifiers for predicting the presence or absence of variable genes occurring in 10%-90% of all publicly available E. coli genomes. Overall, the results show that a large portion of the E. coli variable gene content can be predicted with high accuracy, including genes with functions relating to horizontal gene transfer. This study offers a strategy for predicting gene content using limited input sequence data.

Roseovarius pelagicus sp. nov., a facultatively anaerobic bacterium with potential for degrading polypropylene, isolated from Arctic seawater

A Gram-stain-negative, facultatively anaerobic, non-motile, rod-shaped bacterial strain, designated HL-MP18T, was isolated from Arctic seawater after a prolonged incubation employing polypropylene as the sole carbon source. Phylogenetic analyses of the 16S rRNA gene sequence revealed that strain HL-MP18T was affiliated to the genus Roseovarius with close relatives Roseovarius carneus LXJ103T (96.8 %) and Roseovarius litorisediminis KCTC 32327T (96.5 %). The complete genome sequence of strain HL-MP18T comprised a circular chromosome of 3.86 Mbp and two circular plasmids of 0.17 and 0.24 Mbp. Genomic comparisons based on average nucleotide identity and digital DNA-DNA hybridization showed that strain HL-MP18T was consistently discriminated from its closely related taxa in the genus Roseovarius. Strain HL-MP18T showed optimal growth at 25 °C, pH 7.0 and 2.5 % (w/v) sea salts. The major cellular fatty acids were C18 : 1  ω6c and/or C18 : 1  ω7c (49.6 %), C19 : 0 cyclo ω8c (13.5 %), and C16 : 0 (12.8 %). The major respiratory quinone was ubiquinone-10. The polar lipids consisted of phosphatidylcholine, phosphatidylglycerol, an unidentified aminolipid and three unidentified lipids. The genomic DNA G+C content of the strain was 59.2 mol%. The phylogenetic, genomic, phenotypic and chemotaxonomic results indicate that strain HL-MP18T is distinguishable from the recognized species of the genus Roseovarius. Therefore, we propose that strain HL-MP18T represents a novel species belonging to the genus Roseovarius, for which the name Roseovarius pelagicus sp. nov. is proposed. The type strain is HL-MP18T (=KCCM 90405T=JCM 35639T).

Characterization of ferredoxins from the thermophilic, acetogenic bacterium Thermoanaerobacter kivui

A major electron carrier involved in energy and carbon metabolism in the acetogenic model organism Thermoanaerobacter kivui is ferredoxin, an iron-sulfur-containing, electron-transferring protein. Here, we show that the genome of T. kivui encodes four putative ferredoxin-like proteins (TKV_c09620, TKV_c16450, TKV_c10420 and TKV_c19530). All four genes were cloned, a His-tag encoding sequence was added and the proteins were produced from a plasmid in T. kivui. The purified proteins had an absorption peak at 430 nm typical for ferredoxins. The determined iron-sulfur content is consistent with the presence of two predicted [4Fe4S] clusters in TKV_c09620 and TKV_c19530 or one predicted [4Fe4S] cluster in TKV_c16450 and TKV_c10420 respectively. The reduction potential (Em ) for TKV_c09620, TKV_c16450, TKV_c10420 and TKV_c19530 was determined to be -386 ± 4 mV, -386 ± 2 mV, -559 ± 10 mV and -557 ± 3 mV, respectively. TKV_c09620 and TKV_c16450 served as electron carriers for different oxidoreductases from T. kivui. Deletion of the ferredoxin genes led to only a slight reduction of growth on pyruvate or autotrophically on H2  + CO2 . Transcriptional analysis revealed that TKV_c09620 was upregulated in a ΔTKV_c16450 mutant and vice versa TKV_c16450 in a ΔTKV_c09620 mutant, indicating that TKV_c09620 and TKV_c16450 can replace each other. In sum, our data are consistent with the hypothesis that TKV_c09620 and TKV_c16450 are ferredoxins involved in autotrophic and heterotrophic metabolism of T. kivui.

Identification of two distinct phylogenomic lineages and model strains for the understudied cystic fibrosis lung pathogen Burkholderia multivorans

Burkholderia multivorans is the dominant Burkholderia pathogen recovered from lung infection in people with cystic fibrosis. However, as an understudied pathogen there are knowledge gaps in relation to its population biology, phenotypic traits and useful model strains. A phylogenomic study of B. multivorans was undertaken using a total of 283 genomes, of which 73 were sequenced and 49 phenotypically characterized as part of this study. Average nucleotide identity analysis (ANI) and phylogenetic alignment of core genes demonstrated that the B. multivorans population separated into two distinct evolutionary clades, defined as lineage 1 (n=58 genomes) and lineage 2 (n=221 genomes). To examine the population biology of B. multivorans, a representative subgroup of 77 B. multivorans genomes (28 from the reference databases and the 49 novel short-read genome sequences) were selected based on multilocus sequence typing (MLST), isolation source and phylogenetic placement criteria. Comparative genomics was used to identify B. multivorans lineage-specific genes - ghrB_1 in lineage 1 and glnM_2 in lineage 2 - and diagnostic PCRs targeting them were successfully developed. Phenotypic analysis of 49 representative B. multivorans strains showed considerable inter-strain variance, but the majority of the isolates tested were motile and capable of biofilm formation. A striking absence of B. multivorans protease activity in vitro was observed, but no lineage-specific phenotypic differences were demonstrated. Using phylogenomic and phenotypic criteria, three model B. multivorans CF strains were identified, BCC0084 (lineage 1), BCC1272 (lineage 2a) and BCC0033 lineage 2b, and their complete genome sequences determined. B. multivorans CF strains BCC0033 and BCC0084, and the environmental reference strain, ATCC 17616, were all capable of short-term survival within a murine lung infection model. By mapping the population biology, identifying lineage-specific PCRs and model strains, we provide much needed baseline resources for future studies of B. multivorans.

Evolution of a minimal cell

Possessing only essential genes, a minimal cell can reveal mechanisms and processes that are critical for the persistence and stability of life1,2. Here we report on how an engineered minimal cell3,4 contends with the forces of evolution compared with the Mycoplasma mycoides non-minimal cell from which it was synthetically derived. Mutation rates were the highest among all reported bacteria, but were not affected by genome minimization. Genome streamlining was costly, leading to a decrease in fitness of greater than 50%, but this deficit was regained during 2,000 generations of evolution. Despite selection acting on distinct genetic targets, increases in the maximum growth rate of the synthetic cells were comparable. Moreover, when performance was assessed by relative fitness, the minimal cell evolved 39% faster than the non-minimal cell. The only apparent constraint involved the evolution of cell size. The size of the non-minimal cell increased by 80%, whereas the minimal cell remained the same. This pattern reflected epistatic effects of mutations in ftsZ, which encodes a tubulin-homologue protein that regulates cell division and morphology5,6. Our findings demonstrate that natural selection can rapidly increase the fitness of one of the simplest autonomously growing organisms. Understanding how species with small genomes overcome evolutionary challenges provides critical insights into the persistence of host-associated endosymbionts, the stability of streamlined chassis for biotechnology and the targeted refinement of synthetically engineered cells2,7-9.

High-throughput functional genomics: A (myco)bacterial perspective

Advances in sequencing technologies have enabled unprecedented insights into bacterial genome composition and dynamics. However, the disconnect between the rapid acquisition of genomic data and the (much slower) confirmation of inferred genetic function threatens to widen unless techniques for fast, high-throughput functional validation can be applied at scale. This applies equally to Mycobacterium tuberculosis, the leading infectious cause of death globally and a pathogen whose genome, despite being among the first to be sequenced two decades ago, still contains many genes of unknown function. Here, we summarize the evolution of bacterial high-throughput functional genomics, focusing primarily on transposon (Tn)-based mutagenesis and the construction of arrayed mutant libraries in diverse bacterial systems. We also consider the contributions of CRISPR interference as a transformative technique for probing bacterial gene function at scale. Throughout, we situate our analysis within the context of functional genomics of mycobacteria, focusing specifically on the potential to yield insights into M. tuberculosis pathogenicity and vulnerabilities for new drug and regimen development. Finally, we offer suggestions for future approaches that might be usefully applied in elucidating the complex cellular biology of this major human pathogen.

Repurposing the atypical type I-G CRISPR system for bacterial genome engineering

The CRISPR-Cas system functions as a prokaryotic immune system and is highly diverse, with six major types and numerous sub-types. The most abundant are type I CRISPR systems, which utilize a multi-subunit effector, Cascade, and a CRISPR RNA (crRNA) to detect invading DNA species. Detection leads to DNA loading of the Cas3 helicase-nuclease, leading to long-range deletions in the targeted DNA, thus providing immunity against mobile genetic elements (MGE). Here, we focus on the type I-G system, a streamlined, 4-subunit complex with an atypical Cas3 enzyme. We demonstrate that Cas3 helicase activity is not essential for immunity against MGE in vivo and explore applications of the Thioalkalivibrio sulfidiphilus Cascade effector for genome engineering in Escherichia coli. Long-range, bidirectional deletions were observed when the lacZ gene was targeted. Deactivation of the Cas3 helicase activity dramatically altered the types of deletions observed, with small deletions flanked by direct repeats that are suggestive of microhomology mediated end joining. When donor DNA templates were present, both the wild-type and helicase-deficient systems promoted homology-directed repair (HDR), with the latter system providing improvements in editing efficiency, suggesting that a single nick in the target site may promote HDR in E. coli using the type I-G system. These findings open the way for further application of the type I-G CRISPR systems in genome engineering.

Genomic sequencing of fourteen bacillus thuringiensis isolates: insights into geographic variation and phylogenetic implications

OBJECTIVE

This work was performed in support of a separate study investigating the activity of pesticidal proteins produced by Bacillus thuringiensis against the Asian citrus psyllid, Diaphorina citri. The fourteen Bacillus isolates chosen were selected from a large, geographically diverse collection that was characterized only by biochemical phenotype and morphology of the parasporal crystal, hence, for each isolate it was desired to determine the specific pesticidal proteins produced, assign each to a Bacillus cereus multilocus sequence type (ST), and predict their placement within the classical Bt serotyping system. In addition, phylogenetic distances between the isolates and Bacillus thuringiensis serovar type strains were determined by calculating digital DNA-DNA hybridization (dDDH) values among the isolates.

RESULTS

Based on the assembled sequence data, the isolates were found to be likely representatives of the Bt serovars kurstaki (ST 8), pakistani (ST 550), toumanoffi (ST 240), israelensis (ST 16), thuringiensis (ST 10), entomocidus (ST 239), and finitimus (ST 171). In cases where multiple isolates occurred within a predicted serovar, pesticidal protein profiles were found to be identical, despite the geographic diversity of the isolates. As expected, the dDDH values calculated for pairwise comparisons of the isolates and their apparent corresponding Bt serovar type strains were quite high (> 98%), however dDDH comparisons of the isolates with other serovar type strains were often surprisingly low (< 70%) and suggest unrecognized taxa within Bt and the Bacillus cereus sensu lato.

Core genome multilocus sequence typing scheme for Bacillus cereus group bacteria

Core genome multilocus sequence typing (cgMLST) employs a strategy where the set of orthologous genes common to all members of a group of organisms are used for phylogenetic analysis of the group members. The Bacillus cereus group consists of species with pathogenicity towards insect species as well as warm-blooded animals including humans. While B. cereus is an opportunistic pathogen linked to a range of human disease conditions, including emesis and diarrhoea, Bacillus thuringiensis is an entomopathogenic species with toxicity toward insect larvae, and therefore used as a biological pesticide worldwide. Bacillus anthracis is a classical obligate pathogen causing anthrax, an acute lethal condition in herbivores as well as humans, and which is endemic in many parts of the world. The group also includes a range of additional species, and B. cereus group bacteria have been subject to analysis with a wide variety of phylogenetic typing systems. Here we present, based on analyses of 173 complete genomes from B. cereus group species available in public databases, the identification of a set of 1568 core genes which were used to create a core genome multilocus typing scheme for the group which is implemented in the PubMLST system as an open online database freely available to the community. The new cgMLST system provides unprecedented resolution over existing phylogenetic analysis schemes covering the B. cereus group.

Continuous synthesis of E. coli genome sections and Mb-scale human DNA assembly

Whole-genome synthesis provides a powerful approach for understanding and expanding organism function1-3. To build large genomes rapidly, scalably and in parallel, we need (1) methods for assembling megabases of DNA from shorter precursors and (2) strategies for rapidly and scalably replacing the genomic DNA of organisms with synthetic DNA. Here we develop bacterial artificial chromosome (BAC) stepwise insertion synthesis (BASIS)-a method for megabase-scale assembly of DNA in Escherichia coli episomes. We used BASIS to assemble 1.1 Mb of human DNA containing numerous exons, introns, repetitive sequences, G-quadruplexes, and long and short interspersed nuclear elements (LINEs and SINEs). BASIS provides a powerful platform for building synthetic genomes for diverse organisms. We also developed continuous genome synthesis (CGS)-a method for continuously replacing sequential 100 kb stretches of the E. coli genome with synthetic DNA; CGS minimizes crossovers1,4 between the synthetic DNA and the genome such that the output for each 100 kb replacement provides, without sequencing, the input for the next 100 kb replacement. Using CGS, we synthesized a 0.5 Mb section of the E. coli genome-a key intermediate in its total synthesis1-from five episomes in 10  days. By parallelizing CGS and combining it with rapid oligonucleotide synthesis and episome assembly5,6, along with rapid methods for compiling a single genome from strains bearing distinct synthetic genome sections1,7,8, we anticipate that it will be possible to synthesize entire E. coli genomes from functional designs in less than 2 months.