Fabricated ecosystems (EcoFABs) offer an innovative approach to in situ examination of microbial establishment patterns around plant roots using nondestructive, high-resolution microscopy. Previously high-resolution imaging was challenging because the roots were not constrained to a fixed distance from the objective. Here, we describe a new 'Imaging EcoFAB' and the use of this device to image the entire root system of growing Brachypodium distachyon at high resolutions (20×, 40×) over a 3-week period. The device is capable of investigating root-microbe interactions of multimember communities. We examined nine strains of Pseudomonas simiae with different fluorescent constructs to B. distachyon and individual cells on root hairs were visible. Succession in the rhizosphere using two different strains of P. simiae was examined, where the second addition was shown to be able to establish in the root tissue. The device was suitable for imaging with different solid media at high magnification, allowing for the imaging of fungal establishment in the rhizosphere. Overall, the Imaging EcoFAB could improve our ability to investigate the spatiotemporal dynamics of the rhizosphere, including studies of fluorescently-tagged, multimember, synthetic communities.

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.

Powdery mildew (Golovinomyces cichoracearum), one of the most prolific obligate biotrophic fungal pathogens worldwide, infects its host by penetrating the plant cell wall without activating the plant's innate immune system. The Arabidopsis mutant powdery mildew resistant 5 (pmr5) carries a mutation in a putative pectin acetyltransferase gene that confers enhanced resistance to powdery mildew. Here, we show that heterologously expressed PMR5 protein transfers acetyl groups from [¹⁴ C]-acetyl-CoA to oligogalacturonides. Through site-directed mutagenesis, we show that three amino acids within a highly conserved esterase domain in putative PMR5 orthologs are necessary for PMR5 function. A suppressor screen of mutagenized pmr5 seed selecting for increased powdery mildew susceptibility identified two previously characterized genes affecting the acetylation of plant cell wall polysaccharides, RWA2 and TBR. The rwa2 and tbr mutants also suppress powdery mildew disease resistance in pmr6, a mutant defective in a putative pectate lyase gene. Cell wall analysis of pmr5 and pmr6, and their rwa2 and tbr suppressor mutants, demonstrates minor shifts in cellulose and pectin composition. In direct contrast to their increased powdery mildew resistance, both pmr5 and pmr6 plants are highly susceptibile to multiple strains of the generalist necrotroph Botrytis cinerea, and have decreased camalexin production upon infection with B. cinerea. These results illustrate that cell wall composition is intimately connected to fungal disease resistance and outline a potential route for engineering powdery mildew resistance into susceptible crop species.

Xylan is the second most abundant polysaccharide on Earth, and represents a major component of both dicot wood and the cell walls of grasses. Much knowledge has been gained from studies of xylan biosynthesis in the model plant, Arabidopsis. In particular, the irregular xylem (irx) mutants, named for their collapsed xylem cells, have been essential in gaining a greater understanding of the genes involved in xylan biosynthesis. In contrast, xylan biosynthesis in grass cell walls is poorly understood. We identified three rice genes Os07g49370 (OsIRX9), Os01g48440 (OsIRX9L), and Os06g47340 (OsIRX14), from glycosyltransferase family 43 as putative orthologs to the putative β-1,4-xylan backbone elongating Arabidopsis IRX9, IRX9L, and IRX14 genes, respectively. We demonstrate that the over-expression of the closely related rice genes, in full or partly complement the two well-characterized Arabidopsis irregular xylem (irx) mutants: irx9 and irx14. Complementation was assessed by measuring dwarfed phenotypes, irregular xylem cells in stem cross sections, xylose content of stems, xylosyltransferase (XylT) activity of stems, and stem strength. The expression of OsIRX9 in the irx9 mutant resulted in XylT activity of stems that was over double that of wild type plants, and the stem strength of this line increased to 124% above that of wild type. Taken together, our results suggest that OsIRX9/OsIRX9L, and OsIRX14, have similar functions to the Arabidopsis IRX9 and IRX14 genes, respectively. Furthermore, our expression data indicate that OsIRX9 and OsIRX9L may function in building the xylan backbone in the secondary and primary cell walls, respectively. Our results provide insight into xylan biosynthesis in rice and how expression of a xylan synthesis gene may be modified to increase stem strength.

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.

Developing sustainable agricultural practices will require increasing our understanding of plant-microbe interactions. To study these interactions, new genetic tools for manipulating nonmodel microbes will be needed. To help meet this need, we recently reported development of chassis-independent recombinase-assisted genome engineering (CRAGE). CRAGE relies on cassette exchange between two pairs of mutually exclusive lox sites and allows direct, single-step chromosomal integration of large, complex gene constructs into diverse bacterial species. We then extended CRAGE by introducing a third mutually exclusive lox site, creating CRAGE-Duet, which allows modular integration of two constructs. CRAGE-Duet offers advantages over CRAGE, especially when a cumbersome recloning step is required to build single-integration constructs. To demonstrate the utility of CRAGE-Duet, we created a set of strains from the plant-growth-promoting rhizobacterium Pseudomonas simiae WCS417r that expressed various fluorescence marker genes. We visualized these strains simultaneously under a confocal microscope, demonstrating the usefulness of CRAGE-Duet for creating biological systems to study plant-microbe interactions.

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.

Cell walls of grasses, including cereal crops and biofuel grasses, comprise the majority of plant biomass and intimately influence plant growth, development and physiology. However, the functions of many cell wall synthesis genes, and the relationships among and the functions of cell wall components remain obscure. To better understand the patterns of cell wall accumulation and identify genes that act in grass cell wall biosynthesis, we characterized 30 samples from aerial organs of rice (Oryza sativa cv. Kitaake) at 10 developmental time points, 3-100 d post-germination. Within these samples, we measured 15 cell wall chemical components, enzymatic digestibility and 18 cell wall polysaccharide epitopes/ligands. We also used quantitative reverse transcription-PCR to measure expression of 50 glycosyltransferases, 15 acyltransferases and eight phenylpropanoid genes, many of which had previously been identified as being highly expressed in rice. Most cell wall components vary significantly during development, and correlations among them support current understanding of cell walls. We identified 92 significant correlations between cell wall components and gene expression and establish nine strong hypotheses for genes that synthesize xylans, mixed linkage glucan and pectin components. This work provides an extensive analysis of cell wall composition throughout rice development, identifies genes likely to synthesize grass cell walls, and provides a framework for development of genetically improved grasses for use in lignocellulosic biofuel production and agriculture.