Phylogeny-Aware Analysis of Metagenome Community Ecology Based on Matched Reference Genomes while Bypassing Taxonomy.

We introduce the operational genomic unit (OGU) method, a metagenome analysis strategy that directly exploits sequence alignment hits to individual reference genomes as the minimum unit for assessing the diversity of microbial communities and their relevance to environmental factors. This approach is independent of taxonomic classification, granting the possibility of maximal resolution of community composition, and organizes features into an accurate hierarchy using a phylogenomic tree. The outputs are suitable for contemporary analytical protocols for community ecology, differential abundance, and supervised learning while supporting phylogenetic methods, such as UniFrac and phylofactorization, that are seldom applied to shotgun metagenomics despite being prevalent in 16S rRNA gene amplicon studies. As demonstrated in two real-world case studies, the OGU method produces biologically meaningful patterns from microbiome data sets. Such patterns further remain detectable at very low metagenomic sequencing depths. Compared with taxonomic unit-based analyses implemented in currently adopted metagenomics tools, and the analysis of 16S rRNA gene amplicon sequence variants, this method shows superiority in informing biologically relevant insights, including stronger correlation with body environment and host sex on the Human Microbiome Project data set and more accurate prediction of human age by the gut microbiomes of Finnish individuals included in the FINRISK 2002 cohort. We provide Woltka, a bioinformatics tool to implement this method, with full integration with the QIIME 2 package and the Qiita web platform, to facilitate adoption of the OGU method in future metagenomics studies. IMPORTANCE Shotgun metagenomics is a powerful, yet computationally challenging, technique compared to 16S rRNA gene amplicon sequencing for decoding the composition and structure of microbial communities. Current analyses of metagenomic data are primarily based on taxonomic classification, which is limited in feature resolution. To solve these challenges, we introduce operational genomic units (OGUs), which are the individual reference genomes derived from sequence alignment results, without further assigning them taxonomy. The OGU method advances current read-based metagenomics in two dimensions: (i) providing maximal resolution of community composition and (ii) permitting use of phylogeny-aware tools. Our analysis of real-world data sets shows that it is advantageous over currently adopted metagenomic analysis methods and the finest-grained 16S rRNA analysis methods in predicting biological traits. We thus propose the adoption of OGUs as an effective practice in metagenomic studies.

Agricultural practices drive biological loads, seasonal patterns and potential pathogens in the aerobiome of a mixed-land-use dryland.

Air carries a diverse load of particulate microscopic biological matter in suspension, either aerosolized or aggregated with dust particles, the aerobiome, which is dispersed by winds from sources to sinks. The aerobiome is known to contain microbes, including pathogens, as well as debris or small-sized propagules from plants and animals, but its variability and composition has not been studied comprehensibly. To gain a dynamic insight into the aerobiome existing over a mixed-use dryland setting, we conducted a biologically comprehensive, year-long survey of its composition and dynamics for particles less than 10 µm in diameter based on quantitative analyses of DNA content coupled to genomic sequencing. Airborne biological loads were more dependent on seasonal events than on meteorological conditions and only weakly correlated with dust loads. Core aerobiome species could be understood as a mixture of high elevation (e.g. Microbacteriaceae, Micrococcaceae, Deinococci), and local plant and soil sources (e.g. Sphingomonas, Streptomyces, Acinetobacter). Despite the mixed used of the land surrounding the sampling site, taxa that contributed to high load events were largely traceable to proximal agricultural practices like cotton and livestock farming. This included not only the predominance of specific crop plant signals over those of native vegetation, but also that of their pathogens (bacterial, viral and eukaryotic). Faecal bacterial loads were also seasonally important, possibly sourced in intensive animal husbandry or manure fertilization activity, and this microbial load was enriched in tetracycline resistance genes. The presence of the native opportunistic pathogen, Coccidioides spp., by contrast, was detected only with highly sensitive techniques, and only rarely. We conclude that agricultural activity exerts a much stronger influence that the native vegetation as a mass loss factor to the land system and as an input to dryland aerobiomes, including in the dispersal of plant, animal and human pathogens and their genetic resistance characteristics.

DNA-Stable Isotope Probing Shotgun Metagenomics Reveals the Resilience of Active Microbial Communities to Biochar Amendment in Oxisol Soil.

The functions and interactions of individual microbial populations and their genes in agricultural soils amended with biochar remain elusive but are crucial for a deeper understanding of nutrient cycling and carbon (C) sequestration. In this study, we coupled DNA stable isotope probing (SIP) with shotgun metagenomics in order to target the active community in microcosms which contained soil collected from biochar-amended and control plots under napiergrass cultivation. Our analyses revealed that the active community was composed of high-abundant and low-abundant populations, including Actinobacteria, Proteobacteria, Gemmatimonadetes, and Acidobacteria. Although biochar did not significantly shift the active taxonomic and functional communities, we found that the narG (nitrate reductase) gene was significantly more abundant in the control metagenomes. Interestingly, putative denitrifier genomes generally encoded one gene or a partial denitrification pathway, suggesting denitrification is typically carried out by an assembly of different populations within this Oxisol soil. Altogether, these findings indicate that the impact of biochar on the active soil microbial community are transient in nature. As such, the addition of biochar to soils appears to be a promising strategy for the long-term C sequestration in agricultural soils, does not impart lasting effects on the microbial functional community, and thus mitigates un-intended microbial community shifts that may lead to fertilizer loss through increased N cycling.

Experimental Evolution of Campylobacter jejuni Leads to Loss of Motility, rpoN (σ54) Deletion and Genome Reduction.

Evolution experiments in the laboratory have focused heavily on model organisms, often to the exclusion of clinically relevant pathogens. The foodborne bacterial pathogen Campylobacter jejuni belongs to a genus whose genomes are small compared to those of its closest genomic relative, the free-living genus Sulfurospirillum, suggesting genome reduction during the course of evolution to host association. In an in vitro experiment, C. jejuni serially passaged in rich medium in the laboratory exhibited loss of flagellar motility-an essential function for host colonization. At early time points the motility defect was often reversible, but after 35 days of serial culture, motility was irreversibly lost in most cells in 5 independently evolved populations. Population re-sequencing revealed disruptive mutations to genes in the flagellar transcriptional cascade, rpoN (σ54)-therefore disrupting the expression of the genes σ54 regulates-coupled with deletion of rpoN in all evolved lines. Additional mutations were detected in virulence-related loci. In separate in vivo experiments, we demonstrate that a phase variable (reversible) motility mutant carrying an adenine deletion within a homopolymeric tract resulting in truncation of the flagellar biosynthesis gene fliR was deficient for colonization in a C57BL/6 IL-10-/- mouse disease model. Re-insertion of an adenine residue partially restored motility and ability to colonize mice. Thus, a pathogenic C. jejuni strain was rapidly attenuated by experimental laboratory evolution and demonstrated genomic instability during this evolutionary process. The changes observed suggest C. jejuni is able to evolve in a novel environment through genome reduction as well as transition, transversion, and slip-strand mutations.

Winter warming in Alaska accelerates lignin decomposition contributed by Proteobacteria.

BACKGROUND: In a warmer world, microbial decomposition of previously frozen organic carbon (C) is one of the most likely positive climate feedbacks of permafrost regions to the atmosphere. However, mechanistic understanding of microbial mediation on chemically recalcitrant C instability is limited; thus, it is crucial to identify and evaluate active decomposers of chemically recalcitrant C, which is essential for predicting C-cycle feedbacks and their relative strength of influence on climate change. Using stable isotope probing of the active layer of Arctic tundra soils after depleting soil labile C through a 975-day laboratory incubation, the identity of microbial decomposers of lignin and, their responses to warming were revealed. RESULTS: The ß-Proteobacteria genus Burkholderia accounted for 95.1% of total abundance of potential lignin decomposers. Consistently, Burkholderia isolated from our tundra soils could grow with lignin as the sole C source. A 2.2 °C increase of warming considerably increased total abundance and functional capacities of all potential lignin decomposers. In addition to Burkholderia, α-Proteobacteria capable of lignin decomposition (e.g. Bradyrhizobium and Methylobacterium genera) were stimulated by warming by 82-fold. Those community changes collectively doubled the priming effect, i.e., decomposition of existing C after fresh C input to soil. Consequently, warming aggravates soil C instability, as verified by microbially enabled climate-C modeling. CONCLUSIONS: Our findings are alarming, which demonstrate that accelerated C decomposition under warming conditions will make tundra soils a larger biospheric C source than anticipated. Video Abstract.

Diazotroph Diversity and Nitrogen Fixation in Summer Active Perennial Grasses in a Mediterranean Region Agricultural Soil.

Summer-growing perennial grasses such as Panicum coloratum L. cv. Bambatsi (Bambatsi panic), Chloris gayana Kunth cv. Katambora (Rhodes grass) and Digitaria eriantha Steud. cv. Premier (Premier digit grass) growing in the poor fertility sandy soils in the Mediterranean regions of southern Australia and western Australia mainly depend upon soil N and biological N inputs through diazotrophic (free living or associative) N fixation. We investigated the community composition and diversity (nifH-amplicon sequencing), abundance (qPCR) and functional capacity (15N incubation assay) of the endophytic diazotrophic community in the below and above ground plant parts of field grown and unfertilized grasses. Results showed a diverse and abundant diazotrophic community inside plant both above and below-ground and there was a distinct diazotrophic assemblage in the different plant parts in all the three grasses. There was a limited difference in the diversity between leaves, stems and roots except that Panicum grass roots harbored greater species richness. Nitrogen fixation potentials ranged between 0.24 and 5.9 mg N kg-1 day-1 and N fixation capacity was found in both the above and below ground plant parts. Results confirmed previous reports of plant species-based variation and that Alpha-Proteobacteria were the dominant group of nifH-harboring taxa both in the belowground and aboveground parts of the three grass species. Results also showed a well-structured nifH-harboring community in all plant parts, an example for a functional endophytic community. Overall, the variation in the number and identity of module hubs and connectors among the different plant parts suggests that co-occurrence patterns within the nifH-harboring community specific to individual compartments and local environments of the niches within each plant part may dictate the overall composition of diazotrophs within a plant.

Vineyard Soil Microbiome Composition Related to Rotundone Concentration in Australian Cool Climate 'Peppery' Shiraz Grapes.

Soil microbial communities have an integral association with plants and play an important role in shaping plant nutrition, health, crop productivity and product quality. The influence of bacteria and fungi on wine fermentation is well known. However, little is known about the role of soil microbes, other than microbial pathogens, on grape composition or their role in vintage or site (terroir) impacts on grape composition. In this study, we used an amplicon sequencing approach to investigate the potential relationships between soil microbes and inherent spatial variation in grape metabolite composition - specifically, the concentration of the 'impact aroma compound' rotundone in Shiraz grapes (Vitis vinifera L.) grown in a 6.1 ha vineyard in the Grampians region of Victoria, Australia. Previous work had demonstrated temporal stability in patterns of within-vineyard spatial variation in rotundone concentration, enabling identification of defined 'zones' of inherently 'low' or 'high' concentration of this grape metabolite. 16S rRNA and ITS region-amplicon sequencing analysis of microbial communities in the surface soils collected from these zones indicated marked differences between zones in the genetic diversity and composition of the soil bacterial and fungal microbiome. Soils in the High rotundone zone exhibited higher diversity of bacteria, but lower diversity of fungi, compared to the soils in the Low rotundone zone. In addition, the network analysis of the microbial community in the High rotundone zone soils appeared well structured, especially with respect to the bacterial community, compared to that in the Low rotundone zone soils. The key differences in the microbial community structure between the rotundone zones are obvious for taxa/groups of both bacteria and fungi, particularly for bacteria belonging to Acidobacteria-GP4 and GP7, Rhizobiales, Gaiellaceae, Alphaproteobacteria and the Nectriaceae and Tremellaceae families of fungi. Although mulching in some parts of the vineyard caused changes in bacterial and fungal composition and overall microbial catabolic diversity and activity, its effects did not mask the rotundone zone-based variation. This finding of a systematic rotundone zone-based variation in soil microbiomes suggests an opportunity to bring together understanding of microbial ecology, plant biochemistry, and viticultural management for improved management of grape metabolism, composition and wine flavor.

Comparative Metagenomics Reveals Enhanced Nutrient Cycling Potential after 2 Years of Biochar Amendment in a Tropical Oxisol.

The complex structural and functional responses of agricultural soil microbial communities to the addition of carbonaceous compounds such as biochar remain poorly understood. This severely limits the predictive ability for both the potential enhancement of soil fertility and greenhouse gas mitigation. In this study, we utilized shotgun metagenomics in order to decipher changes in the microbial community in soil microcosms after 14 days of incubation at 23°C, which contained soils from biochar-amended and control plots cultivated with Napier grass. Our analyses revealed that biochar-amended soil microbiomes exhibited significant shifts in both community composition and predicted metabolism. Key metabolic pathways related to carbon turnover, such as the utilization of plant-derived carbohydrates as well as denitrification, were enriched under biochar amendment. These community shifts were in part associated with increased soil carbon, such as labile and aromatic carbon compounds, which was likely stimulated by the increased available nutrients associated with biochar amendment. These findings indicate that the soil microbiome response to the combination of biochar addition and to incubation conditions confers enhanced nutrient cycling and a small decrease in CO2 emissions and potentially mitigates nitrous oxide emissions.IMPORTANCE The incorporation of biochar into soil is a promising management strategy for sustainable agriculture owing to its potential to sequester carbon and improve soil fertility. Expanding the addition of biochar to large-scale agriculture hinges on its lasting beneficial effects on the microbial community. However, there exists a significant knowledge gap regarding the specific role that biochar plays in altering the key biological soil processes that influence plant growth and carbon storage in soil. Previous studies that examined the soil microbiome under biochar amendment principally characterized only how the composition alters in response to biochar amendment. In the present study, we shed light on the functional alterations of the microbial community response 2 years after biochar amendment. Our results show that biochar increased the abundance of genes involved in denitrification and carbon turnover and that biochar-amended soil microcosms had a reduction in cumulative CO2 production.

Size Matters: Assessing Optimum Soil Sample Size for Fungal and Bacterial Community Structure Analyses Using High Throughput Sequencing of rRNA Gene Amplicons.

We examined the effect of different soil sample sizes obtained from an agricultural field, under a single cropping system uniform in soil properties and aboveground crop responses, on bacterial and fungal community structure and microbial diversity indices. DNA extracted from soil sample sizes of 0.25, 1, 5, and 10 g using MoBIO kits and from 10 and 100 g sizes using a bead-beating method (SARDI) were used as templates for high-throughput sequencing of 16S and 28S rRNA gene amplicons for bacteria and fungi, respectively, on the Illumina MiSeq and Roche 454 platforms. Sample size significantly affected overall bacterial and fungal community structure, replicate dispersion and the number of operational taxonomic units (OTUs) retrieved. Richness, evenness and diversity were also significantly affected. The largest diversity estimates were always associated with the 10 g MoBIO extractions with a corresponding reduction in replicate dispersion. For the fungal data, smaller MoBIO extractions identified more unclassified Eukaryota incertae sedis and unclassified glomeromycota while the SARDI method retrieved more abundant OTUs containing unclassified Pleosporales and the fungal genera Alternaria and Cercophora. Overall, these findings indicate that a 10 g soil DNA extraction is most suitable for both soil bacterial and fungal communities for retrieving optimal diversity while still capturing rarer taxa in concert with decreasing replicate variation.