Identical sequences found in distant genomes reveal frequent horizontal transfer across the bacterial domain

Sequencing-based analysis of microbiomes

Microbiomes occupy a range of niches and, in addition to having diverse compositions, they have varied functional roles that have an impact on agriculture, environmental sciences, and human health and disease. The study of microbiomes has been facilitated by recent technological and analytical advances, such as cheaper and higher-throughput DNA and RNA sequencing, improved long-read sequencing and innovative computational analysis methods. These advances are providing a deeper understanding of microbiomes at the genomic, transcriptional and translational level, generating insights into their function and composition at resolutions beyond the species level.

In Silico Analysis of the Phylogenetic and Physiological Characteristics of Sphingobium indicum B90A: A Hexachlorocyclohexane-Degrading Bacterium

The study focuses on the in silico genomic characterization of Sphingobium indicum B90A, revealing a wealth of genes involved in stress response, carbon monoxide oxidation, β-carotene biosynthesis, heavy metal resistance, and aromatic compound degradation, suggesting its potential as a bioremediation agent. Furthermore, genomic adaptations among nine Sphingomonad strains were explored, highlighting shared core genes via pangenome analysis, including those related to the shikimate pathway and heavy metal resistance. The majority of genes associated with aromatic compound degradation, heavy metal resistance, and stress response were found within genomic islands across all strains. Sphingobium indicum UT26S exhibited the highest number of genomic islands, while Sphingopyxis alaskensis RB2256 had the maximum fraction of its genome covered by genomic islands. The distribution of lin genes varied among the strains, indicating diverse genetic responses to environmental pressures. Additionally, in silico evidence of horizontal gene transfer (HGT) between plasmids pSRL3 and pISP3 of the Sphingobium and Sphingomonas genera, respectively, has been provided. The manuscript offers novel insights into strain B90A, highlighting its role in horizontal gene transfer and refining evolutionary relationships among Sphingomonad strains. The discovery of stress response genes and the czcABCD operon emphasizes the potential of Sphingomonads in consortia development, supported by genomic island analysis.

Systematic identification of cargo-mobilizing genetic elements reveals new dimensions of eukaryotic diversity

Cargo-mobilizing mobile elements (CMEs) are genetic entities that faithfully transpose diverse protein coding sequences. Although common in bacteria, we know little about eukaryotic CMEs because no appropriate tools exist for their annotation. For example, Starships are giant fungal CMEs whose functions are largely unknown because they require time-intensive manual curation. To address this knowledge gap, we developed starfish, a computational workflow for high-throughput eukaryotic CME annotation. We applied starfish to 2 899 genomes of 1 649 fungal species and found that starfish recovers known Starships with 95% combined precision and recall while expanding the number of annotated elements ten-fold. Extant Starship diversity is partitioned into 11 families that differ in their enrichment patterns across fungal classes. Starship cargo changes rapidly such that elements from the same family differ substantially in their functional repertoires, which are predicted to contribute to diverse biological processes such as metabolism. Many elements have convergently evolved to insert into 5S rDNA and AT-rich sequence while others integrate into random locations, revealing both specialist and generalist strategies for persistence. Our work establishes a framework for advancing mobile element biology and provides the means to investigate an emerging dimension of eukaryotic genetic diversity, that of genomes within genomes.

Spread of antibiotic resistance genes in drinking water reservoirs: Insights from a deep metagenomic study using a curated database

The exploration of antibiotic resistance genes (ARGs) in drinking water reservoirs is an emerging field. Using a curated database, we enhanced the ARG detection and conducted a comprehensive analysis using 2.2 Tb of deep metagenomic sequencing data to determine the distribution of ARGs across 16 drinking water reservoirs and associated environments. Our findings reveal a greater diversity of ARGs in sediments than in water, underscoring the importance of extensive background surveys. Crucial ARG carriers-specifically Acinetobacter, Pseudomonas, and Mycobacterium were identified in drinking water reservoirs. Extensive analysis of the data uncovered a considerable concern for drinking water safety, particularly in regions reliant on river sources. Mobile genetic elements have been found to contribute markedly to the propagation of ARGs. The results of this research suggest that the establishment of drinking water reservoirs for supplying raw water may be an effective strategy for alleviating the spread of water-mediated ARGs.

A global survey of prokaryotic genomes reveals the eco-evolutionary pressures driving horizontal gene transfer

Horizontal gene transfer, the exchange of genetic material through means other than reproduction, is a fundamental force in prokaryotic genome evolution. Genomic persistence of horizontally transferred genes has been shown to be influenced by both ecological and evolutionary factors. However, there is limited availability of ecological information about species other than the habitats from which they were isolated, which has prevented a deeper exploration of ecological contributions to horizontal gene transfer. Here we focus on transfers detected through comparison of individual gene trees to the species tree, assessing the distribution of gene-exchanging prokaryotes across over a million environmental sequencing samples. By analysing detected horizontal gene transfer events, we show distinct functional profiles for recent versus old events. Although most genes transferred are part of the accessory genome, genes transferred earlier in evolution tend to be more ubiquitous within present-day species. We find that co-occurring, interacting and high-abundance species tend to exchange more genes. Finally, we show that host-associated specialist species are most likely to exchange genes with other host-associated specialist species, whereas species found across different habitats have similar gene exchange rates irrespective of their preferred habitat. Our study covers an unprecedented scale of integrated horizontal gene transfer and environmental information, highlighting broad eco-evolutionary trends.

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.

The defensome of complex bacterial communities

Bacteria have developed various defense mechanisms to avoid infection and killing in response to the fast evolution and turnover of viruses and other genetic parasites. Such pan-immune system (defensome) encompasses a growing number of defense lines that include well-studied innate and adaptive systems such as restriction-modification, CRISPR-Cas and abortive infection, but also newly found ones whose mechanisms are still poorly understood. While the abundance and distribution of defense systems is well-known in complete and culturable genomes, there is a void in our understanding of their diversity and richness in complex microbial communities. Here we performed a large-scale in-depth analysis of the defensomes of 7759 high-quality bacterial population genomes reconstructed from soil, marine, and human gut environments. We observed a wide variation in the frequency and nature of the defensome among large phyla, which correlated with lifestyle, genome size, habitat, and geographic background. The defensome's genetic mobility, its clustering in defense islands, and genetic variability was found to be system-specific and shaped by the bacterial environment. Hence, our results provide a detailed picture of the multiple immune barriers present in environmentally distinct bacterial communities and set the stage for subsequent identification of novel and ingenious strategies of diversification among uncultivated microbes.

Microbial divergence and evolution. The case of anammox bacteria

Species differentiation and the appearance of novel diversity on Earth is a major issue to understand the past and future of microbial evolution. Herein, we propose the analysis of a singular evolutive example, the case of microorganisms carrying out the process of anammox (anaerobic ammonium oxidation). Anammox represents a singular physiology active on Earth from ancient times and, at present, this group is still represented by a relatively limited number of species carrying out a specific metabolism within the Phylum Planctomycetota. The key enzyme on the anammox pathway is hydrazine dehydrogenase (HDH) which has been used as a model in this study. HDH and rRNA (16S subunit) phylogenies are in agreement suggesting a monophyletic origin. The diversity of this singular phylogenetic group is represented by a few enriched bacterial consortia awaiting to be cultured as monospecific taxa. The apparent evolution of the HDH genes in these anammox bacteria is highly related to the diversification of the anammox clades and their genomes as pointed by phylogenomics, their GC content and codon usage profile. This study represents a clear case where bacterial evolution presents a paralleled genome, gene and species diversification through time from a common ancestor; a scenario that most times is masked by a web-like phylogeny and the huge complexity within the prokaryotes. Besides, this contribution suggests that microbial evolution of the anammox bacteria has followed an ordered, vertical diversification through Earth history and will present a potentially similar speciation fate in the future.

Dynamics of bacterial recombination in the human gut microbiome

Horizontal gene transfer (HGT) is a ubiquitous force in microbial evolution. Previous work has shown that the human gut is a hotspot for gene transfer between species, but the more subtle exchange of variation within species-also known as recombination-remains poorly characterized in this ecosystem. Here, we show that the genetic structure of the human gut microbiome provides an opportunity to measure recent recombination events from sequenced fecal samples, enabling quantitative comparisons across diverse commensal species that inhabit a common environment. By analyzing recent recombination events in the core genomes of 29 human gut bacteria, we observed widespread heterogeneities in the rates and lengths of transferred fragments, which are difficult to explain by existing models of ecological isolation or homology-dependent recombination rates. We also show that natural selection helps facilitate the spread of genetic variants across strain backgrounds, both within individual hosts and across the broader population. These results shed light on the dynamics of in situ recombination, which can strongly constrain the adaptability of gut microbial communities.

Modeling the limits of detection for antimicrobial resistance genes in agri-food samples: a comparative analysis of bioinformatics tools

BACKGROUND

Although the spread of antimicrobial resistance (AMR) through food and its production poses a significant concern, there is limited research on the prevalence of AMR bacteria in various agri-food products. Sequencing technologies are increasingly being used to track the spread of AMR genes (ARGs) in bacteria, and metagenomics has the potential to bypass some of the limitations of single isolate characterization by allowing simultaneous analysis of the agri-food product microbiome and associated resistome. However, metagenomics may still be hindered by methodological biases, presence of eukaryotic DNA, and difficulties in detecting low abundance targets within an attainable sequence coverage. The goal of this study was to assess whether limits of detection of ARGs in agri-food metagenomes were influenced by sample type and bioinformatic approaches.

RESULTS

We simulated metagenomes containing different proportions of AMR pathogens and analysed them for taxonomic composition and ARGs using several common bioinformatic tools. Kraken2/Bracken estimates of species abundance were closest to expected values. However, analysis by both Kraken2/Bracken indicated presence of organisms not included in the synthetic metagenomes. Metaphlan3/Metaphlan4 analysis of community composition was more specific but with lower sensitivity than the Kraken2/Bracken analysis. Accurate detection of ARGs dropped drastically below 5X isolate genome coverage. However, it was sometimes possible to detect ARGs and closely related alleles at lower coverage levels if using a lower ARG-target coverage cutoff (< 80%). While KMA and CARD-RGI only predicted presence of expected ARG-targets or closely related gene-alleles, SRST2 (which allows read to map to multiple targets) falsely reported presence of distantly related ARGs at all isolate genome coverage levels. The presence of background microbiota in metagenomes influenced the accuracy of ARG detection by KMA, resulting in mcr-1 detection at 0.1X isolate coverage in the lettuce but not in the beef metagenome.

CONCLUSIONS

This study demonstrates accurate detection of ARGs in synthetic metagenomes using various bioinformatic methods, provided that reads from the ARG-encoding organism exceed approximately 5X isolate coverage (i.e. 0.4% of a 40 million read metagenome). While lowering thresholds for target gene detection improved sensitivity, this led to the identification of alternative ARG-alleles, potentially confounding the identification of critical ARGs in the resistome. Further advancements in sequencing technologies providing increased coverage depth or extended read lengths may improve ARG detection in agri-food metagenomic samples, enabling use of this approach for tracking clinically important ARGs in agri-food samples.

Detecting horizontal gene transfer among microbiota: an innovative pipeline for identifying co-shared genes within the mobilome through advanced comparative analysis

Horizontal gene transfer (HGT) is a key driver in the evolution of bacterial genomes. The acquisition of genes mediated by HGT may enable bacteria to adapt to ever-changing environmental conditions. Long-term application of antibiotics in intensive agriculture is associated with the dissemination of antibiotic resistance genes among bacteria with the consequences causing public health concern. Commensal farm-animal-associated gut microbiota are considered the reservoir of the resistance genes. Therefore, in this study, we identified known and not-yet characterized mobilized genes originating from chicken and porcine fecal samples using our innovative pipeline followed by network analysis to provide appropriate visualization to support proper interpretation.

Interplay Between Antimicrobial Resistance and Global Environmental Change

Antibiotic resistance genes predate the therapeutic uses of antibiotics. However, the current antimicrobial resistance crisis stems from our extensive use of antibiotics and the generation of environmental stressors that impose new selective pressure on microbes and drive the evolution of resistant pathogens that now threaten human health. Similar to climate change, this global threat results from human activities that change habitats and natural microbiomes, which in turn interact with human-associated ecosystems and lead to adverse impacts on human health. Human activities that alter our planet at global scales exacerbate the current resistance crisis and exemplify our central role in large-scale changes in which we are both protagonists and architects of our success but also casualties of unanticipated collateral outcomes. As cognizant participants in this ongoing planetary experiment, we are driven to understand and find strategies to curb the ongoing crises of resistance and climate change.

Intraspecific microdiversity and ecological drivers of lactic acid bacteria in naturally fermented milk ecosystem

Traditional fermented milks are produced by inoculating technique, which selects well-adapted microorganisms that have been passed on through generations. Few reports have used naturally fermented milks as model ecosystems to investigate the mechanism of formation of intra-species microbial diversity. Here, we isolated and whole-genome-sequenced a total of 717 lactic acid bacterial isolates obtained from 12 independent naturally fermented milks collect from 12 regions across five countries. We further analyzed the within-sample intra-species phylogenies of 214 Lactobacillus helveticus isolates, 97 Lactococcus lactis subsp. lactis isolates, and 325 Lactobacillus delbrueckii subsp. bulgaricus isolates. We observed a high degree of intra-species genomic and functional gene diversity within-/between-sample(s). Single nucleotide polymorphism-based phylogenetic reconstruction revealed great within-sample intra-species heterogeneity, evolving from multiple lineages. Further phylogenetic reconstruction (presence-absence gene profile) revealed within-sample inter-clade functional diversity (based on carbohydrate-active enzyme- and peptidase-encoding genes) in all three investigated species/subspecies. By identifying and mapping clade-specific genes of intra-sample clades of the three species/subspecies to the respective fermented milk metagenome, we found extensive potential inter-/intra-species horizontal gene transfer events. Finally, the microbial composition of the samples is closely linked to the nucleotide diversity of the respective species/subspecies. Overall, our results contribute to the conservation of lactic acid bacteria resources, providing ecological insights into the microbial ecosystem of naturally fermented dairy products.

Unipept Desktop 2.0: Construction of Targeted Reference Protein Databases for Metaproteogenomics Analyses

Unipept Desktop 2.0 is the most recent iteration of the Unipept Desktop tool that adds support for the analysis of metaproteogenomics datasets. Unipept Desktop now supports the automatic construction of targeted protein reference databases that only contain proteins (originating from the UniProtKB resource) associated with a predetermined list of taxa. This improves both the taxonomic and functional resolution of a metaproteomic analysis and yields several technical advantages. By limiting the proteins present in a reference database, it is also possible to perform (meta)proteogenomics analyses. Since the protein reference database resides on the user's local machine, they have complete control over the database used during an analysis. Data no longer need to be transmitted over the Internet, decreasing the time required for an analysis and better safeguarding privacy-sensitive data. As a proof of concept, we present a case study in which a human gut metaproteome dataset is analyzed with Unipept Desktop 2.0 using different targeted databases based on matched 16S rRNA gene sequencing data.

Origin and Early Diversification of the Papain Family of Cysteine Peptidases

Peptidases of the papain family play a key role in protein degradation, regulated proteolysis, and the host-pathogen arms race. Although the papain family has been the subject of many studies, knowledge about its diversity, origin, and evolution in Eukaryota, Bacteria, and Archaea is limited; thus, we aimed to address these long-standing knowledge gaps. We traced the origin and expansion of the papain family with a phylogenomic analysis, using sequence data from numerous prokaryotic and eukaryotic proteomes, transcriptomes, and genomes. We identified the full complement of the papain family in all prokaryotic and eukaryotic lineages. Analysis of the papain family provided strong evidence for its early diversification in the ancestor of eukaryotes. We found that the papain family has undergone complex and dynamic evolution through numerous gene duplications, which produced eight eukaryotic ancestral paralogous C1A lineages during eukaryogenesis. Different evolutionary forces operated on C1A peptidases, including gene duplication, horizontal gene transfer, and gene loss. This study challenges the current understanding of the origin and evolution of the papain family and provides valuable insights into their early diversification. The findings of this comprehensive study provide guidelines for future structural and functional studies of the papain family.

The Promise and Pitfalls of Prophages

Phages dominate every ecosystem on the planet. While virulent phages sculpt the microbiome by killing their bacterial hosts, temperate phages provide unique growth advantages to their hosts through lysogenic conversion. Many prophages benefit their host, and prophages are responsible for genotypic and phenotypic differences that separate individual microbial strains. However, the microbes also endure a cost to maintain those phages: additional DNA to replicate and proteins to transcribe and translate. We have never quantified those benefits and costs. Here, we analysed over two and a half million prophages from over half a million bacterial genome assemblies. Analysis of the whole dataset and a representative subset of taxonomically diverse bacterial genomes demonstrated that the normalised prophage density was uniform across all bacterial genomes above 2 Mbp. We identified a constant carrying capacity of phage DNA per bacterial DNA. We estimated that each prophage provides cellular services equivalent to approximately 2.4 % of the cell's energy or 0.9 ATP per bp per hour. We demonstrate analytical, taxonomic, geographic, and temporal disparities in identifying prophages in bacterial genomes that provide novel targets for identifying new phages. We anticipate that the benefits bacteria accrue from the presence of prophages balance the energetics involved in supporting prophages. Furthermore, our data will provide a new framework for identifying phages in environmental datasets, diverse bacterial phyla, and from different locations.

Serratia marcescens antibiotic resistance mechanisms of an opportunistic pathogen: a literature review

Serratia marcescens is a ubiquitous bacterium from order Enterobacterales displaying a high genetic plasticity that allows it to adapt and persist in multiple niches including soil, water, plants, and nosocomial environments. Recently, S. marcescens has gained attention as an emerging pathogen worldwide, provoking infections and outbreaks in debilitated individuals, particularly newborns and patients in intensive care units. S. marcescens isolates recovered from clinical settings are frequently described as multidrug resistant. High levels of antibiotic resistance across Serratia species are a consequence of the combined activity of intrinsic, acquired, and adaptive resistance elements. In this review, we will discuss recent advances in the understanding of mechanisms guiding resistance in this opportunistic pathogen.

Leveraging comparative genomics to uncover alien genes in bacterial genomes

A significant challenge in bacterial genomics is to catalogue genes acquired through the evolutionary process of horizontal gene transfer (HGT). Both comparative genomics and sequence composition-based methods have often been invoked to quantify horizontally acquired genes in bacterial genomes. Comparative genomics methods rely on completely sequenced genomes and therefore the confidence in their predictions increases as the databases become more enriched in completely sequenced genomes. Recent developments including in microbial genome sequencing call for reassessment of alien genes based on information-rich resources currently available. We revisited the comparative genomics approach and developed a new algorithm for alien gene detection. Our algorithm compared favourably with the existing comparative genomics-based methods and is capable of detecting both recent and ancient transfers. It can be used as a standalone tool or in concert with other complementary algorithms for comprehensively cataloguing alien genes in bacterial genomes.

High genomic differentiation and limited gene flow indicate recent cryptic speciation within the genus Laspinema (cyanobacteria)

The sympatric occurrence of closely related lineages displaying conserved morphological and ecological traits is often characteristic of free-living microbes. Gene flow, recombination, selection, and mutations govern the genetic variability between these cryptic lineages and drive their differentiation. However, sequencing conservative molecular markers (e.g., 16S rRNA) coupled with insufficient population-level sampling hindered the study of intra-species genetic diversity and speciation in cyanobacteria. We used phylogenomics and a population genomic approach to investigate the extent of local genomic diversity and the mechanisms underlying sympatric speciation of Laspinema thermale. We found two cryptic lineages of Laspinema. The lineages were highly genetically diverse, with recombination occurring more frequently within than between them. That suggests the existence of a barrier to gene flow, which further maintains divergence. Genomic regions of high population differentiation harbored genes associated with possible adaptations to high/low light conditions and stress stimuli, although with a weak diversifying selection. Overall, the diversification of Laspinema species might have been affected by both genomic and ecological processes.

Assessing Microbial Metabolic and Biological Diversity to Inform Natural Product Library Assembly

The pressing need for novel chemical matter to support bioactive compound discovery has led natural product researchers to explore a wide range of source organisms and environments. One of the implicit guiding principles behind those efforts is the notion that sampling different environments is critical to accessing unique natural products. This idea was tested by comparing fungi from disparate biomes: aquatic sediments from Lake Michigan (USA) and terrestrial samples taken from the surrounding soils. Matched sets of Penicillium brevicompactum, Penicillium expansum, and Penicillium oxalicum from the two source environments were compared, revealing modest differences in physiological performance and chemical output. Analysis of LC-MS/MS-derived molecular feature data showed no source-dependent differences in chemical richness. High levels of scaffold homogeneity were also observed with 78-83% of scaffolds shared among the terrestrial and aquatic Penicillium spp. isolates. A comparison of the culturable fungi from the two biomes indicated that certain genera were more strongly associated with aquatic sediments (e.g., Trichoderma, Pseudeurotium, Cladosporium, and Preussia) versus the surrounding terrestrial environment (e.g., Fusarium, Pseudogymnoascus, Humicola, and Acremonium). Taken together, these results suggest that focusing efforts on sampling the microbial resources that are unique to an environment may have a more pronounced effect on enhancing the sought-after natural product diversity needed for chemical discovery and screening collections.