Root-associated bacterial communities and root metabolite composition are linked to nitrogen use efficiency in sorghum.

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.

Microfabrication of a Chamber for High-Resolution, In Situ Imaging of the Whole Root for Plant-Microbe Interactions.

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.

Feature selection and causal analysis for microbiome studies in the presence of confounding using standardization.

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.

Genome-resolved metagenomics reveals role of iron metabolism in drought-induced rhizosphere microbiome dynamics.

Recent studies have demonstrated that drought leads to dramatic, highly conserved shifts in the root microbiome. At present, the molecular mechanisms underlying these responses remain largely uncharacterized. Here we employ genome-resolved metagenomics and comparative genomics to demonstrate that carbohydrate and secondary metabolite transport functionalities are overrepresented within drought-enriched taxa. These data also reveal that bacterial iron transport and metabolism functionality is highly correlated with drought enrichment. Using time-series root RNA-Seq data, we demonstrate that iron homeostasis within the root is impacted by drought stress, and that loss of a plant phytosiderophore iron transporter impacts microbial community composition, leading to significant increases in the drought-enriched lineage, Actinobacteria. Finally, we show that exogenous application of iron disrupts the drought-induced enrichment of Actinobacteria, as well as their improvement in host phenotype during drought stress. Collectively, our findings implicate iron metabolism in the root microbiome's response to drought and may inform efforts to improve plant drought tolerance to increase food security.

Microbial Community Field Surveys Reveal Abundant Pseudomonas Population in Sorghum Rhizosphere Composed of Many Closely Related Phylotypes.

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.

CRAGE-Duet Facilitates Modular Assembly of Biological Systems for Studying Plant-Microbe Interactions.

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.

PMR5, an acetylation protein at the intersection of pectin biosynthesis and defense against fungal pathogens.

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 [14 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.

Metabolomics of sorghum roots during nitrogen stress reveals compromised metabolic capacity for salicylic acid biosynthesis.

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 Wall Composition and Candidate Biosynthesis Gene Expression During Rice Development.

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.

Three Novel Rice Genes Closely Related to the Arabidopsis IRX9, IRX9L, and IRX14 Genes and Their Roles in Xylan Biosynthesis.

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.

Overexpression of a BAHD acyltransferase, OsAt10, alters rice cell wall hydroxycinnamic acid content and saccharification.

Grass cell wall properties influence food, feed, and biofuel feedstock usage efficiency. The glucuronoarabinoxylan of grass cell walls is esterified with the phenylpropanoid-derived hydroxycinnamic acids ferulic acid (FA) and para-coumaric acid (p-CA). Feruloyl esters undergo oxidative coupling with neighboring phenylpropanoids on glucuronoarabinoxylan and lignin. Examination of rice (Oryza sativa) mutants in a grass-expanded and -diverged clade of BAHD acyl-coenzyme A-utilizing transferases identified four mutants with altered cell wall FA or p-CA contents. Here, we report on the effects of overexpressing one of these genes, OsAt10 (LOC_Os06g39390), in rice. An activation-tagged line, OsAT10-D1, shows a 60% reduction in matrix polysaccharide-bound FA and an approximately 300% increase in p-CA in young leaf tissue but no discernible phenotypic alterations in vegetative development, lignin content, or lignin composition. Two additional independent OsAt10 overexpression lines show similar changes in FA and p-CA content. Cell wall fractionation and liquid chromatography-mass spectrometry experiments isolate the cell wall alterations in the mutant to ester conjugates of a five-carbon sugar with p-CA and FA. These results suggest that OsAT10 is a p-coumaroyl coenzyme A transferase involved in glucuronoarabinoxylan modification. Biomass from OsAT10-D1 exhibits a 20% to 40% increase in saccharification yield depending on the assay. Thus, OsAt10 is an attractive target for improving grass cell wall quality for fuel and animal feed.

XAX1 from glycosyltransferase family 61 mediates xylosyltransfer to rice xylan.

Xylan is the second most abundant polysaccharide on Earth and represents an immense quantity of stored energy for biofuel production. Despite its importance, most of the enzymes that synthesize xylan have yet to be identified. Xylans have a backbone of ß-1,4-linked xylose residues with substitutions that include α-(1→2)-linked glucuronosyl, 4-O-methyl glucuronosyl, and α-1,2- and α-1,3-arabinofuranosyl residues. The substitutions are structurally diverse and vary by taxonomy, with grass xylan representing a unique composition distinct from dicots and other monocots. To date, no enzyme has yet been identified that is specific to grass xylan synthesis. We identified a xylose-deficient loss-of-function rice mutant in Os02g22380, a putative glycosyltransferase in a grass-specific subfamily of family GT61. We designate the mutant xax1 for xylosyl arabinosyl substitution of xylan 1. Enzymatic fingerprinting of xylan showed the specific absence in the mutant of a peak, which was isolated and determined by (1)H-NMR to be (ß-1,4-Xyl)(4) with a ß-Xylp-(1→2)-α-Araf-(1→3). Rice xax1 mutant plants are deficient in ferulic and coumaric acid, aromatic compounds known to be attached to arabinosyl residues in xylan substituted with xylosyl residues. The xax1 mutant plants exhibit an increased extractability of xylan and increased saccharification, probably reflecting a lower degree of diferulic cross-links. Activity assays with microsomes isolated from tobacco plants transiently expressing XAX1 demonstrated xylosyltransferase activity onto endogenous acceptors. Our results provide insight into grass xylan synthesis and how substitutions may be modified for increased saccharification for biofuel generation.

ParA resolvase catalyzes site-specific excision of DNA from the Arabidopsis genome.

The small serine resolvase ParA from bacterial plasmids RK2 and RP4 catalyzes the recombination of two identical 133 bp recombination sites known as MRS. Previously, we reported that ParA is active in the fission yeast Schizosaccharomyces pombe. In this work, the parA recombinase gene was placed under the control of the Arabidopsis OXS3 promoter and introduced into Arabidopsis lines harboring a chromosomally integrated MRS-flanked target. The ParA recombinase excised the MRS-flanked DNA and the excision event was detected in subsequent generations in the absence of ParA, indicating germinal transmission of the excision event. The precise site-specific deletion by the ParA recombination system in planta demonstrates that the ParA recombinase can be used to remove transgenic DNA, such as selectable markers or other introduced transgenes that are no longer desired in the final product.