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.

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.

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.

Systematic analysis of Prx genes in the Brachypodium genus and their expression pattern under abiotic constraints.

Peroxiredoxins (Prx) are ubiquitous peroxidases required for the removal of excess free radicals produced under stress conditions. Peroxiredoxin genes (Prx) in the Brachypodium genus were identified using bioinformatics tools and their expression profiles were determined under abiotic stress using RT-qPCR. The promoter regions of Prx genes contain several cis-acting elements related to stress response. In silico expression analysis showed that B. distachyon Prx genes (BdPrx) are tissue specific. RT-qPCR analysis revealed their differential expression when exposed to salt or PEG-induced dehydration stress. In addition, the upregulation of BdPrx genes was accompanied by accumulation of H2 O2 . Exogenous application of H2 O2 induced expression of almost all BdPrx genes. The identified molecular interaction network indicated that Prx proteins may contribute to abiotic stress tolerance by regulating key enzymes involved in lignin biosynthesis. Overall, our findings suggest the potential role of Prx genes in abiotic stress tolerance and lay the foundation for future functional analyses aiming to engineer genetically improved cereal lines.

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.

Phytochromes transmit photoperiod information via the evening complex in Brachypodium.

BACKGROUND: Daylength is a key seasonal cue for animals and plants. In cereals, photoperiodic responses are a major adaptive trait, and alleles of clock genes such as PHOTOPERIOD1 (PPD1) and EARLY FLOWERING3 (ELF3) have been selected for in adapting barley and wheat to northern latitudes. How monocot plants sense photoperiod and integrate this information into growth and development is not well understood. RESULTS: We find that phytochrome C (PHYC) is essential for flowering in Brachypodium distachyon. Conversely, ELF3 acts as a floral repressor and elf3 mutants display a constitutive long day phenotype and transcriptome. We find that ELF3 and PHYC occur in a common complex. ELF3 associates with the promoters of a number of conserved regulators of flowering, including PPD1 and VRN1. Consistent with observations in barley, we are able to show that PPD1 overexpression accelerates flowering in short days and is necessary for rapid flowering in response to long days. PHYC is in the active Pfr state at the end of the day, but we observe it undergoes dark reversion over the course of the night. CONCLUSIONS: We propose that PHYC acts as a molecular timer and communicates information on night-length to the circadian clock via ELF3.

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.

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.

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.

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.

Analysis of Brachypodium distachyon UVR8 reveals conservation in UV-B receptors.

The Ultraviolet Resistance Locus 8 (UVR8) in plants recognizes ultraviolet-B (UV-B) light and plays a crucial role in regulating plant growth through a series of signal transduction events. However, the UVR8 in monocotyledon crops has not yet been systematically analysed. We identified BdUVR8 (BRADI_3g45740) from the genome of Brachypodium distachyon, a relative of wheat, by analysing the phylogenetic tree, the gene expression pattern, detecting accumulation of UV-B response metabolites, and checking for phenotype recovery. The BdUVR8 protein sequence is similar to the known UVR8 of other species. The phylogenetic tree of UVR8 shows clear divergence between dicotyledons and monocotyledons. Expression analysis revealed that UV-B downregulates BdUVR8 by 70% and upregulates the chalcone synthase (BdCHS) gene 3.4-fold in B. distachyon. The pCAMBIA1300::BdUVR8-mCherry construct introduced into Arabidopsis uvr8 mutants showed that the BdUVR8 protein is localized in the cytoplasm and translocates into the nucleus in response to UV-B irradiation. The introduction of BdUVR8 into uvr8 rescued hypocotyl elongation caused by UV-B and restored expression of HY5, Chalcone synthase, and Flavanone 3-hydroxylase, as well as accumulation of total flavonoids. Together, our results show that BdUVR8 is a photoreceptor that perceives UV-B in B. distachyon.

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.

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.

Identification of Pyrabactin resistance 1-like (PYL) genes in Brachypodium distachyon and functional characterization of BdPYL5.

Abscisic acid (ABA) is an endogenous phytohormone that plays an important role in regulating plant growth, development, and stress response. Pyrabactin resistance 1-like (PYR/PYL) proteins are ABA receptors and core components of ABA signalling in plants. This study identified nine PYL genes in the Brachypodium distachyon genome and they distribute on three chromosomes. Phylogenetical BdPYLs were classified into three clades. 81 protein-protein interactions between 9 BdPYLs and 9 BdPP2C proteins were predicted and 66 pairs were verified by yeast two-hybrid assay previously. Relatively, BdPYL genes are expressed in leaves at high level, and ABA and drought regulate their expression. A homologue of Arabidopsis PYL9, BdPYL5 was selected to overexpress in Arabidopsis to characterize its function. In general, overexpression of BdPYL5 enhanced ABA sensitivity and drought tolerance, implying its conserved function. Our study lays the foundation for further functional elucidation of BdPYL genes.

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.

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.

Targeted and Untargeted Metabolomic Analyses Reveal Organ Specificity of Specialized Metabolites in the Model Grass Brachypodium distachyon.

Brachypodium distachyon, because of its fully sequenced genome, is frequently used as a model grass species. However, its metabolome, which constitutes an indispensable element of complex biological systems, remains poorly characterized. In this study, we conducted comprehensive, liquid chromatography-mass spectrometry (LC-MS)-based metabolomic examination of roots, leaves and spikes of Brachypodium Bd21 and Bd3-1 lines. Our pathway enrichment analysis emphasised the accumulation of specialized metabolites representing the flavonoid biosynthetic pathway in parallel with processes related to nucleotide, sugar and amino acid metabolism. Similarities in metabolite profiles between both lines were relatively high in roots and leaves while spikes showed higher metabolic variance within both accessions. In roots, differences between Bd21 and Bd3-1 lines were manifested primarily in diterpenoid metabolism, while differences within spikes and leaves concerned nucleotide metabolism and nitrogen management. Additionally, sulphate-containing metabolites differentiated Bd21 and Bd3-1 lines in spikes. Structural analysis based on MS fragmentation spectra enabled identification of 93 specialized metabolites. Among them phenylpropanoids and flavonoids derivatives were mainly determined. As compared with closely related barley and wheat species, metabolic profile of Brachypodium is characterized with presence of threonate derivatives of hydroxycinnamic acids.

Mass spectral imaging showing the plant growth-promoting rhizobacteria's effect on the Brachypodium awn.

The plant growth-promoting rhizobacteria (PGPR) on the host plant surface play a key role in biological control and pathogenic response in plant functions and growth. However, it is difficult to elucidate the PGPR effect on plants. Such information is important in biomass production and conversion. Brachypodium distachyon (Brachypodium), a genomics model for bioenergy and native grasses, was selected as a C3 plant model; and the Gram-negative Pseudomonas fluorescens SBW25 (P.) and Gram-positive Arthrobacter chlorophenolicus A6 (A.) were chosen as representative PGPR strains. The PGPRs were introduced to the Brachypodium seed's awn prior to germination, and their possible effects on the seeding and growth were studied using different modes of time-of-flight secondary ion mass spectrometry (ToF-SIMS) measurements, including a high mass-resolution spectral collection and delayed image extraction. We observed key plant metabolic products and biomarkers, such as flavonoids, phenolic compounds, fatty acids, and auxin indole-3-acetic acid in the Brachypodium awns. Furthermore, principal component analysis and two-dimensional imaging analysis reveal that the Brachypodium awns are sensitive to the PGPR, leading to chemical composition and morphology changes on the awn surface. Our results show that ToF-SIMS can be an effective tool to probe cell-to-cell interactions at the biointerface. This work provides a new approach to studying the PGPR effects on awn and shows its potential for the research of plant growth in the future.

Genome-wide identification of YABBY transcription factors in Brachypodium distachyon and functional characterization of Bd DROOPING LEAF.

YABBY transcription factors (TFs) are plant-specific and are characterized by a C2-C2 zinc finger domain at the N-terminus and a YABBY domain at the C-terminus. In this study, eight YABBY genes were identified in the Brachypodium distachyon genome and were unevenly distributed across four chromosomes. Phylogenetic analysis classified BdYABBYs into FIL/YAB3, YAB2, CRC, and INO clades. Sixty-two putative cis-elements were identified in BdYABBY gene putative promoters, among them, CAAT-box, TATA-box, MYB, MYC, ARE, and Box_4 were shared by all. BdYABBY genes are highly expressed in inflorescences, and abiotic stresses regulate their expression. In addition, three transcripts of BdDL were identified. Over-expression in Arabidopsis has shown their different functions in reproductive development, as well as in response to cold stress. Our study lays the foundation for the functional elucidation of BdYABBY genes.