TINY BRANCHED HAIR functions in multicellular trichome development through an ethylene pathway in Cucumis sativus L.

The fruit trichomes of Cucurbitaceae are widely desired in many Asian countries and have been a key determinant of cucumber (Cucumis sativus L.) cultivar selection for commercial production and breeding. However, our understanding of the initiation and development of cucumber trichomes is still limited. Here, we found that the cucumber TINY BRANCHED HAIR (TBH) gene is preferentially expressed in multicellular trichomes. Overexpression of CsTBH in tbh mutants restored the trichome phenotype and increased the percentage of female flowers, whereas silencing of CsTBH in wild-type plants resulted in stunted trichomes with a lower rate of female flowers. Furthermore, we provide evidence that CsTBH can directly bind to the promoters of cucumber 1-Aminocyclopropane-1-Carboxylate Synthase (CsACS) genes and regulate their expression, which affects multicellular trichome development, ethylene accumulation, and sex expression. Two cucumber acs mutants with different trichome morphology and sex morphs compared with their near-isogenic line further support our findings. Collectively, our study provides new information on the molecular mechanism of CsTBH in regulating multicellular trichome development and sex expression through an ethylene pathway.

Multi-omics analysis to visualize the dynamic roles of defense genes in the response of tea plants to gray blight.

Gray blight (GB) is one of the most destructive diseases of tea plants, causing considerable damage and productivity losses; however, the dynamic roles of defense genes during pathogen infection remain largely unclear. To explore the numerous molecular interactions associated with GB stress in tea plants, we employed transcriptome, sRNAome and degradome sequencing from 1 to 13 days post-inoculation (dpi) at 3-day intervals. The transcriptomics results showed that differentially expressed genes (DEGs) related to flavonoid synthesis, such as chalcone synthase (CHS) and phenylalanine ammonia-lyase (PAL), were particularly induced at 4 dpi. Consistent with this, the contents of catechins (especially gallocatechin), which are the dominant flavonoids in tea plants, also increased in the leaves of tea plants infected with GB. Combined analysis of the sRNAome and degradome revealed that microRNAs could mediate tea plant immunity by regulating DEG expression at the post-transcriptional level. Co-expression network analysis demonstrated that miR530b-ethylene responsive factor 96 (ERF96) and miRn211-thaumatin-like protein (TLP) play crucial roles in the response to GB. Accordingly, gene-specific antisense oligonucleotide assays suggested that suppressing ERF96 decreased the levels of reactive oxygen species (ROS), whereas suppressing TLP increased the levels of ROS. Furthermore, ERF96 was induced, but TLP was suppressed, in susceptible tea cultivars. Our results collectively demonstrate that ERF96 is a negative regulator and TLP is a positive regulator in the response of tea plants to GB. Taken together, our comprehensive integrated analysis reveals a dynamic regulatory network linked to GB stress in tea plants and provides candidate genes for improvement of tea plants.

A restorer-of-fertility-like pentatricopeptide repeat protein promotes cytoplasmic male sterility in Arabidopsis thaliana.

Pentatricopeptide repeat (PPR) proteins form a large family of proteins targeted to organelles, where they post-transcriptionally modulate gene expression through binding to specific RNA sequences. Among them, the mitochondria-targeted restorer-of-fertility (Rf) PPRs inhibit peculiar mitochondrial genes that are detrimental to male gametes and cause cytoplasmic male sterility (CMS). Here, we revealed three nuclear loci involved in CMS in a cross between two distant Arabidopsis thaliana strains, Sha and Cvi-0. We identified the causal gene at one of these loci as RFL24, a conserved gene encoding a PPR protein related to known Rf PPRs. By analysing fertile revertants obtained in a male sterile background, we demonstrate that RFL24 promotes pollen abortion, in contrast with the previously described Rf PPRs, which allow pollen to survive in the presence of a sterilizing cytoplasm. We show that the sterility caused by the RFL24 Cvi-0 allele results from higher expression of the gene during early pollen development. Finally, we predict a binding site for RFL24 upstream of two mitochondrial genes, the CMS gene and the important gene cob. These results suggest that the conservation of RFL24 is linked to a primary role of ensuring a proper functioning of mitochondria, and that it was subsequently diverted by the CMS gene to its benefit.

Conservation and diversification of HAIRY MERISTEM gene family in land plants.

The shoot apical meristems (SAMs) of land plants are crucial for plant growth and organ formation. In several angiosperms, the HAIRY MERISTEM (HAM) genes function as key regulators that control meristem development and stem cell homeostasis. To date, the origin and evolutionary history of the HAM family in land plants remains unclear. Potentially shared and divergent functions of HAM family members from angiosperms and non-angiosperms are also not known. In constructing a comprehensive phylogeny of the HAM family, we show that HAM proteins are widely present in land plants and that HAM proteins originated prior to the divergence of bryophytes. The HAM family was duplicated in a common ancestor of angiosperms, leading to two distinct groups: type I and type II. Type-II HAM members are widely present in angiosperms, whereas type-I HAM members were independently lost in different orders of monocots. Furthermore, HAM members from angiosperms and non-angiosperms (including bryophytes, lycophytes, ferns and gymnosperms) are able to replace the role of the type-II HAM genes in Arabidopsis, maintaining established SAMs and promoting the initiation of new stem cell niches. Our results uncover the conserved functions of HAM family members and reveal the conserved regulatory mechanisms underlying HAM expression patterning in meristems, providing insight into the evolution of key stem cell regulators in land plants.

A conservative pathway for coordination of cell wall biosynthesis and cell cycle progression in plants.

The mechanism that coordinates cell growth and cell cycle progression remains poorly understood; in particular, whether the cell cycle and cell wall biosynthesis are coordinated remains unclear. Recently, cell wall biosynthesis and cell cycle progression were reported to respond to wounding. Nonetheless, no genes are reported to synchronize the biosynthesis of the cell wall and the cell cycle. Here, we report that wounding induces the expression of genes associated with cell wall biosynthesis and the cell cycle, and that two genes, AtMYB46 in Arabidopsis thaliana and RrMYB18 in Rosa rugosa, are induced by wounding. We found that AtMYB46 and RrMYB18 promote the biosynthesis of the cell wall by upregulating the expression of cell wall-associated genes, and that both of them also upregulate the expression of a battery of genes associated with cell cycle progression. Ultimately, this response leads to the development of curled leaves of reduced size. We also found that the coordination of cell wall biosynthesis and cell cycle progression by AtMYB46 and RrMYB18 is evolutionarily conservative in multiple species. In accordance with wounding promoting cell regeneration by regulating the cell cycle, these findings also provide novel insight into the coordination between cell growth and cell cycle progression and a method for producing miniature plants.

COPII genes SEC31A/B are essential for gametogenesis and interchangeable in pollen development in Arabidopsis.

In eukaryotes, coat protein complex II (COPII) vesicles mediate anterograde traffic from the endoplasmic reticulum to the Golgi apparatus. Compared to yeasts, plants have multiple COPII coat proteins; however, the functional diversity among them is less well understood. SEC31A and SEC31B are outer coat proteins found in COPII vesicles in Arabidopsis. In this study, we explored the function of SEC31A and compared it with that of SEC31B from various perspectives. SEC31A was widely expressed, but at a significantly lower level than SEC31B. SEC31A-mCherry and SEC31B-GFP exhibited a high co-localization rate in pollen, but a lower rate in growing pollen tubes. The sec31a single mutant exhibited normal growth. SEC31A expression driven by the SEC31B promoter rescued the pollen abortion and infertility observed in sec31b. A sec31asec31b double mutant was unavailable due to lethality of the sec31asec31b gametophyte. Transmission electron microscopy revealed that one quarter of male gametogenesis was arrested at the uninuclear microspore stage, while confocal laser scanning microscopy showed that 1/4 female gametophyte development was suspended at the functional megaspore stage in sec31a-1/+sec31b-3/+ plants. Our study highlights the essential role of SEC31A/B in gametogenesis and their interchangeable functions in pollen development.

The apple gene responsible for columnar tree shape reduces the abundance of biologically active gibberellin.

Ectopic expression of the apple 2-oxoglutarate-dependent dioxygenase (DOX, 2ODD) gene, designated MdDOX-Co, is thought to cause the columnar shape of apple trees. However, the mechanism underlying the formation of such a unique tree shape remains unclear. To solve this problem, we demonstrated that Arabidopsis thaliana overexpressing MdDOX-Co contained reduced levels of biologically active gibberellin (GA) compared with wild type. In summary: (i) with biochemical approaches, the gene product MdDOX-Co was shown to metabolize active GA A4 (GA4 ) to GA58 (12-OH-GA4 ) in vitro. MdDOX-Co also metabolized its precursors GA12 and GA9 to GA111 (12-OH-GA12 ) and GA70 (12-OH-GA9 ), respectively; (ii) Of the three 12-OH-GAs, GA58 was still active physiologically, but not GA70 or GA111 ; (iii) Arabidopsis MdDOX-Co OE transformants converted exogenously applied deuterium-labeled (d2 )-GA12 to d2 -GA111 but not to d2 -GA58 , whereas transformants converted applied d2 -GA9 to d2 -GA58 ; (iv) GA111 is converted poorly to GA70 by GA 20-oxidases in vitro when GA12 is efficiently metabolized to GA9 ; (v) no GA58 was detected endogenously in MdDOX-Co OE transformants. Overall, we conclude that 12-hydroxylation of GA12 by MdDOX-Co prevents the biosynthesis of biologically active GAs in planta, resulting in columnar phenotypes.

Ethylene-activated PpERF105 induces the expression of the repressor-type R2R3-MYB gene PpMYB140 to inhibit anthocyanin biosynthesis in red pear fruit.

Ethylene induces anthocyanin biosynthesis in most fruits, including apple (Malus domestica), strawberry (Fragaria × ananassa), and plum (Prunus spp.). However, ethylene inhibits anthocyanin biosynthesis in pear (Pyrus spp.), but the underlying molecular mechanism has not been characterized. In this study, ethylene induced the expression of PpERF105, which encodes a transcription factor. PpERF105 functioned as a transcriptional activator, but it inhibited anthocyanin biosynthesis in pear. A transcriptome analysis revealed that PpERF105 activated the expression of PpMYB140, which encodes an R2R3-MYB transcriptional repressor. Moreover, PpMYB140 directly inhibited the expression of anthocyanin-related structural genes. It also competed with PpMYB114 for the binding to bHLH3, ultimately resulting in the formation of the MYB140-bHLH-WDR complex rather than the conventional MBW complex, thereby further inhibiting anthocyanin biosynthesis. Furthermore, PpMYB140 prevented the overaccumulation of anthocyanins in the absence of ethylene. Collectively, our study data indicate that ethylene-induced PpERF105 inhibits anthocyanin biosynthesis by upregulating PpMYB140 expression. Our findings may be useful for elucidating the molecular basis of the ethylene-mediated inhibition of anthocyanin biosynthesis in fruit.

APETALA2 control of barley internode elongation.

Many plants dramatically elongate their stems during flowering, yet how this response is coordinated with the reproductive phase is unclear. We demonstrate that microRNA (miRNA) control of APETALA2 (AP2) is required for rapid, complete elongation of stem internodes in barley, especially of the final 'peduncle' internode directly underneath the inflorescence. Disrupted miR172 targeting of AP2 in the Zeo1.b barley mutant caused lower mitotic activity, delayed growth dynamics and premature lignification in the peduncle leading to fewer and shorter cells. Stage- and tissue-specific comparative transcriptomics between Zeo1.b and its parent cultivar showed reduced expression of proliferation-associated genes, ectopic expression of maturation-related genes and persistent, elevated expression of genes associated with jasmonate and stress responses. We further show that applying methyl jasmonate (MeJA) phenocopied the stem elongation of Zeo1.b, and that Zeo1.b itself was hypersensitive to inhibition by MeJA but less responsive to promotion by gibberellin. Taken together, we propose that miR172-mediated restriction of AP2 may modulate the jasmonate pathway to facilitate gibberellin-promoted stem growth during flowering.

Wheat VRN1, FUL2 and FUL3 play critical and redundant roles in spikelet development and spike determinacy.

The spikelet is the basic unit of the grass inflorescence. In this study, we show that wheat MADS-box genes VRN1, FUL2 and FUL3 play critical and redundant roles in spikelet and spike development, and also affect flowering time and plant height. In the vrn1ful2ful3-null triple mutant, the inflorescence meristem formed a normal double-ridge structure, but then the lateral meristems generated vegetative tillers subtended by leaves instead of spikelets. These results suggest an essential role of these three genes in the fate of the upper spikelet ridge and the suppression of the lower leaf ridge. Inflorescence meristems of vrn1ful2ful3-null and vrn1ful2-null remained indeterminate and single vrn1-null and ful2-null mutants showed delayed formation of the terminal spikelet and increased number of spikelets per spike. Moreover, the ful2-null mutant showed more florets per spikelet, which together with a higher number of spikelets, resulted in a significant increase in the number of grains per spike in the field. Our results suggest that a better understanding of the mechanisms underlying wheat spikelet and spike development can inform future strategies to improve grain yield in wheat.

Molecular identification and functional verification of SPL9 and SPL15 of Lilium.

The transformation of plants from juveniles to adults is a key process in plant growth and development, and the main regulatory factors are miR156 and SQUAMOSA promoter binding protein-like (SPL) transcription factors. Lilium is an ornamental bulb, but it has a long maturation time. In this experiment, Lilium bulbs were subjected to a temperature treatment of 15 °C for 4 weeks to initiate vegetative phase change. Transmission electron microscopy indicated the cell wall of bud core tissue undergoing vegetative phase change became thinner, the starch grains were reduced, and the growth of the juvenile stage was accelerated. The key transcription factors LbrSPL9 and LbrSPL15 were cloned, and the phylogenetic analysis showed they possessed high homology with other plant SPLs. Subcellular localization and transcription activation experiments confirmed LbrSPL9 and LbrSPL15 were mainly located in the nucleus and exhibited transcriptional activity. The results of in situ hybridization showed the expression levels of LbrSPL9 and LbrSPL15 were increased after temperature change treatment. The functional verification experiment of the transgenic plants confirmed that the overexpression of LbrSPL9 and LbrSPL15 could shorten maturation time. These findings help elucidate the regulatory mechanisms of phase transition in Lilium and provide a reference for breeding research in other bulbous flowers.

Genome of Crucihimalaya himalaica, a close relative of Arabidopsis, shows ecological adaptation to high altitude.

Crucihimalaya himalaica, a close relative of Arabidopsis and Capsella, grows on the Qinghai-Tibet Plateau (QTP) about 4,000 m above sea level and represents an attractive model system for studying speciation and ecological adaptation in extreme environments. We assembled a draft genome sequence of 234.72 Mb encoding 27,019 genes and investigated its origin and adaptive evolutionary mechanisms. Phylogenomic analyses based on 4,586 single-copy genes revealed that C. himalaica is most closely related to Capsella (estimated divergence 8.8 to 12.2 Mya), whereas both species form a sister clade to Arabidopsis thaliana and Arabidopsis lyrata, from which they diverged between 12.7 and 17.2 Mya. LTR retrotransposons in C. himalaica proliferated shortly after the dramatic uplift and climatic change of the Himalayas from the Late Pliocene to Pleistocene. Compared with closely related species, C. himalaica showed significant contraction and pseudogenization in gene families associated with disease resistance and also significant expansion in gene families associated with ubiquitin-mediated proteolysis and DNA repair. We identified hundreds of genes involved in DNA repair, ubiquitin-mediated proteolysis, and reproductive processes with signs of positive selection. Gene families showing dramatic changes in size and genes showing signs of positive selection are likely candidates for C. himalaica's adaptation to intense radiation, low temperature, and pathogen-depauperate environments in the QTP. Loss of function at the S-locus, the reason for the transition to self-fertilization of C. himalaica, might have enabled its QTP occupation. Overall, the genome sequence of C. himalaica provides insights into the mechanisms of plant adaptation to extreme environments.

Heat stress response in the closest algal relatives of land plants reveals conserved stress signaling circuits.

All land plants (embryophytes) share a common ancestor that likely evolved from a filamentous freshwater alga. Elucidating the transition from algae to embryophytes - and the eventual conquering of Earth's surface - is one of the most fundamental questions in plant evolutionary biology. Here, we investigated one of the organismal properties that might have enabled this transition: resistance to drastic temperature shifts. We explored the effect of heat stress in Mougeotia and Spirogyra, two representatives of Zygnematophyceae - the closest known algal sister lineage to land plants. Heat stress induced pronounced phenotypic alterations in their plastids, and high-performance liquid chromatography-tandem mass spectroscopy-based profiling of 565 transitions for the analysis of main central metabolites revealed significant shifts in 43 compounds. We also analyzed the global differential gene expression responses triggered by heat, generating 92.8 Gbp of sequence data and assembling a combined set of 8905 well-expressed genes. Each organism had its own distinct gene expression profile; less than one-half of their shared genes showed concordant gene expression trends. We nevertheless detected common signature responses to heat such as elevated transcript levels for molecular chaperones, thylakoid components, and - corroborating our metabolomic data - amino acid metabolism. We also uncovered the heat-stress responsiveness of genes for phosphorelay-based signal transduction that links environmental cues, calcium signatures and plastid biology. Our data allow us to infer the molecular heat stress response that the earliest land plants might have used when facing the rapidly shifting temperature conditions of the terrestrial habitat.

The Ficus erecta genome aids Ceratocystis canker resistance breeding in common fig (F. carica).

Ficus erecta, a wild relative of the common fig (F. carica), is a donor of Ceratocystis canker resistance in fig breeding programmes. Interspecific hybridization followed by recurrent backcrossing is an effective method to transfer the resistance trait from wild to cultivated fig. However, this process is time consuming and labour intensive for trees, especially for gynodioecious plants such as fig. In this study, genome resources were developed for F. erecta to facilitate fig breeding programmes. The genome sequence of F. erecta was determined using single-molecule real-time sequencing technology. The resultant assembly spanned 331.6 Mb with 538 contigs and an N50 length of 1.9 Mb, from which 51 806 high-confidence genes were predicted. Pseudomolecule sequences corresponding to the chromosomes of F. erecta were established with a genetic map based on single nucleotide polymorphisms from double-digest restriction-site-associated DNA sequencing. Subsequent linkage analysis and whole-genome resequencing identified a candidate gene for the Ceratocystis canker resistance trait. Genome-wide genotyping analysis enabled the selection of female lines that possessed resistance and effective elimination of the donor genome from the progeny. The genome resources provided in this study will accelerate and enhance disease-resistance breeding programmes in fig.

Combining GAL4 GFP enhancer trap with split luciferase to measure spatiotemporal promoter activity in Arabidopsis.

In multicellular organisms different types of tissues have distinct gene expression profiles associated with specific function or structure of the cell. Quantification of gene expression in whole organs or whole organisms can give misleading information about levels or dynamics of expression in specific cell types. Tissue- or cell-specific analysis of gene expression has potential to enhance our understanding of gene regulation and interactions of cell signalling networks. The Arabidopsis circadian oscillator is a gene network which orchestrates rhythmic expression across the day/night cycle. There is heterogeneity between cell and tissue types of the composition and behaviour of the oscillator. In order to better understand the spatial and temporal patterns of gene expression, flexible tools are required. By combining a Gateway®-compatible split luciferase construct with a GAL4 GFP enhancer trap system, we describe a tissue-specific split luciferase assay for non-invasive detection of spatiotemporal gene expression in Arabidopsis. We demonstrate the utility of this enhancer trap-compatible split luciferase assay (ETSLA) system to investigate tissue-specific dynamics of circadian gene expression. We confirm spatial heterogeneity of circadian gene expression in Arabidopsis leaves and describe the resources available to investigate any gene of interest.

A zinc finger protein BBX19 interacts with ABF3 to affect drought tolerance negatively in chrysanthemum.

Drought is an environmental factor that can severely influence plant development and distribution, and greatly affect the yield and economic value of crops. We characterized CmBBX19, a BBX family subgroup IV member gene, from the transcriptome database of Chrysanthemum morifolium in response to drought stress. Drought stress and ABA treatments downregulated the expression of CmBBX19. We generated CmBBX19-overexpressing (CmBBX19-OX) lines and CmBBX19-suppressing lines (CmBBX19-RNAi), and found that suppressed expression of CmBBX19 led to enhanced drought tolerance compared with the wild-type (WT) controls, while CmBBX19-OX lines exhibited reduced drought tolerance. Downstream gene analysis showed that CmBBX19 modulates drought tolerance mainly through inducing changes in the expression of ABA-dependent pathway genes, including protective protein, redox balance and cell wall biogenesis genes, such as responsive to ABA 18, peroxidase 12, and cellulose synthase-like protein G2. Moreover, CmBBX19 was shown to interact with CmABF3, a master ABA signaling component, to suppress expression of these downstream genes. We conclude that BBX19-ABF3 module functions in the regulation of drought tolerance of chrysanthemum through an ABA-dependent pathway.

Genetic variants and underlying mechanisms influencing variance heterogeneity in maize.

Traditional genetic studies focus on identifying genetic variants associated with the mean difference in a quantitative trait. Because genetic variants also influence phenotypic variation via heterogeneity, we conducted a variance-heterogeneity genome-wide association study to examine the contribution of variance heterogeneity to oil-related quantitative traits. We identified 79 unique variance-controlling single nucleotide polymorphisms (vSNPs) from the sequences of 77 candidate variance-heterogeneity genes for 21 oil-related traits using the Levene test (P < 1.0 × 10-5 ). About 30% of the candidate genes encode enzymes that work in lipid metabolic pathways, most of which define clear expression variance quantitative trait loci. Of the vSNPs specifically associated with the genetic variance heterogeneity of oil concentration, 89% can be explained by additional linked mean-effects genetic variants. Furthermore, we demonstrated that gene × gene interactions play important roles in the formation of variance heterogeneity for fatty acid compositional traits. The interaction pattern was validated for one gene pair (GRMZM2G035341 and GRMZM2G152328) using yeast two-hybrid and bimolecular fluorescent complementation analyses. Our findings have implications for uncovering the genetic basis of hidden additive genetic effects and epistatic interaction effects, and we indicate opportunities to stabilize efficient breeding and selection of high-oil maize (Zea mays L.).

A mutation in the essential and widely conserved DAMAGED DNA BINDING1-Cullin4 ASSOCIATED FACTOR gene OZS3 causes hypersensitivity to zinc excess, cold and UV stress in Arabidopsis thaliana.

The overly zinc sensitive Arabidopsis thaliana mutant ozs3 shows reduced growth of the primary root, which is exacerbated by an excess specifically of Zn ions. In addition, ozs3 plants display various subtle developmental phenotypes, such as longer petioles and early flowering. Also, ozs3 seedlings are completely but reversibly growth-arrested when shifted to 4°C. The causal mutation was mapped to a gene encoding a putative substrate-recognition receptor of cullin4 E3 ligases. OZS3 orthologous genes can be found in almost all eukaryotic genomes. Most species from Schizosaccharomyces pombe to Homo sapiens, and including A. thaliana, possess one ortholog. No functional data are available for these genes in any of the multicellular model systems. CRISPR-Cas9-mediated knockout demonstrated that a complete loss of OZS3 function is embryo-lethal, indicating essentiality of OZS3 and its orthologs. The OZS3 protein interacts with the adaptor protein DAMAGED DNA BINDING1 (DDB1) in the nucleus. Thus, it is indeed a member of the large yet poorly characterized family of DDB1-cullin4 associated factors in plants. Mutant phenotypes of ozs3 plants are apparently caused by the weakened DDB1-OZS3 interaction as a result of the exchange of a conserved amino acid near the conserved WDxR motif.

A quantitative trait locus GW6 controls rice grain size and yield through the gibberellin pathway.

Grain size is one of the essential components determining rice yield and is a target for both domestication and artificial breeding. Gibberellins (GAs) are diterpenoid phytohormones that influence diverse aspects of plant growth and development. Several quantitative trait loci (QTLs) have been identified that control grain size through phytohormone regulation. However, little is known about the role of GAs in the control of grain size. Here we report the cloning and characterization of a QTL, GW6 (GRAIN WIDTH 6), which encodes a GA-regulated GAST family protein and positively regulates grain width and weight. GW6 is highly expressed in the young panicle and increases grain width by promoting cell expansion in the spikelet hull. Knockout of GW6 exhibits reduced grain size and weight, whereas overexpression of GW6 results in increased grain size and weight. GW6 is induced by GA and its knockout downregulates the expression of GA biosynthesis genes and decreases GA content in the young panicle. We found that a natural variation in the cis element CAAT-box in the promoter of GW6 is associated with its expression level and grain width and weight. Furthermore, introduction of GW6 to Oryza indica variety HJX74 can lead to a 10.44% increase in rice grain yield, indicating that GW6 has great potential to improve grain yield in rice.

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.