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.

The H3K4 demethylase JMJ1 is required for proper timing of flowering in Brachypodium distachyon.

Flowering is a key developmental transition in the plant life cycle. In temperate climates, flowering often occurs in response to the perception of seasonal cues such as changes in day-length and temperature. However, the mechanisms that have evolved to control the timing of flowering in temperate grasses are not fully understood. We identified a Brachypodium distachyon mutant whose flowering is delayed under inductive long-day conditions due to a mutation in the JMJ1 gene, which encodes a Jumonji domain-containing protein. JMJ1 is a histone demethylase that mainly demethylates H3K4me2 and H3K4me3 in vitro and in vivo. Analysis of the genome-wide distribution of H3K4me1, H3K4me2, and H3K4me3 in wild-type plants by chromatin immunoprecipitation and sequencing combined with RNA sequencing revealed that H3K4m1 and H3K4me3 are positively associated with gene transcript levels, whereas H3K4me2 is negatively correlated with transcript levels. Furthermore, JMJ1 directly binds to the chromatin of the flowering regulator genes VRN1 and ID1 and affects their transcription by modifying their H3K4me2 and H3K4me3 levels. Genetic analyses indicated that JMJ1 promotes flowering by activating VRN1 expression. Our study reveals a role for JMJ1-mediated chromatin modification in the proper timing of flowering in B. distachyon.

Transposition of HOPPLA in siRNA-deficient plants suggests a limited effect of the environment on retrotransposon mobility in Brachypodium distachyon.

Long terminal repeat retrotransposons (LTR-RTs) are powerful mutagens regarded as a major source of genetic novelty and important drivers of evolution. Yet, the uncontrolled and potentially selfish proliferation of LTR-RTs can lead to deleterious mutations and genome instability, with large fitness costs for their host. While population genomics data suggest that an ongoing LTR-RT mobility is common in many species, the understanding of their dual role in evolution is limited. Here, we harness the genetic diversity of 320 sequenced natural accessions of the Mediterranean grass Brachypodium distachyon to characterize how genetic and environmental factors influence plant LTR-RT dynamics in the wild. When combining a coverage-based approach to estimate global LTR-RT copy number variations with mobilome-sequencing of nine accessions exposed to eight different stresses, we find little evidence for a major role of environmental factors in LTR-RT accumulations in B. distachyon natural accessions. Instead, we show that loss of RNA polymerase IV (Pol IV), which mediates RNA-directed DNA methylation in plants, results in high transcriptional and transpositional activities of RLC_BdisC024 (HOPPLA) LTR-RT family elements, and that these effects are not stress-specific. This work supports findings indicating an ongoing mobility in B. distachyon and reveals that host RNA-directed DNA methylation rather than environmental factors controls their mobility in this wild grass model.

Expanded roles and divergent regulation of FAMA in Brachypodium and Arabidopsis stomatal development.

Stomata, cellular valves found on the surfaces of aerial plant tissues, present a paradigm for studying cell fate and patterning in plants. A highly conserved core set of related basic helix-loop-helix (bHLH) transcription factors regulates stomatal development across diverse species. We characterized BdFAMA in the temperate grass Brachypodium distachyon and found this late-acting transcription factor was necessary and sufficient for specifying stomatal guard cell fate, and unexpectedly, could also induce the recruitment of subsidiary cells in the absence of its paralogue, BdMUTE. The overlap in function is paralleled by an overlap in expression pattern and by unique regulatory relationships between BdMUTE and BdFAMA. To better appreciate the relationships among the Brachypodium stomatal bHLHs, we used in vivo proteomics in developing leaves and found evidence for multiple shared interaction partners. We reexamined the roles of these genes in Arabidopsis thaliana by testing genetic sufficiency within and across species, and found that while BdFAMA and AtFAMA can rescue stomatal production in Arabidopsis fama and mute mutants, only AtFAMA can specify Brassica-specific myrosin idioblasts. Taken together, our findings refine the current models of stomatal bHLH function and regulatory feedback among paralogues within grasses as well as across the monocot/dicot divide.

Scattered differentiation of unlinked loci across the genome underlines ecological divergence of the selfing grass Brachypodium stacei.

Ecological divergence without geographic isolation, as an early speciation process that may lead finally to reproductive isolation through natural selection, remains a captivating topic in evolutionary biology. However, the pattern of genetic divergence underlying this process across the genome may vary between species and mating systems. Here, we present evidence that Brachypodium stacei, an annual and highly selfing grass model species, has undergone sympatric ecological divergence without geographic isolation. Genomic, transcriptomic, and metabolomic analyses together with lab experiments mimicking the two opposite environmental conditions suggest that diploid B. stacei populations have diverged sympatrically in two slopes characterized by distinct biomes at Evolution Canyon I (ECI), Mount Carmel, Israel. Despite ongoing gene flow, primarily facilitated by seed dispersal, the level of gene flow has progressively decreased over time. This local adaptation involves the scattered divergence of many unlinked loci across the total genome that include both coding genes and noncoding regions. Additionally, we have identified significant differential expressions of genes related to the ABA signaling pathway and contrasting metabolome composition between the arid- vs. forest-adapted B. stacei populations in ECI. These results suggest that multiple small loci involved in environmental responses act additively to account for ecological adaptations by this selfing species in contrasting environments.

INDETERMINATE1-mediated expression of FT family genes is required for proper timing of flowering in Brachypodium distachyon.

The transition to flowering is a major developmental switch in plants. In many temperate grasses, perception of indicators of seasonal change, such as changing day-length and temperature, leads to expression of FLOWERING LOCUS T1 (FT1) and FT-Like (FTL) genes that are essential for promoting the transition to flowering. However, little is known about the upstream regulators of FT1 and FTL genes in temperate grasses. Here, we characterize the monocot-specific gene INDETERMINATE1 (BdID1) in Brachypodium distachyon and demonstrate that BdID1 is a regulator of FT family genes. Mutations in ID1 impact the ability of the short-day (SD) vernalization, cold vernalization, and long-day (LD) photoperiod pathways to induce certain FTL genes. BdID1 is required for upregulation of FTL9 (FT-LIKE9) expression by the SD vernalization pathway, and overexpression of FTL9 in an id1 background can partially restore the delayed flowering phenotype of id1. We show that BdID1 binds in vitro to the promoter region of FTL genes suggesting that ID1 directly activates FTL expression. Transcriptome analysis shows that BdID1 is required for FT1, FT2, FTL12, and FTL13 expression under inductive LD photoperiods, indicating that BdID1 is a regulator of the FT gene family. Moreover, overexpression of FT1 in the id1 background results in rapid flowering similar to overexpressing FT1 in the wild type, demonstrating that BdID1 is upstream of FT family genes. Interestingly, ID1 negatively regulates a previously uncharacterized FTL gene, FTL4, and we show that FTL4 is a repressor of flowering. Thus, BdID1 is critical for proper timing of flowering in temperate grasses.

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.

Group VII ethylene response factors forming distinct regulatory loops mediate submergence responses.

Group VII ethylene response factors (ERFVIIs), whose stability is oxygen concentration-dependent, play key roles in regulating hypoxia response genes in Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa) during submergence. To understand the evolution of flooding tolerance in cereal crops, we evaluated whether Brachypodium distachyon ERFVII genes (BdERFVIIs) are related to submergence tolerance. We found that three BdERFVIIs, BdERF108, BdERF018, and BdERF961, form a feedback regulatory loop to mediate downstream responses. BdERF108 and BdERF018 activated the expression of BdERF961 and PHYTOGLOBIN 1 (PGB1), which promoted nitric oxide turnover and preserved ERFVII protein stability. The activation of PGB1 was subsequently counteracted by increased BdERF961 accumulation through negative feedback regulation. Interestingly, we found that OsERF67, the orthologue of BdERF961 in rice, activated PHYTOGLOBIN (OsHB2) expression and formed distinct regulatory loops during submergence. Overall, the divergent regulatory mechanisms exhibited by orthologs collectively offer perspectives for the development of submergence-tolerant crops.

CO2 response screen in grass Brachypodium reveals the key role of a MAP kinase in CO2-triggered stomatal closure.

Plants respond to increased CO2 concentrations through stomatal closure, which can contribute to increased water use efficiency. Grasses display faster stomatal responses than eudicots due to dumbbell-shaped guard cells flanked by subsidiary cells working in opposition. However, forward genetic screening for stomatal CO2 signal transduction mutants in grasses has yet to be reported. The grass model Brachypodium distachyon is closely related to agronomically important cereal crops, sharing largely collinear genomes. To gain insights into CO2 control mechanisms of stomatal movements in grasses, we developed an unbiased forward genetic screen with an EMS-mutagenized B. distachyon M5 generation population using infrared imaging to identify plants with altered leaf temperatures at elevated CO2. Among isolated mutants, a "chill1" mutant exhibited cooler leaf temperatures than wild-type Bd21-3 parent control plants after exposure to increased CO2. chill1 plants showed strongly impaired high CO2-induced stomatal closure despite retaining a robust abscisic acid-induced stomatal closing response. Through bulked segregant whole-genome sequencing analyses followed by analyses of further backcrossed F4 generation plants and generation and characterization of sodium azide and CRISPR-cas9 mutants, chill1 was mapped to a protein kinase, Mitogen-Activated Protein Kinase 5 (BdMPK5). The chill1 mutation impaired BdMPK5 protein-mediated CO2/HCO3- sensing together with the High Temperature 1 (HT1) Raf-like kinase in vitro. Furthermore, AlphaFold2-directed structural modeling predicted that the identified BdMPK5-D90N chill1 mutant residue is located at the interface of BdMPK5 with the BdHT1 Raf-like kinase. BdMPK5 is a key signaling component that mediates CO2-induced stomatal movements and is proposed to function as a component of the primary CO2 sensor in grasses.

Proteomic and metabolic disturbances in lignin-modified Brachypodium distachyon.

Lignin biosynthesis begins with the deamination of phenylalanine and tyrosine (Tyr) as a key branch point between primary and secondary metabolism in land plants. Here, we used a systems biology approach to investigate the global metabolic responses to lignin pathway perturbations in the model grass Brachypodium distachyon. We identified the lignin biosynthetic protein families and found that ammonia-lyases (ALs) are among the most abundant proteins in lignifying tissues in grasses. Integrated metabolomic and proteomic data support a link between lignin biosynthesis and primary metabolism mediated by the ammonia released from ALs that is recycled for the synthesis of amino acids via glutamine. RNA interference knockdown of lignin genes confirmed that the route of the canonical pathway using shikimate ester intermediates is not essential for lignin formation in Brachypodium, and there is an alternative pathway from Tyr via sinapic acid for the synthesis of syringyl lignin involving yet uncharacterized enzymatic steps. Our findings support a model in which plant ALs play a central role in coordinating the allocation of carbon for lignin synthesis and the nitrogen available for plant growth. Collectively, these data also emphasize the value of integrative multiomic analyses to advance our understanding of plant metabolism.

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.

BdRCN4, a Brachypodium distachyon TFL1 homologue, is involved in regulation of apical meristem fate.

In higher plants, the shift from vegetative to reproductive development is governed by complex interplay of internal and external signals. TERMINALFLOWER1 (TFL1) plays a crucial role in the regulation of flowering time and inflorescence architecture in Arabidopsis thaliana. This study aimed to explore the function of BdRCN4, a homolog of TFL1 in Brachypodium distachyon, through functional analyses in mutant and transgenic plants. The results revealed that overexpression of BdRCN4 in B. distachyon leads to an extended vegetative phase and reduced production of spikelets. Similar results were found in A. thaliana, where constitutive expression of BdRCN4 promoted a delay in flowering time, followed by the development of hypervegetative shoots, with no flowers or siliques produced. Our results suggest that BdRCN4 acts as a flowering repressor analogous to TFL1, negatively regulating AP1, but no LFY expression. To further validate this hypothesis, a 35S::LFY-GR co-transformation approach on 35::BdRCN4 lines was performed. Remarkably, AP1 expression levels and flower formation were restored to normal in co-transformed plants when treated with dexamethasone. Although further molecular studies will be necessary, the evidence in B. distachyon support the idea that a balance between LFY and BdRCN4/TFL1 seems to be essential for activating AP1 expression and initiating floral organ identity gene expression. This study also demonstrates interesting conservation through the molecular pathways that regulate flowering meristem transition and identity across the evolution of monocot and dicot plants.

Subgenomic Stability of Progenitor Genomes During Repeated Allotetraploid Origins of the Same Grass Brachypodium hybridum.

Both homeologous exchanges and homeologous expression bias are generally found in most allopolyploid species. Whether homeologous exchanges and homeologous expression bias differ between repeated allopolyploid speciation events from the same progenitor species remains unknown. Here, we detected a third independent and recent allotetraploid origin for the model grass Brachypodium hybridum. Our homeologous exchange with replacement analyses indicated the absence of significant homeologous exchanges in any of the three types of wild allotetraploids, supporting the integrity of their progenitor subgenomes and the immediate creation of the amphidiploids. Further homeologous expression bias tests did not uncover significant subgenomic dominance in different tissues and conditions of the allotetraploids. This suggests a balanced expression of homeologs under similar or dissimilar ecological conditions in their natural habitats. We observed that the density of transposons around genes was not associated with the initial establishment of subgenome dominance; rather, this feature is inherited from the progenitor genome. We found that drought response genes were highly induced in the two subgenomes, likely contributing to the local adaptation of this species to arid habitats in the third allotetraploid event. These findings provide evidence for the consistency of subgenomic stability of parental genomes across multiple allopolyploidization events that led to the same species at different periods. Our study emphasizes the importance of selecting closely related progenitor species genomes to accurately assess homeologous exchange with replacement in allopolyploids, thereby avoiding the detection of false homeologous exchanges when using less related progenitor species genomes.

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.

Does time matter? Intraspecific diversity of ribosomal RNA genes in lineages of the allopolyploid model grass Brachypodium hybridum with different evolutionary ages.

BACKGROUND: Polyploidisation often results in genome rearrangements that may involve changes in both the single-copy sequences and the repetitive genome fraction. In this study, we performed a comprehensive comparative analysis of repetitive DNA, with a particular focus on ribosomal DNA (rDNA), in Brachypodium hybridum (2n = 4x = 30, subgenome composition DDSS), an allotetraploid resulting from a natural cross between two diploid species that resemble the modern B. distachyon (2n = 10; DD) and B. stacei (2n = 20; SS). Taking advantage of the recurrent origin of B. hybridum, we investigated two genotypes, Bhyb26 and ABR113, differing markedly in their evolutionary age (1.4 and 0.14 Mya, respectively) and which resulted from opposite cross directions. To identify the origin of rDNA loci we employed cytogenetic and molecular methods (FISH, gCAPS and Southern hybridisation), phylogenetic and genomic approaches. RESULTS: Unlike the general maintenance of doubled gene dosage in B. hybridum, the rRNA genes showed a remarkable tendency towards diploidisation at both locus and unit levels. While the partial elimination of 35S rDNA units occurred in the younger ABR113 lineage, unidirectional elimination of the entire locus was observed in the older Bhyb26 lineage. Additionally, a novel 5S rDNA family was amplified in Bhyb26 replacing the parental units. The 35S and 5S rDNA units were preferentially eliminated from the S- and D-subgenome, respectively. Thus, in the more ancient B. hybridum lineage, Bhyb26, 5S and 35S rRNA genes are likely expressed from different subgenomes, highlighting the complexity of polyploid regulatory networks. CONCLUSION: Comparative analyses between two B. hybridum lineages of distinct evolutionary ages revealed that although the recent lineage ABR113 exhibited an additive pattern of rDNA loci distribution, the ancient lineage Bhyb26 demonstrated a pronounced tendency toward diploidisation manifested by the reduction in the number of both 35S and 5S loci. In conclusion, the age of the allopolyploid appears to be a decisive factor in rDNA turnover in B. hybridum.

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.

The auxin efflux carrier PIN1a regulates vascular patterning in cereal roots.

Barley (Hordeum vulgare) is an important global cereal crop and a model in genetic studies. Despite advances in characterising barley genomic resources, few mutant studies have identified genes controlling root architecture and anatomy, which plays a critical role in capturing soil resources. Our phenotypic screening of a TILLING mutant collection identified line TM5992 exhibiting a short-root phenotype compared with wild-type (WT) Morex background. Outcrossing TM5992 with barley variety Proctor and subsequent SNP array-based bulk segregant analysis, fine mapped the mutation to a cM scale. Exome sequencing pinpointed a mutation in the candidate gene HvPIN1a, further confirming this by analysing independent mutant alleles. Detailed analysis of root growth and anatomy in Hvpin1a mutant alleles exhibited a slower growth rate, shorter apical meristem and striking vascular patterning defects compared to WT. Expression and mutant analyses of PIN1 members in the closely related cereal brachypodium (Brachypodium distachyon) revealed that BdPIN1a and BdPIN1b were redundantly expressed in root vascular tissues but only Bdpin1a mutant allele displayed root vascular defects similar to Hvpin1a. We conclude that barley PIN1 genes have sub-functionalised in cereals, compared to Arabidopsis (Arabidopsis thaliana), where PIN1a sequences control root vascular patterning.

Repeated migration, interbreeding and bottlenecking shaped the phylogeography of the selfing grass Brachypodium stacei.

Brachypodium stacei is the most ancestral lineage in the genus Brachypodium, a model system for grass functional genomics. B. stacei shows striking and sometimes contradictory biological and evolutionary features, including a high selfing rate yet extensive admixture, an ancient Miocene origin yet with recent evolutionary radiation, and adaptation to different dry climate conditions in its narrow distribution range. Therefore, it constitutes an ideal system to study these life history traits. We studied the phylogeography of 17 native circum-Mediterranean B. stacei populations (39 individuals) using genome-wide RADseq SNP data and complete plastome sequences. Nuclear SNP data revealed the existence of six distinct genetic clusters, low levels of intra-population genetic diversity and high selfing rates, albeit with signatures of admixture. Coalescence-based dating analysis detected a recent split between crown lineages in the Late Quaternary. Plastome sequences showed incongruent evolutionary relationships with those recovered by the nuclear data, suggesting interbreeding and chloroplast capture events between genetically distant populations. Demographic and population dispersal coalescent models identified an ancestral origin of B. stacei in the western-central Mediterranean islands, followed by an early colonization of the Canary Islands and two independent colonization events of the eastern Mediterranean region through long-distance dispersal and bottleneck events as the most likely evolutionary history. Climate niche data identified three arid niches of B. stacei in the southern Mediterranean region. Our findings indicate that the phylogeography of B. stacei populations was shaped by recent radiations, frequent extinctions, long-distance dispersal events, occasional interbreeding, and adaptation to local climates.