Ecological selection for small microbial genomes along a temperate-to-thermal soil gradient

Acidic proteomes are linked to microbial alkaline preference in African lakes

Microbial amino acid composition (AA) reflects adaptive strategies of cellular and molecular regulations such as a high proportion of acidic AAs, including glutamic and aspartic acids in alkaliphiles. It remains understudied how microbial AA content is linked to their pH adaptation especially in natural environments. Here we examined prokaryotic communities and their AA composition of genes with metagenomics for 39 water and sediments of East African lakes along a gradient of pH spanning from 7.2 to 10.1. We found that Shannon diversity declined with the increasing pH and that species abundance were either positively or negatively associated with pH, indicating their distinct habitat preference in lakes. Microbial communities showed higher acidic proteomes in alkaline than neutral lakes. Species acidic proteomes were also positively correlated with their pH preference, which was consistent across major bacterial lineages. These results suggest selective pressure associated with high pH likely shape microbial amino acid composition both at the species and community levels. Comparative genome analyses further revealed that alkaliphilic microbes contained more functional genes with higher acidic AAs when compared to those in neutral conditions. These traits included genes encoding diverse classes of cation transmembrane transporters, antiporters, and compatible solute transporters, which are involved in cytoplasmic pH homeostasis and osmotic stress defense under high pH conditions. Our results provide the field evidence for the strong relationship between prokaryotic AA composition and their habitat preference and highlight amino acid optimization as strategies for environmental adaptation.

Comparative genomic analysis of thermophilic fungi reveals convergent evolutionary adaptations and gene losses

Thermophily is a trait scattered across the fungal tree of life, with its highest prevalence within three fungal families (Chaetomiaceae, Thermoascaceae, and Trichocomaceae), as well as some members of the phylum Mucoromycota. We examined 37 thermophilic and thermotolerant species and 42 mesophilic species for this study and identified thermophily as the ancestral state of all three prominent families of thermophilic fungi. Thermophilic fungal genomes were found to encode various thermostable enzymes, including carbohydrate-active enzymes such as endoxylanases, which are useful for many industrial applications. At the same time, the overall gene counts, especially in gene families responsible for microbial defense such as secondary metabolism, are reduced in thermophiles compared to mesophiles. We also found a reduction in the core genome size of thermophiles in both the Chaetomiaceae family and the Eurotiomycetes class. The Gene Ontology terms lost in thermophilic fungi include primary metabolism, transporters, UV response, and O-methyltransferases. Comparative genomics analysis also revealed higher GC content in the third base of codons (GC3) and a lower effective number of codons in fungal thermophiles than in both thermotolerant and mesophilic fungi. Furthermore, using the Support Vector Machine classifier, we identified several Pfam domains capable of discriminating between genomes of thermophiles and mesophiles with 94% accuracy. Using AlphaFold2 to predict protein structures of endoxylanases (GH10), we built a similarity network based on the structures. We found that the number of disulfide bonds appears important for protein structure, and the network clusters based on protein structures correlate with the optimal activity temperature. Thus, comparative genomics offers new insights into the biology, adaptation, and evolutionary history of thermophilic fungi while providing a parts list for bioengineering applications.

Two RNA Folds from One Sequence: A Ribozyme with Versatile Substrate Processing Abilities

We report the design of a single RNA sequence capable of adopting one of two ribozyme folds and catalyzing the cleavage and/or ligation of the respective substrates. The RNA is able to change its conformation in response to its environment, hence it is called chameleon ribozyme (CHR). Efficient RNA cleavage of two different substrates as well as RNA ligation by CHR is demonstrated in separate experiments and in a one pot reaction. Our study shows that sequence variants of the hairpin ribozyme intersect with the hammerhead ribozyme and that rather short RNA molecules can have comprehensive conformational flexibility, which is an important feature for the emergence of new functional folds in early evolution.

Eco-evolutionary strategies for relieving carbon limitation under salt stress differ across microbial clades

With the continuous expansion of saline soils under climate change, understanding the eco-evolutionary tradeoff between the microbial mitigation of carbon limitation and the maintenance of functional traits in saline soils represents a significant knowledge gap in predicting future soil health and ecological function. Through shotgun metagenomic sequencing of coastal soils along a salinity gradient, we show contrasting eco-evolutionary directions of soil bacteria and archaea that manifest in changes to genome size and the functional potential of the soil microbiome. In salt environments with high carbon requirements, bacteria exhibit reduced genome sizes associated with a depletion of metabolic genes, while archaea display larger genomes and enrichment of salt-resistance, metabolic, and carbon-acquisition genes. This suggests that bacteria conserve energy through genome streamlining when facing salt stress, while archaea invest in carbon-acquisition pathways to broaden their resource usage. These findings suggest divergent directions in eco-evolutionary adaptations to soil saline stress amongst microbial clades and serve as a foundation for understanding the response of soil microbiomes to escalating climate change.

Deciphering the microbial communities of alkaline hot spring in Panamik, Ladakh, India using a high-throughput sequencing approach

Due to their distinctive physicochemical characteristics, hot springs are extremely important. The whole genome metagenomic sequencing technology can be utilized to analyze the diverse microbial community that thrives in this habitat due to the particular selection pressure that prevails there. The current investigation emphasizes on culture-independent metagenomic study of the Panamik hot spring and its nearby areas from Ladakh, India. Based on different diversity indices, sequence analysis of the soil reservoir showed higher species richness and diversity in comparison to water and sediment samples. The mineral content and various physicochemical pameters like temperature, pH had an impact on the composition of the microbial community of the geothermal springs. The phyla Proteobacteria, Cyanobacteria, Bacteroidetes, Actinobacter, Firmicutes, and Verrucomicrobia in bacterial domain dominate the thermos-alkaline spring at Panamik in different concentrations. Economically significant microbes from the genera Actinobacter, Thermosynechoccus, Candidatus Solibacter, Chthoniobacter, Synechoccus, Pseudomonas and Sphingomonas, were prevalent in hot spring. In the archaeal domain, the most dominant phylum and genera were Euryarchaeota and Thermococcus in all the samples. Further, the most abundant species were Methanosarcina barkeri, Nitrospumilus maritimus and Methanosarcina acetivorans. The present study which only examined one of the several thermal springs present in the Himalayan geothermal area, should be regarded as a preliminary investigation of the microbiota that live in the hot springs on these remote areas. These findings suggest that further investigations should be undertaken to characterize the ecosystems of the Panamik hot spring, which serve as a repository for unidentified microbial lineages.

Metagenomic Thermometer

Various microorganisms exist in environments, and each of them has its optimal growth temperature (OGT). The relationship between genomic information and OGT of each species has long been studied, and one such study revealed that OGT of prokaryotes can be accurately predicted based on the fraction of seven amino acids (IVYWREL) among all encoded amino-acid sequences in its genome. Extending this discovery, we developed a 'Metagenomic Thermometer' as a means of predicting environmental temperature based on metagenomic sequences. Temperature prediction of diverse environments using publicly available metagenomic data revealed that the Metagenomic Thermometer can predict environmental temperatures with small temperature changes and little influx of microorganisms from other environments. The accuracy of the Metagenomic Thermometer was also confirmed by a demonstration experiment using an artificial hot water canal. The Metagenomic Thermometer was also applied to human gut metagenomic samples, yielding a reasonably accurate value for human body temperature. The result further suggests that deep body temperature determines the dominant lineage of the gut community. Metagenomic Thermometer provides a new insight into temperature-driven community assembly based on amino-acid composition rather than microbial taxa.

Thermal adaptation of soil microbial growth traits in response to chronic warming

IMPORTANCE

Soils are the largest terrestrial carbon sink and the foundation of our food, fiber, and fuel systems. Healthy soils are carbon sinks, storing more carbon than they release. This reduces the amount of carbon dioxide released into the atmosphere and buffers against climate change. Soil microbes drive biogeochemical cycling and contribute to soil health through organic matter breakdown, plant growth promotion, and nutrient distribution. In this study, we determined how soil microbial growth traits respond to long-term soil warming. We found that bacterial isolates from warmed plots showed evidence of adaptation of optimum growth temperature. This suggests that increased microbial biomass and growth in a warming world could result in greater carbon storage. As temperatures increase, greater microbial activity may help reduce the soil carbon feedback loop. Our results provide insight on how atmospheric carbon cycling and soil health may respond in a warming world.

Life history strategies of soil bacterial communities across global terrestrial biomes

The life history strategies of soil microbes determine their metabolic potential and their response to environmental changes. Yet these strategies remain poorly understood. Here we use shotgun metagenomes from terrestrial biomes to characterize overarching covariations of the genomic traits that capture dominant life history strategies in bacterial communities. The emerging patterns show a triangle of life history strategies shaped by two trait dimensions, supporting previous theoretical and isolate-based studies. The first dimension ranges from streamlined genomes with simple metabolisms to larger genomes and expanded metabolic capacities. As metabolic capacities expand, bacterial communities increasingly differentiate along a second dimension that reflects a trade-off between increasing capacities for environmental responsiveness or for nutrient recycling. Random forest analyses show that soil pH, C:N ratio and precipitation patterns together drive the dominant life history strategy of soil bacterial communities and their biogeographic distribution. Our findings provide a trait-based framework to compare life history strategies of soil bacteria.

Environmental selection and evolutionary process jointly shape genomic and functional profiles of mangrove rhizosphere microbiomes

Mangrove reforestation with introduced species has been an important strategy to restore mangrove ecosystem functioning. However, how such activities affect microbially driven methane (CH4), nitrogen (N), and sulfur (S) cycling of rhizosphere microbiomes remains unclear. To understand the effect of environmental selection and the evolutionary process on microbially driven biogeochemical cycles in native and introduced mangrove rhizospheres, we analyzed key genomic and functional profiles of rhizosphere microbiomes from native and introduced mangrove species by metagenome sequencing technologies. Compared with the native mangrove (Kandelia obovata, KO), the introduced mangrove (Sonneratia apetala, SA) rhizosphere microbiome had significantly (p < 0.05) higher average genome size (AGS) (5.8 vs. 5.5 Mb), average 16S ribosomal RNA gene copy number (3.5 vs. 3.1), relative abundances of mobile genetic elements, and functional diversity in terms of the Shannon index (7.88 vs. 7.84) but lower functional potentials involved in CH4 cycling (e.g., mcrABCDG and pmoABC), N2 fixation (nifHDK), and inorganic S cycling (dsrAB, dsrC, dsrMKJOP, soxB, sqr, and fccAB). Similar results were also observed from the recovered Proteobacterial metagenome-assembled genomes with a higher AGS and distinct functions in the introduced mangrove rhizosphere. Additionally, salinity and ammonium were identified as the main environmental drivers of functional profiles of mangrove rhizosphere microbiomes through deterministic processes. This study advances our understanding of microbially mediated biogeochemical cycling of CH4, N, and S in the mangrove rhizosphere and provides novel insights into the influence of environmental selection and evolutionary processes on ecosystem functions, which has important implications for future mangrove reforestation.

Hot spring distribution and survival mechanisms of thermophilic comammox Nitrospira

The recent discovery of Nitrospira species capable of complete ammonia oxidation (comammox) in non-marine natural and engineered ecosystems under mesothermal conditions has changed our understanding of microbial nitrification. However, little is known about the occurrence of comammox bacteria or their ability to survive in moderately thermal and/or hyperthermal habitats. Here, we report the wide distribution of comammox Nitrospira in five terrestrial hot springs at temperatures ranging from 36 to 80°C and provide metagenome-assembled genomes of 11 new comammox strains. Interestingly, the identification of dissimilatory nitrate reduction to ammonium (DNRA) in thermophilic comammox Nitrospira lineages suggests that they have versatile ecological functions as both sinks and sources of ammonia, in contrast to the described mesophilic comammox lineages, which lack the DNRA pathway. Furthermore, the in situ expression of key genes associated with nitrogen metabolism, thermal adaptation, and oxidative stress confirmed their ability to survive in the studied hot springs and their contribution to nitrification in these environments. Additionally, the smaller genome size and higher GC content, less polar and more charged amino acids in usage profiles, and the expression of a large number of heat shock proteins compared to mesophilic comammox strains presumably confer tolerance to thermal stress. These novel insights into the occurrence, metabolic activity, and adaptation of comammox Nitrospira in thermal habitats further expand our understanding of the global distribution of comammox Nitrospira and have significant implications for how these unique microorganisms have evolved thermal tolerance strategies.

High speciation rate of niche specialists in hot springs

Ecological and evolutionary processes simultaneously regulate microbial diversity, but the evolutionary processes and their driving forces remain largely unexplored. Here we investigated the ecological and evolutionary characteristics of microbiota in hot springs spanning a broad temperature range (54.8-80 °C) by sequencing the 16S rRNA genes. Our results demonstrated that niche specialists and niche generalists are embedded in a complex interaction of ecological and evolutionary dynamics. On the thermal tolerance niche axis, thermal (T) sensitive (at a specific temperature) versus T-resistant (at least in five temperatures) species were characterized by different niche breadth, community abundance and dispersal potential, consequently differing in potential evolutionary trajectory. The niche-specialized T-sensitive species experienced strong temperature barriers, leading to completely species shift and high fitness but low abundant communities at each temperature ("home niche"), and such trade-offs thus reinforced peak performance, as evidenced by high speciation across temperatures and increasing diversification potential with temperature. In contrast, T-resistant species are advantageous of niche expansion but with poor local performance, as shown by wide niche breadth with high extinction, indicating these niche generalists are "jack-of-all-trades, master-of-none". Despite of such differences, the T-sensitive and T-resistant species are evolutionarily interacted. Specifically, the continuous transition from T-sensitive to T-resistant species insured the exclusion probability of T-resistant species at a relatively constant level across temperatures. The co-evolution and co-adaptation of T-sensitive and T-resistant species were in line with the red queen theory. Collectively, our findings demonstrate that high speciation of niche specialists could alleviate the environmental-filtering-induced negative effect on diversity.

Prokaryotic Life Associated with Coal-Fire Gas Vents Revealed by Metagenomics

The natural combustion of underground coal seams leads to the formation of gas, which contains molecular hydrogen and carbon monoxide. In places where hot coal gases are released to the surface, specific thermal ecosystems are formed. Here, 16S rRNA gene profiling and shotgun metagenome sequencing were employed to characterize the taxonomic diversity and genetic potential of prokaryotic communities of the near-surface ground layer near hot gas vents in an open quarry heated by a subsurface coal fire. The communities were dominated by only a few groups of spore-forming Firmicutes, namely the aerobic heterotroph Candidatus Carbobacillus altaicus, the aerobic chemolitoautotrophs Kyrpidia tusciae and Hydrogenibacillus schlegelii, and the anaerobic chemolithoautotroph Brockia lithotrophica. Genome analysis predicted that these species can obtain energy from the oxidation of hydrogen and/or carbon monoxide in coal gases. We assembled the first complete closed genome of a member of uncultured class-level division DTU015 in the phylum Firmicutes. This bacterium, 'Candidatus Fermentithermobacillus carboniphilus' Bu02, was predicted to be rod-shaped and capable of flagellar motility and sporulation. Genome analysis showed the absence of aerobic and anaerobic respiration and suggested chemoheterotrophic lifestyle with the ability to ferment peptides, amino acids, N-acetylglucosamine, and tricarboxylic acid cycle intermediates. Bu02 bacterium probably plays the role of a scavenger, performing the fermentation of organics formed by autotrophic Firmicutes supported by coal gases. A comparative genome analysis of the DTU015 division revealed that most of its members have a similar lifestyle.

Warmer and drier ecosystems select for smaller bacterial genomes in global soils

Bacterial genome size reflects bacterial evolutionary processes and metabolic lifestyles, with implications for microbial community assembly and ecosystem functions. However, to understand the extent of genome-mediated microbial responses to environmental selections, we require studies that observe genome size distributions along environmental gradients representing different conditions that soil bacteria normally encounter. In this study, we used surface soils collected from 237 sites across the globe and analyzed how environmental conditions (e.g., soil carbon and nutrients, aridity, pH, and temperature) affect soil bacterial occurrences and genome size at the community level using bacterial community profiling. We used a joint species distribution model to quantify the effects of environments on species occurrences and found that aridity was a major regulator of genome size with warmer and drier environments selecting bacteria with smaller genomes. Drought-induced physiological constraints on bacterial growth (e.g., water scarcity for cell metabolisms) may have led to these correlations. This finding suggests that increasing cover by warmer and drier ecosystems may result in bacterial genome simplifications by a reduction of genome size.

Updated Virophage Taxonomy and Distinction from Polinton-like Viruses

Virophages are small dsDNA viruses that hijack the machinery of giant viruses during the co-infection of a protist (i.e., microeukaryotic) host and represent an exceptional case of "hyperparasitism" in the viral world. While only a handful of virophages have been isolated, a vast diversity of virophage-like sequences have been uncovered from diverse metagenomes. Their wide ecological distribution, idiosyncratic infection and replication strategy, ability to integrate into protist and giant virus genomes and potential role in antiviral defense have made virophages a topic of broad interest. However, one limitation for further studies is the lack of clarity regarding the nomenclature and taxonomy of this group of viruses. Specifically, virophages have been linked in the literature to other "virophage-like" mobile genetic elements and viruses, including polinton-like viruses (PLVs), but there are no formal demarcation criteria and proper nomenclature for either group, i.e., virophage or PLVs. Here, as part of the ICTV Virophage Study Group, we leverage a large set of genomes gathered from published datasets as well as newly generated protist genomes to propose delineation criteria and classification methods at multiple taxonomic ranks for virophages 'sensu stricto', i.e., genomes related to the prototype isolates Sputnik and mavirus. Based on a combination of comparative genomics and phylogenetic analyses, we show that this group of virophages forms a cohesive taxon that we propose to establish at the class level and suggest a subdivision into four orders and seven families with distinctive ecogenomic features. Finally, to illustrate how the proposed delineation criteria and classification method would be used, we apply these to two recently published datasets, which we show include both virophages and other virophage-related elements. Overall, we see this proposed classification as a necessary first step to provide a robust taxonomic framework in this area of the virosphere, which will need to be expanded in the future to cover other virophage-related viruses such as PLVs.

Identifying candidate structured RNAs in CRISPR operons

Noncoding RNAs with secondary structures play important roles in CRISPR-Cas systems. Many of these structures likely remain undiscovered. We used a large-scale comparative genomics approach to predict 156 novel candidate structured RNAs from 36,111 CRISPR-Cas systems. A number of these were found to overlap with coding genes, including palindromic candidates that overlapped with a variety of Cas genes in type I and III systems. Among these 156 candidates, we identified 46 new models of CRISPR direct repeats and 1 tracrRNA. This tracrRNA model occasionally overlapped with predicted cas9 coding regions, emphasizing the importance of expanding our search windows for novel structure RNAs in coding regions. We also demonstrated that the antirepeat sequence in this tracrRNA model can be used to accurately assign thousands of predicted CRISPR arrays to type II-C systems. This study highlights the importance of unbiased identification of candidate structured RNAs across CRISPR-Cas systems.

Discovery of the non-cosmopolitan lineages in Candidatus Thermoprofundales

The recently proposed order Candidatus Thermoprofundales, currently containing only one family-level lineage Marine Benthic Group-D (MBG-D), is distributed in global subsurface ecosystems and ecologically important, but its diversity, evolution and metabolism remain largely unknown. Here we described two novel family-level specialized lineages in Ca. Thermoprofundales, JdFR-43 and HyVt, which are restricted to specific biotopes (primarily in marine hydrothermal vents and occasionally in oil reservoirs and hot springs) in contrast to the cosmopolitan lineage MBG-D. The comparative genomics revealed that the specialized lineages have streamlined genomes, higher GC contents, enriched genes associated with nucleotide biosynthesis, ribosome biogenesis and DNA repair and additional thermostable aminopeptidases, enabling them to adapt to high-temperature habitats such as marine hydrothermal vents, deep subsurface oil reservoirs and hot springs. On the contrary, the unique metabolic traits of the cosmopolitan MBG-D, motility, glycolysis, butanoate metabolism, secondary metabolites production and additional genes for specific peptides and carbohydrates degradation potentially enhance its response to environmental change. Substrate preference is found for most MAGs across all lineages with the ability to utilize both polysaccharides (chitin and starch) and proteinaceous substances, whereas JdFR-43 members from oil reservoirs can only utilize proteins. These results expand the diversity of Ca. Thermoprofundales significantly and further improve our understandings of the adaptations of Ca. Thermoprofundales to various environments. This article is protected by copyright. All rights reserved.

Ecological selection of bacterial taxa with larger genome sizes in response to polycyclic aromatic hydrocarbons stress

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous priority pollutants that cause great damage to the natural environment and health. Average genome size in a community is critical for shedding light on microbiome's functional response to pollution stress within an environment. Here, microcosms under different concentrations were performed to evaluate the selection of PAHs stress on the average genome size in a community. We found the distinct communities of significantly larger genome size with the increase of PAHs concentration gradients in soils, and consistent trends were discovered in soils at different latitudes. The abundance of Proteobacteria and Deinococcus-Thermus with relatively larger genomes increased along with PAHs stress and well adapted to polluted environments. In contrast, the abundance of Patescibacteria with a highly streamlined and smaller genome decreased, implying complex interactions between environmental selection and functional fitness resulted in bacteria with larger genomes becoming more abundant. Moreover, we confirmed the increased capacity for horizontal transfer of degrading genes between communities by showing an increased connection number per node positively related to the nidA gene along the concentration gradients in the co-occurrence network. Our findings suggest PAHs tend to select bacterial taxa with larger genome sizes, with significant consequences for community stability and potential biodegradation strategies.

Variation in genomic traits of microbial communities among ecosystems

Free-living bacteria in nutrient limited environments often exhibit traits which may reduce the cost of reproduction, such as smaller genome size, low GC content and fewer sigma (σ) factor and 16S rRNA gene copies. Despite the potential utility of these traits to detect relationships between microbial communities and ecosystem-scale properties, few studies have assessed these traits on a community-scale. Here, we analysed these traits from publicly available metagenomes derived from marine, soil, host-associated and thermophilic communities. In marine and thermophilic communities, genome size and GC content declined in parallel, consistent with genomic streamlining, with GC content in thermophilic communities generally higher than in marine systems. In contrast, soil communities averaging smaller genomes featured higher GC content and were often from low-carbon environments, suggesting unique selection pressures in soil bacteria. The abundance of specific σ-factors varied with average genome size and ecosystem type. In oceans, abundance of fliA, a σ-factor controlling flagella biosynthesis, was positively correlated with community average genome size—reflecting known trade-offs between nutrient conservation and chemotaxis. In soils, a high abundance of the stress response σ-factor gene rpoS was associated with smaller average genome size and often located in harsh and/or carbon-limited environments—a result which tracks features observed in culture and indicates an increased capacity for stress response in nutrient-poor soils. This work shows how ecosystem-specific constraints are associated with trade-offs which are embedded in the genomic features of bacteria in microbial communities, and which can be detected at the community level, highlighting the importance of genomic features in microbial community analysis.

Minnesota peat viromes reveal terrestrial and aquatic niche partitioning for local and global viral populations

BACKGROUND

Peatlands are expected to experience sustained yet fluctuating higher temperatures due to climate change, leading to increased microbial activity and greenhouse gas emissions. Despite mounting evidence for viral contributions to these processes in peatlands underlain with permafrost, little is known about viruses in other peatlands. More generally, soil viral biogeography and its potential drivers are poorly understood at both local and global scales. Here, 87 metagenomes and five viral size-fraction metagenomes (viromes) from a boreal peatland in northern Minnesota (the SPRUCE whole-ecosystem warming experiment and surrounding bog) were analyzed for dsDNA viral community ecological patterns, and the recovered viral populations (vOTUs) were compared with our curated PIGEON database of 266,125 vOTUs from diverse ecosystems.

RESULTS

Within the SPRUCE experiment, viral community composition was significantly correlated with peat depth, water content, and carbon chemistry, including CH₄ and CO₂ concentrations, but not with temperature during the first 2 years of warming treatments. Peat vOTUs with aquatic-like signatures (shared predicted protein content with marine and/or freshwater vOTUs) were significantly enriched in more waterlogged surface peat depths. Predicted host ranges for SPRUCE vOTUs were relatively narrow, generally within a single bacterial genus. Of the 4326 SPRUCE vOTUs, 164 were previously detected in other soils, mostly peatlands. None of the previously identified 202,371 marine and freshwater vOTUs in our PIGEON database were detected in SPRUCE peat, but 0.4% of 80,714 viral clusters (VCs, grouped by predicted protein content) were shared between soil and aquatic environments. On a per-sample basis, vOTU recovery was 32 times higher from viromes compared with total metagenomes.

CONCLUSIONS

Results suggest strong viral "species" boundaries between terrestrial and aquatic ecosystems and to some extent between peat and other soils, with differences less pronounced at higher taxonomic levels. The significant enrichment of aquatic-like vOTUs in more waterlogged peat suggests that viruses may also exhibit niche partitioning on more local scales. These patterns are presumably driven in part by host ecology, consistent with the predicted narrow host ranges. Although more samples and increased sequencing depth improved vOTU recovery from total metagenomes, the substantially higher per-sample vOTU recovery after viral particle enrichment highlights the utility of soil viromics. Video abstract The importance of Minnesota peat viromes in revealing terrestrial and aquatic niche partitioning for viral populations.

Mechanistic insights into the success of xenobiotic degraders resolved from metagenomes of microbial enrichment cultures

Even though microbial communities can be more effective at degrading xenobiotics than cultured micro-organisms, yet little is known about the microbial strategies that underpin xenobiotic biodegradation by microbial communities. Here, we employ metagenomic community sequencing to explore the mechanisms that drive the development of 49 xenobiotic-degrading microbial communities, which were enriched from 7 contaminated soils or sediments with a range of xenobiotic compounds. We show that multiple microbial strategies likely drive the development of xenobiotic degrading communities, notably (i) presence of genes encoding catabolic enzymes to degrade xenobiotics; (ii) presence of genes encoding efflux pumps; (iii) auxiliary catabolic genes on plasmids; and (iv) positive interactions dominate microbial communities with efficient degradation. Overall, the integrated analyses of microbial ecological strategies advance our understanding of microbial processes driving the biodegradation of xenobiotics and promote the design of bioremediation systems.

Metagenomic Analysis of the Microbial Community in the Underground Coal Fire Area (Kemerovo Region, Russia) Revealed Predominance of Thermophilic Members of the Phyla Deinococcus-Thermus, Aquificae, and Firmicutes

Abstract

Underground burning of coal seams accompanied by release of gases leads to development of local thermal ecosystems. We investigated the microbial community of the ground heated to 72°C in the release area of hot gases resulting from underground combustion of coal mining waste at the Bungurskiy-Severny coal deposit in the Kemerovo region of Russia. Analysis of the composition of the microbial community by 16S rRNA gene profiling revealed predominance of thermophilic bacteria of the phyla Deinococcus-Thermus, Aquificae, and Firmicutes. As a result of metagenomic analysis, 18 genomes of the main members of the microbial community were assembled, including the complete genomes of Hydrogenobacter thermophiles, a member of the candidate genus UBA11096 of the phylum Aquificae (RBS10-58), Thermoflexus hugenholtzii, and Thermus antranikianii. Analysis of the RBS10-58 genome indicates that this bacterium can autotrophically fix carbon in the reductive tricarboxylic acid cycle and obtain energy via oxidation of hydrogen and sulfur compounds with oxygen or nitrate as electron acceptors. Genome analysis of the two dominant Firmicutes species, Hydrogenibacillus schlegelii and an uncultured member of the class Thermaerobacteria, showed that these bacteria could grow aerobically by oxidizing hydrogen and carbon monoxide. Overall, the community was dominated by aerobic bacteria capable of growing autotrophically and obtaining energy via oxidation of the main components of coal gases, hydrogen and carbon monoxide. Thermus antranikianii, which makes up about a half of the microbial community, probably uses organic matter produced by autotrophic members of Firmicutes and Aquificae.

Cell size, genome size, and maximum growth rate are near-independent dimensions of ecological variation across bacteria and archaea

Among bacteria and archaea, maximum relative growth rate, cell diameter, and genome size are widely regarded as important influences on ecological strategy. Via the most extensive data compilation so far for these traits across all clades and habitats, we ask whether they are correlated and if so how. Overall, we found little correlation among them, indicating they should be considered as independent dimensions of ecological variation. Nor was correlation evident within particular habitat types. A weak nonlinearity (6% of variance) was found whereby high maximum growth rates (temperature-adjusted) tended to occur in the midrange of cell diameters. Species identified in the literature as oligotrophs or copiotrophs were clearly separated on the dimension of maximum growth rate, but not on the dimensions of genome size or cell diameter.

Thermophilic Chloroflexi Dominate in the Microbial Community Associated with Coal-Fire Gas Vents in the Kuznetsk Coal Basin, Russia

Thermal ecosystems associated with areas of underground burning coal seams are rare and poorly understood in comparison with geothermal objects. We studied the microbial communities associated with gas vents from the coal-fire in the mining wastes in the Kemerovo region of the Russian Federation. The temperature of the ground heated by the hot coal gases and steam coming out to the surface was 58 °C. Analysis of the composition of microbial communities revealed the dominance of Ktedonobacteria (the phylum Chloroflexi), known to be capable of oxidizing hydrogen and carbon monoxide. Thermophilic hydrogenotrophic Firmicutes constituted a minor part of the community. Among the well-known thermophiles, members of the phyla Aquificae, Deinococcus-Thermus and Bacteroidetes were also found. In the upper ground layer, Acidobacteria, Verrucomicrobia, Actinobacteria, Planctomycetes, as well as Proteobacteria of the alpha and gamma classes, typical of soils, were detected; their relative abundancies decreased with depth. The phylum Verrucomicrobia was dominated by Candidatus Udaeobacter, aerobic heterotrophs capable of generating energy through the oxidation of hydrogen present in the atmosphere in trace amounts. Archaea made up a small part of the communities and were represented by thermophilic ammonium-oxidizers. Overall, the community was dominated by bacteria, whose cultivated relatives are able to obtain energy through the oxidation of the main components of coal gases, hydrogen and carbon monoxide, under aerobic conditions.

Trait dimensions in bacteria and archaea compared to vascular plants

Bacteria and archaea have very different ecology compared to plants. One similarity, though, is that much discussion of their ecological strategies has invoked concepts such as oligotrophy or stress tolerance. For plants, so-called 'trait ecology'-strategy description reframed along measurable trait dimensions-has made global syntheses possible. Among widely measured trait dimensions for bacteria and archaea three main axes are evident. Maximum growth rate in association with rRNA operon copy number expresses a rate-yield trade-off that is analogous to the acquisitive-conservative spectrum in plants, though underpinned by different trade-offs. Genome size in association with signal transduction expresses versatility. Cell size has influence on diffusive uptake and on relative wall costs. These trait dimensions, and potentially others, offer promise for interpreting ecology. At the same time, there are very substantial differences from plant trait ecology. Traits and their underpinning trade-offs are different. Also, bacteria and archaea use a variety of different substrates. Bacterial strategies can be viewed both through the facet of substrate-use pathways, and also through the facet of quantitative traits such as maximum growth rate. Preliminary evidence shows the quantitative traits vary widely within substrate-use pathways. This indicates they convey information complementary to substrate use.

Metabolic capabilities mute positive response to direct and indirect impacts of warming throughout the soil profile

Increasing global temperatures are predicted to stimulate soil microbial respiration. The direct and indirect impacts of warming on soil microbes, nevertheless, remain unclear. This is particularly true for understudied subsoil microbes. Here, we show that 4.5 years of whole-profile soil warming in a temperate mixed forest results in altered microbial community composition and metabolism in surface soils, partly due to carbon limitation. However, microbial communities in the subsoil responded differently to warming than in the surface. Throughout the soil profile-but to a greater extent in the subsoil-physiologic and genomic measurements show that phylogenetically different microbes could utilize complex organic compounds, dampening the effect of altered resource availability induced by warming. We find subsoil microbes had 20% lower carbon use efficiencies and 47% lower growth rates compared to surface soils, which constrain microbial communities. Collectively, our results show that unlike in surface soils, elevated microbial respiration in subsoils may continue without microbial community change in the near-term.

Microbial tropicalization driven by a strengthening western ocean boundary current

Western boundary currents (WBCs) redistribute heat and oligotrophic seawater from the tropics to temperate latitudes, with several displaying substantial climate change-driven intensification over the last century. Strengthening WBCs have been implicated in the poleward range expansion of marine macroflora and fauna, however, the impacts on the structure and function of temperate microbial communities are largely unknown. Here we show that the major subtropical WBC of the South Pacific Ocean, the East Australian Current (EAC), transports microbial assemblages that maintain tropical and oligotrophic (k-strategist) signatures, to seasonally displace more copiotrophic (r-strategist) temperate microbial populations within temperate latitudes of the Tasman Sea. We identified specific characteristics of EAC microbial assemblages compared with non-EAC assemblages, including strain transitions within the SAR11 clade, enrichment of Prochlorococcus, predicted smaller genome sizes and shifts in the importance of several functional genes, including those associated with cyanobacterial photosynthesis, secondary metabolism and fatty acid and lipid transport. At a temperate time-series site in the Tasman Sea, we observed significant reductions in standing stocks of total carbon and chlorophyll a, and a shift towards smaller phytoplankton and carnivorous copepods, associated with the seasonal impact of the EAC microbial assemblage. In light of the substantial shifts in microbial assemblage structure and function associated with the EAC, we conclude that climate-driven expansions of WBCs will expand the range of tropical oligotrophic microbes, and potentially profoundly impact the trophic status of temperate waters.

Meta-Analysis of Microbial Communities in Hot Springs: Recurrent Taxa and Complex Shaping Factors beyond pH and Temperature

The study of microbial communities from extreme environments is a fascinating topic. With every study, biologists and ecologists reveal interesting facts and questions that dispel the old belief that these are inhospitable environments. In this work, we assess the microbial diversity of three hot springs from Neuquén, Argentina, using high-throughput amplicon sequencing. We predicted a distinct metabolic profile in the acidic and the circumneutral samples, with the first ones being dominated by chemolithotrophs and the second ones by chemoheterotrophs. Then, we collected data of the microbial communities of hot springs around the world in an effort to comprehend the roles of pH and temperature as shaping factors. Interestingly, there was a covariation between both parameters and the phylogenetic distance between communities; however, neither of them could explain much of the microbial profile in an ordination model. Moreover, there was no correlation between alpha diversity and these parameters. Therefore, the microbial communities' profile seemed to have complex shaping factors beyond pH and temperature. Lastly, we looked for taxa associated with different environmental conditions. Several such taxa were found. For example, Hydrogenobaculum was frequently present in acidic springs, as was the Sulfolobaceae family; on the other hand, Candidatus Hydrothermae phylum was strongly associated with circumneutral conditions. Interestingly, some singularities related to sites featuring certain taxa were also observed.

Dormancy dynamics and dispersal contribute to soil microbiome resilience

In disturbance ecology, stability is composed of resistance to change and resilience towards recovery after the disturbance subsides. Two key microbial mechanisms that can support microbiome stability include dormancy and dispersal. Specifically, microbial populations that are sensitive to disturbance can be re-seeded by local dormant pools of viable and reactivated cells, or by immigrants dispersed from regional metacommunities. However, it is difficult to quantify the contributions of these mechanisms to stability without, first, distinguishing the active from inactive membership, and, second, distinguishing the populations recovered by local resuscitation from those recovered by dispersed immigrants. Here, we investigate the contributions of dormancy dynamics (activation and inactivation), and dispersal to soil microbial community resistance and resilience. We designed a replicated, 45-week time-series experiment to quantify the responses of the active soil microbial community to a thermal press disturbance, including unwarmed control mesocosms, disturbed mesocosms without dispersal, and disturbed mesocosms with dispersal after the release of the stressor. Communities changed in structure within one week of warming. Though the disturbed mesocosms did not fully recover within 29 weeks, resuscitation of thermotolerant taxa was key for community transition during the press, and both resuscitation of opportunistic taxa and immigration contributed to community resilience. Also, mesocosms with dispersal were more resilient than mesocosms without. This work advances the mechanistic understanding of how microbiomes respond to disturbances in their environment.

This article is part of the theme issue ‘Conceptual challenges in microbial community ecology’.

Influence of Substrate Concentration on the Culturability of Heterotrophic Soil Microbes Isolated by High-Throughput Dilution-to-Extinction Cultivation

The vast majority of microbes inhabiting oligotrophic shallow subsurface soil environments have not been isolated or studied under controlled laboratory conditions. In part, the challenges associated with isolating shallow subsurface microbes may persist because microbes in deeper soils are adapted to low nutrient availability or quality. Here, we use high-throughput dilution-to-extinction culturing to isolate shallow subsurface microbes from a conifer forest in Arizona, USA. We hypothesized that the concentration of heterotrophic substrates in microbiological growth medium would affect which microbial taxa were culturable from these soils. To test this, we diluted cells extracted from soil into one of two custom-designed defined growth media that differed by 100-fold in the concentration of amino acids and organic carbon. Across the two media, we isolated a total of 133 pure cultures, all of which were classified as Actinobacteria or Alphaproteobacteria The substrate availability dictated which actinobacterial phylotypes were culturable but had no significant effect on the culturability of Alphaproteobacteria We isolated cultures that were representative of the most abundant phylotype in the soil microbial community (Bradyrhizobium spp.) and representatives of five of the top 10 most abundant Actinobacteria phylotypes, including Nocardioides spp., Mycobacterium spp., and several other phylogenetically divergent lineages. Flow cytometry of nucleic acid-stained cells showed that cultures isolated on low-substrate medium had significantly lower nucleic acid fluorescence than those isolated on high-substrate medium. These results show that dilution-to-extinction is an effective method to isolate abundant soil microbes and that the concentration of substrates in culture medium influences the culturability of specific microbial lineages.IMPORTANCE Isolating environmental microbes and studying their physiology under controlled conditions are essential aspects of understanding their ecology. Subsurface ecosystems are typically nutrient-poor environments that harbor diverse microbial communities-the majority of which are thus far uncultured. In this study, we use modified high-throughput cultivation methods to isolate subsurface soil microbes. We show that a component of whether a microbe is culturable from subsurface soils is the concentration of growth substrates in the culture medium. Our results offer new insight into technical approaches and growth medium design that can be used to access the uncultured diversity of soil microbes.

Charting the Complexity of the Marine Microbiome through Single-Cell Genomics

Marine bacteria and archaea play key roles in global biogeochemistry. To improve our understanding of this complex microbiome, we employed single-cell genomics and a randomized, hypothesis-agnostic cell selection strategy to recover 12,715 partial genomes from the tropical and subtropical euphotic ocean. A substantial fraction of known prokaryoplankton coding potential was recovered from a single, 0.4 mL ocean sample, which indicates that genomic information disperses effectively across the globe. Yet, we found each genome to be unique, implying limited clonality within prokaryoplankton populations. Light harvesting and secondary metabolite biosynthetic pathways were numerous across lineages, highlighting the value of single-cell genomics to advance the identification of ecological roles and biotechnology potential of uncultured microbial groups. This genome collection enabled functional annotation and genus-level taxonomic assignments for >80% of individual metagenome reads from the tropical and subtropical surface ocean, thus offering a model to improve reference genome databases for complex microbiomes.

Predictive genomic traits for bacterial growth in culture versus actual growth in soil

Relationships between microbial genes and performance are often evaluated in the laboratory in pure cultures, with little validation in nature. Here, we show that genomic traits related to laboratory measurements of maximum growth potential failed to predict the growth rates of bacteria in unamended soil, but successfully predicted growth responses to resource pulses: growth increased with 16S rRNA gene copy number and declined with genome size after substrate addition to soils, responses that were repeated in four different ecosystems. Genome size best predicted growth rate in response to addition of glucose alone; adding ammonium with glucose weakened the relationship, and the relationship was absent in nutrient-replete pure cultures, consistent with the idea that reduced genome size is a mechanism of nutrient conservation. Our findings demonstrate that genomic traits of soil bacteria can map to their ecological performance in nature, but the mapping is poor under native soil conditions, where genomic traits related to stress tolerance may prove more predictive. These results remind that phenotype depends on environmental context, underscoring the importance of verifying proposed schemes of trait-based strategies through direct measurement of performance in nature, an important and currently missing foundation for translating microbial processes from genes to ecosystems.

Predictive genomic traits for bacterial growth in culture versus actual growth in soil

Relationships between microbial genes and performance are often evaluated in the laboratory in pure cultures, with little validation in nature. Here, we show that genomic traits related to laboratory measurements of maximum growth potential failed to predict the growth rates of bacteria in unamended soil, but successfully predicted growth responses to resource pulses: growth increased with 16S rRNA gene copy number and declined with genome size after substrate addition to soils, responses that were repeated in four different ecosystems. Genome size best predicted growth rate in response to addition of glucose alone; adding ammonium with glucose weakened the relationship, and the relationship was absent in nutrient-replete pure cultures, consistent with the idea that reduced genome size is a mechanism of nutrient conservation. Our findings demonstrate that genomic traits of soil bacteria can map to their ecological performance in nature, but the mapping is poor under native soil conditions, where genomic traits related to stress tolerance may prove more predictive. These results remind that phenotype depends on environmental context, underscoring the importance of verifying proposed schemes of trait-based strategies through direct measurement of performance in nature, an important and currently missing foundation for translating microbial processes from genes to ecosystems.

Niche-directed evolution modulates genome architecture in freshwater Planctomycetes

Freshwater environments teem with microbes that do not have counterparts in culture collections or genetic data available in genomic repositories. Currently, our apprehension of evolutionary ecology of freshwater bacteria is hampered by the difficulty to establish organism models for the most representative clades. To circumvent the bottlenecks inherent to the cultivation-based techniques, we applied ecogenomics approaches in order to unravel the evolutionary history and the processes that drive genome architecture in hallmark freshwater lineages from the phylum Planctomycetes. The evolutionary history inferences showed that sediment/soil Planctomycetes transitioned to aquatic environments, where they gave rise to new freshwater-specific clades. The most abundant lineage was found to have the most specialised lifestyle (increased regulatory genetic circuits, metabolism tuned for mineralization of proteinaceous sinking aggregates, psychrotrophic behaviour) within the analysed clades and to harbour the smallest freshwater Planctomycetes genomes, highlighting a genomic architecture shaped by niche-directed evolution (through loss of functions and pathways not needed in the newly acquired freshwater niche).

RefSoil+: a Reference Database for Genes and Traits of Soil Plasmids

Soil-associated plasmids have the potential to transfer antibiotic resistance genes from environmental to clinical microbial strains, which is a public health concern. A specific resource is needed to aggregate the knowledge of soil plasmid characteristics so that the content, host associations, and dynamics of antibiotic resistance genes can be assessed and then tracked between the environment and the clinic. Here, we present RefSoil+, a database of soil-associated plasmids. RefSoil+ presents a contemporary snapshot of antibiotic resistance genes in soil that can serve as a reference as novel plasmids and transferred antibiotic resistances are discovered. Our study broadens our understanding of plasmids in soil and provides a community resource of important plasmid-associated genes, including antibiotic resistance genes.

Plasmids harbor transferable genes that contribute to the functional repertoire of microbial communities, yet their contributions to metagenomes are often overlooked. Environmental plasmids have the potential to spread antibiotic resistance to clinical microbial strains. In soils, high microbiome diversity and high variability in plasmid characteristics present a challenge for studying plasmids. To improve the understanding of soil plasmids, we present RefSoil+, a database containing plasmid sequences from 922 soil microorganisms. Soil plasmids were larger than other described plasmids, which is a trait associated with plasmid mobility. There was a weak relationship between chromosome size and plasmid size and no relationship between chromosome size and plasmid number, suggesting that these genomic traits are independent in soil. We used RefSoil+ to inform the distributions of antibiotic resistance genes among soil microorganisms compared to those among nonsoil microorganisms. Soil-associated plasmids, but not chromosomes, had fewer antibiotic resistance genes than other microorganisms. These data suggest that soils may offer limited opportunity for plasmid-mediated transfer of described antibiotic resistance genes. RefSoil+ can serve as a reference for the diversity, composition, and host associations of plasmid-borne functional genes in soil, a utility that will be enhanced as the database expands. Our study improves the understanding of soil plasmids and provides a resource for assessing the dynamics of the genes that they carry, especially genes conferring antibiotic resistances.

IMPORTANCE Soil-associated plasmids have the potential to transfer antibiotic resistance genes from environmental to clinical microbial strains, which is a public health concern. A specific resource is needed to aggregate the knowledge of soil plasmid characteristics so that the content, host associations, and dynamics of antibiotic resistance genes can be assessed and then tracked between the environment and the clinic. Here, we present RefSoil+, a database of soil-associated plasmids. RefSoil+ presents a contemporary snapshot of antibiotic resistance genes in soil that can serve as a reference as novel plasmids and transferred antibiotic resistances are discovered. Our study broadens our understanding of plasmids in soil and provides a community resource of important plasmid-associated genes, including antibiotic resistance genes.