Plant and soil microbiome studies are becoming increasingly important for understanding the roles microorganisms play in agricultural productivity. The purpose of this manuscript is to provide detail on how to rapidly sample soil, rhizosphere, and endosphere of replicated field trials and analyze changes that may occur in the microbial communities due to sample type, treatment, and plant genotype. The experiment used to demonstrate these methods consists of replicated field plots containing two, pure, warm-season grasses (Panicum virgatum and Andropogon gerardii) and a low-diversity grass mixture (A. gerardii, Sorghastrum nutans, and Bouteloua curtipendula). Briefly, plants are excavated, a variety of roots are cut and placed in phosphate buffer, and then shaken to collect the rhizosphere. Roots are brought to the laboratory on ice and surface sterilized with bleach and ethanol (EtOH). The rhizosphere is filtered and concentrated by centrifugation. Excavated soil from around the root ball is placed into plastic bags and brought to the lab where a small amount of soil is taken for DNA extractions. DNA is extracted from roots, soil, and rhizosphere and then amplified with primers for the V4 region of the 16S rRNA gene. Amplicons are sequenced, then analyzed with open access bioinformatics tools. These methods allow researchers to test how the microbial community diversity and composition varies due to sample type, treatment, and plant genotype. Using these methods along with statistical models, the representative results demonstrate there are significant differences in microbial communities of roots, rhizosphere, and soil. Methods presented here provide a complete set of steps for how to collect field samples, isolate, extract, quantify, amplify, and sequence DNA, and analyze microbial community diversity and composition in replicated field trials.

Rising atmospheric concentrations of CO₂ and O₃ are key features of global environmental change. To investigate changes in the belowground bacterial community composition in response to elevated CO₂ and O₃ (eCO₂ and eO₃) the endosphere, rhizosphere and soil were sampled from soybeans under eCO₂ and maize under eO₃. The maize rhizosphere and endosphere α-diversity was higher than soybean, which may be due to a high relative abundance of Rhizobiales. Only the rhizosphere microbiome composition of the soybeans changed in response to eCO₂, associated with an increased abundance of nitrogen fixing microbes. In maize, the microbiome composition was altered by the genotype and linked to differences in root exudate profiles. The eO₃ treatment did not change the microbial communities in the rhizosphere, but altered the soil communities where hybrid maize was grown. In contrast to previous studies that focused exclusively on the soil, this study provides new insights into the effects of plant root exudates on the composition of the belowground microbiome in response to changing atmospheric conditions. Our results demonstrate that plant species and plant genotype were key factors driving the changes in the belowground bacterial community composition in agroecosystems that experience rising levels of atmospheric CO₂ and O₃.

Motivation

High-throughput sequencing technologies have greatly facilitated microbiome research and have generated a large volume of microbiome data with the potential to answer key questions regarding microbiome assembly, structure and function. Cluster analysis aims to group features that behave similarly across treatments, and such grouping helps to highlight the functional relationships among features and may provide biological insights into microbiome networks. However, clustering microbiome data are challenging due to the sparsity and high dimensionality.

Results

We propose a model-based clustering method based on Poisson hurdle models for sparse microbiome count data. We describe an expectation-maximization algorithm and a modified version using simulated annealing to conduct the cluster analysis. Moreover, we provide algorithms for initialization and choosing the number of clusters. Simulation results demonstrate that our proposed methods provide better clustering results than alternative methods under a variety of settings. We also apply the proposed method to a sorghum rhizosphere microbiome dataset that results in interesting biological findings.

Availability and implementation

R package is freely available for download at https://cran.r-project.org/package=PHclust.

Supplementary information

Supplementary data are available at Bioinformatics online.

Background

Microbiome studies have uncovered associations between microbes and human, animal, and plant health outcomes. This has led to an interest in developing microbial interventions for treatment of disease and optimization of crop yields which requires identification of microbiome features that impact the outcome in the population of interest. That task is challenging because of the high dimensionality of microbiome data and the confounding that results from the complex and dynamic interactions among host, environment, and microbiome. In the presence of such confounding, variable selection and estimation procedures may have unsatisfactory performance in identifying microbial features with an effect on the outcome.

Results

In this manuscript, we aim to estimate population-level effects of individual microbiome features while controlling for confounding by a categorical variable. Due to the high dimensionality and confounding-induced correlation between features, we propose feature screening, selection, and estimation conditional on each stratum of the confounder followed by a standardization approach to estimation of population-level effects of individual features. Comprehensive simulation studies demonstrate the advantages of our approach in recovering relevant features. Utilizing a potential-outcomes framework, we outline assumptions required to ascribe causal, rather than associational, interpretations to the identified microbiome effects. We conducted an agricultural study of the rhizosphere microbiome of sorghum in which nitrogen fertilizer application is a confounding variable. In this study, the proposed approach identified microbial taxa that are consistent with biological understanding of potential plant-microbe interactions.

Conclusions

Standardization enables more accurate identification of individual microbiome features with an effect on the outcome of interest compared to other variable selection and estimation procedures when there is confounding by a categorical variable.

The development of cereal crops with high nitrogen use efficiency (NUE) is a priority for worldwide agriculture. In addition to conventional plant breeding and genetic engineering, the use of the plant microbiome offers another approach to improving crop NUE. To gain insight into the bacterial communities associated with sorghum lines that differ in NUE, a field experiment was designed comparing 24 diverse Sorghum bicolor lines under sufficient and deficient nitrogen (N). Amplicon sequencing and untargeted gas chromatography-mass spectrometry were used to characterize the bacterial communities and the root metabolome associated with sorghum genotypes varying in sensitivity to low N. We demonstrated that N stress and sorghum type (energy, sweet, and grain sorghum) significantly impacted the root-associated bacterial communities and root metabolite composition of sorghum. We found a positive correlation between sorghum NUE and bacterial richness and diversity in the rhizosphere. The greater alpha diversity in high NUE lines was associated with the decreased abundance of a dominant bacterial taxon, Pseudomonas. Multiple strong correlations were detected between root metabolites and rhizosphere bacterial communities in response to low N stress. This indicates that the shift in the sorghum microbiome due to low N is associated with the root metabolites of the host plant. Taken together, our findings suggest that host genetic regulation of root metabolites plays a role in defining the root-associated microbiome of sorghum genotypes differing in NUE and tolerance to low N stress.IMPORTANCEThe development of crops that are more nitrogen use-efficient (NUE) is critical for the future of the enhanced sustainability of agriculture worldwide. This objective has been pursued mainly through plant breeding and plant molecular engineering, but these approaches have had only limited success. Therefore, a different strategy that leverages soil microbes needs to be fully explored because it is known that soil microbes improve plant growth through multiple mechanisms. To design approaches that use the soil microbiome to increase NUE, it will first be essential to understand the relationship among soil microbes, root metabolites, and crop productivity. Using this approach, we demonstrated that certain key metabolites and specific microbes are associated with high and low sorghum NUE in a field study. This important information provides a new path forward for developing crop genotypes that have increased NUE through the positive contribution of soil microbes.

While the root-associated microbiome is typically less diverse than the surrounding soil due to both plant selection and microbial competition for plant derived resources, it typically retains considerable complexity, harboring many hundreds of distinct bacterial species. Here, we report a time-dependent deviation from this trend in the rhizospheres of field grown sorghum. In this study, 16S rRNA amplicon sequencing was used to determine the impact of nitrogen fertilization on the development of the root-associated microbiomes of 10 sorghum genotypes grown in eastern Nebraska. We observed that early rhizosphere samples exhibit a significant reduction in overall diversity due to a high abundance of the bacterial genus Pseudomonas that occurred independent of host genotype in both high and low nitrogen fields and was not observed in the surrounding soil or associated root endosphere samples. When clustered at 97% identity, nearly all the Pseudomonas reads in this dataset were assigned to a single operational taxonomic unit (OTU); however, exact sequence variant (ESV)-level resolution demonstrated that this population comprised a large number of distinct Pseudomonas lineages. Furthermore, single-molecule long-read sequencing enabled high-resolution taxonomic profiling revealing further heterogeneity in the Pseudomonas lineages that was further confirmed using shotgun metagenomic sequencing. Finally, field soil enriched with specific carbon compounds recapitulated the increase in Pseudomonas, suggesting a possible connection between the enrichment of these Pseudomonas species and a plant-driven exudate profile.

Sorghum (Sorghum bicolor [L.] Moench) is the fifth most productive cereal crop worldwide with some hybrids having high biomass yield traits making it promising for sustainable, economical biofuel production. To maximize biofuel feedstock yields, a more complete understanding of metabolic responses to low nitrogen (N) will be useful for incorporation in crop improvement efforts. In this study, 10 diverse sorghum entries (including inbreds and hybrids) were field-grown under low and full N conditions and roots were sampled at two time points for metabolomics and 16S amplicon sequencing. Roots of plants grown under low N showed altered metabolic profiles at both sampling dates including metabolites important in N storage and synthesis of aromatic amino acids. Complementary investigation of the rhizosphere microbiome revealed dominance by a single operational taxonomic unit (OTU) in an early sampling that was taxonomically assigned to the genus Pseudomonas. Abundance of this Pseudomonas OTU was significantly greater under low N in July and was decreased dramatically in September. Correlation of Pseudomonas abundance with root metabolites revealed a strong negative association with the defense hormone salicylic acid (SA) under full N but not under low N, suggesting reduced defense response. Roots from plants with N stress also contained reduced phenylalanine, a precursor for SA, providing further evidence for compromised metabolic capacity for defense response under low N conditions. Our findings suggest that interactions between biotic and abiotic stresses may affect metabolic capacity for plant defense and need to be concurrently prioritized as breeding programs become established for biofuels production on marginal soils.

Correction to: The ISME Journal, published online 6 May 2022 In this article the author name Kira M. Veley was incorrectly written as Kira W. Veley. The original article has been corrected.

Drought is a major abiotic stress limiting agricultural productivity. Previous field-level experiments have demonstrated that drought decreases microbiome diversity in the root and rhizosphere. How these changes ultimately affect plant health remains elusive. Toward this end, we combined reductionist, transitional and ecological approaches, applied to the staple cereal crop sorghum to identify key root-associated microbes that robustly affect drought-stressed plant phenotypes. Fifty-three Arabidopsis-associated bacteria were applied to sorghum seeds and their effect on root growth was monitored. Two Arthrobacter strains caused root growth inhibition (RGI) in Arabidopsis and sorghum. In the context of synthetic communities, Variovorax strains were able to protect plants from Arthrobacter-caused RGI. As a transitional system, high-throughput phenotyping was used to test the synthetic communities. During drought stress, plants colonized by Arthrobacter had reduced growth and leaf water content. Plants colonized by both Arthrobacter and Variovorax performed as well or better than control plants. In parallel, we performed a field trial wherein sorghum was evaluated across drought conditions. By incorporating data on soil properties into the microbiome analysis, we accounted for experimental noise with a novel method and were able to observe the negative correlation between the abundance of Arthrobacter and plant growth. Having validated this approach, we cross-referenced datasets from the high-throughput phenotyping and field experiments and report a list of bacteria with high confidence that positively associated with plant growth under drought stress. In conclusion, a three-tiered experimental system successfully spanned the lab-to-field gap and identified beneficial and deleterious bacterial strains for sorghum under drought.