Relating genomic characteristics to environmental preferences and ubiquity in different microbial taxa

Cross-biome microbial networks reveal functional redundancy and suggest genome reduction through functional complementarity

The structure of microbial communities arises from a multitude of factors, including the interactions of microorganisms with each other and with the environment. In this work, we sought to disentangle those drivers by performing a cross-study, cross-biome meta-analysis of microbial occurrence data in more than 5000 samples, applying a novel network clustering algorithm aimed to capture conditional taxa co-occurrences. We then examined the phylogenetic and functional composition of the resulting clusters, and searched for global patterns of assembly both at the community level and in the presence/absence of individual metabolic pathways.Our analysis highlighted the prevalence of functional redundancy in microbial communities, particularly between taxa that co-occur in more than one environment, pointing to a relationship between functional redundancy and environmental adaptation. In spite of this, certain pathways were observed in fewer taxa than expected by chance, suggesting the presence of auxotrophy, and presumably cooperation among community members. This hypothetical cooperation may play a role in genome reduction, since we observed a negative relationship between the size of bacterial genomes and the size of the community they belong to.Overall, our results suggest the microbial community assembly is driven by universal principles that operate consistently across different biomes and taxonomic groups.

A systematic comparison of natural product potential, with an emphasis on RiPPs, by mining of bacteria of three large ecosystems

The implementation of several global microbiome studies has yielded extensive insights into the biosynthetic potential of natural microbial communities. However, studies on the distribution of several classes of ribosomally synthesized and post-translationally modified peptides (RiPPs), non-ribosomal peptides (NRPs) and polyketides (PKs) in different large microbial ecosystems have been very limited. Here, we collected a large set of metagenome-assembled bacterial genomes from marine, freshwater and terrestrial ecosystems to investigate the biosynthetic potential of these bacteria. We demonstrate the utility of public dataset collections for revealing the different secondary metabolite biosynthetic potentials among these different living environments. We show that there is a higher occurrence of RiPPs in terrestrial systems, while in marine systems, we found relatively more terpene-, NRP-, and PK encoding gene clusters. Among the many new biosynthetic gene clusters (BGCs) identified, we analyzed various Nif-11-like and nitrile hydratase leader peptide (NHLP) containing gene clusters that would merit further study, including promising products, such as mersacidin-, LAP- and proteusin analogs. This research highlights the significance of public datasets in elucidating the biosynthetic potential of microbes in different living environments and underscores the wide bioengineering opportunities within the RiPP family.

A global survey of host, aquatic, and soil microbiomes reveals shared abundance and genomic features between bacterial and fungal generalists

Environmental change, coupled with alteration in human lifestyles, is profoundly impacting the microbial communities critical to the health of the Earth and its inhabitants. To identify bacteria and fungi that are resistant and susceptible to habitat change, we analyze thousands of genera detected in 1,580 host, soil, and aquatic samples. This large-scale analysis identifies 48 bacterial and 4 fungal genera that are abundant across the three biomes, demonstrating fitness in diverse environmental conditions. Samples containing these generalists have significantly higher alpha diversity. These generalists play a significant role in shaping cross-kingdom community structure, boasting larger genomes with more secondary metabolism and antimicrobial resistance genes. Conversely, 30 bacterial and 19 fungal genera are only found in a single habitat, suggesting a limited ability to adapt to different and changing environments. These findings contribute to our understanding of microbial niche breadth and its consequences for global biodiversity loss.

Taxonomic and environmental distribution of bacterial amino acid auxotrophies

Many microorganisms are auxotrophic-unable to synthesize the compounds they require for growth. With this work, we quantify the prevalence of amino acid auxotrophies across a broad diversity of bacteria and habitats. We predicted the amino acid biosynthetic capabilities of 26,277 unique bacterial genomes spanning 12 phyla using a metabolic pathway model validated with empirical data. Amino acid auxotrophy is widespread across bacterial phyla, but we conservatively estimate that the majority of taxa (78.4%) are able to synthesize all amino acids. Our estimates indicate that amino acid auxotrophies are more prevalent among obligate intracellular parasites and in free-living taxa with genomic attributes characteristic of 'streamlined' life history strategies. We predicted the amino acid biosynthetic capabilities of bacterial communities found in 12 unique habitats to investigate environmental associations with auxotrophy, using data compiled from 3813 samples spanning major aquatic, terrestrial, and engineered environments. Auxotrophic taxa were more abundant in host-associated environments (including the human oral cavity and gut) and in fermented food products, with auxotrophic taxa being relatively rare in soil and aquatic systems. Overall, this work contributes to a more complete understanding of amino acid auxotrophy across the bacterial tree of life and the ecological contexts in which auxotrophy can be a successful strategy.

Ribosomal RNA gene operon copy number, a functional trait indicating the hydrocarbon degradation level of bacterial communities

The lack of universal indicators for predicting microbial biodegradation potential and assessing remediation effects limits the generalization of bioremediation. The community-level ribosomal RNA gene operon (rrn) copy number, an important functional trait, has the potential to serve as a key indicator of the bioremediation of organic pollutants. A meta-analysis based on 1275 samples from 26 hydrocarbon-related studies revealed a positive relationship between the microbial hydrocarbon biodegradation level and the community-level rrn copy number in soil, seawater and culture. Subsequently, a microcosm experiment was performed to decipher the community-level rrn copy number response mechanism during total petroleum hydrocarbon (TPH) biodegradation. The treatment combining straw with resuscitation-promoting factor (Rpf) exhibited the highest community-level rrn copy number and the most effective biodegradation compared with other treatments, and the initial TPH content (20,000 mg kg-1) was reduced by 67.67% after 77 days of incubation. TPH biodegradation rate was positively correlated with the average community-level rrn copy number (p = 0.001, R2 = 0.5781). Both meta and community analyses showed that rrn copy number may reflect the potential of hydrocarbon degradation and microbial dormancy. Our findings provide insight into the applicability of the community-level rrn copy number to assess bacterial biodegradation for pollution remediation.

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.

Genomic and phenotypic features of Acinetobacter baumannii isolated from oil reservoirs reveal a novel subspecies specialized in degrading hazardous hydrocarbons

The genus Acinetobacter encompasses biotechnologically relevant species and nosocomial pathogens. In this study, nine isolates recovered from different oil reservoir samples showed the ability to grow with petroleum as the only carbon source and possessed the ability to emulsify kerosene. The whole genomes of the nine strains were sequenced and analyzed. The average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) values of all strains were compared to the reference strains, and the results were below the reference values (<97.88 and 82, respectively), suggesting that the isolates belong to a new subspecies of Acinetobacter baumannii. The name Acinetobacter baumannii oleum ficedula is proposed. A comparison of the whole genome repertoire of 290 Acinetobacter species indicated that the strains in this study resemble non-pathogenic Acinetobacter strains. However, the new isolates resemble A. baumannii when comparing virulence factors. The isolates in this study carry many genes involved in hydrocarbon degradation, indicating the potential to degrade most toxic compounds listed by environmental regulatory agencies such as ATSDR, EPA, and CONAMA. In addition, despite the absence of known biosurfactant or bioemulsifier genes, the strains showed emulsifying activity, suggesting the presence of new pathways or genes related to this process. This study investigated the genomic, phenotypic, and biochemical features of the novel environmental subspecies A. baumannii oleum ficedula, revealing their potential to degrade hydrocarbons and to produce biosurfactants or bioemulsifiers. Applying these environmental subspecies in bioaugmentation strategies sheds light on future approaches to bioremediation. The study shows the importance of genomic analysis of environmental strains and their inclusion in metabolic pathways databases, highlighting unique enzymes/alternative pathways for consuming hazardous hydrocarbons.

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.

Exploring the Interspecific Interactions and the Metabolome of the Soil Isolate Hylemonella gracilis

Microbial community analysis of aquatic environments showed that an important component of its microbial diversity consists of bacteria with cell sizes of ~0.1 μm. Such small bacteria can show genomic reductions and metabolic dependencies with other bacteria. However, so far, no study has investigated if such bacteria exist in terrestrial environments like soil. Here, we isolated soil bacteria that passed through a 0.1-μm filter. The complete genome of one of the isolates was sequenced and the bacterium was identified as Hylemonella gracilis. A set of coculture assays with phylogenetically distant soil bacteria with different cell and genome sizes was performed. The coculture assays revealed that H. gracilis grows better when interacting with other soil bacteria like Paenibacillus sp. AD87 and Serratia plymuthica. Transcriptomics and metabolomics showed that H. gracilis was able to change gene expression, behavior, and biochemistry of the interacting bacteria without direct cell-cell contact. Our study indicates that in soil there are bacteria that can pass through a 0.1-μm filter. These bacteria may have been overlooked in previous research on soil microbial communities. Such small bacteria, exemplified here by H. gracilis, can induce transcriptional and metabolomic changes in other bacteria upon their interactions in soil. In vitro, the studied interspecific interactions allowed utilization of growth substrates that could not be utilized by monocultures, suggesting that biochemical interactions between substantially different sized soil bacteria may contribute to the symbiosis of soil bacterial communities. IMPORTANCE Analysis of aquatic microbial communities revealed that parts of its diversity consist of bacteria with cell sizes of ~0.1 μm. Such bacteria can show genomic reductions and metabolic dependencies with other bacteria. So far, no study investigated if such bacteria exist in terrestrial environments such as soil. Here, we show that such bacteria also exist in soil. The isolated bacteria were identified as Hylemonella gracilis. Coculture assays with phylogenetically different soil bacteria revealed that H. gracilis grows better when cocultured with other soil bacteria. Transcriptomics and metabolomics showed that H. gracilis was able to change gene expression, behavior, and biochemistry of the interacting bacteria without direct contact. Our study revealed that bacteria are present in soil that can pass through 0.1-μm filters. Such bacteria may have been overlooked in previous research on soil microbial communities and may contribute to the symbiosis of soil bacterial communities.

Mixotrophy broadens the ecological niche range of the iron oxidizer Sideroxydans sp. CL21 isolated from an iron-rich peatland

Sideroxydans sp. CL21 is a microaerobic, acid-tolerant Fe(II)-oxidizer, isolated from the Schlöppnerbrunnen fen. Since the genome size of Sideroxydans sp. CL21 is 21% larger than that of the neutrophilic Sideroxydans lithotrophicus ES-1, we hypothesized that strain CL21 contains additional metabolic traits to thrive in the fen. The common genomic content of both strains contains homologs of the putative Fe(II) oxidation genes, mtoAB and cyc2. A large part of the accessory genome in strain CL21 contains genes linked to utilization of alternative electron donors, including NiFe uptake hydrogenases, and genes encoding lactate uptake and utilization proteins, motility and biofilm formation, transposable elements, and pH homeostasis mechanisms. Next, we incubated the strain in different combinations of electron donors and characterized the fen microbial communities. Sideroxydans spp. comprised 3.33% and 3.94% of the total relative abundance in the peatland soil and peatland water, respectively. Incubation results indicate Sideroxydans sp. CL21 uses H2 and thiosulfate, while lactate only enhances growth when combined with Fe, H2, or thiosulfate. Rates of H2 utilization were highest in combination with other substrates. Thus, Sideroxydans sp. CL21 is a mixotroph, growing best by simultaneously using substrate combinations, which helps to thrive in dynamic and complex habitats.

Trait-trait relationships and tradeoffs vary with genome size in prokaryotes

We report genomic traits that have been associated with the life history of prokaryotes and highlight conflicting findings concerning earlier observed trait correlations and tradeoffs. In order to address possible explanations for these contradictions we examined trait-trait variations of 11 genomic traits from ~18,000 sequenced genomes. The studied trait-trait variations suggested: (i) the predominance of two resistance and resilience-related orthogonal axes and (ii) at least in free living species with large effective population sizes whose evolution is little affected by genetic drift an overlap between a resilience axis and an oligotrophic-copiotrophic axis. These findings imply that resistance associated traits of prokaryotes are globally decoupled from resilience related traits and in the case of free-living communities also from traits associated with resource availability. However, further inspection of pairwise scatterplots showed that resistance and resilience traits tended to be positively related for genomes up to roughly five million base pairs and negatively for larger genomes. Genome size distributions differ across habitats and our findings therefore point to habitat dependent tradeoffs between resistance and resilience. This in turn may preclude a globally consistent assignment of prokaryote genomic traits to the competitor - stress-tolerator - ruderal (CSR) schema that sorts species depending on their location along disturbance and productivity gradients into three ecological strategies and may serve as an explanation for conflicting findings from earlier studies. All reviewed genomic traits featured significant phylogenetic signals and we propose that our trait table can be applied to extrapolate genomic traits from taxonomic marker genes. This will enable to empirically evaluate the assembly of these genomic traits in prokaryotic communities from different habitats and under different productivity and disturbance scenarios as predicted via the resistance-resilience framework formulated here.

Macroecological distributions of gene variants highlight the functional organization of soil microbial systems

The recent application of macroecological tools and concepts has made it possible to identify consistent patterns in the distribution of microbial biodiversity, which greatly improved our understanding of the microbial world at large scales. However, the distribution of microbial functions remains largely uncharted from the macroecological point of view. Here, we used macroecological models to examine how the genes encoding the functional capabilities of microorganisms are distributed within and across soil systems. Models built using functional gene array data from 818 soil microbial communities showed that the occupancy-frequency distributions of genes were bimodal in every studied site, and that their rank-abundance distributions were best described by a lognormal model. In addition, the relationships between gene occupancy and abundance were positive in all sites. This allowed us to identify genes with high abundance and ubiquitous distribution (core) and genes with low abundance and limited spatial distribution (satellites), and to show that they encode different sets of microbial traits. Common genes encode microbial traits related to the main biogeochemical cycles (C, N, P and S) while rare genes encode traits related to adaptation to environmental stresses, such as nutrient limitation, resistance to heavy metals and degradation of xenobiotics. Overall, this study characterized for the first time the distribution of microbial functional genes within soil systems, and highlight the interest of macroecological models for understanding the functional organization of microbial systems across spatial scales.

Complete Genome Sequencing and Comparative Genomics of Three Potential Probiotic Strains, Lacticaseibacillus casei FBL6, Lacticaseibacillus chiayiensis FBL7, and Lacticaseibacillus zeae FBL8

Lacticaseibacillus casei, Lacticaseibacillus chiayiensis, and Lacticaseibacillus zeae are very closely related Lacticaseibacillus species. L. casei has long been proposed as a probiotic, whereas studies on functional characterization for L. chiayiensis and L. zeae are some compared to L. casei. In this study, L. casei FBL6, L. chiayiensis FBL7, and L. zeae FBL8 were isolated from raw milk, and their probiotic properties were investigated. Genomic analysis demonstrated the role of L. chiayiensis and L. zeae as probiotic candidates. The three strains were tolerant to acid and bile salt, with inhibitory action against pathogenic bacterial strains and capacity of antioxidants. Complete genome sequences of the three strains were analyzed to highlight the probiotic properties at the genetic level, which results in the discovery of genes corresponding to phenotypic characterization. Moreover, genes known to confer probiotic characteristics were identified, including genes related to biosynthesis, defense machinery, adhesion, and stress adaptation. The comparative genomic analysis with other available genomes revealed 256, 214, and 32 unique genes for FBL6, FBL7, and FBL8, respectively. These genomes contained individual genes encoding proteins that are putatively involved in carbohydrate transport and metabolism, prokaryotic immune system for antiviral defense, and physiological control processes. In particular, L. casei FBL6 had a bacteriocin gene cluster that was not present in other genomes of L. casei, resulting in this strain may exhibit a wide range of antimicrobial activity compared to other L. casei strains. Our data can help us understand the probiotic functionalities of the three strains and suggest that L. chiayiensis and L. zeae species, which are closely related to L. casei, can also be considered as novel potential probiotic candidate strains.

The genus Caulobacter and its role in plant microbiomes

Recent omics approaches have revealed the prevalent microbial taxa that constitute the microbiome of various plant species. Across global scales and environmental conditions, strains belonging to the bacterial genus Caulobacter have consistently been found in association with various plant species. Aligned with agroecological relevance and biotechnological advances, many scientific communications have demonstrated that several Caulobacter strains (spanning several Caulobacter species) harbor the potential to enhance plant biomass for various plant species ranging from Arabidopsis to Citrullus and Zea mays. In the past several years, co-occurrence data have driven mechanistically resolved communications about select Caulobacter-plant interactions. Given the long-standing history of Caulobacter as a model organism for cell cycle regulation, genetic studies, and the prevalence of Caulobacter species in various plant microbiomes, the genus Caulobacter offers researchers a unique opportunity to leverage for investigating plant-microbe interactions and realizing targeted biotechnological applications. In this review, recent developments regarding Caulobacter-plant interactions are presented in terms of model utility for future biotechnological investigations.

A Genomic Perspective Across Earth's Microbiomes Reveals That Genome Size in Archaea and Bacteria Is Linked to Ecosystem Type and Trophic Strategy

Our view of genome size in Archaea and Bacteria has remained skewed as the data has been dominated by genomes of microorganisms that have been cultivated under laboratory settings. However, the continuous effort to catalog Earth's microbiomes, specifically propelled by recent extensive work on uncultivated microorganisms, provides an opportunity to revise our perspective on genome size distribution. We present a meta-analysis that includes 26,101 representative genomes from 3 published genomic databases; metagenomic assembled genomes (MAGs) from GEMs and stratfreshDB, and isolates from GTDB. Aquatic and host-associated microbial genomes present on average the smallest estimated genome sizes (3.1 and 3.0 Mbp, respectively). These are followed by terrestrial microbial genomes (average 3.7 Mbp), and genomes from isolated microorganisms (average 4.3 Mbp). On the one hand, aquatic and host-associated ecosystems present smaller genomes sizes in genera of phyla with genome sizes above 3 Mbp. On the other hand, estimated genome size in phyla with genomes under 3 Mbp showed no difference between ecosystems. Moreover, we observed that when using 95% average nucleotide identity (ANI) as an estimator for genetic units, only 3% of MAGs cluster together with genomes from isolated microorganisms. Although there are potential methodological limitations when assembling and binning MAGs, we found that in genome clusters containing both environmental MAGs and isolate genomes, MAGs were estimated only an average 3.7% smaller than isolate genomes. Even when assembly and binning methods introduce biases, estimated genome size of MAGs and isolates are very similar. Finally, to better understand the ecological drivers of genome size, we discuss on the known and the overlooked factors that influence genome size in different ecosystems, phylogenetic groups, and trophic 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.

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.

Exploring Microbiome Functional Dynamics through Space and Time with Trait-Based Theory

Microbiomes play essential roles in the health and function of animal and plant hosts and drive nutrient cycling across ecosystems. Integrating novel trait-based approaches with ecological theory can facilitate the prediction of microbial functional traits important for ecosystem functioning and health. In particular, the yield-acquisition-stress (Y-A-S) framework considers dominant microbial life history strategies across gradients of resource availability and stress. However, microbiomes are dynamic, and spatial and temporal shifts in taxonomic and trait composition can affect ecosystem functions. We posit that extending the Y-A-S framework to microbiomes during succession and across biogeographic gradients can lead to generalizable rules for how microbiomes and their functions respond to resources and stress across space, time, and diverse ecosystems. We demonstrate the potential of this framework by applying it to the microbiomes hosted by the carnivorous pitcher plant Sarracenia purpurea, which have clear successional trajectories and are distributed across a broad climatic gradient.

Exploring the functional composition of the human microbiome using a hand-curated microbial trait database

BACKGROUND

Even when microbial communities vary wildly in their taxonomic composition, their functional composition is often surprisingly stable. This suggests that a functional perspective could provide much deeper insight into the principles governing microbiome assembly. Much work to date analyzing the functional composition of microbial communities, however, relies heavily on inference from genomic features. Unfortunately, output from these methods can be hard to interpret and often suffers from relatively high error rates.

RESULTS

We built and analyzed a domain-specific microbial trait database from known microbe-trait pairs recorded in the literature to better understand the functional composition of the human microbiome. Using a combination of phylogentically conscious machine learning tools and a network science approach, we were able to link particular traits to areas of the human body, discover traits that determine the range of body areas a microbe can inhabit, and uncover drivers of metabolic breadth.

CONCLUSIONS

Domain-specific trait databases are an effective compromise between noisy methods to infer complex traits from genomic data and exhaustive, expensive attempts at database curation from the literature that do not focus on any one subset of taxa. They provide an accurate account of microbial traits and, by limiting the number of taxa considered, are feasible to build within a reasonable time-frame. We present a database specific for the human microbiome, in the hopes that this will prove useful for research into the functional composition of human-associated microbial communities.

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.

Verrucomicrobial methanotrophs: ecophysiology of metabolically versatile acidophiles

Methanotrophs are an important group of microorganisms that counteract methane emissions to the atmosphere. Methane-oxidising bacteria of the Alpha- and Gammaproteobacteria have been studied for over a century, while methanotrophs of the phylum Verrucomicrobia are a more recent discovery. Verrucomicrobial methanotrophs are extremophiles that live in very acidic geothermal ecosystems. Currently, more than a dozen strains have been isolated, belonging to the genera Methylacidiphilum and Methylacidimicrobium. Initially, these methanotrophs were thought to be metabolically confined. However, genomic analyses and physiological and biochemical experiments over the past years revealed that verrucomicrobial methanotrophs, as well as proteobacterial methanotrophs, are much more metabolically versatile than previously assumed. Several inorganic gases and other molecules present in acidic geothermal ecosystems can be utilised, such as methane, hydrogen gas, carbon dioxide, ammonium, nitrogen gas and perhaps also hydrogen sulfide. Verrucomicrobial methanotrophs could therefore represent key players in multiple volcanic nutrient cycles and in the mitigation of greenhouse gas emissions from geothermal ecosystems. Here, we summarise the current knowledge on verrucomicrobial methanotrophs with respect to their metabolic versatility and discuss the factors that determine their diversity in their natural environment. In addition, key metabolic, morphological and ecological characteristics of verrucomicrobial and proteobacterial methanotrophs are reviewed.

Plasmids encode niche-specific traits in Lactobacillaceae

Species belonging to the family Lactobacillaceae are found in highly diverse environments and play an important role in fermented foods and probiotic products. Many of these species have been individually reported to harbour plasmids that encode important genes. In this study, we performed comparative genomic analysis of publicly available data for 512 plasmids from 282 strains represented by 51 species of this family and correlated the genomic features of plasmids with the ecological niches in which these species are found. Two-thirds of the species had at least one plasmid-harbouring strain. Plasmid abundance and GC content were significantly lower in vertebrate-adapted species as compared to nomadic and free-living species. Hierarchical clustering highlighted the distinct nature of plasmids from the nomadic and free-living species than those from the vertebrate-adapted species. EggNOG-assisted functional annotation revealed that genes associated with transposition, conjugation, DNA repair and recombination, exopolysaccharide production, metal ion transport, toxin–antitoxin system, and stress tolerance were significantly enriched on the plasmids of the nomadic and in some cases nomadic and free-living species. On the other hand, genes related to anaerobic metabolism, ABC transporters and the major facilitator superfamily were overrepresented on the plasmids of the vertebrate-adapted species. These genomic signatures correlate with the comparatively nutrient-depleted, stressful and dynamic environments of nomadic and free-living species and nutrient-rich and anaerobic environments of vertebrate-adapted species. Thus, these results indicate the contribution of the plasmids in the adaptation of lactobacilli to their respective habitats. This study also underlines the potential application of these plasmids in improving the technological and probiotic properties of lactic acid bacteria.

Gene duplication drives genome expansion in a major lineage of Thaumarchaeota

Ammonia-oxidising archaea of the phylum Thaumarchaeota are important organisms in the nitrogen cycle, but the mechanisms driving their radiation into diverse ecosystems remain underexplored. Here, existing thaumarchaeotal genomes are complemented with 12 genomes belonging to the previously under-sampled Nitrososphaerales to investigate the impact of lateral gene transfer (LGT), gene duplication and loss across thaumarchaeotal evolution. We reveal a major role for gene duplication in driving genome expansion subsequent to early LGT. In particular, two large LGT events are identified into Nitrososphaerales and the fate of these gene families is highly lineage-specific, being lost in some descendant lineages, but undergoing extensive duplication in others, suggesting niche-specific roles. Notably, some genes involved in carbohydrate transport or coenzyme metabolism were duplicated, likely facilitating niche specialisation in soils and sediments. Overall, our results suggest that LGT followed by gene duplication drives Nitrososphaerales evolution, highlighting a previously under-appreciated mechanism of genome expansion in archaea.

The Ingenuity of Bacterial Genomes

The genomes of bacteria contain fewer genes and substantially less noncoding DNA than those of eukaryotes, and as a result, they have much less raw material to invent new traits. Yet, bacteria are vastly more taxonomically diverse, numerically abundant, and globally successful in colonizing new habitats compared to eukaryotes. Although bacterial genomes are generally considered to be optimized for efficient growth and rapid adaptation, nonadaptive processes have played a major role in shaping the size, contents, and compact organization of bacterial genomes and have allowed the establishment of deleterious traits that serve as the raw materials for genetic innovation.

Globally Abundant "Candidatus Udaeobacter" Benefits from Release of Antibiotics in Soil and Potentially Performs Trace Gas Scavenging

Verrucomicrobia affiliated with "Candidatus Udaeobacter" belong to the most abundant soil bacteria worldwide. Although the synthesis of antibiotics presumably evolved in soil, and environmental pollution with antimicrobials increases, the impact of these complex molecules on "Ca Udaeobacter" remains to be elucidated. In this study, we demonstrate that "Ca. Udaeobacter" representatives residing in grassland as well as forest soil ecosystems show multidrug resistance and even take advantage of antibiotics release. Soils treated with up to six different antibiotics exhibited a higher "Ca. Udaeobacter" abundance than corresponding controls after 3, 8, and 20 days of incubation. In this context, we provide evidence that "Ca. Udaeobacter" representatives may utilize nutrients which are released due to antibiotic-driven lysis of other soil microbes and thereby reduce energetically expensive synthesis of required biomolecules. Moreover, genomic analysis revealed the presence of genes conferring resistance to multiple classes of antibiotics and indicated that "Ca. Udaeobacter" representatives most likely oxidize the trace gas H₂ to generate energy. This energy might be required for long-term persistence in terrestrial habitats, as already suggested for other dominant soil bacteria. Our study illustrates, for the first time, that globally abundant "Ca. Udaeobacter" benefits from release of antibiotics, which confers advantages over other soil bacteria and represents a so-far overlooked fundamental lifestyle feature of this poorly characterized verrucomicrobial genus. Furthermore, our study suggests that "Ca. Udaeobacter" representatives can utilize H₂ as an alternative electron donor.IMPORTANCE Soil bacteria have been investigated for more than a century, but one of the most dominant terrestrial groups on Earth, "Candidatus Udaeobacter," remains elusive and largely unexplored. Its natural habitat is considered a major reservoir of antibiotics, which directly or indirectly impact phylogenetically diverse microorganisms. Here, we found that "Ca. Udaeobacter" representatives exhibit multidrug resistance and not only evade harmful effects of antimicrobials but even benefit from antibiotic pressure in soil. Therefore, "Ca. Udaeobacter" evidently affects the composition of soil resistomes worldwide and might represent a winner of rising environmental pollution with antimicrobials. In addition, our study indicates that "Ca. Udaeobacter" representatives utilize H₂ and thereby contribute to global hydrogen cycling. The here-reported findings provide insights into elementary lifestyle features of "Ca. Udaeobacter," potentially contributing to its successful global dissemination.

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

Microbial organisms inhabit virtually all environments and encompass a vast biological diversity. The pangenome concept aims to facilitate an understanding of diversity within defined phylogenetic groups. Hence, pangenomes are increasingly used to characterize the strain diversity of prokaryotic species. To understand the interdependence of pangenome features (such as the number of core and accessory genes) and to study the impact of environmental and phylogenetic constraints on the evolution of conspecific strains, we computed pangenomes for 155 phylogenetically diverse species (from ten phyla) using 7,000 high-quality genomes to each of which the respective habitats were assigned. Species habitat ubiquity was associated with several pangenome features. In particular, core-genome size was more important for ubiquity than accessory genome size. In general, environmental preferences had a stronger impact on pangenome evolution than phylogenetic inertia. Environmental preferences explained up to 49% of the variance for pangenome features, compared with 18% by phylogenetic inertia. This observation was robust when the dataset was extended to 10,100 species (59 phyla). The importance of environmental preferences was further accentuated by convergent evolution of pangenome features in a given habitat type across different phylogenetic clades. For example, the soil environment promotes expansion of pangenome size, while host-associated habitats lead to its reduction. Taken together, we explored the global principles of pangenome evolution, quantified the influence of habitat, and phylogenetic inertia on the evolution of pangenomes and identified criteria governing species ubiquity and habitat specificity.

RNA-based spatial community analysis revealed intra-reactor variation and expanded collection of direct interspecies electron transfer microorganisms in anaerobic digestion

Granular activated carbon (GAC) has been shown to mediate direct interspecies electron transfer (DIET) in anaerobic digestion. Adding GAC to up-flow anaerobic sludge bed reactor increased the total biomass slightly from 20.0 to 26.6 gVSS/reactor, and maximum organic removal capacity remarkably from 285 to 1660 mgCOD/L/d. Since GAC occupied 7% of reactor volume (denser than suspended sludge, settled to the reactor bottom), we used a spatial sampling strategy (sludge bed top/mid/bottom layers, and tightly attached GAC-biofilm) and DNA- and RNA-based community analyses. RNA-based analysis demonstrated significant community differences between the non-GAC and GAC-amended reactors (p < 0.05) based on ANOSIM statistical analysis. In comparison, DNA-based analysis showed little community difference between these reactors (p > 0.05). RNA-based analysis revealed active enrichments in GAC-biofilm, including bacteria Geobacter, Syntrophus, Desulfovibrio and Blvii28, and archaea Methanosaeta and Methanospirillum. These are potential electro-active syntrophic microorganisms related with DIET, which expand the previously defined list of DIET microorganisms.

Oligotrophic bacterium Hymenobacter latericoloratus CGMCC 16346 degrades the neonicotinoid imidacloprid in surface water

The intensive and extensive application of imidacloprid in agriculture has resulted in water pollution and risks to aquatic invertebrates. However, pure bacteria remediation of imidacloprid in surface water environments has not been studied. Here, we isolated an imidacloprid-degrading bacterium from a water environment, examined its imidacloprid degradation in pure culture and surface water, sequenced its genome, and compared its Clusters of Orthologous Groups (COG) protein categorization with that for another imidacloprid-degrading bacterium. The isolate was an obligate oligotrophic bacterium, Hymenobacter latericoloratus CGMCC 16346, which degraded imidacloprid via hydroxylation by co-metabolism in pure culture. Resting cells degraded 64.4% of 100 mg/L imidacloprid in 6 days in the presence of co-substrate maltose, and growing culture degraded 40.8% of imidacloprid in 10 days. H. latericoloratus CGMCC 16346 degraded imidacloprid in surface water without co-substrate supplementation and retained imidacloprid-degrading activity after 30 days. The half-life of imidacloprid in surface water was decreased from 173.3 days in the control to 57.8 days by CGMCC 16346 inoculation. Genome sequencing and COG analysis indicated that carbohydrate metabolism and transport, cell wall/membrane biogenesis, and defense mechanisms are enriched in H. latericoloratus CGMCC 16346 compared with the copiotrophic imidacloprid-degrading Pseudoxanthomonas indica CGMCC 6648, indicating that H. latericoloratus CGMCC 16346 is adapted to live in oligotrophic water environments and biofilms. H. latericoloratus CGMCC 16346 is a promising bioremediation agent for elimination of imidacloprid contamination from surface water.

Evolutionary relationships among bifidobacteria and their hosts and environments

BACKGROUND

The assembly of animal microbiomes is influenced by multiple environmental factors and host genetics, although the relative importance of these factors remains unclear. Bifidobacteria (genus Bifidobacterium, phylum Actinobacteria) are common first colonizers of gut microbiomes in humans and inhabit other mammals, social insects, food, and sewages. In humans, the presence of bifidobacteria in the gut has been correlated with health-promoting benefits. Here, we compared the genome sequences of a subset of the over 400 Bifidobacterium strains publicly available to investigate the adaptation of bifidobacteria diversity. We tested 1) whether bifidobacteria show a phylogenetic signal with their isolation sources (hosts and environments) and 2) whether key traits encoded by the bifidobacteria genomes depend on the host or environment from which they were isolated. We analyzed Bifidobacterium genomes available in the PATRIC and NCBI repositories and identified the hosts and/or environment from which they were isolated. A multilocus phylogenetic analysis was conducted to compare the genetic relatedness the strains harbored by different hosts and environments. Furthermore, we examined differences in genomic traits and genes related to amino acid biosynthesis and degradation of carbohydrates.

RESULTS

We found that bifidobacteria diversity appears to have evolved with their hosts as strains isolated from the same host were non-randomly associated with their phylogenetic relatedness. Moreover, bifidobacteria isolated from different sources displayed differences in genomic traits such as genome size and accessory gene composition and on particular traits related to amino acid production and degradation of carbohydrates. In contrast, when analyzing diversity within human-derived bifidobacteria, we observed no phylogenetic signal or differences on specific traits (amino acid biosynthesis genes and CAZymes).

CONCLUSIONS

Overall, our study shows that bifidobacteria diversity is strongly adapted to specific hosts and environments and that several genomic traits were associated with their isolation sources. However, this signal is not observed in human-derived strains alone. Looking into the genomic signatures of bifidobacteria strains in different environments can give insights into how this bacterial group adapts to their environment and what types of traits are important for these adaptations.

Microdiversity ensures the maintenance of functional microbial communities under changing environmental conditions

Microdiversity can lead to different ecotypes within the same species. These are assumed to provide stability in time and space to those species. However, the role of microdiversity in the stability of whole microbial communities remains underexplored. Understanding the drivers of microbial community stability is necessary to predict community response to future disturbances. Here, we analyzed 16S rRNA gene amplicons from eight different temperate bog lakes at the 97% OTU and amplicon sequence variant (ASV) levels and found ecotypes within the same OTU with different distribution patterns in space and time. We observed that these ecotypes are adapted to different values of environmental factors such as water temperature and oxygen concentration. Our results showed that the existence of several ASVs within a OTU favored its persistence across changing environmental conditions. We propose that microdiversity aids the stability of microbial communities in the face of fluctuations in environmental factors.

Chromids Aid Genome Expansion and Functional Diversification in the Family Burkholderiaceae

Multipartite genomes, containing at least two large replicons, are found in diverse bacteria; however, the advantage of this genome structure remains incompletely understood. Here, we perform comparative genomics of hundreds of finished β-proteobacterial genomes to gain insights into the role and emergence of multipartite genomes. Almost all essential secondary replicons (chromids) of the β-proteobacteria are found in the family Burkholderiaceae. These replicons arose from just two plasmid acquisition events, and they were likely stabilized early in their evolution by the presence of core genes. On average, Burkholderiaceae genera with multipartite genomes had a larger total genome size, but smaller chromosome, than genera without secondary replicons. Pangenome-level functional enrichment analyses suggested that interreplicon functional biases are partially driven by the enrichment of secondary replicons in the accessory pangenome fraction. Nevertheless, the small overlap in orthologous groups present in each replicon's pangenome indicated a clear functional separation of the replicons. Chromids appeared biased to environmental adaptation, as the functional categories enriched on chromids were also overrepresented on the chromosomes of the environmental genera (Paraburkholderia and Cupriavidus) compared with the pathogenic genera (Burkholderia and Ralstonia). Using ancestral state reconstruction, it was predicted that the rate of accumulation of modern-day genes by chromids was more rapid than the rate of gene accumulation by the chromosomes. Overall, the data are consistent with a model where the primary advantage of secondary replicons is in facilitating increased rates of gene acquisition through horizontal gene transfer, consequently resulting in replicons enriched in genes associated with adaptation to novel environments.

Wastewater microbial community structure and functional traits change over short timescales

Wastewater contains microorganisms coming from various sources, e.g. feces discharges, soil infiltrations and sewer biofilms and sediments. The primary objective of this work was to determine if end-of-pipe wastewater microbial community structures exhibits short-timescale variation, and assess possible microbial origins. To this end, we measured hourly physicochemical characteristics of wastewater influent for 2 days and analyzed the microbial community at 4-h intervals using 16S rRNA gene amplicon sequencing. Results showed large variations in the microbial community composition at phylum and genus levels, i.e. Proteobacteria ranged from 44 to 63% of the total relative abundance and Arcobacter ranged from 11 to 22%. Diurnal patterns were observed in the alpha-diversity, beta-diversity and the prevalence of several taxa. Wastewater physicochemical characteristics explained 61% of the total microbial community variance by Canonical Correspondence Analysis (CCA), with flow rate being the main explanatory variable exhibiting a clear diurnal profile. Comparison with public databases using closed reference OTUs revealed that only 7.3% of the sequences were shared with human gut microbiota and 21.7% with soil microbiota, the majority being from the sewer biofilms and sediments. The functional trait, weighted average ribosomal RNA operon (rrn) copy number per genome, was found to be relatively high in the wastewater microbiota (average 3.6, soil 2.1, and human gut 2.6) and significantly correlated with flow, inferring active microbial enrichments in the sewer. The prevalence of Methylophilaceae, methanol oxidation genes and denitrification genes were related to high influent methanol and NO₃^- concentration in the influent wastewater. These functional organisms and genes indicate important carbon and nutrient removal related functions in the sewer. Together, the observed temporal patterns of the microbial community and functional traits suggest that high wastewater flow causes greater transport of active sewer microorganisms which are functionally important.

Co-occurrence patterns of litter decomposing communities in mangroves indicate a robust community resistant to disturbances

Background

Mangroves are important coastal ecosystems known for high photosynthetic productivity and the ability to support marine food chains through supply of dissolved carbon or particular organic matter. Most of the carbon found in mangroves is produced by its vegetation and is decomposed in root associated sediment. This process involves a tight interaction between microbial populations, litter chemical composition, and environmental parameters. Here, we study the complex interactions found during litter decomposition in mangroves by applying network analysis to metagenomic data.

Methods

Leaves of three species of mangrove trees typically found in the southeast of Brazil (Rhizophora mangle, Laguncularia racemosa, and Avicennia schaueriana) were collected in separate litter bags and left on three different mangroves for 60 days. These leaves were subsequently used for metagenome sequencing using Ion Torrent technology. Sequences were annotated in MG-RAST and used for network construction using MENAp.

Results

The most common phyla were Proteobacteria (classes Gamma and Alphaproteobacteria) followed by Firmicutes (Clostridia and Bacilli). The most abundant protein clusters were associated with the metabolism of carbohydrates, amino acids, and proteins. Non-metric multidimensional scaling of the metagenomic data indicated that substrate (i.e., tree species) did not significantly select for a specific community. Both networks exhibited scale-free characteristics and small world structure due to the low mean shortest path length and high average clustering coefficient. These networks also had a low number of hub nodes most of which were module hubs.

Discussion

This study demonstrates that under different environmental pressures (i.e., plant species or mangrove location) the microbial community associated with the decaying material forms a robust and stable network.

Functional Microbial Features Driving Community Assembly During Seed Germination and Emergence

Microbial interactions occurring on and around seeds are especially important for plant fitness since seed-borne microorganisms are the initial source of inoculum for the plant microbiota. In this study, we analyze structural and functional changes occurring within the plant microbiota at these early stages of the plant cycle, namely germination and emergence. To this purpose, we performed shotgun DNA sequencing of microbial assemblages associated to seeds, germinating seeds and seedlings of two plant species: bean and radish. We observed an enrichment of Enterobacteriales and Pseudomonadales during emergence and a set of functional traits linked to copiotrophy that could be responsible for this selection as a result of an increase of nutrient availability after germination. Representative bacterial isolates of taxa that are selected in seedlings showed indeed faster bacterial growth rate in comparison to seed-associated bacteria isolates. Finally, binning of metagenomics contigs results in the reconstruction of population genomes of the major bacterial taxa associated to the samples. Together, our results demonstrate that, although seed microbiota varied across plant species, nutrient availability during germination elicits changes of the composition of microbial communities by potentially selecting microbial groups with functional traits linked to copiotrophy. The data presented here represents the first attempts to empirically assess changes in the microbial community during plant emergence and moves us toward a more holistic understanding of the plant microbiome.

Variation in range size and dispersal capabilities of microbial taxa

Geographic range size can span orders of magnitude for plant and animal species, with the study of why range sizes vary having preoccupied biogeographers for decades. In contrast, there have been few comparable studies of how range size varies across microbial taxa and what traits may be associated with this variation. We determined the range sizes of 74,134 bacterial and archaeal taxa found in settled dust collected from 1,065 locations across the United States. We found that most microorganisms have small ranges and few have large ranges, a pattern similar to the range size distributions commonly observed for macrobes. However, contrary to expectations, those microbial taxa that were locally abundant did not necessarily have larger range sizes. The observed differences in microbial range sizes were generally predictable from taxonomic identity, phenotypic traits, genomic attributes, and habitat preferences, findings that provide insight into the factors shaping patterns of microbial biogeography.