RiboDB Database: A Comprehensive Resource for Prokaryotic Systematics

Multi-proteins similarity-based sampling to select representative genomes from large databases

BACKGROUND

Genome sequence databases are growing exponentially, but with high redundancy and uneven data quality. For these reasons, selecting representative subsets of genomes is an essential step for almost all studies. However, most current sampling approaches are biased and unable to process large datasets in a reasonable time.

METHODS

Here we present MPS-Sampling (Multiple-Protein Similarity-based Sampling), a fast, scalable, and efficient method for selecting reliable and representative samples of genomes from very large datasets. Using families of homologous proteins as input, MPS-Sampling delineates homogeneous groups of genomes through two successive clustering steps. Representative genomes are then selected within these groups according to predefined or user-defined priority criteria.

RESULTS

MPS-Sampling was applied to a dataset of 48 ribosomal protein families from 178,203 bacterial genomes to generate representative genome sets of various size, corresponding to a sampling of 32.17% down to 0.3% of the complete dataset. An in-depth analysis shows that the selected genomes are both taxonomically and phylogenetically representative of the complete dataset, demonstrating the relevance of the approach.

CONCLUSION

MPS-Sampling provides an efficient, fast and scalable way to sample large collections of genomes in an acceptable computational time. MPS-Sampling does not rely on taxonomic information and does not require the inference of phylogenetic trees, thus avoiding the biases inherent in these approaches. As such, MPS-Sampling meets the needs of a growing number of users.

Lineage-dependent partitioning of activities in chemoclines defines Woesearchaeota ecotypes in an extreme aquatic ecosystem

BACKGROUND

DPANN archaea, including Woesearchaeota, encompass a large fraction of the archaeal diversity, yet their genomic diversity, lifestyle, and role in natural microbiomes remain elusive. With an archaeal assemblage naturally enriched in Woesearchaeota and steep vertical geochemical gradients, Lake Dziani Dzaha (Mayotte) provides an ideal model to decipher their in-situ activity and ecology.

RESULTS

Using genome-resolved metagenomics and phylogenomics, we identified highly diversified Woesearchaeota populations and defined novel halophilic clades. Depth distribution of these populations in the water column showed an unusual double peak of abundance, located at two distinct chemoclines that are hotspots of microbial diversity in the water column. Genome-centric metatranscriptomics confirmed this vertical distribution and revealed a fermentative activity, with acetate and lactate as end products, and active cell-to-cell processes, supporting strong interactions with other community members at chemoclines. Our results also revealed distinct Woesearchaeota ecotypes, with different transcriptional patterns, contrasted lifestyles, and ecological strategies, depending on environmental/host conditions.

CONCLUSIONS

This work provides novel insights into Woesearchaeota in situ activity and metabolism, revealing invariant, bimodal, and adaptative lifestyles among halophilic Woesearchaeota. This challenges our precepts of an invariable host-dependent metabolism for all the members of this taxa and revises our understanding of their contributions to ecosystem functioning and microbiome assemblage. Video Abstract.

Phosphorothioate DNA modification by BREX Type 4 systems in the human gut microbiome

Among dozens of microbial DNA modifications regulating gene expression and host defense, phosphorothioation (PT) is the only known backbone modification, with sulfur inserted at a non-bridging oxygen by dnd and ssp gene families. Here we explored the distribution of PT genes in 13,663 human gut microbiome genomes, finding that 6.3% possessed dnd or ssp genes predominantly in Bacillota, Bacteroidota, and Pseudomonadota. This analysis uncovered several putative new PT synthesis systems, including Type 4 Bacteriophage Exclusion (BREX) brx genes, which were genetically validated in Bacteroides salyersiae. Mass spectrometric analysis of DNA from 226 gut microbiome isolates possessing dnd, ssp, and brx genes revealed 8 PT dinucleotide settings confirmed in 6 consensus sequences by PT-specific DNA sequencing. Genomic analysis showed PT enrichment in rRNA genes and depletion at gene boundaries. These results illustrate the power of the microbiome for discovering prokaryotic epigenetics and the widespread distribution of oxidation-sensitive PTs in gut microbes.

Biochemical and Structural Characterization of Chi-Class Glutathione Transferases: A Snapshot on the Glutathione Transferase Encoded by sll0067 Gene in the Cyanobacterium Synechocystis sp. Strain PCC 6803

Glutathione transferases (GSTs) constitute a widespread superfamily of enzymes notably involved in detoxification processes and/or in specialized metabolism. In the cyanobacterium Synechocsytis sp. PCC 6803, SynGSTC1, a chi-class GST (GSTC), is thought to participate in the detoxification process of methylglyoxal, a toxic by-product of cellular metabolism. A comparative genomic analysis showed that GSTCs were present in all orders of cyanobacteria with the exception of the basal order Gloeobacterales. These enzymes were also detected in some marine and freshwater noncyanobacterial bacteria, probably as a result of horizontal gene transfer events. GSTCs were shorter of about 30 residues compared to most cytosolic GSTs and had a well-conserved SRAS motif in the active site (¹⁰SRAS¹³ in SynGSTC1). The crystal structure of SynGSTC1 in complex with glutathione adopted the canonical GST fold with a very open active site because the α4 and α5 helices were exceptionally short. A transferred multipolar electron-density analysis allowed a fine description of the solved structure. Unexpectedly, Ser10 did not have an electrostatic influence on glutathione as usually observed in serinyl-GSTs. The S10A variant was only slightly less efficient than the wild-type and molecular dynamics simulations suggested that S10 was a stabilizer of the protein backbone rather than an anchor site for glutathione.

Contamination detection in genomic data: more is not enough

The decreasing cost of sequencing and concomitant augmentation of publicly available genomes have created an acute need for automated software to assess genomic contamination. During the last 6 years, 18 programs have been published, each with its own strengths and weaknesses. Deciding which tools to use becomes more and more difficult without an understanding of the underlying algorithms. We review these programs, benchmarking six of them, and present their main operating principles. This article is intended to guide researchers in the selection of appropriate tools for specific applications. Finally, we present future challenges in the developing field of contamination detection.

A divide-and-conquer phylogenomic approach based on character supermatrices resolves early steps in the evolution of the Archaea

Background

The recent rise in cultivation-independent genome sequencing has provided key material to explore uncharted branches of the Tree of Life. This has been particularly spectacular concerning the Archaea, projecting them at the center stage as prominently relevant to understand early stages in evolution and the emergence of fundamental metabolisms as well as the origin of eukaryotes. Yet, resolving deep divergences remains a challenging task due to well-known tree-reconstruction artefacts and biases in extracting robust ancient phylogenetic signal, notably when analyzing data sets including the three Domains of Life. Among the various strategies aimed at mitigating these problems, divide-and-conquer approaches remain poorly explored, and have been primarily based on reconciliation among single gene trees which however notoriously lack ancient phylogenetic signal.

Results

We analyzed sub-sets of full supermatrices covering the whole Tree of Life with specific taxonomic sampling to robustly resolve different parts of the archaeal phylogeny in light of their current diversity. Our results strongly support the existence and early emergence of two main clades, Cluster I and Cluster II, which we name Ouranosarchaea and Gaiarchaea, and we clarify the placement of important novel archaeal lineages within these two clades. However, the monophyly and branching of the fast evolving nanosized DPANN members remains unclear and worth of further study.

Conclusions

We inferred a well resolved rooted phylogeny of the Archaea that includes all recently described phyla of high taxonomic rank. This phylogeny represents a valuable reference to study the evolutionary events associated to the early steps of the diversification of the archaeal domain. Beyond the specifics of archaeal phylogeny, our results demonstrate the power of divide-and-conquer approaches to resolve deep phylogenetic relationships, which should be applied to progressively resolve the entire Tree of Life.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01952-0.

ORPER: A Workflow for Constrained SSU rRNA Phylogenies

The continuous increase in sequenced genomes in public repositories makes the choice of interesting bacterial strains for future sequencing projects ever more complicated, as it is difficult to estimate the redundancy between these strains and the already available genomes. Therefore, we developed the Nextflow workflow “ORPER”, for “ORganism PlacER”, containerized in Singularity, which allows the determination the phylogenetic position of a collection of organisms in the genomic landscape. ORPER constrains the phylogenetic placement of SSU (16S) rRNA sequences in a multilocus reference tree based on ribosomal protein genes extracted from public genomes. We demonstrate the utility of ORPER on the Cyanobacteria phylum, by placing 152 strains of the BCCM/ULC collection.

A Comprehensive Evolutionary Scenario of Cell Division and Associated Processes in the Firmicutes

The cell cycle is a fundamental process that has been extensively studied in bacteria. However, many of its components and their interactions with machineries involved in other cellular processes are poorly understood. Furthermore, most knowledge relies on the study of a few models, but the real diversity of the cell division apparatus and its evolution are largely unknown. Here, we present a massive in-silico analysis of cell division and associated processes in around 1,000 genomes of the Firmicutes, a major bacterial phylum encompassing models (i.e. Bacillus subtilis, Streptococcus pneumoniae, and Staphylococcus aureus), as well as many important pathogens. We analyzed over 160 proteins by using an original approach combining phylogenetic reconciliation, phylogenetic profiles, and gene cluster survey. Our results reveal the presence of substantial differences among clades and pinpoints a number of evolutionary hotspots. In particular, the emergence of Bacilli coincides with an expansion of the gene repertoires involved in cell wall synthesis and remodeling. We also highlight major genomic rearrangements at the emergence of Streptococcaceae. We establish a functional network in Firmicutes that allows identifying new functional links inside one same process such as between FtsW (peptidoglycan polymerase) and a previously undescribed Penicilin-Binding Protein or between different processes, such as replication and cell wall synthesis. Finally, we identify new candidates involved in sporulation and cell wall synthesis. Our results provide a previously undescribed view on the diversity of the bacterial cell cycle, testable hypotheses for further experimental studies, and a methodological framework for the analysis of any other biological system.

ToRQuEMaDA: tool for retrieving queried Eubacteria, metadata and dereplicating assemblies

TQMD is a tool for high-performance computing clusters which downloads, stores and produces lists of dereplicated prokaryotic genomes. It has been developed to counter the ever-growing number of prokaryotic genomes and their uneven taxonomic distribution. It is based on word-based alignment-free methods (k-mers), an iterative single-linkage approach and a divide-and-conquer strategy to remain both efficient and scalable. We studied the performance of TQMD by verifying the influence of its parameters and heuristics on the clustering outcome. We further compared TQMD to two other dereplication tools (dRep and Assembly-Dereplicator). Our results showed that TQMD is primarily optimized to dereplicate at higher taxonomic levels (phylum/class), as opposed to the other dereplication tools, but also works at lower taxonomic levels (species/strain) like the other dereplication tools. TQMD is available from source and as a Singularity container at [https://bitbucket.org/phylogeno/tqmd ].

Elevated rates of horizontal gene transfer in the industrialized human microbiome

Industrialization has impacted the human gut ecosystem, resulting in altered microbiome composition and diversity. Whether bacterial genomes may also adapt to the industrialization of their host populations remains largely unexplored. Here, we investigate the extent to which the rates and targets of horizontal gene transfer (HGT) vary across thousands of bacterial strains from 15 human populations spanning a range of industrialization. We show that HGTs have accumulated in the microbiome over recent host generations and that HGT occurs at high frequency within individuals. Comparison across human populations reveals that industrialized lifestyles are associated with higher HGT rates and that the functions of HGTs are related to the level of host industrialization. Our results suggest that gut bacteria continuously acquire new functionality based on host lifestyle and that high rates of HGT may be a recent development in human history linked to industrialization.

A library of human gut bacterial isolates paired with longitudinal multiomics data enables mechanistic microbiome research

Our understanding of how the gut microbiome interacts with its human host has been restrained by limited access to longitudinal datasets to examine stability and dynamics, and by having only a few isolates to test mechanistic hypotheses. Here, we present the Broad Institute-OpenBiome Microbiome Library (BIO-ML), a comprehensive collection of 7,758 gut bacterial isolates paired with 3,632 genome sequences and longitudinal multi-omics data. We show that microbial species maintain stable population sizes within and across humans and that commonly used 'omics' survey methods are more reliable when using averages over multiple days of sampling. Variation of gut metabolites within people over time is associated with amino acid levels, and differences across people are associated with differences in bile acids. Finally, we show that genomic diversification can be used to infer eco-evolutionary dynamics and in vivo selection pressures for strains within individuals. The BIO-ML is a unique resource designed to enable hypothesis-driven microbiome research.

Taxonomic assignment of uncultured prokaryotes with long range PCR targeting the spectinomycin operon

The taxonomic assignment of uncultured prokaryotes to known taxa is a major challenge in microbial systematics. This relies usually on the phylogenetic analysis of the ribosomal small subunit RNA or a few housekeeping genes. Recent works have disclosed ribosomal proteins as valuable markers for systematics and, due to the boom in complete genome sequencing, their use has become widespread. Yet, in the case of uncultured strains, for which complete genome sequences cannot be easily obtained, sequencing many markers is complicated and time consuming. Taking the advantage of the organization of ribosomal protein coding genes in large gene clusters, we amplified a 32 kb conserved region encompassing the spectinomycin (spc) operon using long range PCR from isolated and from uncultured nodular endophytic Frankia strains. The phylogenetic analysis of the 27 ribosomal protein genes contained in this region provided a robust phylogenetic tree consistent with phylogenies based on larger set of markers, indicating that this subset of ribosomal proteins contains enough phylogenetic signal to address systematic issues. This work shows that using long range PCR could break down the barrier preventing the use of ribosomal proteins as phylogenetic markers when complete genome sequences cannot be easily obtained.

The phytopathogenic nature of Dickeya aquatica 174/2 and the dynamic early evolution of Dickeya pathogenicity

Dickeya is a genus of phytopathogenic enterobacterales causing soft rot in a variety of plants (e.g. potato, chicory, maize). Among the species affiliated to this genus, Dickeya aquatica, described in 2014, remained particularly mysterious because it had no known host. Furthermore, while D. aquatica was proposed to represent a deep-branching species among Dickeya genus, its precise phylogenetic position remained elusive. Here, we report the complete genome sequence of the D. aquatica type strain 174/2. We demonstrate the affinity of D. aquatica strain 174/2 for acidic fruits such as tomato and cucumber and show that exposure of this bacterium to acidic pH induces twitching motility. An in-depth phylogenomic analysis of all available Dickeya proteomes pinpoints D. aquatica as the second deepest branching lineage within this genus and reclassifies two lineages that likely correspond to new genomospecies (gs.): Dickeya gs. poaceaephila (Dickeya sp NCPPB 569) and Dickeya gs. undicola (Dickeya sp 2B12), together with a new putative genus, tentatively named Prodigiosinella. Finally, from comparative analyses of Dickeya proteomes, we infer the complex evolutionary history of this genus, paving the way to study the adaptive patterns and processes of Dickeya to different environmental niches and hosts. In particular, we hypothesize that the lack of xylanases and xylose degradation pathways in D. aquatica could reflect adaptation to aquatic charophyte hosts which, in contrast to land plants, do not contain xyloglucans.

A Novel Subfamily of Bacterial AAT-Fold Basic Amino Acid Decarboxylases and Functional Characterization of Its First Representative: Pseudomonas aeruginosa LdcA

Polyamines are small amino-acid derived polycations capable of binding negatively charged macromolecules. Bacterial polyamines are structurally and functionally diverse, and are mainly produced biosynthetically by pyridoxal-5-phosphate-dependent amino acid decarboxylases referred to as Lysine-Arginine-Ornithine decarboxylases (LAOdcs). In a phylogenetically limited group of bacteria, LAOdcs are also induced in response to acid stress. Here, we performed an exhaustive phylogenetic analysis of the AAT-fold LAOdcs which showcased the ancient nature of their short forms in Cyanobacteria and Firmicutes, and emergence of distinct subfamilies of long LAOdcs in Proteobacteria. We identified a novel subfamily of lysine decarboxylases, LdcA, ancestral in Betaproteobacteria and Pseudomonadaceae. We analyzed the expression of LdcA from Pseudomonas aeruginosa, and uncovered its role, intimately linked to cadaverine (Cad) production, in promoting growth and reducing persistence of this multidrug resistant human pathogen during carbenicillin treatment. Finally, we documented a certain redundancy in the function of the three main polyamines—Cad, putrescine (Put), and spermidine (Spd)—in P. aeruginosa by demonstrating the link between their intracellular level, as well as the capacity of Put and Spd to complement the growth phenotype of the ldcA mutant.

Metagenomic assembly of new (sub)polar Cyanobacteria and their associated microbiome from non-axenic cultures

Cyanobacteria form one of the most diversified phyla of Bacteria. They are important ecologically as primary producers, for Earth evolution and biotechnological applications. Yet, Cyanobacteria are notably difficult to purify and grow axenically, and most strains in culture collections contain heterotrophic bacteria that were probably associated with Cyanobacteria in the environment. Obtaining cyanobacterial DNA without contaminant sequences is thus a challenging and time-consuming task. Here, we describe a metagenomic pipeline that enables the easy recovery of genomes from non-axenic cultures. We tested this pipeline on 17 cyanobacterial cultures from the BCCM/ULC public collection and generated novel genome sequences for 12 polar or subpolar strains and three temperate ones, including three early-branching organisms that will be useful for phylogenomics. In parallel, we assembled 31 co-cultivated bacteria (12 nearly complete) from the same cultures and showed that they mostly belong to Bacteroidetes and Proteobacteria, some of them being very closely related in spite of geographically distant sampling sites.

Consensus assessment of the contamination level of publicly available cyanobacterial genomes

Publicly available genomes are crucial for phylogenetic and metagenomic studies, in which contaminating sequences can be the cause of major problems. This issue is expected to be especially important for Cyanobacteria because axenic strains are notoriously difficult to obtain and keep in culture. Yet, despite their great scientific interest, no data are currently available concerning the quality of publicly available cyanobacterial genomes. As reliably detecting contaminants is a complex task, we designed a pipeline combining six methods in a consensus strategy to assess the contamination level of 440 genome assemblies of Cyanobacteria. Two methods are based on published reference databases of ribosomal genes (SSU rRNA 16S and ribosomal proteins), one is indirectly based on a reference database of marker genes (CheckM), and three are based on complete genome analysis. Among those genome-wide methods, Kraken and DIAMOND blastx share the same reference database that we derived from Ensembl Bacteria, whereas CONCOCT does not require any reference database, instead relying on differences in DNA tetramer frequencies. Given that all the six methods appear to have their own strengths and limitations, we used the consensus of their rankings to infer that >5% of cyanobacterial genome assemblies are highly contaminated by foreign DNA (i.e., contaminants were detected by 5 or 6 methods). Our results will help researchers to check the quality of publicly available genomic data before use in their own analyses. Moreover, we argue that journals should make mandatory the submission of raw read data along with genome assemblies in order to facilitate the detection of contaminants in sequence databases.

PASTA repeats of the protein kinase StkP interconnect cell constriction and separation of Streptococcus pneumoniae

Eukaryotic-like serine/threonine kinases (eSTKs) with extracellular PASTA repeats are key membrane regulators of bacterial cell division. How PASTA repeats govern eSTK activation and function remains elusive. Using evolution- and structural-guided approaches combined with cell imaging, we disentangle the role of each PASTA repeat of the eSTK StkP from Streptococcus pneumoniae. While the three membrane-proximal PASTA repeats behave as interchangeable modules required for the activation of StkP independently of cell wall binding, they also control the septal cell wall thickness. In contrast, the fourth and membrane-distal PASTA repeat directs StkP localization at the division septum and encompasses a specific motif that is critical for final cell separation through interaction with the cell wall hydrolase LytB. We propose a model in which the extracellular four-PASTA domain of StkP plays a dual function in interconnecting the phosphorylation of StkP endogenous targets along with septal cell wall remodelling to allow cell division of the pneumococcus.