Dual role of BdMUTE during stomatal development in the model grass Brachypodium distachyon.

Grasses form morphologically derived, four-celled stomata, where two dumbbell-shaped guard cells (GCs) are flanked by two lateral subsidiary cells (SCs). This innovative form enables rapid opening and closing kinetics and efficient plant-atmosphere gas exchange. The mobile bHLH transcription factor MUTE is required for SC formation in grasses. Yet whether and how MUTE also regulates GC development and whether MUTE mobility is required for SC recruitment is unclear. Here, we transgenically impaired BdMUTE mobility from GC to SC precursors in the emerging model grass Brachypodium distachyon. Our data indicate that reduced BdMUTE mobility severely affected the spatiotemporal coordination of GC and SC development. Furthermore, although BdMUTE has a cell-autonomous role in GC division orientation, complete dumbbell morphogenesis of GCs required SC recruitment. Finally, leaf-level gas exchange measurements showed that dosage-dependent complementation of the four-celled grass morphology was mirrored in a gradual physiological complementation of stomatal kinetics. Together, our work revealed a dual role of grass MUTE in regulating GC division orientation and SC recruitment, which in turn is required for GC morphogenesis and the rapid kinetics of grass stomata.

PHYTOCHROME C regulation of photoperiodic flowering via PHOTOPERIOD1 is mediated by EARLY FLOWERING 3 in Brachypodium distachyon.

Daylength sensing in many plants is critical for coordinating the timing of flowering with the appropriate season. Temperate climate-adapted grasses such as Brachypodium distachyon flower during the spring when days are becoming longer. The photoreceptor PHYTOCHROME C is essential for long-day (LD) flowering in B. distachyon. PHYC is required for the LD activation of a suite of genes in the photoperiod pathway including PHOTOPERIOD1 (PPD1) that, in turn, result in the activation of FLOWERING LOCUS T (FT1)/FLORIGEN, which causes flowering. Thus, B. distachyon phyC mutants are extremely delayed in flowering. Here we show that PHYC-mediated activation of PPD1 occurs via EARLY FLOWERING 3 (ELF3), a component of the evening complex in the circadian clock. The extreme delay of flowering of the phyC mutant disappears when combined with an elf3 loss-of-function mutation. Moreover, the dampened PPD1 expression in phyC mutant plants is elevated in phyC/elf3 mutant plants consistent with the rapid flowering of the double mutant. We show that loss of PPD1 function also results in reduced FT1 expression and extremely delayed flowering consistent with results from wheat and barley. Additionally, elf3 mutant plants have elevated expression levels of PPD1, and we show that overexpression of ELF3 results in delayed flowering associated with a reduction of PPD1 and FT1 expression, indicating that ELF3 represses PPD1 transcription consistent with previous studies showing that ELF3 binds to the PPD1 promoter. Indeed, PPD1 is the main target of ELF3-mediated flowering as elf3/ppd1 double mutant plants are delayed flowering. Our results indicate that ELF3 operates downstream from PHYC and acts as a repressor of PPD1 in the photoperiod flowering pathway of B. distachyon.

Natural variation in Brachypodium distachyon responses to combined abiotic stresses.

The demand for agricultural production is becoming more challenging as climate change increases global temperature and the frequency of extreme weather events. This study examines the phenotypic variation of 149 accessions of Brachypodium distachyon under drought, heat, and the combination of stresses. Heat alone causes the largest amounts of tissue damage while the combination of stresses causes the largest decrease in biomass compared to other treatments. Notably, Bd21-0, the reference line for B. distachyon, did not have robust growth under stress conditions, especially the heat and combined drought and heat treatments. The climate of origin was significantly associated with B. distachyon responses to the assessed stress conditions. Additionally, a GWAS found loci associated with changes in plant height and the amount of damaged tissue under stress. Some of these SNPs were closely located to genes known to be involved in responses to abiotic stresses and point to potential causative loci in plant stress response. However, SNPs found to be significantly associated with a response to heat or drought individually are not also significantly associated with the combination of stresses. This, with the phenotypic data, suggests that the effects of these abiotic stresses are not simply additive, and the responses to the combined stresses differ from drought and heat alone.

Comparative and functional analysis unveils the contribution of photoperiod to DNA methylation, sRNA accumulation, and gene expression variations in short-day and long-day grasses.

Photoperiod employs complicated networks to regulate various developmental processes in plants, including flowering transition. However, the specific mechanisms by which photoperiod affects epigenetic modifications and gene expression variations in plants remain elusive. In this study, we conducted a comprehensive analysis of DNA methylation, small RNA (sRNA) accumulation, and gene expressions under different daylengths in facultative long-day (LD) grass Brachypodium distachyon and short-day (SD) grass rice. Our results showed that while overall DNA methylation levels were minimally affected by different photoperiods, CHH methylation levels were repressed under their favorable light conditions, particularly in rice. We identified numerous differentially methylated regions (DMRs) that were influenced by photoperiod in both plant species. Apart from differential sRNA clusters, we observed alterations in the expression of key components of the RNA-directed DNA methylation pathway, DNA methyltransferases, and demethylases, which may contribute to the identified photoperiod-influenced CHH DMRs. Furthermore, we identified many differentially expressed genes in response to different daylengths, some of which were associated with the DMRs. Notably, we discovered a photoperiod-responsive gene MYB11 in the transcriptome of B. distachyon, and further demonstrated its role as a flowering inhibitor by repressing FT1 transcription. Together, our comparative and functional analysis sheds light on the effects of daylength on DNA methylation, sRNA accumulation, and gene expression variations in LD and SD plants, thereby facilitating better designing breeding programs aimed at developing high-yield crops that can adapt to local growing seasons.

Information theory and machine learning illuminate large-scale metabolomic responses of Brachypodium distachyon to environmental change.

Plant responses to environmental change are mediated via changes in cellular metabolomes. However, <5% of signals obtained from liquid chromatography tandem mass spectrometry (LC-MS/MS) can be identified, limiting our understanding of how metabolomes change under biotic/abiotic stress. To address this challenge, we performed untargeted LC-MS/MS of leaves, roots, and other organs of Brachypodium distachyon (Poaceae) under 17 organ-condition combinations, including copper deficiency, heat stress, low phosphate, and arbuscular mycorrhizal symbiosis. We found that both leaf and root metabolomes were significantly affected by the growth medium. Leaf metabolomes were more diverse than root metabolomes, but the latter were more specialized and more responsive to environmental change. We found that 1 week of copper deficiency shielded the root, but not the leaf metabolome, from perturbation due to heat stress. Machine learning (ML)-based analysis annotated approximately 81% of the fragmented peaks versus approximately 6% using spectral matches alone. We performed one of the most extensive validations of ML-based peak annotations in plants using thousands of authentic standards, and analyzed approximately 37% of the annotated peaks based on these assessments. Analyzing responsiveness of each predicted metabolite class to environmental change revealed significant perturbations of glycerophospholipids, sphingolipids, and flavonoids. Co-accumulation analysis further identified condition-specific biomarkers. To make these results accessible, we developed a visualization platform on the Bio-Analytic Resource for Plant Biology website (https://bar.utoronto.ca/efp_brachypodium_metabolites/cgi-bin/efpWeb.cgi), where perturbed metabolite classes can be readily visualized. Overall, our study illustrates how emerging chemoinformatic methods can be applied to reveal novel insights into the dynamic plant metabolome and stress adaptation.

Duplicated antagonistic EPF peptides optimize grass stomatal initiation.

Peptide signaling has emerged as a key component of plant growth and development, including stomatal patterning, which is crucial for plant productivity and survival. Although exciting progress has been made in understanding EPIDERMAL PATTERNING FACTOR (EPF) signaling in Arabidopsis, the mechanisms by which EPF peptides control different stomatal patterns and morphologies in grasses are poorly understood. Here, by examining expression patterns, overexpression transgenics and cross-species complementation, the antagonistic stomatal ligands orthologous to Arabidopsis AtEPF2 and AtSTOMAGEN/AtEPFL9 peptides were identified in Triticum aestivum (wheat) and the grass model organism Brachypodium distachyon. Application of bioactive BdEPF2 peptides inhibited stomatal initiation, but not the progression or differentiation of stomatal precursors in Brachypodium. Additionally, the inhibitory roles of these EPF peptides during grass stomatal development were suppressed by the contrasting positive action of the BdSTOMAGEN peptide in a dose-dependent manner. These results not only demonstrate how conserved EPF peptides that control different stomatal patterns exist in nature, but also suggest new strategies to improve crop yield through the use of plant-derived antagonistic peptides that optimize stomatal density on the plant epidermis.

Altering the substitution and cross-linking of glucuronoarabinoxylans affects cell wall architecture in Brachypodium distachyon.

The Poaceae family of plants provides cereal crops that are critical for human and animal nutrition, and also, they are an important source of biomass. Interacting plant cell wall components give rise to recalcitrance to digestion; thus, understanding the wall molecular architecture is important to improve biomass properties. Xylan is the main hemicellulose in grass cell walls. Recently, we reported structural variation in grass xylans, suggesting functional specialisation and distinct interactions with cellulose and lignin. Here, we investigated the functions of these xylans by perturbing the biosynthesis of specific xylan types. We generated CRISPR/Cas9 knockout mutants in Brachypodium distachyon XAX1 and GUX2 genes involved in xylan substitution. Using carbohydrate gel electrophoresis, we identified biochemical changes in different xylan types. Saccharification, cryo-SEM, subcritical water extraction and ssNMR were used to study wall architecture. BdXAX1A and BdGUX2 enzymes modify different types of grass xylan. Brachypodium mutant walls are likely more porous, suggesting the xylan substitutions directed by both BdXAX1A and GUX2 enzymes influence xylan-xylan and/or xylan-lignin interactions. Since xylan substitutions influence wall architecture and digestibility, our findings open new avenues to improve cereals for food and to use grass biomass for feed and the production of bioenergy and biomaterials.

A stress-inducible protein regulates drought tolerance and flowering time in Brachypodium and Arabidopsis.

To cope with environmental stresses and ensure maximal reproductive success, plants have developed strategies to adjust the timing of their transition to reproductive growth. This has a substantial impact on the stress resilience of crops and ultimately on agricultural productivity. Here, we report a previously uncharacterized, plant-specific gene family designated as Regulator of Flowering and Stress (RFS). Overexpression of the BdRFS gene in Brachypodium distachyon delayed flowering, increased biomass accumulation, and promoted drought tolerance, whereas clustered regularly interspaced short palindromic repeats/CRISPR-associated protein 9 (CRISPR/Cas9)-mediated knockout mutants exhibited opposite phenotypes. A double T-DNA insertional mutant in the two Arabidopsis (Arabidopsis thaliana) homologs replicated the effects on flowering and water deprivation seen in the B. distachyon CRISPR knockout lines, highlighting the functional conservation of the family between monocots and dicots. Lipid analysis of B. distachyon and Arabidopsis revealed that digalactosyldiacylglycerol (DGDG) and phosphatidylcholine (PC) contents were significantly, and reciprocally, altered in overexpressor and knockout mutants. Importantly, alteration of C16:0-containing PC, a Flowering Locus T-interacting lipid, associated with flowering phenotype, with elevated levels corresponding to earlier flowering. Co-immunoprecipitation analysis suggested that BdRFS interacts with phospholipase Dα1 as well as several other abscisic acid-related proteins. Furthermore, reduction of C18:3 fatty acids in DGDG corresponded with reduced jasmonic acid metabolites in CRISPR mutants. Collectively, we suggest that stress-inducible RFS proteins represent a regulatory component of lipid metabolism that impacts several agronomic traits of biotechnological importance.

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.

MiR398-regulated antioxidants contribute to Bamboo mosaic virus accumulation and symptom manifestation.

Virus infections that cause mosaic or mottling in leaves commonly also induce increased levels of reactive oxygen species (ROS). However, how ROS contributes to symptoms is less well documented. Bamboo mosaic virus (BaMV) causes chlorotic mosaic symptoms in both Brachypodium distachyon and Nicotiana benthamiana. The BaMV △CPN35 mutant with an N-terminal deletion of its coat protein gene exhibits asymptomatic infection independently of virus titer. Histochemical staining of ROS in mock-, BaMV-, and BaMV△CPN35-infected leaves revealed that hydrogen peroxide (H2O2) accumulated solely in BaMV-induced chlorotic spots. Moreover, exogenous H2O2 treatment enhanced yellowish chlorosis in BaMV-infected leaves. Both BaMV and BaMV△CPN35 infection could induce the expression of Cu/Zu superoxide dismutase (CSD) antioxidants at messenger RNA and protein level. However, BaMV triggered the abundant accumulation of full-length NbCSD2 preprotein (prNbCSD2, without transit peptide cleavage), whereas BaMV△CPN35 induced a truncated prNbCSD2. Confocal microscopy showed that majority of NbCSD2-green fluorescent protein (GFP) predominantly localized in the cytosol upon BaMV infection, but BaMV△CPN35 infection tended to cause NbCSD2-GFP to remain in chloroplasts. By 5'-RNA ligase-mediated rapid amplification of cDNA ends, we validated CSDs are the targets of miR398 in vivo. Furthermore, BaMV infection increased the level of miR398, while the level of BaMV titer was regulated positively by miR398 but negatively by CSD2. In contrast, overexpression of cytosolic form NbCSD2, impairing the transport into chloroplasts, greatly enhanced BaMV accumulation. Taken together, our results indicate that induction of miR398 by BaMV infection may facilitate viral titer accumulation, and cytosolic prNbCSD2 induction may contribute to H2O2 accumulation, resulting in the development of BaMV chlorotic symptoms in plants.

REGULATOR OF FLOWERING AND STRESS manipulates stomatal density and size in Brachypodium.

Stomata are crucial for gas exchange and water evaporation, and environmental stimuli influence their density (SD) and size (SS). Although genes and mechanisms underlying stomatal development have been elucidated, stress-responsive regulators of SD and SS are less well-known. Previous studies have shown that the stress-inducible Brachypodium RFS (REGULATOR OF FLOWERING AND STRESS, BdRFS) gene affects heading time and enhances drought tolerance by reducing leaf water loss. Here, we report that overexpression lines (OXs) of BdRFS have reduced SD and increased SS, regardless of soil water status. Furthermore, biomass and plant water content of OXs were significantly increased compared to wild type. CRISPR/Cas9-mediated BdRFS knockout mutant (KO) exhibited the opposite stomatal characteristics and biomass changes. Reverse transcription-quantitative polymerase chain reaction analysis revealed that expression of BdICE1 was reversely altered in OXs and KO, pointing to a potential cause for the observed changes in stomatal phenotypes. Stomatal and transcriptional changes were not observed in the Arabidopsis rfs double mutant. Taken together, RFS is a novel regulator of SD and SS and is a promising candidate for genetic engineering of climate-resilient crops.

The Aux/IAA protein TaIAA15-1A confers drought tolerance in Brachypodium by regulating abscisic acid signal pathway.

KEY MESSAGE: Overexpression of the Aux/IAA protein TaIAA15-1A from wheat improves drought tolerance by regulating the ABA signalling pathway in transgenic Brachypodium. Drought is a major abiotic stress that causes severe crop yield loss. Aux/IAA genes have been shown to be involved in drought stress responses. However, to the best of our knowledge, there has been little research on the molecular mechanism of the wheat Aux/IAA gene in the context of drought tolerance. In this study, we found that expression of the wheat Aux/IAA gene TaIAA15-1A was upregulated by PEG6000, NaCl, SA, JA, IAA and ABA. Transgenic plants overexpressing TaIAA15-1A showed higher drought tolerance than wild-type (WT) plants. The physiological analyses showed that the transgenic lines exhibited a higher survival rate, shoot length, and relative water content than the WT plants. The activities of superoxide dismutase (SOD), catalase (CAT) and peroxidase (POD) were enhanced in transgenic lines, causing a reduction in the hydrogen peroxide (H2O2) and superoxide anion radical (O2-) contents. Transcriptome analysis showed that TaIAA15-1A overexpression alters the expression of these genes involved in the auxin signalling pathway, ABA signalling pathway, phenolamides and antioxidant pathways. The results of exogenous ABA treatment suggested that TaIAA15-1A overexpression increased sensitivity to ABA at the germination and postgermination stages compared to WT plants. These results indicate that TaIAA15-1A plays a positive role in plant drought tolerance by regulating ABA-related genes and improving antioxidative stress ability and has potential application in genetically modified crops.

Characterization of diterpene synthase genes in Brachypodium distachyon, a monocotyledonous model plant, provides evolutionary insight into their multiple homologs in cereals.

Gibberellins are diterpenoid phytohormones that regulate plant growth, and are biosynthesized from a diterpene intermediate, ent-kaurene, which is produced from geranylgeranyl diphosphate via ent-copalyl diphosphate (ent-CDP). The successive 2 cyclization reactions are catalyzed by 2 distinct diterpene synthases, ent-CDP synthase (ent-CPS) and ent-kaurene synthase (KS). Various diterpene synthase genes involved in specialized metabolism were likely created through duplication and neofunctionalization of gibberellin-biosynthetic ent-CPS and KS genes in crops. Brachypodium distachyon is a monocotyledonous species that is a model plant in grasses. We herein found 1 ent-CPS gene homolog BdCPS and 4 tandemly arrayed KS-like genes BdKS1, KSL2, KSL3, and KSL4 in the B. distachyon genome, a simpler collection of paralogs than in crops. Phylogenetic and biochemical analyses showed that BdCPS and BdKS1 are responsible for gibberellin biosynthesis. BdKSL2 and BdKSL3 are suggested to be involved in specialized diterpenoid metabolism. Moreover, we restored KS activity of BdKSL2 through amino acid substitution.

Submergence Stress Alters the Expression of Clock Genes and Configures New Zeniths and Expression of Outputs in Brachypodium distachyon.

Plant networks of oscillating genes coordinate internal processes with external cues, contributing to increased fitness. We hypothesized that the response to submergence stress may dynamically change during different times of the day. In this work, we determined the transcriptome (RNA sequencing) of the model monocotyledonous plant, Brachypodium distachyon, during a day of submergence stress, low light, and normal growth. Two ecotypes of differential tolerance, Bd21 (sensitive) and Bd21-3 (tolerant), were included. We submerged 15-day-old plants under a long-day diurnal cycle (16 h light/8 h dark) and collected samples after 8 h of submergence at ZT0 (dawn), ZT8 (midday), ZT16 (dusk), ZT20 (midnight), and ZT24 (dawn). Rhythmic processes were enriched both with up- and down-regulated genes, and clustering highlighted that the morning and daytime oscillator components (PRRs) show peak expression in the night, and a decrease in the amplitude of the clock genes (GI, LHY, RVE) was observed. Outputs included photosynthesis-related genes losing their known rhythmic expression. Up-regulated genes included oscillating suppressors of growth, hormone-related genes with new late zeniths (e.g., JAZ1, ZEP), and mitochondrial and carbohydrate signaling genes with shifted zeniths. The results highlighted genes up-regulated in the tolerant ecotype such as METALLOTHIONEIN3 and ATPase INHIBITOR FACTOR. Finally, we show by luciferase assays that Arabidopsis thaliana clock genes are also altered by submergence changing their amplitude and phase. This study can guide the research of chronocultural strategies and diurnal-associated tolerance mechanisms.

Mixed-Linkage Glucan Is the Main Carbohydrate Source and Starch Is an Alternative Source during Brachypodium Grain Germination.

Seeds of the model grass Brachypodium distachyon are unusual because they contain very little starch and high levels of mixed-linkage glucan (MLG) accumulated in thick cell walls. It was suggested that MLG might supplement starch as a storage carbohydrate and may be mobilised during germination. In this work, we observed massive degradation of MLG during germination in both endosperm and nucellar epidermis. The enzymes responsible for the MLG degradation were identified in germinated grains and characterized using heterologous expression. By using mutants targeting MLG biosynthesis genes, we showed that the expression level of genes coding for MLG and starch-degrading enzymes was modified in the germinated grains of knocked-out cslf6 mutants depleted in MLG but with higher starch content. Our results suggest a substrate-dependent regulation of the storage sugars during germination. These overall results demonstrated the function of MLG as the main carbohydrate source during germination of Brachypodium grain. More astonishingly, cslf6 Brachypodium mutants are able to adapt their metabolism to the lack of MLG by modifying the energy source for germination and the expression of genes dedicated for its use.

A genetically encoded biosensor reveals spatiotemporal variation in cellular phosphate content in Brachypodium distachyon mycorrhizal roots.

Arbuscular mycorrhizal (AM) symbiosis is accompanied by alterations to root cell metabolism and physiology, and to the pathways of orthophosphate (Pi) entry into the root, which increase with Pi delivery to cortical cells via arbuscules. How AM symbiosis influences the Pi content and Pi response dynamics of cells in the root cortex and epidermis is unknown. Using fluorescence resonance energy transfer (FRET)-based Pi biosensors, we mapped the relative cytosolic and plastidic Pi content of Brachypodium distachyon mycorrhizal root cells, analyzed responses to extracellular Pi and traced extraradical hyphae-mediated Pi transfer to colonized cells. Colonized cortical cells had a higher cytosolic Pi content relative to noncolonized cortical and epidermal cells, while plastidic Pi content was highest in cells at the infection front. Pi application to the entire mycorrhizal root resulted in transient changes in cytosolic Pi that differed in direction and magnitude depending on cell type and arbuscule status; cells with mature arbuscules showed a substantial transient increase in cytosolic Pi while those with collapsed arbuscules showed a decrease. Directed Pi application to extraradical hyphae resulted in measurable changes in cytosolic Pi of colonized cells 18 h after application. Our experiments reveal that cells within a mycorrhizal root vary in Pi content and Pi response dynamics.

N-dependent dynamics of root growth and nitrate and ammonium uptake are altered by the bacterium Herbaspirillum seropedicae in the cereal model Brachypodium distachyon.

Nitrogen (N) fixation in cereals by root-associated bacteria is a promising solution for reducing use of chemical N fertilizers in agriculture. However, plant and bacterial responses are unpredictable across environments. We hypothesized that cereal responses to N-fixing bacteria are dynamic, depending on N supply and time. To quantify the dynamics, a gnotobiotic, fabricated ecosystem (EcoFAB) was adapted to analyse N mass balance, to image shoot and root growth, and to measure gene expression of Brachypodium distachyon inoculated with the N-fixing bacterium Herbaspirillum seropedicae. Phenotyping throughput of EcoFAB-N was 25-30 plants h-1 with open software and imaging systems. Herbaspirillum seropedicae inoculation of B. distachyon shifted root and shoot growth, nitrate versus ammonium uptake, and gene expression with time; directions and magnitude depended on N availability. Primary roots were longer and root hairs shorter regardless of N, with stronger changes at low N. At higher N, H. seropedicae provided 11% of the total plant N that came from sources other than the seed or the nutrient solution. The time-resolved phenotypic and molecular data point to distinct modes of action: at 5 mM NH4NO3 the benefit appears through N fixation, while at 0.5 mM NH4NO3 the mechanism appears to be plant physiological, with H. seropedicae promoting uptake of N from the root medium.Future work could fine-tune plant and root-associated microorganisms to growth and nutrient dynamics.

Genes involved in mRNA surveillance are induced in Brachypodium distachyon under cadmium toxicity.

BACKGROUND: Cd accumulation in plant cells results in dramatic problems including oxidative stress and inhibition of vital enzymes. It also affects mineral uptakes by disrupting membrane permeability. Interaction among Cd and other plant nutrient elements changes the nutritional contents of crops and reduces their yield. METHODS AND RESULTS: In the present study, Cd stress in Brachypodium distachyon led to the upregulation of some heavy metal transport genes (influx or efflux) encoding cation-efflux proteins, heavy metal-associated proteins and NRAMP proteins. The Arabidopsis orthologs of the differentially expressed B. distachyon genes (DEGs) under Cd toxicity were identified, which exhibited Bradi4g26905 was an ortholog of AtALY1-2. Detailed co-expression network and gene ontology analyses found the potential involvement of the mRNA surveillance pathway in Cd tolerance in B. distachyon. These genes were shown to be downregulated by sulfur (S) deficiency. CONCLUSIONS: This is the first transcriptomic study investigating the effect of Cd toxicity in B. distachyon, a model plant for genomic studies in Poaceae (Gramineae) species. The results are expected to provide valuable information for more comprehensive research related to heavy metal toxicity in plants.

Phylogeny and evolution of plant Phytochrome Interacting Factors (PIFs) gene family and functional analyses of PIFs in Brachypodium distachyon.

KEY MESSAGES: Plant PIFs have been characterized, WGDs contributed to the expansion of class II PIFs; BdPIFs localized in the nucleus; BdPIF4/5C most likely response to high temperature and light stress. Phytochrome interacting factors (PIFs) belong to a small subset of basic helix-loop-helix (bHLH) transcription factors (TFs). As cellular signaling hubs, PIFs integrate multiple external and internal signals to orchestrate the regulation of the transcriptional network, thereby actuating the pleiotropic aspects of downstream morphogenesis. Nevertheless, the origin, phylogeny and function of plant PIFs are not well understood. To elucidate their evolution history and biological function, the comprehensive genomic analysis of the PIF genes was conducted using 40 land plant genomes plus additionally four alga lineages and also performed their gene organizations, sequence features and expression patterns in different subfamilies. In this study, phylogenetic analysis displayed that 246 PIF gene members retrieved from all embryophytes could be divided into three main clades, which were further felled into five distinct classes (Class I-V). The duplications of Class II PIFs were associated specially with whole genome duplication (WGD) events during the plant evolution process. Sequence analysis showed that PIF proteins had a conserved APB motif, and its crucial amino acid residues were relatively high proportion in the average abundance. As expected, subcellular localization analysis revealed that all BdPIF proteins were localized to the nucleus. Especially, BdPIF4/5C showed the highest expression level at high temperature, and the most significant hypocotyl elongation phenotype of overexpression of BdPIFs in Arabidopsis, which was consistent with the function and phenotype of AtPIF4. In brief, our findings provide a novel perspective on the origin and evolutionary history of plant PIFs, and lays a foundation for further investigation on its functions in plant growth and development.

BdGUCD1 and Cyclic GMP Are Required for Responses of Brachypodium distachyon to Fusarium pseudograminearum in the Mechanism Involving Jasmonate.

Guanosine 3',5'-cyclic monophosphate (cGMP) is an important signaling molecule in plants. cGMP and guanylyl cyclases (GCs), enzymes that catalyze the synthesis of cGMP from GTP, are involved in several physiological processes and responses to environmental factors, including pathogen infections. Using in vitro analysis, we demonstrated that recombinant BdGUCD1 is a protein with high guanylyl cyclase activity and lower adenylyl cyclase activity. In Brachypodium distachyon, infection by Fusarium pseudograminearum leads to changes in BdGUCD1 mRNA levels, as well as differences in endogenous cGMP levels. These observed changes may be related to alarm reactions induced by pathogen infection. As fluctuations in stress phytohormones after infection have been previously described, we performed experiments to determine the relationship between cyclic nucleotides and phytohormones. The results revealed that inhibition of cellular cGMP changes disrupts stress phytohormone content and responses to pathogen. The observations made here allow us to conclude that cGMP is an important element involved in the processes triggered as a result of infection and changes in its levels affect jasmonic acid. Therefore, stimuli-induced transient elevation of cGMP in plants may play beneficial roles in priming an optimized response, likely by triggering the mechanisms of feedback control.