Exclusivity offers a sound yet practical species criterion for bacteria despite abundant gene flow

Accurately clustering biological sequences in linear time by relatedness sorting

Clustering biological sequences into similar groups is an increasingly important task as the number of available sequences continues to grow exponentially. Search-based approaches to clustering scale super-linearly with the number of input sequences, making it impractical to cluster very large sets of sequences. Approaches to clustering sequences in linear time currently lack the accuracy of super-linear approaches. Here, I set out to develop and characterize a strategy for clustering with linear time complexity that retains the accuracy of less scalable approaches. The resulting algorithm, named Clusterize, sorts sequences by relatedness to linearize the clustering problem. Clusterize produces clusters with accuracy rivaling popular programs (CD-HIT, MMseqs2, and UCLUST) but exhibits linear asymptotic scalability. Clusterize generates higher accuracy and oftentimes much larger clusters than Linclust, a fast linear time clustering algorithm. I demonstrate the utility of Clusterize by accurately solving different clustering problems involving millions of nucleotide or protein sequences.

Clostridium strain FAM25158, a unique endospore-forming bacterium related to Clostridium tyrobutyricum and isolated from Emmental cheese shows low tolerance to salt

The genus Clostridium is a large and diverse group of species that can cause food spoilage, including late blowing defect (LBD) in cheese. In this study, we investigated the taxonomic status of strain FAM25158 isolated from Emmental cheese with LBD using a polyphasic taxonomic and comparative genomic approach. A 16S rRNA gene sequence phylogeny suggested affiliation to the Clostridium sensu stricto cluster, with Clostridium tyrobutyricum DSM 2637T being the closest related type strain (99.16% sequence similarity). Average Nucleotide Identity (ANI) analysis revealed that strain FAM25158 is at the species threshold with C. tyrobutyricum, with ANI values ranging from 94.70 to 95.26%, while the digital DNA-DNA hybridization values were below the recommended threshold, suggesting that FAM25158 is significantly different from C. tyrobutyricum at the genomic level. Moreover, comparative genomic analysis between FAM25158 and its four closest C. tyrobutyricum relatives revealed a diversity of metabolic pathways, with FAM25158 differing from other C. tyrobutyricum strains by the presence of genes such as scrA, srcB, and scrK, responsible for sucrose utilization, and the absence of many important functional genes associated with cold and osmolality adaptation, which was further supported by phenotypic analyses. Surprisingly, strain FAM25158 exhibited unique physiologic traits, such as an optimal growth temperature of 30°C, in contrast to its closest relatives, C. tyrobutyricum species with an optimal growth temperature of 37°C. Additionally, the growth of FAM25158 was inhibited at NaCl concentrations higher than 0.5%, a remarkable observation considering its origin from cheese. While the results of this study provide novel information on the genetic content of strain FAM25158, the relationship between its genetic content and the observed phenotype remains a topic requiring further investigation.

A universal and constant rate of gene content change traces pangenome flux to LUCA

Prokaryotic genomes constantly undergo gene flux via lateral gene transfer, generating a pangenome structure consisting of a conserved core genome surrounded by a more variable accessory genome shell. Over time, flux generates change in genome content. Here, we measure and compare the rate of genome flux for 5655 prokaryotic genomes as a function of amino acid sequence divergence in 36 universally distributed proteins of the informational core (IC). We find a clock of gene content change. The long-term average rate of gene content flux is remarkably constant across all higher prokaryotic taxa sampled, whereby the size of the accessory genome-the proportion of the genome harboring gene content difference for genome pairs-varies across taxa. The proportion of species-level accessory genes per genome, varies from 0% (Chlamydia) to 30%-33% (Alphaproteobacteria, Gammaproteobacteria, and Clostridia). A clock-like rate of gene content change across all prokaryotic taxa sampled suggest that pangenome structure is a general feature of prokaryotic genomes and that it has been in existence since the divergence of bacteria and archaea.

Phenotypic and Genomic Characterization of Streptomyces pakalii sp. nov., a Novel Species with Anti-Biofilm and Anti-Quorum Sensing Activity in ESKAPE Bacteria

The increasing number of infections caused by antimicrobial multi-resistant microorganisms has led to the search for new microorganisms capable of producing novel antibiotics. This work proposes Streptomyces pakalii sp. nov. as a new member of the Streptomycetaceae family. The strain ENCB-J15 was isolated from the jungle soil in Palenque National Park, Chiapas, Mexico. The strain formed pale brown, dry, tough, and buried colonies in the agar with no diffusible pigment in GAE (glucose-asparagine-yeast extract) medium. Scanning electron micrographs showed typical mycelium with long chains of smooth and oval-shaped spores (3-10 m). The strain grew in all of the International Streptomyces Project (ISP)'s media at 28-37 °C with a pH of 6-9 and 0-10% NaCl. S. pakalii ENCB-J15 assimilated diverse carbon as well as organic and inorganic nitrogen sources. The strain also exhibited significant inhibitory activity against the prodigiosin synthesis of Serratia marcescens and the inhibition of the formation and destruction of biofilms of ESKAPE strains of Acinetobacter baumannii and Klebsiella pneumoniae. The draft genome sequencing of ENCB-J15 revealed a 7.6 Mb genome with a high G + C content (71.6%), 6833 total genes, and 6746 genes encoding putative proteins. A total of 26 accessory clusters of proteins associated with carbon sources and amino acid catabolism, DNA modification, and the antibiotic biosynthetic process were annotated. The 16S rRNA gene phylogeny, core-proteome phylogenomic tree, and virtual genome fingerprints support that S. pakalii ENCB-J15 is a new species related to Streptomyces badius and Streptomyces globisporus. Similarly, its average nucleotide identity (ANI) (96.4%), average amino acid identity (AAI) (96.06%), and virtual DNA-DNA hybridization (67.3%) provide evidence to recognize it as a new species. Comparative genomics revealed that S. pakalli and its closest related species maintain a well-conserved genomic synteny. This work proposes Streptomyces pakalii sp. nov. as a novel species that expresses anti-biofilm and anti-quorum sensing activities.

Exploring the universal healthy human gut microbiota around the World

The human gut holds a special place in the study of different microbial environments due to growing evidence that the gut microbiota is related to host health. However, despite extensive research, there is still a lack of knowledge about the core taxa forming the gut microbiota and, moreover, available information is biased towards western microbiomes in both genome databases and most core taxa studies. To tackle these limitations, we tested a database enrichment strategy and analyzed public datasets of whole-genome shotgun data, generated from 545 fecal samples, comprising three gradients of westernization. The NT database was selected as a baseline of biological diversity, subsequently being combined with various studies of interest related to the human microbiota. This enrichment strategy made it possible to improve classification capacity, compared to the original unenriched database, regarding the various lifestyles and populations studied. The effects of incomplete-taxonomy metagenome-assembled genomes on genome database enrichment were also examined, revealing that, while they are helpful, they should be used with caution depending on the taxonomic level of interest. Moreover, in terms of high prevalence, the core analysis revealed a conserved set of bacterial taxa in the healthy human gut microbiota worldwide, despite apparent lifestyle differences. Such taxa show a set of traits, metabolic roles, and ancestral status, making them suitable candidates for a hypothetical phylogenetic core of mutualistic microorganisms co-evolving with the human species.

Benchmarking DNA Extraction Methods for Phylogenomic Analysis of Sub-Antarctic Rhodococcus and Williamsia Species

Bacteria containing mycolic acids in their cell envelope are often recalcitrant to cell lysis, so extracting DNA of sufficient quality for third-generation sequencing and high-fidelity genome assembly requires optimization, even when using commercial kits with protocols for hard-to-lyse bacteria. We benchmarked three spin-column-based kits against a classical DNA extraction method employing lysozyme, proteinase K and SDS for six lysozyme-resistant, sub-Antarctic strains of Corynebaceriales. Prior cultivation in broths containing glycine at highly growth-inhibitory concentrations (4.0–4.5%) improved cell lysis using both classical and kit methods. The classical method produced DNA with average fragment sizes of 27–59 Kbp and tight fragment size ranges, meeting quality standards for genome sequencing, assembly and phylogenomic analyses. By 16S rRNA gene sequencing, we classified two strains as Williamsia and four strains as Rhodococcus species. Pairwise comparison of average nucleotide identity (ANI) and alignment fraction (AF), plus genome clustering analysis, confirmed Rhodococcus sp. 1163 and 1168 and Williamsia sp. 1135 and 1138 as novel species. Phylogenetic, lipidomic and biochemical analyses classified psychrotrophic strains 1139 and 1159 as R. qingshengii and R. erythropolis, respectively, using ANI similarity of >98% and AF >60% for species delineation. On this basis, some members of the R. erythropolis genome cluster groups, including strains currently named as R. enclensis, R. baikonurensis, R. opacus and R. rhodochrous, would be reclassified either as R. erythropolis or R. qingshengii.

A spectrum of verticality across genes

Lateral gene transfer (LGT) has impacted prokaryotic genome evolution, yet the extent to which LGT compromises vertical evolution across individual genes and individual phyla is unknown, as are the factors that govern LGT frequency across genes. Estimating LGT frequency from tree comparisons is problematic when thousands of genomes are compared, because LGT becomes difficult to distinguish from phylogenetic artefacts. Here we report quantitative estimates for verticality across all genes and genomes, leveraging a well-known property of phylogenetic inference: phylogeny works best at the tips of trees. From terminal (tip) phylum level relationships, we calculate the verticality for 19,050,992 genes from 101,422 clusters in 5,655 prokaryotic genomes and rank them by their verticality. Among functional classes, translation, followed by nucleotide and cofactor biosynthesis, and DNA replication and repair are the most vertical. The most vertically evolving lineages are those rich in ecological specialists such as Acidithiobacilli, Chlamydiae, Chlorobi and Methanococcales. Lineages most affected by LGT are the α-, β-, γ-, and δ- classes of Proteobacteria and the Firmicutes. The 2,587 eukaryotic clusters in our sample having prokaryotic homologues fail to reject eukaryotic monophyly using the likelihood ratio test. The low verticality of α-proteobacterial and cyanobacterial genomes requires only three partners-an archaeal host, a mitochondrial symbiont, and a plastid ancestor-each with mosaic chromosomes, to directly account for the prokaryotic origin of eukaryotic genes. In terms of phylogeny, the 100 most vertically evolving prokaryotic genes are neither representative nor predictive for the remaining 97% of an average genome. In search of factors that govern LGT frequency, we find a simple but natural principle: Verticality correlates strongly with gene distribution density, LGT being least likely for intruding genes that must replace a preexisting homologue in recipient chromosomes. LGT is most likely for novel genetic material, intruding genes that encounter no competing copy.

Accurate reconstruction of bacterial pan- and core genomes with PEPPAN

Bacterial genomes can contain traces of a complex evolutionary history, including extensive homologous recombination, gene loss, gene duplications, and horizontal gene transfer. To reconstruct the phylogenetic and population history of a set of multiple bacteria, it is necessary to examine their pangenome, the composite of all the genes in the set. Here we introduce PEPPAN, a novel pipeline that can reliably construct pangenomes from thousands of genetically diverse bacterial genomes that represent the diversity of an entire genus. PEPPAN outperforms existing pangenome methods by providing consistent gene and pseudogene annotations extended by similarity-based gene predictions, and identifying and excluding paralogs by combining tree- and synteny-based approaches. The PEPPAN package additionally includes PEPPAN_parser, which implements additional downstream analyses, including the calculation of trees based on accessory gene content or allelic differences between core genes. To test the accuracy of PEPPAN, we implemented SimPan, a novel pipeline for simulating the evolution of bacterial pangenomes. We compared the accuracy and speed of PEPPAN with four state-of-the-art pangenome pipelines using both empirical and simulated data sets. PEPPAN was more accurate and more specific than any of the other pipelines and was almost as fast as any of them. As a case study, we used PEPPAN to construct a pangenome of approximately 40,000 genes from 3052 representative genomes spanning at least 80 species of Streptococcus The resulting gene and allelic trees provide an unprecedented overview of the genomic diversity of the entire Streptococcus genus.

A Genome-Based Species Taxonomy of the Lactobacillus Genus Complex

The Lactobacillus genus complex is a group of bacteria that constitutes an important source of strains with medical and food applications. The number of bacterial whole-genome sequences available for this taxon has been increasing rapidly in recent years. Despite this wealth of information, the species within this group are still largely defined by older techniques. Here, we constructed a completely new species-level taxonomy for the Lactobacillus genus complex based on ∼2,500 whole-genome sequences. As a result of this effort, we found that many genomes are not classified to their correct species, and we were able to correct these. In addition, we found that some published species are abnormally large, while others are too small. Finally, we discovered at least eight completely novel species that have not been published before. Our work will help the field to evolve toward a more meaningful and complete taxonomy, based on whole-genome sequences.

There are more than 200 published species within the Lactobacillus genus complex (LGC), the majority of which have sequenced type strain genomes available. Although genome-based species delimitation cutoffs are accepted as the gold standard by the community, these are seldom actually checked for new or already published species. In addition, the availability of genome data is revealing inconsistencies in the species-level classification of many strains. We constructed a de novo species taxonomy for the LGC based on 2,459 publicly available genomes, using a 94% core nucleotide identity cutoff. We reconciled these de novo species with published species and subspecies names by (i) identifying genomes of type strains and (ii) comparing 16S rRNA genes of the genomes with 16S rRNA genes of type strains. We found that genomes within the LGC could be divided into 239 de novo species that were discontinuous and exclusive. Comparison of these de novo species to published species led to the identification of nine sets of published species that can be merged and one species that can be split. Further, we found at least eight de novo species that constitute new, unpublished species. Finally, we reclassified 74 genomes on the species level and identified for the first time the species of 98 genomes. Overall, the current state of LGC species taxonomy is largely consistent with genome-based species delimitation cutoffs. There are, however, exceptions that should be resolved to evolve toward a taxonomy where species share a consistent diversity in terms of sequence divergence.

IMPORTANCE The Lactobacillus genus complex is a group of bacteria that constitutes an important source of strains with medical and food applications. The number of bacterial whole-genome sequences available for this taxon has been increasing rapidly in recent years. Despite this wealth of information, the species within this group are still largely defined by older techniques. Here, we constructed a completely new species-level taxonomy for the Lactobacillus genus complex based on ∼2,500 whole-genome sequences. As a result of this effort, we found that many genomes are not classified to their correct species, and we were able to correct these. In addition, we found that some published species are abnormally large, while others are too small. Finally, we discovered at least eight completely novel species that have not been published before. Our work will help the field to evolve toward a more meaningful and complete taxonomy, based on whole-genome sequences.