Conservation of beneficial microbes between the rhizosphere and the cyanosphere

Biocrusts are phototroph-driven communities inhabiting arid soil surfaces. Like plants, most photoautotrophs (largely cyanobacteria) in biocrusts are thought to exchange fixed carbon for essential nutrients like nitrogen with cyanosphere bacteria. Here, we aim to compare beneficial interactions in rhizosphere and cyanosphere environments, including finding growth-promoting strains for hosts from both environments. To examine this, we performed a retrospective analysis of 16S rRNA gene sequencing datasets, host-microbe co-culture experiments between biocrust communities/biocrust isolates and a model grass (Brachypodium distachyon) or a dominant biocrust cyanobacterium (Microcoleus vaginatus), and metabolomic analysis. All 18 microbial phyla in the cyanosphere were also present in the rhizosphere, with additional 17 phyla uniquely found in the rhizosphere. The biocrust microbes promoted the growth of the model grass, and three biocrust isolates (Bosea sp._L1B56, Pseudarthrobacter sp._L1D14 and Pseudarthrobacter picheli_L1D33) significantly promoted the growth of both hosts. Moreover, pantothenic acid was produced by Pseudarthrobacter sp._L1D14 when grown on B. distachyon exudates, and supplementation of plant growth medium with this metabolite increased B. distachyon biomass by over 60%. These findings suggest that cyanobacteria and other diverse photoautotrophic hosts can be a source for new plant growth-promoting microbes and metabolites.

Grass xylan structural variation suggests functional specialization and distinctive interaction with cellulose and lignin

Xylan is the most abundant non-cellulosic polysaccharide in grass cell walls, and it has important structural roles. The name glucuronoarabinoxylan (GAX) is used to describe this variable hemicellulose. It has a linear backbone of β-1,4-xylose (Xyl) residues that may be substituted with α-1,2-linked (4-O-methyl)-glucuronic acid (GlcA), α-1,3-linked arabinofuranose (Araf), and sometimes acetylation at the O-2 and/or O-3 positions. The role of these substitutions remains unclear, although there is increasing evidence that they affect the way xylan interacts with other cell wall components, particularly cellulose and lignin. Here, we used substitution-dependent endo-xylanase enzymes to investigate the variability of xylan substitution in grass culm cell walls. We show that there are at least three different types of xylan: (i) an arabinoxylan with evenly distributed Araf substitutions without GlcA (AXe); (ii) a glucuronoarabinoxylan with clustered GlcA modifications (GAXc); and (iii) a highly substituted glucuronoarabinoxylan (hsGAX). Immunolocalization of AXe and GAXc in Brachypodium distachyon culms revealed that these xylan types are not restricted to a few cell types but are instead widely detected in Brachypodium cell walls. We hypothesize that there are functionally specialized xylan types within the grass cell wall. The even substitutions of AXe may permit folding and binding on the surface of cellulose fibrils, whereas the more complex substitutions of the other xylans may support a role in the matrix and interaction with other cell wall components.

Immune gene activation by NPR and TGA transcriptional regulators in the model monocot Brachypodium distachyon

The nonexpressor of pathogenesis-related (NPR) gene family is well known to play a crucial role in transactivation of TGA transcription factors for salicylic acid (SA)-responsive genes, including pathogenesis-related protein 1 (PR1), during plants' immune response after pathogen attack in the model dicot Arabidopsis thaliana. However, little is known about NPR gene functions in monocots. We therefore explored the functions of NPRs in SA signaling in the model monocot Brachypodium distachyon. BdNPR1 and BdNPR2/3 share structural similarities with A. thaliana AtNPR1/2 and AtNPR3/4 subfamilies, respectively. The transcript level of BdNPR2 but not BdNPR1/3 appeared to be positively regulated in leaves in response to methyl salicylate. Reporter assays in protoplasts showed that BdNPR2 positively regulated BdTGA1-mediated activation of PR1. This transactivation occurred in an SA-dependent manner through SA binding at Arg468 of BdNPR2. In contrast, BdNPR1 functioned as a suppressor of BdNPR2/BdTGA1-mediated transcription of PR1. Collectively, our findings reveal that the TGA-promoted transcription of SA-inducible PR1 is orchestrated by the activator BdNPR2 and the repressor BdNPR1, which function competitively in B. distachyon.

Disruption of Brachypodium lichenase alters metabolism of mixed-linkage glucan and starch

Mixed-linkage glucan, which is widely distributed in grasses, is a polysaccharide highly abundant in cell walls of grass endosperm and young vegetative tissues. Lichenases are enzymes that hydrolyze mixed-linkage glucan first identified in mixed-linkage glucan-rich lichens. In this study, we identify a gene encoding a lichenase we name Brachypodium distachyon LICHENASE 1 (BdLCH1), which is highly expressed in the endosperm of germinating seeds and coleoptiles and at lower amounts in mature shoots. RNA in situ hybridization showed that BdLCH1 is primarily expressed in chlorenchyma cells of mature leaves and internodes. Disruption of BdLCH1 resulted in an eight-fold increase in mixed-linkage glucan content in senesced leaves. Consistent with the in situ hybridization data, immunolocalization results showed that mixed-linkage glucan was not removed in chlorenchyma cells of lch1 mutants as it was in wild type and implicate the BdLCH1 enzyme in removing mixed-linkage glucan in chlorenchyma cells in mature vegetative tissues. We also show that mixed-linkage glucan accumulation in lch1 mutants was resistant to dark-induced degradation, and 8-week-old lch1 plants showed a faster rate of starch breakdown than wild type in darkness. Our results suggest a role for BdLCH1 in modifying the cell wall to support highly metabolically active cells.

A phosphorus-limitation induced, functionally conserved DUF506 protein is a repressor of root hair elongation in plants

Root hairs (RHs) function in nutrient and water acquisition, root metabolite exudation, soil anchorage and plant-microbe interactions. Longer or more abundant RHs are potential breeding traits for developing crops that are more resource-use efficient and can improve soil health. While many genes are known to promote RH elongation, relatively little is known about genes and mechanisms that constrain RH growth. Here we demonstrate that a DOMAIN OF UNKNOWN FUNCTION 506 (DUF506) protein, AT3G25240, negatively regulates Arabidopsis thaliana RH growth. The AT3G25240 gene is strongly and specifically induced during phosphorus (P)-limitation. Mutants of this gene, which we call REPRESSOR OF EXCESSIVE ROOT HAIR ELONGATION 1 (RXR1), have much longer RHs, higher phosphate content and seedling biomass, while overexpression of the gene exhibits opposite phenotypes. Co-immunoprecipitation, pull-down and bimolecular fluorescence complementation (BiFC) analyses reveal that RXR1 physically interacts with a RabD2c GTPase in nucleus, and a rabd2c mutant phenocopies the rxr1 mutant. Furthermore, N-terminal variable region of RXR1 is crucial for inhibiting RH growth. Overexpression of a Brachypodium distachyon RXR1 homolog results in repression of RH elongation in Brachypodium. Taken together, our results reveal a novel DUF506-GTPase module with a prominent role in repression of plant RH elongation especially under P stress.

An aptamer highly specific to cellulose enables the analysis of the association of cellulose with matrix cell wall polymers in vitro and in muro

The current toolbox of cell wall-directed molecular probes has been pivotal for advancing basic and application-oriented plant carbohydrate research; however, it still exhibits limitations regarding target diversity and specificity. Scarcity of probes targeting intramolecular associations between cell wall polymers particularly hinders our understanding of the cell wall microstructure and affects the development of effective means for its efficient deconstruction for bioconversion. Here we report a detailed characterization of a cellulose-binding DNA aptamer CELAPT MINI using a combination of various in vitro biochemical, biophysical, and molecular biology techniques. Our results show evidence for its high specificity towards long non-substituted β-(1-4)-glucan chains in both crystalline and amorphous forms. Fluorescent conjugates of CELAPT MINI are applicable as in situ cellulose probes and are well suited for various microscopy techniques, including super-resolution imaging. Compatibility of fluorescent CELAPT MINI variants with immunodetection of cell wall matrix polymers enabled them simultaneously to resolve the fibrillar organization of complex cellulose-enriched pulp material and to quantify the level of cellulose masking by xyloglucan and xylan. Using enzymatically, chemically, or genetically modulated Brachypodium internode sections we showed the diversity in cell wall packing among various cell types and even cell wall microdomains. We showed that xylan is the most prominent, but not the only, cellulose-masking agent in Brachypodium internode tissues. These results collectively highlight the hitherto unexplored potential to expand the cell wall probing toolbox with highly specific and versatile in vitro generated polynucleotide probes.

BdFAR4, a root-specific fatty acyl-coenzyme A reductase, is involved in fatty alcohol synthesis of root suberin polyester in Brachypodium distachyon

Suberin is a complex hydrophobic polymer of aliphatic and phenolic compounds which controls the movement of gases, water, and solutes and protects plants from environmental stresses and pathogenic infection. The synthesis and regulatory pathways of suberin remain unknown in Brachypodium distachyon. Here we describe the identification of a B. distachyon gene, BdFAR4, encoding a fatty acyl-coenzyme A reductase (FAR) by a reverse genetic approach, and investigate the molecular relevance of BdFAR4 in the root suberin synthesis of B. distachyon. BdFAR4 is specifically expressed throughout root development. Heterologous expression of BdFAR4 in yeast (Saccharomyces cerevisiae) afforded the production of C20:0 and C22:0 fatty alcohols. The loss-of-function knockout of BdFAR4 by CRISPR/Cas9-mediated gene editing significantly reduced the content of C20:0 and C22:0 fatty alcohols associated with root suberin. In contrast, overexpression of BdFAR4 in B. distachyon and tomato (Solanum lycopersicum) resulted in the accumulation of root suberin-associated C20:0 and C22:0 fatty alcohols, suggesting that BdFAR4 preferentially accepts C20:0 and C22:0 fatty acyl-CoAs as substrates. The BdFAR4 protein was localized to the endoplasmic reticulum in Arabidopsis thaliana protoplasts and Nicotiana benthamiana leaf epidermal cells. BdFAR4 transcript levels can be increased by abiotic stresses and abscisic acid treatment. Furthermore, yeast one-hybrid, dual-luciferase activity, and electrophoretic mobility shift assays indicated that the R2R3-MYB transcription factor BdMYB41 directly binds to the promoter of BdFAR4. Taken together, these results imply that BdFAR4 is essential for the production of root suberin-associated fatty alcohols, especially under stress conditions, and that its activity is transcriptionally regulated by the BdMYB41 transcription factor.

Transcriptional memories mediate the plasticity of cold stress responses to enable morphological acclimation in Brachypodium distachyon

Plants that successfully acclimate to stress can resume growth under stressful conditions. The grass Brachypodium distachyon can grow a cold-adaptive morphology during cold acclimation. Studies on transcriptional memory (TM) have revealed that plants can be primed for stress by adjusting their transcriptional responses, but the function of TM in stress acclimation is not well understood. We investigated the function of TM during cold acclimation in B. distachyon. Quantitative polymerase chain reaction (qPCR), RNA-seq and chromatin immunoprecipitation qPCR analyses were performed on plants exposed to repeated episodes of cold to characterize the presence and stability of TM during the stress and growth responses of cold acclimation. Transcriptional memory mainly dampened stress responses as growth resumed and as B. distachyon became habituated to cold stress. Although permanent on vernalization gene VRN1, TMs were short-term and reversible on cold-stress genes. Growing under cold conditions also coincided with the acquisition of new and targeted cold-induced transcriptional responses. Overall, TM provided plasticity to cold stress responses during cold acclimation in B. distachyon, leading to stress habituation, acquired stress responses, and resumed growth. Our study shows that chromatin-associated TMs are involved in tuning plant responses to environmental change and, as such, regulate both stress and developmental components that characterize cold-climate adaptation in B. distachyon.

BdWRKY38 is required for the incompatible interaction of Brachypodium distachyon with the necrotrophic fungus Rhizoctonia solani

Rhizoctonia solani is a soil-borne necrotrophic fungus that causes sheath blight in grasses. The basal resistance of compatible interactions between R. solani and rice is known to be modulated by some WRKY transcription factors (TFs). However, genes and defense responses involved in incompatible interaction with R. solani remain unexplored, because no such interactions are known in any host plants. Recently, we demonstrated that Bd3-1, an accession of the model grass Brachypodium distachyon, is resistant to R. solani and, upon inoculation with the fungus, undergoes rapid induction of genes responsive to the phytohormone salicylic acid (SA) that encode the WRKY TFs BdWRKY38 and BdWRKY44. Here, we show that endogenous SA and these WRKY TFs positively regulate this accession-specific R. solani resistance. In contrast to a susceptible accession (Bd21), the infection process in the resistant accessions Bd3-1 and Tek-3 was suppressed at early stages before the development of fungal biomass and infection machinery. A comparative transcriptome analysis during pathogen infection revealed that putative WRKY-dependent defense genes were induced faster in the resistant accessions than in Bd21. A gene regulatory network (GRN) analysis based on the transcriptome dataset demonstrated that BdWRKY38 was a GRN hub connected to many target genes specifically in resistant accessions, whereas BdWRKY44 was shared in the GRNs of all three accessions. Moreover, overexpression of BdWRKY38 increased R. solani resistance in Bd21. Our findings demonstrate that these resistant accessions can activate an incompatible host response to R. solani, and BdWRKY38 regulates this response by mediating SA signaling.

Mutations in the predicted DNA polymerase subunit POLD3 result in more rapid flowering of Brachypodium distachyon

The timing of reproduction is a critical developmental decision in the life cycle of many plant species. Fine mapping of a rapid-flowering mutant was done using whole-genome sequence data from bulked DNA from a segregating F2 mapping populations. The causative mutation maps to a gene orthologous with the third subunit of DNA polymerase δ (POLD3), a previously uncharacterized gene in plants. Expression analyses of POLD3 were conducted via real time qPCR to determine when and in what tissues the gene is expressed. To better understand the molecular basis of the rapid-flowering phenotype, transcriptomic analyses were conducted in the mutant vs wild-type. Consistent with the rapid-flowering mutant phenotype, a range of genes involved in floral induction and flower development are upregulated in the mutant. Our results provide the first characterization of the developmental and gene expression phenotypes that result from a lesion in POLD3 in plants.

Changes in ambient temperature are the prevailing cue in determining Brachypodium distachyon diurnal gene regulation

Plants are continuously exposed to diurnal fluctuations in light and temperature, and spontaneous changes in their physical or biotic environment. The circadian clock coordinates regulation of gene expression with a 24 h period, enabling the anticipation of these events. We used RNA sequencing to characterize the Brachypodium distachyon transcriptome under light and temperature cycles, as well as under constant conditions. Approximately 3% of the transcriptome was regulated by the circadian clock, a smaller proportion than reported in most other species. For most transcripts that were rhythmic under all conditions, including many known clock genes, the period of gene expression lengthened from 24 to 27 h in the absence of external cues. To functionally characterize the cyclic transcriptome in B. distachyon, we used Gene Ontology enrichment analysis, and found several terms significantly associated with peak expression at particular times of the day. Furthermore, we identified sequence motifs enriched in the promoters of similarly phased genes, some potentially associated with transcription factors. When considering the overlap in rhythmic gene expression and specific pathway behavior, thermocycles was the prevailing cue that controlled diurnal gene regulation. Taken together, our characterization of the rhythmic B. distachyon transcriptome represents a foundational resource with implications in other grass species.

Grass secondary cell walls, Brachypodium distachyon as a model for discovery

A key aspect of plant growth is the synthesis and deposition of cell walls. In specific tissues and cell types including xylem and fibre, a thick secondary wall comprised of cellulose, hemicellulose and lignin is deposited. Secondary cell walls provide a physical barrier that protects plants from pathogens, promotes tolerance to abiotic stresses and fortifies cells to withstand the forces associated with water transport and the physical weight of plant structures. Grasses have numerous cell wall features that are distinct from eudicots and other plants. Study of the model species Brachypodium distachyon as well as other grasses has revealed numerous features of the grass cell wall. These include the characterisation of xylosyl and arabinosyltransferases, a mixed-linkage glucan synthase and hydroxycinnamate acyltransferases. Perhaps the most fertile area for discovery has been the formation of lignins, including the identification of novel substrates and enzyme activities towards the synthesis of monolignols. Other enzymes function as polymerising agents or transferases that modify lignins and facilitate interactions with polysaccharides. The regulatory aspects of cell wall biosynthesis are largely overlapping with those of eudicots, but salient differences among species have been resolved that begin to identify the determinants that define grass cell walls.

Brachypodium distachyon MAP20 functions in metaxylem pit development and contributes to drought recovery

Pits are regions in the cell walls of plant tracheary elements that lack secondary walls. Each pit consists of a space within the secondary wall called a pit chamber, and a modified primary wall called the pit membrane. The pit membrane facilitates transport of solutions between vessel cells and restricts embolisms during drought. Here we analyzed the role of an angiosperm-specific TPX2-like microtubule protein MAP20 in pit formation using Brachypodium distachyon as a model system. Live cell imaging was used to analyze the interaction of MAP20 with microtubules and the impact of MAP20 on microtubule dynamics. MAP20-specific antibody was used to study expression and localization of MAP20 in different cell types during vascular bundle development. We used an artificial microRNAs (amiRNA) knockdown approach to determine the function of MAP20. MAP20 is expressed during the late stages of vascular bundle development and localizes around forming pits and under secondary cell wall thickenings in metaxylem cells. MAP20 suppresses microtubule depolymerization; however, unlike the animal TPX2 counterpart, MAP20 does not cooperate with the γ-tubulin ring complex in microtubule nucleation. Knockdown of MAP20 causes bigger pits, thinner pit membranes, perturbed vasculature development, lower reproductive potential and higher drought susceptibility. We conclude that MAP20 may contribute to drought adaptation by modulating pit size and pit membrane thickness in metaxylem.

Roles of BdUNICULME4 and BdLAXATUM-A in the non-domesticated grass Brachypodium distachyon

In cultivated grasses, tillering, spike architecture and seed shattering represent major agronomical traits. In barley, maize and rice, the NOOT-BOP-COCH-LIKE (NBCL) genes play important roles in development, especially in ligule development, tillering and flower identity. However, compared with dicots, the role of grass NBCL genes is underinvestigated. To better understand the role of grass NBCLs and to overcome any effects of domestication that might conceal their original functions, we studied TILLING nbcl mutants in the non-domesticated grass Brachypodium distachyon. In B. distachyon, the NBCL genes BdUNICULME4 (CUL4) and BdLAXATUM-A (LAXA) are orthologous, respectively, to the barley HvUniculme4 and HvLaxatum-a, to the maize Zmtassels replace upper ears1 and Zmtassels replace upper ears2 and to the rice OsBLADE-ON-PETIOLE1 and OsBLADE-ON-PETIOLE2/3. In B. distachyon, our reverse genetics study shows that CUL4 is not essential for the establishment of the blade-sheath boundary but is necessary for the development of the ligule and auricles. We report that CUL4 also exerts a positive role in tillering and a negative role in spikelet meristem activity. On the other hand, we demonstrate that LAXA plays a negative role in tillering, positively participates in spikelet development and contributes to the control of floral organ number and identity. In this work, we functionally characterized two new NBCL genes in a context of non-domesticated grass and highlighted original roles for grass NBCL genes that are related to important agronomical traits.

Multilab EcoFAB study shows highly reproducible physiology and depletion of soil metabolites by a model grass

There is a dynamic reciprocity between plants and their environment: soil physiochemical properties influence plant morphology and metabolism, and root morphology and exudates shape the environment surrounding roots. Here, we investigate the reproducibility of plant trait changes in response to three growth environments. We utilized fabricated ecosystem (EcoFAB) devices to grow the model grass Brachypodium distachyon in three distinct media across four laboratories: phosphate-sufficient and -deficient mineral media allowed assessment of the effects of phosphate starvation, and a complex, sterile soil extract represented a more natural environment with yet uncharacterized effects on plant growth and metabolism. Tissue weight and phosphate content, total root length, and root tissue and exudate metabolic profiles were consistent across laboratories and distinct between experimental treatments. Plants grown in soil extract were morphologically and metabolically distinct, with root hairs four times longer than with other growth conditions. Further, plants depleted half of the metabolites investigated from the soil extract. To interact with their environment, plants not only adapt morphology and release complex metabolite mixtures, but also selectively deplete a range of soil-derived metabolites. The EcoFABs utilized here generated high interlaboratory reproducibility, demonstrating their value in standardized investigations of plant traits.

Multilab EcoFAB study shows highly reproducible physiology and depletion of soil metabolites by a model grass

There is a dynamic reciprocity between plants and their environment: soil physiochemical properties influence plant morphology and metabolism, and root morphology and exudates shape the environment surrounding roots. Here, we investigate the reproducibility of plant trait changes in response to three growth environments. We utilized fabricated ecosystem (EcoFAB) devices to grow the model grass Brachypodium distachyon in three distinct media across four laboratories: phosphate-sufficient and -deficient mineral media allowed assessment of the effects of phosphate starvation, and a complex, sterile soil extract represented a more natural environment with yet uncharacterized effects on plant growth and metabolism. Tissue weight and phosphate content, total root length, and root tissue and exudate metabolic profiles were consistent across laboratories and distinct between experimental treatments. Plants grown in soil extract were morphologically and metabolically distinct, with root hairs four times longer than with other growth conditions. Further, plants depleted half of the metabolites investigated from the soil extract. To interact with their environment, plants not only adapt morphology and release complex metabolite mixtures, but also selectively deplete a range of soil-derived metabolites. The EcoFABs utilized here generated high interlaboratory reproducibility, demonstrating their value in standardized investigations of plant traits.

SECONDARY WALL ASSOCIATED MYB1 is a positive regulator of secondary cell wall thickening in Brachypodium distachyon and is not found in the Brassicaceae

Grass biomass is comprised chiefly of secondary walls that surround fiber and xylem cells. A regulatory network of interacting transcription factors in part regulates cell wall thickening. We identified Brachypodium distachyon SECONDARY WALL ASSOCIATED MYB1 (SWAM1) as a potential regulator of secondary cell wall biosynthesis based on gene expression, phylogeny, and transgenic plant phenotypes. SWAM1 interacts with cellulose and lignin gene promoters with preferential binding to AC-rich sequence motifs commonly found in the promoters of cell wall-related genes. SWAM1 overexpression (SWAM-OE) lines had greater above-ground biomass with only a slight change in flowering time while SWAM1 dominant repressor (SWAM1-DR) plants were severely dwarfed with a striking reduction in lignin of sclerenchyma fibers and stem epidermal cell length. Cellulose, hemicellulose, and lignin genes were significantly down-regulated in SWAM1-DR plants and up-regulated in SWAM1-OE plants. There was no reduction in bioconversion yield in SWAM1-OE lines; however, it was significantly increased for SWAM1-DR samples. Phylogenetic and syntenic analyses strongly suggest that the SWAM1 clade was present in the last common ancestor between eudicots and grasses, but is not in the Brassicaceae. Collectively, these data suggest that SWAM1 is a transcriptional activator of secondary cell wall thickening and biomass accumulation in B. distachyon.

Genome-wide scans of selection highlight the impact of biotic and abiotic constraints in natural populations of the model grass Brachypodium distachyon

Grasses are essential plants for ecosystem functioning. Quantifying the selective pressures that act on natural variation in grass species is therefore essential regarding biodiversity maintenance. In this study, we investigate the selection pressures that act on two distinct populations of the grass model Brachypodium distachyon without prior knowledge about the traits under selection. We took advantage of whole-genome sequencing data produced for 44 natural accessions of B. distachyon and used complementary genome-wide selection scans (GWSS) methods to detect genomic regions under balancing and positive selection. We show that selection is shaping genetic diversity at multiple temporal and spatial scales in this species, and affects different genomic regions across the two populations. Gene ontology annotation of candidate genes reveals that pathogens may constitute important factors of positive and balancing selection in B. distachyon. We eventually cross-validated our results with quantitative trait locus data available for leaf-rust resistance in this species and demonstrate that, when paired with classical trait mapping, GWSS can help pinpointing candidate genes for further molecular validation. Thanks to a near base-perfect reference genome and the large collection of freely available natural accessions collected across its natural range, B. distachyon appears as a prime system for studies in ecology, population genomics and evolutionary biology.

Global analysis of H3K4me3/H3K27me3 in Brachypodium distachyon reveals VRN3 as critical epigenetic regulation point in vernalization and provides insights into epigenetic memory

Vernalization, the requirement of plants for long-term exposure to low environmental temperature for flowering, is an epigenetic phenomenon. Histone modification regulation has been revealed in vernalization, but is limited to key genes. Now, we know that VRN1 is epigenetically critical for monocots. Genome-wide analysis is still unavailable, however. We performed chromatin immunoprecipitation-sequencing for H3K4me3/H3K27me3 in Brachypodium distachyon to obtain a global view of histone modifications in vernalization on a genome-wide scale and for different pathways/genes. Our data showed that H3K4me3 and H3K27me3 play distinct roles in vernalization. Unlike H3K4me3, H3K27me3 exhibited regional regulation, showed main regulation targets in vernalization and contributed to epigenetic memory. For genes in four flowering regulation pathways, only FT2 (functional ortholog of VRN3 in B. distachyon) and VRN1 showed coordinated changes in H3K4me3/H3K27me3. The epigenetic response at VRN3 was weaker under short-day than under long-day conditions. VRN3 was revealed as an epigenetic regulation point integrating vernalization and day length signals. We globally identified genes maintaining vernalization-induced epigenetic changes. Most of these genes showed dose-dependent vernalization responses, revealing a quantitative 'recording system' for vernalization. Our studies shed light on the epigenetic role of VRN3 and H3K4me3/H3K27me3 in vernalization and reveal genes underlying epigenetic memory, laying the foundation for further study.

In the grass species Brachypodium distachyon, the production of mixed-linkage (1,3;1,4)-β-glucan (MLG) occurs in the Golgi apparatus

Mixed-linkage (1,3;1,4)-β-glucan (MLG) is a glucose polymer with beneficial effects on human health and high potential for the agricultural industry. MLG is present predominantly in the cell wall of grasses and is synthesized by cellulose synthase-like F or H families of proteins, with CSLF6 being the best-characterized MLG synthase. Although the function of this enzyme in MLG production has been established, the site of MLG synthesis in the cell is debated. It has been proposed that MLG is synthesized at the plasma membrane, as occurs for cellulose and callose; in contrast, it has also been proposed that MLG is synthesized in the Golgi apparatus, as occurs for other matrix polysaccharides of the cell wall. Testing these conflicting possibilities is fundamentally important in the general understanding of the biosynthesis of the plant cell wall. Using immuno-localization analyses with MLG-specific antibody in Brachypodium and in barley, we found MLG present in the Golgi, in post-Golgi structures and in the cell wall. Accordingly, analyses of a functional fluorescent protein fusion of CSLF6 stably expressed in Brachypodium demonstrated that the enzyme is localized in the Golgi. We also established that overproduction of MLG causes developmental and growth defects in Brachypodium as also occur in barley. Our results indicated that MLG production occurs in the Golgi similarly to other cell wall matrix polysaccharides, and supports the broadly applicable model in grasses that tight mechanisms control optimal MLG accumulation in the cell wall during development and growth. This work addresses the fundamental question of where mixed linkage (1,3;1,4)-β-glucan (MLG) is synthesized in plant cells. By analyzing the subcellular localization of MLG and MLG synthase in an endogenous system, we demonstrated that MLG synthesis occurs at the Golgi in Brachypodium and barley. A growth inhibition due to overproduced MLG in Brachypodium supports the general applicability of the model that a tight control of the cell wall polysaccharides accumulation is needed to maintain growth homeostasis during development.

Salicylic acid-dependent immunity contributes to resistance against Rhizoctonia solani, a necrotrophic fungal agent of sheath blight, in rice and Brachypodium distachyon

Rhizoctonia solani is a soil-borne fungus causing sheath blight. In consistent with its necrotrophic life style, no rice cultivars fully resistant to R. solani are known, and agrochemical plant defense activators used for rice blast, which upregulate a phytohormonal salicylic acid (SA)-dependent pathway, are ineffective towards this pathogen. As a result of the unavailability of genetics, the infection process of R. solani remains unclear. We used the model monocotyledonous plants Brachypodium distachyon and rice, and evaluated the effects of phytohormone-induced resistance to R. solani by pharmacological, genetic and microscopic approaches to understand fungal pathogenicity. Pretreatment with SA, but not with plant defense activators used in agriculture, can unexpectedly induce sheath blight resistance in plants. SA treatment inhibits the advancement of R. solani to the point in the infection process in which fungal biomass shows remarkable expansion and specific infection machinery is developed. The involvement of SA in R. solani resistance is demonstrated by SA-deficient NahG transgenic rice and the sheath blight-resistant B. distachyon accessions, Bd3-1 and Gaz-4, which activate SA-dependent signaling on inoculation. Our findings suggest a hemi-biotrophic nature of R. solani, which can be targeted by SA-dependent plant immunity. Furthermore, B. distachyon provides a genetic resource that can confer disease resistance against R. solani to plants.

Non-redundant functions of the dimeric ABA receptor BdPYL1 in the grass Brachypodium

Abiotic stresses have severe detrimental effects on agricultural productivity worldwide. Abscisic acid (ABA) levels rise in response to abiotic stresses, and play a role in coordinating physiological responses. ABA elicits its effects by binding a family of soluble receptors, increasing affinity of the receptors to type 2C phosphatases (PP2Cs) leading to phosphatase inhibition. In the current study, we conducted a comprehensive analysis of the ABA signaling pathway in the cereal model grass Brachypodium distachyon. The Brachypodium genome encodes a family of 10 functionally conserved ABA receptors. The 10th in the series, BdPYL10, encodes a defective receptor and is likely a pseudogene. Combinatorial protein interaction assay further validated computational analysis, which grouped Brachypodium ABA receptors into three subfamilies, similarly to Arabidopsis classification. Brachypodium subfamily III receptors inhibited PP2C activity in vitro and complemented Arabidopsis quadruple (pyr1/pyl1/pyl2/pyl4) mutant. BdPYL1 T-DNA mutant exhibited clear ABA hyposensitivity phenotypes during seedling establishment and in mature plants. Single receptor predominance is in agreement with high transcriptional abundance of only a small Brachypodium ABA receptors subset, harboring the higher marginal significance of BdPYL1. Our findings suggest that unlike the highly redundant ABA core signaling components of Arabidopsis, Brachypodium encompasses a more compact and specialized ABA receptor apparatus. This organization may contribute to plant adaptations to ecological niches. These results lay the groundwork for targeting the prominent ABA receptors for stress perception in grasses, and reveal functional differences and commonalities between monocots and eudicots.

Overexpression of AtPAD4 in transgenic Brachypodium distachyon enhances resistance to Puccinia brachypodii

Brachypodium distachyon (L.) has recently emerged as a model for temperate grasses for investigating the molecular basis of plant-pathogen interactions. Phytoalexin deficient 4 (PAD4) plays a regulatory role in mediating expression of genes involved in plant defence. In this research, we generated transgenic B. distachyon plants constitutively overexpressing AtPAD4. Two transgenic B. distachyon lines were verified using PCR and GUS phenotype. Constitutive expression of AtPAD4 in B. distachyon enhanced resistance to Puccinia brachypodii. P. brachypodii generated less urediniospores on transgenic than on wild-type plants. AtPAD4 overexpression enhanced salicylic acid (SA) levels in B. distachyon-infected tissues. qRT-PCR showed that expression of pathogenesis-related 1 (PR1) and other defence-related genes were up-regulated in transformed B. distachyon following infection with P. brachypodii. Our results indicate that AtPAD4 overexpression in B. distachyon plants led to SA accumulation and induced PR gene expression that reduced the rate of colonisation by P. brachypodii.

Expression atlas and comparative coexpression network analyses reveal important genes involved in the formation of lignified cell wall in Brachypodium distachyon

While Brachypodium distachyon (Brachypodium) is an emerging model for grasses, no expression atlas or gene coexpression network is available. Such tools are of high importance to provide insights into the function of Brachypodium genes. We present a detailed Brachypodium expression atlas, capturing gene expression in its major organs at different developmental stages. The data were integrated into a large-scale coexpression database ( www.gene2function.de), enabling identification of duplicated pathways and conserved processes across 10 plant species, thus allowing genome-wide inference of gene function. We highlight the importance of the atlas and the platform through the identification of duplicated cell wall modules, and show that a lignin biosynthesis module is conserved across angiosperms. We identified and functionally characterised a putative ferulate 5-hydroxylase gene through overexpression of it in Brachypodium, which resulted in an increase in lignin syringyl units and reduced lignin content of mature stems, and led to improved saccharification of the stem biomass. Our Brachypodium expression atlas thus provides a powerful resource to reveal functionally related genes, which may advance our understanding of important biological processes in grasses.

Sequencing and functional validation of the JGI Brachypodium distachyon T-DNA collection

Due to a large and growing collection of genomic and experimental resources, Brachypodium distachyon has emerged as a powerful experimental model for the grasses. To add to these resources we sequenced 21 165 T-DNA lines, 15 569 of which were produced in this study. This increased the number of unique insertion sites in the T-DNA collection by 21 078, bringing the overall total to 26 112. Thirty-seven per cent (9754) of these insertion sites are within genes (including untranslated regions and introns) and 28% (7217) are within 500 bp of a gene. Approximately 31% of the genes in the v.2.1 annotation have been tagged in this population. To demonstrate the utility of this collection, we phenotypically characterized six T-DNA lines with insertions in genes previously shown in other systems to be involved in cellulose biosynthesis, hemicellulose biosynthesis, secondary cell wall development, DNA damage repair, wax biosynthesis and chloroplast synthesis. In all cases, the phenotypes observed supported previous studies, demonstrating the utility of this collection for plant functional genomics. The Brachypodium T-DNA collection can be accessed at http://jgi.doe.gov/our-science/science-programs/plant-genomics/brachypodium/brachypodium-t-dna-collection/.

Interactive effects of water limitation and elevated temperature on the physiology, development and fitness of diverse accessions of Brachypodium distachyon

An enduring question in plant physiology and evolution is how single genotypes of plants optimize performance in diverse, often highly variable, environments. We grew 35 natural accessions of the grass Brachypodium distachyon in four environments in the glasshouse, contrasting soil water deficit, elevated temperature and their interaction. We modeled treatment, genotype and interactive effects on leaf-level and whole-plant traits, including fecundity. We also assessed the relationship between glasshouse-measured traits and parameters related to climate at the place of origin. We found abundant genetic variation in both constitutive and induced traits related to plant-water relations. Most traits showed strong interaction between temperature and water availability, and we observed genotype-by-environment interaction for several traits. Notably, leaf free proline abundance showed a strong effect of genotype × temperature × water. We found strong associations between phenology, biomass and water use efficiency (WUE) with parameters describing climate of origin. Plants respond to multiple stressors in ways not directly predictable from single stressors, underscoring the complex and trait-specific mechanisms of environmental response. Climate-trait correlations support a role for WUE and phenology in local adaptation to climate in B. distachyon.

Plant growth responses to elevated atmospheric CO2 are increased by phosphorus sufficiency but not by arbuscular mycorrhizas

Capturing the full growth potential in crops under future elevated CO2 (eCO2) concentrations would be facilitated by improved understanding of eCO2 effects on uptake and use of mineral nutrients. This study investigates interactions of eCO2, soil phosphorus (P), and arbuscular mycorrhizal (AM) symbiosis in Medicago truncatula and Brachypodium distachyon grown under the same conditions. The focus was on eCO2 effects on vegetative growth, efficiency in acquisition and use of P, and expression of phosphate transporter (PT) genes. Growth responses to eCO2 were positive at P sufficiency, but under low-P conditions they ranged from non-significant in M. truncatula to highly significant in B. distachyon Growth of M. truncatula was increased by AM at low P conditions at both CO2 levels and eCO2×AM interactions were sparse. Elevated CO2 had small effects on P acquisition, but enhanced conversion of tissue P into biomass. Expression of PT genes was influenced by eCO2, but effects were inconsistent across genes and species. The ability of eCO2 to partly mitigate P limitation-induced growth reductions in B. distachyon was associated with enhanced P use efficiency, and requirements for P fertilizers may not increase in such species in future CO2-rich climates.

The arrangement of Brachypodium distachyon chromosomes in interphase nuclei

The spatial organization of chromatin within the interphase nucleus and the interactions between chromosome territories (CTs) are essential for various biological processes, such as DNA replication, transcription, and repair. However, detailed data about the CT arrangement in monocotyledonous plants are scarce. In this study, chromosome painting was used to analyse the distribution and associations of individual chromosomes in the 3-D preserved nuclei of Brachypodium distachyon root cells in order to determine the factors that may have an impact on the homologous CT arrangement. It was shown that the frequency of CT association is linked to the steric constraints imposed by the limited space within the nucleus and may depend on chromosome size and morphology as well as on the nuclear shape. Furthermore, in order to assess whether the distribution of interphase chromosomes is random or is subject to certain patterns, a comparison between the experimental data and the results of a computer simulation (ChroTeMo), which was based on a fully probabilistic distribution of the CTs, was performed. This comparison revealed that homologous chromosome arm CTs associate more often than if they were randomly arranged inside the interphase nucleus.