Resting for viability: Gordonia polyisoprenivorans ZM27, a robust generalist for petroleum bioremediation under hypersaline stress

The intrinsic issue associated with the application of microbes for practical pollution remediation involves maintaining the expected activity of engaged strains or consortiums as effectively as that noted under laboratory conditions. Faced with various stress factors, degraders with dormancy ability are more likely to survive and exhibit degradation activity. In this study, a hydrocarbonoclastic and halotolerant strain, Gordonia polyisoprenivorans ZM27, was isolated via stimulation with resuscitation-promoting factor (Rpf). Long-term exposure to dual stresses of 10% NaCl and starvation induced ZM27 to enter a viable but nonculturable (VBNC)-like state, and ZM27 cells could be resuscitated upon Rpf stimulation. Notable changes in both morphological and physiological characteristics between VBNC-like ZM27 cells and resuscitated cells confirmed the response to Rpf and their robust resistance against harsh environments. Whole-genome sequencing and analysis indicated ZM27 could be a generalist degrader with dormancy ability. Subsequently, VBNC-like ZM27 was applied in a soil microcosm experiment to investigate the practical application potential under harsh conditions. VBNC-like ZM27 combined with Rpf stimulation exhibited the most effective biodegradation performance, and the initial n-hexadecane content (1000 mg kg-1) decreased by 63.29% after 14-day incubation. Based on 16S rRNA amplicon sequencing and analysis, Gordonia exhibited a positive response to Rpf stimulation. The relative abundance of genus Gordonia was negatively correlated with that of Alcanivorax, a genus of obligate hydrocarbon degrader with the greatest abundance during soil incubation. Based on the degradation profile and community analysis, generalist Gordonia may be more efficient in hydrocarbon degradation than specialist Alcanivorax under harsh conditions. The characteristics of ZM27, including its sustainable culturability under long-term stress, response to Rpf and robust performance in soil microcosms, are valuable for the remediation of petroleum pollution under stressful conditions. Our work validated the importance of dormancy and highlighted the underestimated role of low-activity degraders in petroleum remediation.

Sterile sentinels and MinION sequencing capture active soil microbial communities that differentiate crop rotations

BACKGROUND

Soil microbial communities are difficult to measure and critical to soil processes. The bulk soil microbiome is highly diverse and spatially heterogeneous, which can make it difficult to detect and monitor the responses of microbial communities to differences or changes in management, such as different crop rotations in agricultural research. Sampling a subset of actively growing microbes should promote monitoring how soil microbial communities respond to management by reducing the variation contributed by high microbial spatial and temporal heterogeneity and less active microbes. We tested an in-growth bag method using sterilized soil in root-excluding mesh, "sterile sentinels," for the capacity to differentiate between crop rotations. We assessed the utility of different incubation times and compared colonized sentinels to concurrently sampled bulk soils for the statistical power to differentiate microbial community composition in low and high diversity crop rotations. We paired this method with Oxford Nanopore MinION sequencing to assess sterile sentinels as a standardized, fast turn-around monitoring method.

RESULTS

Compared to bulk soil, sentinels provided greater statistical power to distinguish between crop rotations for bacterial communities and equivalent power for fungal communities. The incubation time did not affect the statistical power to detect treatment differences in community composition, although longer incubation time increased total biomass. Bulk and sentinel soil samples contained shared and unique microbial taxa that were differentially abundant between crop rotations.

CONCLUSIONS

Overall, compared to bulk soils, the sentinels captured taxa with copiotrophic or ruderal traits, and plant-associated taxa. The sentinels show promise as a sensitive, scalable method to monitor soil microbial communities and provide information complementary to traditional soil sampling.

r/K selection of GC content in prokaryotes

The guanine/cytosine (GC) content of prokaryotic genomes is species-specific, taking values from 16% to 77%. This diversity of selection for GC content remains contentious. We analyse the correlations between GC content and a range of phenotypic and genotypic data in thousands of prokaryotes. GC content integrates well with these traits into r/K selection theory when phenotypic plasticity is considered. High GC-content prokaryotes are r-strategists with cheaper descendants thanks to a lower average amino acid metabolic cost, colonize unstable environments thanks to flagella and a bacillus form and are generalists in terms of resource opportunism and their defence mechanisms. Low GC content prokaryotes are K-strategists specialized for stable environments that maintain homeostasis via a high-cost outer cell membrane and endospore formation as a response to nutrient deprivation, and attain a higher nutrient-to-biomass yield. The lower proteome cost of high GC content prokaryotes is driven by the association between GC-rich codons and cheaper amino acids in the genetic code, while the correlation between GC content and genome size may be partly due to functional diversity driven by r/K selection. In all, molecular diversity in the GC content of prokaryotes may be a consequence of ecological r/K selection.

Scaling of Protein Function across the Tree of Life

Scaling laws are a powerful way to compare genomes because they put all organisms onto a single curve and reveal nontrivial generalities as genomes change in size. The abundance of functional categories across genomes has previously been found to show power law scaling with respect to the total number of functional categories, suggesting that universal constraints shape genomic category abundance. Here, we look across the tree of life to understand how genome evolution may be related to functional scaling. We revisit previous observations of functional genome scaling with an expanded taxonomy by analyzing 3,726 bacterial, 220 archaeal, and 79 unicellular eukaryotic genomes. We find that for some functional classes, scaling is best described by multiple exponents, revealing previously unobserved shifts in scaling as genome-encoded protein annotations increase or decrease. Furthermore, we find that scaling varies between phyletic groups at both the domain and phyla levels and is less universal than previously thought. This variability in functional scaling is not related to taxonomic phylogeny resolved at the phyla level, suggesting that differences in cell plan or physiology outweigh broad patterns of taxonomic evolution. Since genomes are maintained and replicated by the functional proteins encoded by them, these results point to functional degeneracy between taxonomic groups and unique evolutionary trajectories toward these. We also find that individual phyla frequently span scaling exponents of functional classes, revealing that individual clades can move across scaling exponents. Together, our results reveal unique shifts in functions across the tree of life and highlight that as genomes grow or shrink, proteins of various functions may be added or lost.

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.

Bacterial enzymatic degradation of recalcitrant organic pollutants: catabolic pathways and genetic regulations

Contamination of soil and natural water bodies driven by increased organic pollutants remains a universal concern. Naturally, organic pollutants contain carcinogenic and toxic properties threatening all known life forms. The conventional physical and chemical methods employed to remove these organic pollutants ironically produce toxic and non-ecofriendly end-products. Whereas microbial-based degradation of organic pollutants provides an edge, they are usually cost-effective and take an eco-friendly approach towards remediation. Bacterial species, including Pseudomonas, Comamonas, Burkholderia, and Xanthomonas, have the unique genetic makeup to metabolically degrade toxic pollutants, conferring their survival in toxic environments. Several catabolic genes, such as alkB, xylE, catA, and nahAc, that encode enzymes and allow bacteria to degrade organic pollutants have been identified, characterized, and even engineered for better efficacy. Aerobic and anaerobic processes are followed by bacteria to metabolize aliphatic saturated and unsaturated hydrocarbons such as alkanes, cycloalkanes, aldehydes, and ethers. Bacteria use a variety of degrading pathways, including catechol, protocatechuate, gentisate, benzoate, and biphenyl, to remove aromatic organic contaminants such as polychlorinated biphenyls, polycyclic aromatic hydrocarbons, and pesticides from the environment. A better understanding of the principle, mechanisms, and genetics would be beneficial for improving the metabolic efficacy of bacteria to such ends. With a focus on comprehending the mechanisms involved in various catabolic pathways and the genetics of the biotransformation of these xenobiotic compounds, the present review offers insight into the various sources and types of known organic pollutants and their toxic effects on health and the environment.

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.

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.

Ecological insights into soil health according to the genomic traits and environment-wide associations of bacteria in agricultural soils

Soil microbiomes are sensitive to current and previous soil conditions, and bacterial 'bioindicators' of biological, physical, and chemical soil properties have considerable potential for soil health assessment. However, the lack of ecological or physiological information for most soil microorganisms limits our ability to interpret the associations of bioindicators and, thus, their utility for guiding management. We identified bioindicators of tillage intensity and twelve soil properties used to rate soil health using a 16S rRNA gene-based survey of farmland across North America. We then inferred the genomic traits of bioindicators and evaluated their environment-wide associations (EWAS) with respect to agricultural management practice, disturbance, and plant associations with 89 studies from agroecosystems. Most bioindicators were either positively correlated with biological properties (e.g., organic matter) or negatively correlated with physical and chemical properties. Higher soil health ratings corresponded with smaller genome size and higher coding density, while lower ratings corresponded with larger genomes and higher rrn copy number. Community-weighted genome size explained most variation in health ratings. EWAS linked prominent bioindicators with the impacts of environmental disturbances. Our findings provide ecological insights into bioindicators of soil properties relevant to soil health management, illustrating the tight coupling of microbiome and soil function.

Genome size distributions in bacteria and archaea are strongly linked to evolutionary history at broad phylogenetic scales

The evolutionary forces that determine genome size in bacteria and archaea have been the subject of intense debate over the last few decades. Although the preferential loss of genes observed in prokaryotes is explained through the deletional bias, factors promoting and preventing the fixation of such gene losses often remain unclear. Importantly, statistical analyses on this topic typically do not consider the potential bias introduced by the shared ancestry of many lineages, which is critical when using species as data points because of the potential dependence on residuals. In this study, we investigated the genome size distributions across a broad diversity of bacteria and archaea to evaluate if this trait is phylogenetically conserved at broad phylogenetic scales. After model fit, Pagel’s lambda indicated a strong phylogenetic signal in genome size data, suggesting that the diversification of this trait is influenced by shared evolutionary histories. We used a phylogenetic generalized least-squares analysis (PGLS) to test whether phylogeny influences the predictability of genome size from dN/dS ratios and 16S copy number, two variables that have been previously linked to genome size. These results confirm that failure to account for evolutionary history can lead to biased interpretations of genome size predictors. Overall, our results indicate that although bacteria and archaea can rapidly gain and lose genetic material through gene transfers and deletions, respectively, phylogenetic signal for genome size distributions can still be recovered at broad phylogenetic scales that should be taken into account when inferring the drivers of genome size evolution.

The evolutionary forces driving genome size in bacteria and archaea have been subject to debate during the last decades. Typically, independent comparative analyses have suggested that unique variables, such as the strength of selection, environmental complexity, and mutation rate, are the main drivers of this trait, without considering for potential biases derived from shared ancestry. Here, we applied a phylogeny-based statistical approach to assess how tightly genome size in bacteria and archaea is linked to evolutionary history. Moreover, we also evaluated the predictability of genome size from the strength of purifying selection and ecological strategy on a broad diversity of bacteria and archaea genomes under a phylogenetic comparative framework. Our approach indicates that despite the ability of bacteria and archaea to rapidly exchange genes, a strong phylogenetic signal to genome size distributions can be recovered at broad phylogenetic scales.

Phylogeny, distribution and potential metabolism of candidate bacterial phylum KSB1

Candidate phylum KSB1 is composed of uncultured bacteria and has been reported across various environments. However, the phylogeny and metabolic potential of KSB1 have not been studied comprehensively. In this study, phylogenomic analysis of KSB1 genomes from public databases and eleven metagenome-assembled genomes (MAGs) from marine and hydrothermal sediments revealed that those genomes were clustered into four clades. Isolation source and relative abundance of KSB1 genomes showed that clade I was particularly abundant in bioreactor sludge. Genes related to dissimilatory reduction of nitrate to ammonia (DNRA), the last step of denitrification converting nitrous oxide to nitrogen and assimilatory sulfur reduction were observed in the expanded genomes of clade I, which may due to horizontal gene transfer that frequently occurred in bioreactor. Annotation and metabolic reconstruction of clades II and IV showed flagellum assembly and chemotaxis genes in the genomes, which may indicate that exploration and sensing for nutrients and chemical gradients are critical for the two clades in deep-sea and hydrothermal sediment. Metabolic potentials of fatty acids and short-chain hydrocarbons utilization were predicted in clades I and IV of KSB1. Collectively, phylogenomic and metabolic analyses of KSB1 clades provide insight into their anaerobic heterotrophic lifestyle and differentiation in potential ecological roles.

A Large-Scale Genome-Based Survey of Acidophilic Bacteria Suggests That Genome Streamlining Is an Adaption for Life at Low pH

The genome streamlining theory suggests that reduction of microbial genome size optimizes energy utilization in stressful environments. Although this hypothesis has been explored in several cases of low-nutrient (oligotrophic) and high-temperature environments, little work has been carried out on microorganisms from low-pH environments, and what has been reported is inconclusive. In this study, we performed a large-scale comparative genomics investigation of more than 260 bacterial high-quality genome sequences of acidophiles, together with genomes of their closest phylogenetic relatives that live at circum-neutral pH. A statistically supported correlation is reported between reduction of genome size and decreasing pH that we demonstrate is due to gene loss and reduced gene sizes. This trend is independent from other genome size constraints such as temperature and G + C content. Genome streamlining in the evolution of acidophilic bacteria is thus supported by our results. The analyses of predicted Clusters of Orthologous Genes (COG) categories and subcellular location predictions indicate that acidophiles have a lower representation of genes encoding extracellular proteins, signal transduction mechanisms, and proteins with unknown function but are enriched in inner membrane proteins, chaperones, basic metabolism, and core cellular functions. Contrary to other reports for genome streamlining, there was no significant change in paralog frequencies across pH. However, a detailed analysis of COG categories revealed a higher proportion of genes in acidophiles in the following categories: "replication and repair," "amino acid transport," and "intracellular trafficking". This study brings increasing clarity regarding the genomic adaptations of acidophiles to life at low pH while putting elements, such as the reduction of average gene size, under the spotlight of streamlining theory.

Strategic traits of bacteria and archaea vary widely within substrate-use groups

Quantitative traits such as maximum growth rate and cell radial diameter are one facet of ecological strategy variation across bacteria and archaea. Another facet is substrate-use pathways, such as iron reduction or methylotrophy. Here, we ask how these two facets intersect, using a large compilation of data for culturable species and examining seven quantitative traits (genome size, signal transduction protein count, histidine kinase count, growth temperature, temperature-adjusted maximum growth rate, cell radial diameter and 16S rRNA operon copy number). Overall, quantitative trait variation within groups of organisms possessing a particular substrate-use pathway was very broad, outweighing differences between substrate-use groups. Although some substrate-use groups had significantly different means for some quantitative traits, standard deviation of quantitative trait values within each substrate-use pathway mostly averaged between 1.6 and 1.8 times larger than standard deviation across group means. Most likely, this wide variation reflects ecological strategy: for example, fast maximum growth rate is likely to express an early successional or copiotrophic strategy, and maximum growth varies widely within most substrate-use pathways. In general, it appears that these quantitative traits express different and complementary information about ecological strategy, compared with substrate use.

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.

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.