Disentangling the impact of environmental and phylogenetic constraints on prokaryotic within-species diversity

Genomics-based identification of a cold adapted clade in Deinococcus

BACKGROUND

Microbes in the cold polar and alpine environments play a critical role in feedbacks that amplify the effects of climate change. Defining the cold adapted ecotype is one of the prerequisites for understanding the response of polar and alpine microbes to climate change.

RESULTS

Here, we analysed 85 high-quality, de-duplicated genomes of Deinococcus, which can survive in a variety of harsh environments. By leveraging genomic and phenotypic traits with reverse ecology, we defined a cold adapted clade from eight Deinococcus strains isolated from Arctic, Antarctic and high alpine environments. Genome-wide optimization in amino acid composition and regulation and signalling enable the cold adapted clade to produce CO2 from organic matter and boost the bioavailability of mineral nitrogen.

CONCLUSIONS

Based primarily on in silico genomic analysis, we defined a potential cold adapted clade in Deinococcus and provided an updated view of the genomic traits and metabolic potential of Deinococcus. Our study would facilitate the understanding of microbial processes in the cold polar and alpine environments.

Strain dynamics of contaminating bacteria modulate the yield of ethanol biorefineries

Bioethanol is a sustainable energy alternative and can contribute to global greenhouse-gas emission reductions by over 60%. Its industrial production faces various bottlenecks, including sub-optimal efficiency resulting from bacteria. Broad-spectrum removal of these contaminants results in negligible gains, suggesting that the process is shaped by ecological interactions within the microbial community. Here, we survey the microbiome across all process steps at two biorefineries, over three timepoints in a production season. Leveraging shotgun metagenomics and cultivation-based approaches, we identify beneficial bacteria and find improved outcome when yeast-to-bacteria ratios increase during fermentation. We provide a microbial gene catalogue which reveals bacteria-specific pathways associated with performance. We also show that Limosilactobacillus fermentum overgrowth lowers production, with one strain reducing yield by ~5% in laboratory fermentations, potentially due to its metabolite profile. Temperature is found to be a major driver for strain-level dynamics. Improved microbial management strategies could unlock environmental and economic gains in this US $ 60 billion industry enabling its wider adoption.

Bacterial lifestyle shapes pangenomes

Pangenomes vary across bacteria. Some species have fluid pangenomes, with a high proportion of genes varying between individual genomes. Other species have less fluid pangenomes, with different genomes tending to contain the same genes. Two main hypotheses have been suggested to explain this variation: differences in species' bacterial lifestyle and effective population size. However, previous studies have not been able to test between these hypotheses because the different features of lifestyle and effective population size are highly correlated with each other, and phylogenetically conserved, making it hard to disentangle their relative importance. We used phylogeny-based analyses, across 126 bacterial species, to tease apart the causal role of different factors. We found that pangenome fluidity was lower in i) host-associated compared with free-living species and ii) host-associated species that are obligately dependent on a host, live inside cells, and are more pathogenic and less motile. In contrast, we found no support for the competing hypothesis that larger effective population sizes lead to more fluid pangenomes. Effective population size appears to correlate with pangenome variation because it is also driven by bacterial lifestyle, rather than because of a causal relationship.

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.

Genomic analysis of Ralstonia pickettii reveals the genetic features for potential pathogenicity and adaptive evolution in drinking water

Ralstonia pickettii, the most critical clinical pathogen of the genus Ralstonia, has been identified as a causative agent of numerous harmful infections. Additionally, Ralstonia pickettii demonstrates adaptability to extreme environmental conditions, such as those found in drinking water. In this study, we conducted a comprehensive genomic analysis to investigate the genomic characteristics related to potential pathogenicity and adaptive evolution in drinking water environments of Ralstonia pickettii. Through phylogenetic analysis and population genetic analysis, we divided Ralstonia pickettii into five Groups, two of which were associated with drinking water environments. The open pan-genome with a large and flexible gene repertoire indicated a high genetic plasticity. Significant differences in functional enrichment were observed between the core- and pan-genome of different groups. Diverse mobile genetic elements (MGEs), extensive genomic rearrangements, and horizontal gene transfer (HGT) events played a crucial role in generating genetic diversity. In drinking water environments, Ralstonia pickettii exhibited strong adaptability, and the acquisition of specific adaptive genes was potentially facilitated by genomic islands (GIs) and HGT. Furthermore, environmental pressures drove the adaptive evolution of Ralstonia pickettii, leading to the accumulation of unique mutations in key genes. These mutations may have a significant impact on various physiological functions, particularly carbon metabolism and energy metabolism. The presence of virulence-related elements associated with macromolecular secretion systems, virulence factors, and antimicrobial resistance indicated the potential pathogenicity of Ralstonia pickettii, making it capable of causing multiple nosocomial infections. This study provides comprehensive insights into the potential pathogenicity and adaptive evolution of Ralstonia pickettii in drinking water environments from a genomic perspective.

Changing fitness effects of mutations through long-term bacterial evolution

The distribution of fitness effects of new mutations shapes evolution, but it is challenging to observe how it changes as organisms adapt. Using Escherichia coli lineages spanning 50,000 generations of evolution, we quantify the fitness effects of insertion mutations in every gene. Macroscopically, the fraction of deleterious mutations changed little over time whereas the beneficial tail declined sharply, approaching an exponential distribution. Microscopically, changes in individual gene essentiality and deleterious effects often occurred in parallel; altered essentiality is only partly explained by structural variation. The identity and effect sizes of beneficial mutations changed rapidly over time, but many targets of selection remained predictable because of the importance of loss-of-function mutations. Taken together, these results reveal the dynamic-but statistically predictable-nature of mutational fitness effects.

A bacterial sensor taxonomy across earth ecosystems for machine learning applications

Microbial communities have evolved to colonize all ecosystems of the planet, from the deep sea to the human gut. Microbes survive by sensing, responding, and adapting to immediate environmental cues. This process is driven by signal transduction proteins such as histidine kinases, which use their sensing domains to bind or otherwise detect environmental cues and "transduce" signals to adjust internal processes. We hypothesized that an ecosystem's unique stimuli leave a sensor "fingerprint," able to identify and shed insight on ecosystem conditions. To test this, we collected 20,712 publicly available metagenomes from Host-associated, Environmental, and Engineered ecosystems across the globe. We extracted and clustered the collection's nearly 18M unique sensory domains into 113,712 similar groupings with MMseqs2. We built gradient-boosted decision tree machine learning models and found we could classify the ecosystem type (accuracy: 87%) and predict the levels of different physical parameters (R2 score: 83%) using the sensor cluster abundance as features. Feature importance enables identification of the most predictive sensors to differentiate between ecosystems which can lead to mechanistic interpretations if the sensor domains are well annotated. To demonstrate this, a machine learning model was trained to predict patient's disease state and used to identify domains related to oxygen sensing present in a healthy gut but missing in patients with abnormal conditions. Moreover, since 98.7% of identified sensor domains are uncharacterized, importance ranking can be used to prioritize sensors to determine what ecosystem function they may be sensing. Furthermore, these new predictive sensors can function as targets for novel sensor engineering with applications in biotechnology, ecosystem maintenance, and medicine.IMPORTANCEMicrobes infect, colonize, and proliferate due to their ability to sense and respond quickly to their surroundings. In this research, we extract the sensory proteins from a diverse range of environmental, engineered, and host-associated metagenomes. We trained machine learning classifiers using sensors as features such that it is possible to predict the ecosystem for a metagenome from its sensor profile. We use the optimized model's feature importance to identify the most impactful and predictive sensors in different environments. We next use the sensor profile from human gut metagenomes to classify their disease states and explore which sensors can explain differences between diseases. The sensors most predictive of environmental labels here, most of which correspond to uncharacterized proteins, are a useful starting point for the discovery of important environment signals and the development of possible diagnostic interventions.

Depth shapes microbiome assembly and network stability in the Mariana Trench

IMPORTANCE

Exploring microbial interactions and their stability/resilience from the surface to the hadal ocean is critical for further understanding of the microbiome structure and ecosystem function in the Mariana Trench. Vertical gradients did not destabilize microbial communities after long-term evolution and adaption. The uniform niche breadth, diversity, community complexity, and stability of microbiomes in both upper bathypelagic and hadal waters suggest the consistent roles of microbiomes in elemental cycling and adaptive strategies to overcome extreme environmental conditions. Compared with microeukaryotes, bacteria and archaea play a pivotal role in shaping the stability of the hadal microbiome. The consistent co-occurrence stability of microbiomes across vertical gradients was observed in the Mariana Trench. These results illuminate a key principle of microbiomes inhabiting the deepest trench: although distinct microbial communities occupy specific habitats, the interactions within microbial communities remain consistently stable from the upper bathypelagic to the hadal waters.

Community-scale models of microbiomes: Articulating metabolic modelling and metagenome sequencing

Building models is essential for understanding the functions and dynamics of microbial communities. Metabolic models built on genome-scale metabolic network reconstructions (GENREs) are especially relevant as a means to decipher the complex interactions occurring among species. Model reconstruction increasingly relies on metagenomics, which permits direct characterisation of naturally occurring communities that may contain organisms that cannot be isolated or cultured. In this review, we provide an overview of the field of metabolic modelling and its increasing reliance on and synergy with metagenomics and bioinformatics. We survey the means of assigning functions and reconstructing metabolic networks from (meta-)genomes, and present the variety and mathematical fundamentals of metabolic models that foster the understanding of microbial dynamics. We emphasise the characterisation of interactions and the scaling of model construction to large communities, two important bottlenecks in the applicability of these models. We give an overview of the current state of the art in metagenome sequencing and bioinformatics analysis, focusing on the reconstruction of genomes in microbial communities. Metagenomics benefits tremendously from third-generation sequencing, and we discuss the opportunities of long-read sequencing, strain-level characterisation and eukaryotic metagenomics. We aim at providing algorithmic and mathematical support, together with tool and application resources, that permit bridging the gap between metagenomics and metabolic modelling.

Microbial interactions strengthen deterministic processes during community assembly in a subtropical estuary

Systematic studies on the assembly process and driving mechanisms of microbial communities in estuaries with diverse seasonal and spatial scales are still limited. In this study, high-throughput sequencing, and microbial network analysis were combined to decipher the impact of environmental changes and biological interactions on the maintenance of microbial diversity patterns in the Jiulong River Estuary (JRE). The results showed that overall, stochastic processes dominated the bacterioplankton community assembly in the estuary, accounting for 49.66-74.78 % of the total. Additionally, bacterioplankton community diversity varied significantly across seasons and subzones. Specifically, the concentration of soluble reactive phosphorus (SRP) in the estuary steadily reduced from winter to summer, and the corresponding bacterioplankton community interactions gradually shifted from the weakest interaction in winter to the strongest in summer. The deterministic processes contributed more than half (50.34 %) to microbial assembly in the summer, but only 25.22 % in winter. Deterministic processes prevailed in the seaward with low SRP concentrations and strong bacterioplankton community interactions, while stochastic processes contributed 70.14 % to the assembly of microbial communities riverward. Biotic and abiotic factors, such as nutrients and microbial interactions, jointly drove the seasonal and spatial patterns of bacterioplankton community assembly, but overall, nutrients played a dominant role. Nevertheless, the contributions of nutrients and microbial interactions were equivalent in spatial assembly processes, albeit nutrients were the primary seasonal driver of the bacterioplankton community assembly process. This study emphasizes the significance of microbial interactions in the bacterioplankton community assemblage. These findings provide new and comprehensive insights into the microbial communities' organization in estuaries.

Mosaic environment-driven evolution of the deep-sea mussel Gigantidas platifrons bacterial endosymbiont

BACKGROUND

The within-species diversity of symbiotic bacteria represents an important genetic resource for their environmental adaptation, especially for horizontally transmitted endosymbionts. Although strain-level intraspecies variation has recently been detected in many deep-sea endosymbionts, their ecological role in environmental adaptation, their genome evolution pattern under heterogeneous geochemical environments, and the underlying molecular forces remain unclear.

RESULTS

Here, we conducted a fine-scale metagenomic analysis of the deep-sea mussel Gigantidas platifrons bacterial endosymbiont collected from distinct habitats: hydrothermal vent and methane seep. Endosymbiont genomes were assembled using a pipeline that distinguishes within-species variation and revealed highly heterogeneous compositions in mussels from different habitats. Phylogenetic analysis separated the assemblies into three distinct environment-linked clades. Their functional differentiation follows a mosaic evolutionary pattern. Core genes, essential for central metabolic function and symbiosis, were conserved across all clades. Clade-specific genes associated with heavy metal resistance, pH homeostasis, and nitrate utilization exhibited signals of accelerated evolution. Notably, transposable elements and plasmids contributed to the genetic reshuffling of the symbiont genomes and likely accelerated adaptive evolution through pseudogenization and the introduction of new genes.

CONCLUSIONS

The current study uncovers the environment-driven evolution of deep-sea symbionts mediated by mobile genetic elements. Its findings highlight a potentially common and critical role of within-species diversity in animal-microbiome symbioses. Video Abstract.

A genomic catalogue of soil microbiomes boosts mining of biodiversity and genetic resources

Soil harbors a vast expanse of unidentified microbes, termed as microbial dark matter, presenting an untapped reservo)ir of microbial biodiversity and genetic resources, but has yet to be fully explored. In this study, we conduct a large-scale excavation of soil microbial dark matter by reconstructing 40,039 metagenome-assembled genome bins (the SMAG catalogue) from 3304 soil metagenomes. We identify 16,530 of 21,077 species-level genome bins (SGBs) as unknown SGBs (uSGBs), which expand archaeal and bacterial diversity across the tree of life. We also illustrate the pivotal role of uSGBs in augmenting soil microbiome's functional landscape and intra-species genome diversity, providing large proportions of the 43,169 biosynthetic gene clusters and 8545 CRISPR-Cas genes. Additionally, we determine that uSGBs contributed 84.6% of previously unexplored viral-host associations from the SMAG catalogue. The SMAG catalogue provides an useful genomic resource for further studies investigating soil microbial biodiversity and genetic resources.

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.

Comparison of gene clustering criteria reveals intrinsic uncertainty in pangenome analyses

BACKGROUND

A key step for comparative genomics is to group open reading frames into functionally and evolutionarily meaningful gene clusters. Gene clustering is complicated by intraspecific duplications and horizontal gene transfers that are frequent in prokaryotes. In consequence, gene clustering methods must deal with a trade-off between identifying vertically transmitted representatives of multicopy gene families, which are recognizable by synteny conservation, and retrieving complete sets of species-level orthologs. We studied the implications of adopting homology, orthology, or synteny conservation as formal criteria for gene clustering by performing comparative analyses of 125 prokaryotic pangenomes.

RESULTS

Clustering criteria affect pangenome functional characterization, core genome inference, and reconstruction of ancestral gene content to different extents. Species-wise estimates of pangenome and core genome sizes change by the same factor when using different clustering criteria, allowing robust cross-species comparisons regardless of the clustering criterion. However, cross-species comparisons of genome plasticity and functional profiles are substantially affected by inconsistencies among clustering criteria. Such inconsistencies are driven not only by mobile genetic elements, but also by genes involved in defense, secondary metabolism, and other accessory functions. In some pangenome features, the variability attributed to methodological inconsistencies can even exceed the effect sizes of ecological and phylogenetic variables.

CONCLUSIONS

Choosing an appropriate criterion for gene clustering is critical to conduct unbiased pangenome analyses. We provide practical guidelines to choose the right method depending on the research goals and the quality of genome assemblies, and a benchmarking dataset to assess the robustness and reproducibility of future comparative studies.

Origins and Evolution of Novel Bacteroides in Captive Apes

Bacterial strains evolve in response to the gut environment of their hosts, with genomic changes that influence their interactions with hosts as well as with other members of the gut community. Great apes in captivity have acquired strains of Bacteroides xylanisolvens, which are common within gut microbiome of humans but not typically found other apes, thereby enabling characterization of strain evolution following colonization. Here, we isolate, sequence and reconstruct the history of gene gain and loss events in numerous captive-ape-associated strains since their divergence from their closest human-associated strains. We show that multiple captive-ape-associated B. xylanisolvens lineages have independently acquired gene complexes that encode functions related to host mucin metabolism. Our results support the finding of high genome fluidity in Bacteroides, in that several strains, in moving from humans to captive apes, have rapidly gained large genomic regions that augment metabolic properties not previously present in their relatives.

Insights into early evolutionary adaptations of the Akkermansia genus to the vertebrate gut

Akkermansia, a relevant mucin degrader from the vertebrate gut microbiota, is a member of the deeply branched Verrucomicrobiota, as well as the only known member of this phylum to be described as inhabitants of the gut. Only a few Akkermansia species have been officially described so far, although there is genomic evidence addressing the existence of more species-level variants for this genus. This niche specialization makes Akkermansia an interesting model for studying the evolution of microorganisms to their adaptation to the gastrointestinal tract environment, including which kind of functions were gained when the Akkermansia genus originated or how the evolutionary pressure functions over those genes. In order to gain more insight into Akkermansia adaptations to the gastrointestinal tract niche, we performed a phylogenomic analysis of 367 high-quality Akkermansia isolates and metagenome-assembled genomes, in addition to other members of Verrucomicrobiota. This work was focused on three aspects: the definition of Akkermansia genomic species clusters and the calculation and functional characterization of the pangenome for the most represented species; the evolutionary relationship between Akkermansia and their closest relatives from Verrucomicrobiota, defining the gene families which were gained or lost during the emergence of the last Akkermansia common ancestor (LAkkCA) and; the evaluation of the evolutionary pressure metrics for each relevant gene family of main Akkermansia species. This analysis found 25 Akkermansia genomic species clusters distributed in two main clades, divergent from their non-Akkermansia relatives. Pangenome analyses suggest that Akkermansia species have open pangenomes, and the gene gain/loss model indicates that genes associated with mucin degradation (both glycoside hydrolases and peptidases), (micro)aerobic metabolism, surface interaction, and adhesion were part of LAkkCA. Specifically, mucin degradation is a very ancestral innovation involved in the origin of Akkermansia. Horizontal gene transfer detection suggests that Akkermansia could receive genes mostly from unknown sources or from other Gram-negative gut bacteria. Evolutionary metrics suggest that Akkemansia species evolved differently, and even some conserved genes suffered different evolutionary pressures among clades. These results suggest a complex evolutionary landscape of the genus and indicate that mucin degradation could be an essential feature in Akkermansia evolution as a symbiotic species.

Human Pangenomics: Promises and Challenges of a Distributed Genomic Reference

A pangenome is a collection of the common and unique genomes that are present in a given species. It combines the genetic information of all the genomes sampled, resulting in a large and diverse range of genetic material. Pangenomic analysis offers several advantages compared to traditional genomic research. For example, a pangenome is not bound by the physical constraints of a single genome, so it can capture more genetic variability. Thanks to the introduction of the concept of pangenome, it is possible to use exceedingly detailed sequence data to study the evolutionary history of two different species, or how populations within a species differ genetically. In the wake of the Human Pangenome Project, this review aims at discussing the advantages of the pangenome around human genetic variation, which are then framed around how pangenomic data can inform population genetics, phylogenetics, and public health policy by providing insights into the genetic basis of diseases or determining personalized treatments, targeting the specific genetic profile of an individual. Moreover, technical limitations, ethical concerns, and legal considerations are discussed.

Comparative genomics and phylogenomics of Campylobacter unveil potential novel species and provide insights into niche segregation

Campylobacter is a bacterial genus associated with community outbreaks and gastrointestinal symptoms. Studies on Campylobacter generally focus on specific pathogenic species such as C. coli and C. jejuni. Currently, there are thousands of publicly available Campylobacter genomes, allowing a more complete assessment of the genus diversity. In this work, we report a network-based analysis of all available Campylobacter genomes to explore the genus structure and diversity, revealing potentially new species and elucidating genus features. We also hypothesize that the previously established Clade III of C. coli is in fact a novel species (referred here as Campylobacter spp12). Finally, we found a negative correlation between pangenome fluidity and saturation coefficient, with potential implications to the lifestyles of distinct Campylobacter species. Since pangenome analysis depends on the number of available genomes, this correlation could help estimate pangenome metrics of Campylobacter species with less sequenced genomes, helping understand their lifestyle and niche adaptation. Together, our results indicate that the Campylobacter genus should be re-evaluated, with particular attention to the interplay between genome structure and niche segregation.

A social niche breadth score reveals niche range strategies of generalists and specialists

Generalists can survive in many environments, whereas specialists are restricted to a single environment. Although a classical concept in ecology, niche breadth has remained challenging to quantify for microorganisms because it depends on an objective definition of the environment. Here, by defining the environment of a microorganism as the community it resides in, we integrated information from over 22,000 environmental sequencing samples to derive a quantitative measure of the niche, which we call social niche breadth. At the level of genera, we explored niche range strategies throughout the prokaryotic tree of life. We found that social generalists include opportunists that stochastically dominate local communities, whereas social specialists are stable but low in abundance. Social generalists have a more diverse and open pan-genome than social specialists, but we found no global correlation between social niche breadth and genome size. Instead, we observed two distinct evolutionary strategies, whereby specialists have relatively small genomes in habitats with low local diversity, but relatively large genomes in habitats with high local diversity. Together, our analysis shines data-driven light on microbial niche range strategies.

Linking prokaryotic genome size variation to metabolic potential and environment

While theories and models have appeared to explain genome size as a result of evolutionary processes, little work has shown that genome sizes carry ecological signatures. Our work delves into the ecological implications of microbial genome size variation in benthic and pelagic habitats across environmental gradients of the brackish Baltic Sea. While depth is significantly associated with genome size in benthic and pelagic brackish metagenomes, salinity is only correlated to genome size in benthic metagenomes. Overall, we confirm that prokaryotic genome sizes in Baltic sediments (3.47 Mbp) are significantly bigger than in the water column (2.96 Mbp). While benthic genomes have a higher number of functions than pelagic genomes, the smallest genomes coded for a higher number of module steps per Mbp for most of the functions irrespective of their environment. Some examples of this functions are amino acid metabolism and central carbohydrate metabolism. However, we observed that nitrogen metabolism was almost absent in pelagic genomes and was mostly present in benthic genomes. Finally, we also show that Bacteria inhabiting Baltic sediments and water column not only differ in taxonomy, but also in their metabolic potential, such as the Wood-Ljungdahl pathway or the presence of different hydrogenases. Our work shows how microbial genome size is linked to abiotic factors in the environment, metabolic potential and taxonomic identity of Bacteria and Archaea within aquatic ecosystems.

About the dark corners in the gene function space of Escherichia coli remaining without illumination by scientific literature

BACKGROUND

Although Escherichia coli (E. coli) is the most studied prokaryote organism in the history of life sciences, many molecular mechanisms and gene functions encoded in its genome remain to be discovered. This work aims at quantifying the illumination of the E. coli gene function space by the scientific literature and how close we are towards the goal of a complete list of E. coli gene functions.

RESULTS

The scientific literature about E. coli protein-coding genes has been mapped onto the genome via the mentioning of names for genomic regions in scientific articles both for the case of the strain K-12 MG1655 as well as for the 95%-threshold softcore genome of 1324 E. coli strains with known complete genome. The article match was quantified with the ratio of a given gene name's occurrence to the mentioning of any gene names in the paper. The various genome regions have an extremely uneven literature coverage. A group of elite genes with ≥ 100 full publication equivalents (FPEs, FPE = 1 is an idealized publication devoted to just a single gene) attracts the lion share of the papers. For K-12, ~ 65% of the literature covers just 342 elite genes; for the softcore genome, ~ 68% of the FPEs is about only 342 elite gene families (GFs). We also find that most genes/GFs have at least one mentioning in a dedicated scientific article (with the exception of at least 137 protein-coding transcripts for K-12 and 26 GFs from the softcore genome). Whereas the literature growth rates were highest for uncharacterized or understudied genes until 2005-2010 compared with other groups of genes, they became negative thereafter. At the same time, literature for anyhow well-studied genes started to grow explosively with threshold T10 (≥ 10 FPEs). Typically, a body of ~ 20 actual articles generated over ~ 15 years of research effort was necessary to reach T10. Lineage-specific co-occurrence analysis of genes belonging to the accessory genome of E. coli together with genomic co-localization and sequence-analytic exploration hints previously completely uncharacterized genes yahV and yddL being associated with osmotic stress response/motility mechanisms.

CONCLUSION

If the numbers of scientific articles about uncharacterized and understudied genes remain at least at present levels, full gene function lists for the strain K-12 MG1655 and the E. coli softcore genome are in reach within the next 25-30 years. Once the literature body for a gene crosses 10 FPEs, most of the critical fundamental research risk appears overcome and steady incremental research becomes possible.

Elevated temperature and CO2 strongly affect the growth strategies of soil bacteria

The trait-based strategies of microorganisms appear to be phylogenetically conserved, but acclimation to climate change may complicate the scenario. To study the roles of phylogeny and environment on bacterial responses to sudden moisture increases, we determine bacterial population-specific growth rates by ¹⁸O-DNA quantitative stable isotope probing (¹⁸O-qSIP) in soils subjected to a free-air CO₂ enrichment (FACE) combined with warming. We find that three growth strategies of bacterial taxa - rapid, intermediate and slow responders, defined by the timing of the peak growth rates - are phylogenetically conserved, even at the sub-phylum level. For example, members of class Bacilli and Sphingobacteriia are mainly rapid responders. Climate regimes, however, modify the growth strategies of over 90% of species, partly confounding the initial phylogenetic pattern. The growth of rapid bacterial responders is more influenced by phylogeny, whereas the variance for slow responders is primarily explained by environmental conditions. Overall, these results highlight the role of phylogenetic and environmental constraints in understanding and predicting the growth strategies of soil microorganisms under global change scenarios.

Machine learning enables prediction of metabolic system evolution in bacteria

Evolution prediction is a long-standing goal in evolutionary biology, with potential impacts on strategic pathogen control, genome engineering, and synthetic biology. While laboratory evolution studies have shown the predictability of short-term and sequence-level evolution, that of long-term and system-level evolution has not been systematically examined. Here, we show that the gene content evolution of metabolic systems is generally predictable by applying ancestral gene content reconstruction and machine learning techniques to ~3000 bacterial genomes. Our framework, Evodictor, successfully predicted gene gain and loss evolution at the branches of the reference phylogenetic tree, suggesting that evolutionary pressures and constraints on metabolic systems are universally shared. Investigation of pathway architectures and meta-analysis of metagenomic datasets confirmed that these evolutionary patterns have physiological and ecological bases as functional dependencies among metabolic reactions and bacterial habitat changes. Last, pan-genomic analysis of intraspecies gene content variations proved that even "ongoing" evolution in extant bacterial species is predictable in our framework.

proGenomes3: approaching one million accurately and consistently annotated high-quality prokaryotic genomes

The interpretation of genomic, transcriptomic and other microbial 'omics data is highly dependent on the availability of well-annotated genomes. As the number of publicly available microbial genomes continues to increase exponentially, the need for quality control and consistent annotation is becoming critical. We present proGenomes3, a database of 907 388 high-quality genomes containing 4 billion genes that passed stringent criteria and have been consistently annotated using multiple functional and taxonomic databases including mobile genetic elements and biosynthetic gene clusters. proGenomes3 encompasses 41 171 species-level clusters, defined based on universal single copy marker genes, for which pan-genomes and contextual habitat annotations are provided. The database is available at http://progenomes.embl.de/.

Comparison of prokaryotes between Mount Everest and the Mariana Trench

BACKGROUND

Mount Everest and the Mariana Trench represent the highest and deepest places on Earth, respectively. They are geographically separated, with distinct extreme environmental parameters that provide unique habitats for prokaryotes. Comparison of prokaryotes between Mount Everest and the Mariana Trench will provide a unique perspective to understanding the composition and distribution of environmental microbiomes on Earth.

RESULTS

Here, we compared prokaryotic communities between Mount Everest and the Mariana Trench based on shotgun metagenomic analysis. Analyzing 25 metagenomes and 1176 metagenome-assembled genomes showed distinct taxonomic compositions between Mount Everest and the Mariana Trench, with little taxa overlap, and significant differences in genome size, GC content, and predicted optimal growth temperature. However, community metabolic capabilities exhibited striking commonality, with > 90% of metabolic modules overlapping among samples of Mount Everest and the Mariana Trench, with the only exception for CO2 fixations (photoautotrophy in Mount Everest but chemoautotrophy in the Mariana Trench). Most metabolic pathways were common but performed by distinct taxa in the two extreme habitats, even including some specialized metabolic pathways, such as the versatile degradation of various refractory organic matters, heavy metal metabolism (e.g., As and Se), stress resistance, and antioxidation. The metabolic commonality indicated the overall consistent roles of prokaryotes in elemental cycling and common adaptation strategies to overcome the distinct stress conditions despite the intuitively huge differences in Mount Everest and the Mariana Trench.

CONCLUSION

Our results, the first comparison between prokaryotes in the highest and the deepest habitats on Earth, may highlight the principles of prokaryotic diversity: although taxa are habitat-specific, primary metabolic functions could be always conserved. Video abstract.

Alternative nitrogenase of Paenibacillus sonchi genomovar Riograndensis: An insight in the origin of Fe-nitrogenase in the Paenibacillaceae family

Paenibacillus sonchi genomovar Riograndensis is a nitrogen-fixing bacteria isolated from wheat that displays diverse plant growth-promoting abilities. Beyond conventional Mo-nitrogenase, this organism also harbors an alternative Fe-nitrogenase, whose many aspects related to regulation, physiology, and evolution remain to be elucidated. In this work, the origins of this alternative system were investigated, exploring the distribution and diversification of nitrogenases in the Panibacillaceae family. Our analysis showed that diazotrophs represent 17% of Paenibacillaceae genomes, of these, only 14.4% (2.5% of all Paenibacillaceae genomes) also contained Fe or V- nitrogenases. Diverse nif-like sequences were also described, occurring mainly in genomes that also harbor the alternative systems. The analysis of genomes containing Fe-nitrogenase showed a conserved cluster of nifEN anfHDGK across three genera: Gorillibacterium, Fontibacillus, and Paenibacillus. A phylogeny of anfHDGK separated the Fe-nitrogenases into three main groups. Our analysis suggested that Fe-nitrogenase was acquired by the ancestral lineage of Fontibacillus, Gorillibacterium, and Paenibacillus genera via horizontal gene transfer (HGT), and further events of transfer and gene loss marked the evolution of this alternative nitrogenase in these groups. The species phylogeny of N-fixing Paenibacillaceae separated the diazotrophs into five clades, one of these containing all occurrences of strains harboring alternative nitrogenases in the Paenibacillus genus. The pangenome of this clade is open and composed of more than 96% of accessory genes. Diverse functional categories were enriched in the flexible genome, including functions related to replication and repair. The latter involved diverse genes related to HGT, suggesting that such events may have an important role in the evolution of diazotrophic Paenibacillus. This study provided an insight into the organization, distribution, and evolution of alternative nitrogenase genes in Paenibacillaceae, considering different genomic aspects.

The newest COI molecular detection of Asian redtail catfish Hemibagrus nemurus (Valenciennes, 1840) in Progo River, Magelang, Central Java, Indonesia

Objective

This study describes the newest deoxyribonucleic acid (DNA) barcoding Asian redtail catfish (Hemibagrus nemurus) in the Progo River, Magelang, Central Java, Indonesia.

Materials and Methods

Ten fish were caught in the Progo River, Magelang, Central Java, Indonesia. The polymerase chain reaction was the molecular diagnosis to detect the sequences of DNA of Cytochrome Oxidase I compared to National Center for Biotechnology Information data (GenBank).

Results

The results showed that the percent identity was not 100% with H. nemurus data from other locations (GenBank), including Indonesia. The closest percent identity was H. nemurus from Java Island (Accession ID: MK312566.1) with 97.6% similarity. The genetic mutation that happened might be due to environmental change (pollution) in the Progo River recently.

Conclusions

This study showed a genetic mutation in H. nemurus from Progo River may be caused by environmental change. Low pollution exposure levels may not be detrimental (lethal) to fish. However, it can affect fish fertility, which leads to population degradation (gene variation). Attention must be increased for fish survival in the future.

Distinct evolutionary trajectories in the Escherichia coli pangenome occur within sequence types

The Escherichia coli species contains a diverse set of sequence types and there remain important questions regarding differences in genetic content within this population that need to be addressed. Pangenomes are useful vehicles for studying gene content within sequence types. Here, we analyse 21 E. coli sequence type pangenomes using comparative pangenomics to identify variance in both pangenome structure and content. We present functional breakdowns of sequence type core genomes and identify sequence types that are enriched in metabolism, transcription and cell membrane biogenesis genes. We also uncover metabolism genes that have variable core classification, depending on which allele is present. Our comparative pangenomics approach allows for detailed exploration of sequence type pangenomes within the context of the species. We show that ongoing gene gain and loss in the E. coli pangenome is sequence type-specific, which may be a consequence of distinct sequence type-specific evolutionary drivers.

Nutrition or nature: using elementary flux modes to disentangle the complex forces shaping prokaryote pan-genomes

BACKGROUND

Microbial pan-genomes are shaped by a complex combination of stochastic and deterministic forces. Even closely related genomes exhibit extensive variation in their gene content. Understanding what drives this variation requires exploring the interactions of gene products with each other and with the organism's external environment. However, to date, conceptual models of pan-genome dynamics often represent genes as independent units and provide limited information about their mechanistic interactions.

RESULTS

We simulated the stochastic process of gene-loss using the pooled genome-scale metabolic reaction networks of 46 taxonomically diverse bacterial and archaeal families as proxies for their pan-genomes. The frequency by which reactions are retained in functional networks when stochastic gene loss is simulated in diverse environments allowed us to disentangle the metabolic reactions whose presence depends on the metabolite composition of the external environment (constrained by "nutrition") from those that are independent of the environment (constrained by "nature"). By comparing the frequency of reactions from the first group with their observed frequencies in bacterial and archaeal families, we predicted the metabolic niches that shaped the genomic composition of these lineages. Moreover, we found that the lineages that were shaped by a more diverse metabolic niche also occur in more diverse biomes as assessed by global environmental sequencing datasets.

CONCLUSION

We introduce a computational framework for analyzing and interpreting pan-reactomes that provides novel insights into the ecological and evolutionary drivers of pan-genome dynamics.

Comparative Genomics Reveal the High Conservation and Scarce Distribution of Nitrogen Fixation nif Genes in the Plant-Associated Genus Herbaspirillum

The genus Herbaspirillum gained the spotlight due to the several reports of diazotrophic strains and promising results in plant-growth field assays. However, as diversity exploration of Herbaspirillum species gained momentum, it became clearer that the plant beneficial lifestyle was not the only form of ecological interaction in this genus, due to reports of phytopathogenesis and nosocomial infections. Here we performed a deep search across all publicly available Herbaspirillum genomes. Using a robust core genome phylogeny, we have found that all described species are well delineated, being the only exception H. aquaticum and H. huttiense clade. We also uncovered that the nif genes are only highly prevalent in H. rubrisubalbicans; however, irrespective to the species, all nif genes share the same gene arrangement with high protein identity, and are present in only two main types, in inverted strands. By means of a NifHDKENB phylogenetic tree, we have further revealed that the Herbaspirillum nif sequences may have been acquired from the same last common ancestor belonging to the Nitrosomonadales order.

Pangenome Evolution in Environmentally Transmitted Symbionts of Deep-Sea Mussels Is Governed by Vertical Inheritance

Microbial pangenomes vary across species; their size and structure are determined by genetic diversity within the population and by gene loss and horizontal gene transfer (HGT). Many bacteria are associated with eukaryotic hosts where the host colonization dynamics may impact bacterial genome evolution. Host-associated lifestyle has been recognized as a barrier to HGT in parentally transmitted bacteria. However, pangenome evolution of environmentally acquired symbionts remains understudied, often due to limitations in symbiont cultivation. Using high-resolution metagenomics, here we study pangenome evolution of two co-occurring endosymbionts inhabiting Bathymodiolus brooksi mussels from a single cold seep. The symbionts, sulfur-oxidizing (SOX) and methane-oxidizing (MOX) gamma-proteobacteria, are environmentally acquired at an early developmental stage and individual mussels may harbor multiple strains of each symbiont species. We found differences in the accessory gene content of both symbionts across individual mussels, which are reflected by differences in symbiont strain composition. Compared with core genes, accessory genes are enriched in genome plasticity functions. We found no evidence for recent HGT between both symbionts. A comparison between the symbiont pangenomes revealed that the MOX population is less diverged and contains fewer accessory genes, supporting that the MOX association with B. brooksi is more recent in comparison to that of SOX. Our results show that the pangenomes of both symbionts evolved mainly by vertical inheritance. We conclude that genome evolution of environmentally transmitted symbionts that associate with individual hosts over their lifetime is affected by a narrow symbiosis where the frequency of HGT is constrained.

To kill or to be killed: pangenome analysis of Escherichia coli strains reveals a tailocin specific for pandemic ST131

BACKGROUND

Escherichia coli (E. coli) has been one of the most studied model organisms in the history of life sciences. Initially thought just to be commensal bacteria, E. coli has shown wide phenotypic diversity including pathogenic isolates with great relevance to public health. Though pangenome analysis has been attempted several times, there is no systematic functional characterization of the E. coli subgroups according to the gene profile.

RESULTS

Systematically scanning for optimal parametrization, we have built the E. coli pangenome from 1324 complete genomes. The pangenome size is estimated to be ~25,000 gene families (GFs). Whereas the core genome diminishes as more genomes are added, the softcore genome (≥95% of strains) is stable with ~3000 GFs regardless of the total number of genomes. Apparently, the softcore genome (with a 92% or 95% generation threshold) can define the genome of a bacterial species listing the critically relevant, evolutionarily most conserved or important classes of GFs. Unsupervised clustering of common E. coli sequence types using the presence/absence GF matrix reveals distinct characteristics of E. coli phylogroups B1, B2, and E. We highlight the bi-lineage nature of B1, the variation of the secretion and of the iron acquisition systems in ST11 (E), and the incorporation of a highly conserved prophage into the genome of ST131 (B2). The tail structure of the prophage is evolutionarily related to R2-pyocin (a tailocin) from Pseudomonas aeruginosa PAO1. We hypothesize that this molecular machinery is highly likely to play an important role in protecting its own colonies; thus, contributing towards the rapid rise of pandemic E. coli ST131.

CONCLUSIONS

This study has explored the optimized pangenome development in E. coli. We provide complete GF lists and the pangenome matrix as supplementary data for further studies. We identified biological characteristics of different E. coli subtypes, specifically for phylogroups B1, B2, and E. We found an operon-like genome region coding for a tailocin specific for ST131 strains. The latter is a potential killer weapon providing pandemic E. coli ST131 with an advantage in inter-bacterial competition and, suggestively, explains their dominance as human pathogen among E. coli strains.

Pangenome analyses of Bacillus pumilus, Bacillus safensis, and Priestia megaterium exploring the plant-associated features of bacilli strains isolated from canola

Previous genome mining of the strains Bacillus pumilus 7PB, Bacillus safensis 1TAz, 8Taz, and 32PB, and Priestia megaterium 16PB isolated from canola revealed differences in the profile of antimicrobial biosynthetic genes when compared to the species type strains. To evaluate not only the similarities among B. pumilus, B. safensis, and P. megaterium genomes but also the specificities found in the canola bacilli, we performed comparative genomic analyses through the pangenome evaluation of each species. Besides that, other genome features were explored, especially focusing on plant-associated and biotechnological characteristics. The combination of the genome metrics Average Nucleotide Identity and digital DNA-DNA hybridization formulas 1 and 3 adopting the universal thresholds of 95 and 70%, respectively, was suitable to verify the identification of strains from these groups. On average, core genes corresponded to 45%, 52%, and 34% of B. pumilus, B. safensis, and P. megaterium open pangenomes, respectively. Many genes related to adaptations to plant-associated lifestyles were predicted, especially in the Bacillus genomes. These included genes for acetoin production, polyamines utilization, root exudate chemoreceptors, biofilm formation, and plant cell-wall degrading enzymes. Overall, we could observe that strains of these species exhibit many features in common, whereas most of their variable genome portions have features yet to be uncovered. The observed antifungal activity of canola bacilli might be a result of the synergistic action of secondary metabolites, siderophores, and chitinases. Genome analysis confirmed that these species and strains have biotechnological potential to be used both as agricultural inoculants or hydrolases producers. Up to our knowledge, this is the first work that evaluates the pangenome features of P. megaterium.

Giant Starship Elements Mobilize Accessory Genes in Fungal Genomes

Accessory genes are variably present among members of a species and are a reservoir of adaptive functions. In bacteria, differences in gene distributions among individuals largely result from mobile elements that acquire and disperse accessory genes as cargo. In contrast, the impact of cargo-carrying elements on eukaryotic evolution remains largely unknown. Here, we show that variation in genome content within multiple fungal species is facilitated by Starships, a newly discovered group of massive mobile elements that are 110 kb long on average, share conserved components, and carry diverse arrays of accessory genes. We identified hundreds of Starship-like regions across every major class of filamentous Ascomycetes, including 28 distinct Starships that range from 27 to 393 kb and last shared a common ancestor ca. 400 Ma. Using new long-read assemblies of the plant pathogen Macrophomina phaseolina, we characterize four additional Starships whose activities contribute to standing variation in genome structure and content. One of these elements, Voyager, inserts into 5S rDNA and contains a candidate virulence factor whose increasing copy number has contrasting associations with pathogenic and saprophytic growth, suggesting Voyager's activity underlies an ecological trade-off. We propose that Starships are eukaryotic analogs of bacterial integrative and conjugative elements based on parallels between their conserved components and may therefore represent the first dedicated agents of active gene transfer in eukaryotes. Our results suggest that Starships have shaped the content and structure of fungal genomes for millions of years and reveal a new concerted route for evolution throughout an entire eukaryotic phylum.

Toward quantifying the adaptive role of bacterial pangenomes during environmental perturbations

Metagenomic surveys have revealed that natural microbial communities are predominantly composed of sequence-discrete, species-like populations but the genetic and/or ecological processes that maintain such populations remain speculative, limiting our understanding of population speciation and adaptation to perturbations. To address this knowledge gap, we sequenced 112 Salinibacter ruber isolates and 12 companion metagenomes from four adjacent saltern ponds in Mallorca, Spain that were experimentally manipulated to dramatically alter salinity and light intensity, the two major drivers of this ecosystem. Our analyses showed that the pangenome of the local Sal. ruber population is open and similar in size (~15,000 genes) to that of randomly sampled Escherichia coli genomes. While most of the accessory (noncore) genes were isolate-specific and showed low in situ abundances based on the metagenomes compared to the core genes, indicating that they were functionally unimportant and/or transient, 3.5% of them became abundant when salinity (but not light) conditions changed and encoded for functions related to osmoregulation. Nonetheless, the ecological advantage of these genes, while significant, was apparently not strong enough to purge diversity within the population. Collectively, our results provide an explanation for how this immense intrapopulation gene diversity is maintained, which has implications for the prokaryotic species concept.

A comparative whole-genome approach identifies bacterial traits for marine microbial interactions

Microbial interactions shape the structure and function of microbial communities with profound consequences for biogeochemical cycles and ecosystem health. Yet, most interaction mechanisms are studied only in model systems and their prevalence is unknown. To systematically explore the functional and interaction potential of sequenced marine bacteria, we developed a trait-based approach, and applied it to 473 complete genomes (248 genera), representing a substantial fraction of marine microbial communities. We identified genome functional clusters (GFCs) which group bacterial taxa with common ecology and life history. Most GFCs revealed unique combinations of interaction traits, including the production of siderophores (10% of genomes), phytohormones (3-8%) and different B vitamins (57-70%). Specific GFCs, comprising Alpha- and Gammaproteobacteria, displayed more interaction traits than expected by chance, and are thus predicted to preferentially interact synergistically and/or antagonistically with bacteria and phytoplankton. Linked trait clusters (LTCs) identify traits that may have evolved to act together (e.g., secretion systems, nitrogen metabolism regulation and B vitamin transporters), providing testable hypotheses for complex mechanisms of microbial interactions. Our approach translates multidimensional genomic information into an atlas of marine bacteria and their putative functions, relevant for understanding the fundamental rules that govern community assembly and dynamics.

Integrative Assessments on Molecular Taxonomy of Acidiferrobacter thiooxydans ZJ and Its Environmental Adaptation Based on Mobile Genetic Elements

Acidiferrobacter spp. are facultatively anaerobic acidophiles that belong to a distinctive Acidiferrobacteraceae family, which are similar to Ectothiorhodospiraceae phylogenetically, and are closely related to Acidithiobacillia class/subdivision physiologically. The limited genome information has kept them from being studied on molecular taxonomy and environmental adaptation in depth. Herein, Af. thiooxydans ZJ was isolated from acid mine drainage (AMD), and the complete genome sequence was reported to scan its genetic constitution for taxonomic and adaptative feature exploration. The genome has a single chromosome of 3,302,271 base pairs (bp), with a GC content of 63.61%. The phylogenetic tree based on OrthoANI highlighted the unique position of Af. thiooxydans ZJ, which harbored more unique genes among the strains from Ectothiorhodospiraceae and Acidithiobacillaceae by pan-genome analysis. The diverse mobile genetic elements (MGEs), such as insertion sequence (IS), clustered regularly interspaced short palindromic repeat (CRISPR), prophage, and genomic island (GI), have been identified and characterized in Af. thiooxydans ZJ. The results showed that Af. thiooxydans ZJ may effectively resist the infection of foreign viruses and gain functional gene fragments or clusters to shape its own genome advantageously. This study will offer more evidence of the genomic plasticity and improve our understanding of evolutionary adaptation mechanisms to extreme AMD environment, which could expand the potential utilization of Af. thiooxydans ZJ as an iron and sulfur oxidizer in industrial bioleaching.

Pangenomes as a Resource to Accelerate Breeding of Under-Utilised Crop Species

Pangenomes are a rich resource to examine the genomic variation observed within a species or genera, supporting population genetics studies, with applications for the improvement of crop traits. Major crop species such as maize (Zea mays), rice (Oryza sativa), Brassica (Brassica spp.), and soybean (Glycine max) have had pangenomes constructed and released, and this has led to the discovery of valuable genes associated with disease resistance and yield components. However, pangenome data are not available for many less prominent crop species that are currently under-utilised. Despite many under-utilised species being important food sources in regional populations, the scarcity of genomic data for these species hinders their improvement. Here, we assess several under-utilised crops and review the pangenome approaches that could be used to build resources for their improvement. Many of these under-utilised crops are cultivated in arid or semi-arid environments, suggesting that novel genes related to drought tolerance may be identified and used for introgression into related major crop species. In addition, we discuss how previously collected data could be used to enrich pangenome functional analysis in genome-wide association studies (GWAS) based on studies in major crops. Considering the technological advances in genome sequencing, pangenome references for under-utilised species are becoming more obtainable, offering the opportunity to identify novel genes related to agro-morphological traits in these species.

Experimental Evolution Expands the Breadth of Adaptation to an Environmental Gradient Correlated With Genome Reduction

Whether and how adaptive evolution adjusts the breadth of adaptation in coordination with the genome are essential issues for connecting evolution with ecology. To address these questions, experimental evolution in five Escherichia coli strains carrying either the wild-type genome or a reduced genome was performed in a defined minimal medium (C0). The ancestral and evolved populations were subsequently subjected to fitness and chemical niche analyses across an environmental gradient with 29 combinations of eight chemical components of the minimal medium. The results showed that adaptation was achieved not only specific to the evolutionary condition (C0), but also generally, to the environmental gradient; that is, the breadth of adaptation to the eight chemical niches was expanded. The magnitudes of the adaptive improvement and the breadth increase were both correlated with genome reduction and were highly significant in two out of eight niches (i.e., glucose and sulfate). The direct adaptation-induced correlated adaptation to the environmental gradient was determined by only a few genome mutations. An additive increase in fitness associated with the stepwise fixation of mutations was consistently observed in the reduced genomes. In summary, this preliminary survey demonstrated that evolution finely tuned the breadth of adaptation correlated with genome reduction.

Prokaryote pangenomes are dynamic entities

Prokaryote pangenomes are influenced heavily by environmental factors and the opportunity for gene gain and loss events. As the field of pangenome analysis has expanded, so has the need to fully understand the complexity of how eco-evolutionary dynamics shape pangenomes. Here, we describe current models of pangenome evolution and discuss their suitability and accuracy. We suggest that pangenomes are dynamic entities under constant flux, highlighting the influence of two-way interactions between pangenome and environment. New classifications of core and accessory genes are also considered, underscoring the need for continuous evaluation of nomenclature in a fast-moving field. We conclude that future models of pangenome evolution should incorporate eco-evolutionary dynamics to fully encompass their dynamic, changeable nature.

Comparative pangenomics: analysis of 12 microbial pathogen pangenomes reveals conserved global structures of genetic and functional diversity

BACKGROUND

With the exponential growth of publicly available genome sequences, pangenome analyses have provided increasingly complete pictures of genetic diversity for many microbial species. However, relatively few studies have scaled beyond single pangenomes to compare global genetic diversity both within and across different species. We present here several methods for "comparative pangenomics" that can be used to contextualize multi-pangenome scale genetic diversity with gene function for multiple species at multiple resolutions: pangenome shape, genes, sequence variants, and positions within variants.

RESULTS

Applied to 12,676 genomes across 12 microbial pathogenic species, we observed several shared resolution-specific patterns of genetic diversity: First, pangenome openness is associated with species' phylogenetic placement. Second, relationships between gene function and frequency are conserved across species, with core genomes enriched for metabolic and ribosomal genes and accessory genomes for trafficking, secretion, and defense-associated genes. Third, genes in core genomes with the highest sequence diversity are functionally diverse. Finally, certain protein domains are consistently mutation enriched across multiple species, especially among aminoacyl-tRNA synthetases where the extent of a domain's mutation enrichment is strongly function-dependent.

CONCLUSIONS

These results illustrate the value of each resolution at uncovering distinct aspects in the relationship between genetic and functional diversity across multiple species. With the continued growth of the number of sequenced genomes, these methods will reveal additional universal patterns of genetic diversity at the pangenome scale.

Towards the biogeography of prokaryotic genes

Microbial genes encode the majority of the functional repertoire of life on earth. However, despite increasing efforts in metagenomic sequencing of various habitats^1-3, little is known about the distribution of genes across the global biosphere, with implications for human and planetary health. Here we constructed a non-redundant gene catalogue of 303 million species-level genes (clustered at 95% nucleotide identity) from 13,174 publicly available metagenomes across 14 major habitats and use it to show that most genes are specific to a single habitat. The small fraction of genes found in multiple habitats is enriched in antibiotic-resistance genes and markers for mobile genetic elements. By further clustering these species-level genes into 32 million protein families, we observed that a small fraction of these families contain the majority of the genes (0.6% of families account for 50% of the genes). The majority of species-level genes and protein families are rare. Furthermore, species-level genes, and in particular the rare ones, show low rates of positive (adaptive) selection, supporting a model in which most genetic variability observed within each protein family is neutral or nearly neutral.

Ultra-resolution Metagenomics: When Enough Is Not Enough

Technological advances in community sequencing have steadily increased the taxonomic resolution at which microbes can be delineated. In high-resolution metagenomics, bacterial strains can now be resolved, enhancing medical microbiology and the description of microbial evolution in vivo. In the Hildebrand lab, we are researching novel approaches to further increase the phylogenetic resolution of metagenomics. I propose that ultra-resolution metagenomics will be the next qualitative level of community sequencing, classified by the accurate resolution of ultra-rare genetic events, such as subclonal mutations present in all populations of evolving cells. This will be used to quantify evolutionary processes at ecologically relevant scales, monitor the progress of infections within a patient, and accurately track pathogens in food and infection chains. However, to develop this next metagenomic generation, we first need to understand the currently imposed limits of sequencing technologies, metagenomic strain delineation, and genome reconstructions.

Metapangenomics reveals depth-dependent shifts in metabolic potential for the ubiquitous marine bacterial SAR324 lineage

BACKGROUND

Oceanic microbiomes play a pivotal role in the global carbon cycle and are central to the transformation and recycling of carbon and energy in the ocean's interior. SAR324 is a ubiquitous but poorly understood uncultivated clade of Deltaproteobacteria that inhabits the entire water column, from ocean surface waters to its deep interior. Although some progress has been made in elucidating potential metabolic traits of SAR324 in the dark ocean, very little is known about the ecology and the metabolic capabilities of this group in the euphotic and twilight zones. To investigate the comparative genomics, ecology, and physiological potential of the SAR324 clade, we examined the distribution and variability of key genomic features and metabolic pathways in this group from surface waters to the abyss in the North Pacific Subtropical Gyre, one of the largest biomes on Earth.

RESULTS

We leveraged a pangenomic ecological approach, combining spatio-temporally resolved single-amplified genome, metagenomic, and metatranscriptomic datasets. The data revealed substantial genomic diversity throughout the SAR324 clade, with distinct depth and temporal distributions that clearly differentiated ecotypes. Phylogenomic subclade delineation, environmental distributions, genomic feature similarities, and metabolic capacities revealed strong congruence. The four SAR324 ecotypes delineated in this study revealed striking divergence from one another with respect to their habitat-specific metabolic potentials. The ecotypes living in the dark or twilight oceans shared genomic features and metabolic capabilities consistent with a sulfur-based chemolithoautotrophic lifestyle. In contrast, those inhabiting the sunlit ocean displayed higher plasticity energy-related metabolic pathways, supporting a presumptive photoheterotrophic lifestyle. In epipelagic SAR324 ecotypes, we observed the presence of two types of proton-pumping rhodopsins, as well as genomic, transcriptomic, and ecological evidence for active photoheterotrophy, based on xanthorhodopsin-like light-harvesting proteins.

CONCLUSIONS

Combining pangenomic and both metagenomic and metatranscriptomic profiling revealed a striking divergence in the vertical distribution, genomic composition, metabolic potential, and predicted lifestyle strategies of geographically co-located members of the SAR324 bacterial clade. The results highlight the utility of metapangenomic approaches employed across environmental gradients, to decipher the properties and variation in function and ecological traits of specific phylogenetic clades within complex microbiomes. Video abstract.

Genomic Analysis of Family UBA6911 (Group 18 Acidobacteria) Expands the Metabolic Capacities of the Phylum and Highlights Adaptations to Terrestrial Habitats

Approaches for recovering and analyzing genomes belonging to novel, hitherto-unexplored bacterial lineages have provided invaluable insights into the metabolic capabilities and ecological roles of yet-uncultured taxa. The phylum Acidobacteria is one of the most prevalent and ecologically successful lineages on Earth, yet currently, multiple lineages within this phylum remain unexplored. Here, we utilize genomes recovered from Zodletone Spring, an anaerobic sulfide and sulfur-rich spring in southwestern Oklahoma, as well as from multiple disparate soil and nonsoil habitats, to examine the metabolic capabilities and ecological role of members of family UBA6911 (group 18) Acidobacteria. The analyzed genomes clustered into five distinct genera, with genera Gp18_AA60 and QHZH01 recovered from soils, genus Ga0209509 from anaerobic digestors, and genera Ga0212092 and UBA6911 from freshwater habitats. All genomes analyzed suggested that members of Acidobacteria group 18 are metabolically versatile heterotrophs capable of utilizing a wide range of proteins, amino acids, and sugars as carbon sources, possess respiratory and fermentative capacities, and display few auxotrophies. Soil-dwelling genera were characterized by larger genome sizes, higher numbers of CRISPR loci, an expanded carbohydrate active enzyme (CAZyme) machinery enabling debranching of specific sugars from polymers, possession of a C1 (methanol and methylamine) degradation machinery, and a sole dependence on aerobic respiration. In contrast, nonsoil genomes encoded a more versatile respiratory capacity for oxygen, nitrite, sulfate, and trimethylamine N-oxide (TMAO) respiration, as well as the potential for utilizing the Wood-Ljungdahl (WL) pathway as an electron sink during heterotrophic growth. Our results not only expand our knowledge of the metabolism of a yet-uncultured bacterial lineage but also provide interesting clues on how terrestrialization and niche adaptation drive metabolic specialization within the Acidobacteria. IMPORTANCE Members of the Acidobacteria are important players in global biogeochemical cycles, especially in soils. A wide range of acidobacterial lineages remain currently unexplored. We present a detailed genomic characterization of genomes belonging to family UBA6911 (also known as group 18) within the phylum Acidobacteria. The genomes belong to different genera and were obtained from soil (genera Gp18_AA60 and QHZH01), freshwater habitats (genera Ga0212092 and UBA6911), and an anaerobic digestor (genus Ga0209509). While all members of the family shared common metabolic features, e.g., heterotrophic respiratory abilities, broad substrate utilization capacities, and few auxotrophies, distinct differences between soil and nonsoil genera were observed. Soil genera were characterized by expanded genomes, higher numbers of CRISPR loci, a larger carbohydrate active enzyme (CAZyme) repertoire enabling monomer extractions from polymer side chains, and methylotrophic (methanol and methylamine) degradation capacities. In contrast, nonsoil genera encoded more versatile respiratory capacities for utilizing nitrite, sulfate, TMAO, and the WL pathway, in addition to oxygen as electron acceptors. Our results not only broaden our understanding of the metabolic capacities within the Acidobacteria but also provide interesting clues on how terrestrialization shaped Acidobacteria evolution and niche adaptation.

Pangenomics reveals alternative environmental lifestyles among chlamydiae

Chlamydiae are highly successful strictly intracellular bacteria associated with diverse eukaryotic hosts. Here we analyzed metagenome-assembled genomes of the "Genomes from Earth's Microbiomes" initiative from diverse environmental samples, which almost double the known phylogenetic diversity of the phylum and facilitate a highly resolved view at the chlamydial pangenome. Chlamydiae are defined by a relatively large core genome indicative of an intracellular lifestyle, and a highly dynamic accessory genome of environmental lineages. We observe chlamydial lineages that encode enzymes of the reductive tricarboxylic acid cycle and for light-driven ATP synthesis. We show a widespread potential for anaerobic energy generation through pyruvate fermentation or the arginine deiminase pathway, and we add lineages capable of molecular hydrogen production. Genome-informed analysis of environmental distribution revealed lineage-specific niches and a high abundance of chlamydiae in some habitats. Together, our data provide an extended perspective of the variability of chlamydial biology and the ecology of this phylum of intracellular microbes.

CoreCruncher: Fast and Robust Construction of Core Genomes in Large Prokaryotic Data Sets

The core genome represents the set of genes shared by all, or nearly all, strains of a given population or species of prokaryotes. Inferring the core genome is integral to many genomic analyses, however, most methods rely on the comparison of all the pairs of genomes; a step that is becoming increasingly difficult given the massive accumulation of genomic data. Here, we present CoreCruncher; a program that robustly and rapidly constructs core genomes across hundreds or thousands of genomes. CoreCruncher does not compute all pairwise genome comparisons and uses a heuristic based on the distributions of identity scores to classify sequences as orthologs or paralogs/xenologs. Although it is much faster than current methods, our results indicate that our approach is more conservative than other tools and less sensitive to the presence of paralogs and xenologs. CoreCruncher is freely available from: https://github.com/lbobay/CoreCruncher. CoreCruncher is written in Python 3.7 and can also run on Python 2.7 without modification. It requires the python library Numpy and either Usearch or Blast. Certain options require the programs muscle or mafft.

Bacillus pumilus Group Comparative Genomics: Toward Pangenome Features, Diversity, and Marine Environmental Adaptation

Background

Members of the Bacillus pumilus group (abbreviated as the Bp group) are quite diverse and ubiquitous in marine environments, but little is known about correlation with their terrestrial counterparts. In this study, 16 marine strains that we had isolated before were sequenced and comparative genome analyses were performed with a total of 52 Bp group strains. The analyses included 20 marine isolates (which included the 16 new strains) and 32 terrestrial isolates, and their evolutionary relationships, differentiation, and environmental adaptation.

Results

Phylogenomic analysis revealed that the marine Bp group strains were grouped into three species: B. pumilus, B. altitudinis and B. safensis. All the three share a common ancestor. However, members of B. altitudinis were observed to cluster independently, separating from the other two, thus diverging from the others. Consistent with the universal nature of genes involved in the functioning of the translational machinery, the genes related to translation were enriched in the core genome. Functional genomic analyses revealed that the marine-derived and the terrestrial strains showed differences in certain hypothetical proteins, transcriptional regulators, K⁺ transporter (TrK) and ABC transporters. However, species differences showed the precedence of environmental adaptation discrepancies. In each species, land specific genes were found with possible functions that likely facilitate survival in diverse terrestrial niches, while marine bacteria were enriched with genes of unknown functions and those related to transcription, phage defense, DNA recombination and repair.

Conclusion

Our results indicated that the Bp isolates show distinct genomic features even as they share a common core. The marine and land isolates did not evolve independently; the transition between marine and non-marine habitats might have occurred multiple times. The lineage exhibited a priority effect over the niche in driving their dispersal. Certain intra-species niche specific genes could be related to a strains adaptation to its respective marine or terrestrial environment(s). In summary, this report describes the systematic evolution of 52 Bp group strains and will facilitate future studies toward understanding their ecological role and adaptation to marine and/or terrestrial environments.

The past, present and future of the tree of life

The advent of comparative genomics in the late 1990s led to the discovery of extensive lateral gene transfer in prokaryotes. The resulting debate over whether life as a whole is best represented as a tree or a network has since given way to a general consensus in which trees and networks co-exist rather than stand in opposition. Embracing this consensus allows us to move beyond the question of which is true or false. The future of the tree of life debate lies in asking what trees and networks can, and should, do for science.

A genomic catalog of Earth's microbiomes

The reconstruction of bacterial and archaeal genomes from shotgun metagenomes has enabled insights into the ecology and evolution of environmental and host-associated microbiomes. Here we applied this approach to >10,000 metagenomes collected from diverse habitats covering all of Earth's continents and oceans, including metagenomes from human and animal hosts, engineered environments, and natural and agricultural soils, to capture extant microbial, metabolic and functional potential. This comprehensive catalog includes 52,515 metagenome-assembled genomes representing 12,556 novel candidate species-level operational taxonomic units spanning 135 phyla. The catalog expands the known phylogenetic diversity of bacteria and archaea by 44% and is broadly available for streamlined comparative analyses, interactive exploration, metabolic modeling and bulk download. We demonstrate the utility of this collection for understanding secondary-metabolite biosynthetic potential and for resolving thousands of new host linkages to uncultivated viruses. This resource underscores the value of genome-centric approaches for revealing genomic properties of uncultivated microorganisms that affect ecosystem processes.