Brachypodium: A Monocot Grass Model Genus for Plant Biology

Investigation of wheat cold response pathway regulated by TaICE41 and TaCBFⅣd-B9 through Brachypodium distachyon transformation

Wheat (Triticum aestivum L.), a major global crop, is vulnerable to freezing stress, particularly during late spring frosts. Enhancing freezing tolerance through cold acclimation, primarily via the ICE-CBF-COR pathway, is crucial for improving wheat productivity. This study focuses on identifying genes regulated by the ICE-CBF pathway and those that function independently in response to freezing stress. TaICE41 and TaCBFⅣd-B9, two key genes associated with cold tolerance, were cloned and analyzed for their phylogenetic characteristics and subcellular localization. Transgenic Brachypodium distachyon overexpressing these genes demonstrated enhanced freezing tolerance, with increased survival rates and proline content, compared to wild-type plants. RNA-seq analysis revealed distinct gene expression profiles under cold stress, highlighting both shared and unique pathways regulated by ICE41 and CBF. Notably, the TaICE41-overexpressing lines exhibited upregulation of genes involved in phenylpropanoid biosynthesis and starch-sucrose metabolism, contributing to stress response. This study provides new insights into the ICE-CBF pathway and its role in cold tolerance, emphasizing the importance of understanding both ICE-CBF-regulated and independent cold-responsive genes for improving freezing tolerance in crops.

Natural variation in the adjustment of primary metabolism determines ammonium tolerance in the model grass Brachypodium distachyon

Nitrogen (N) fertilization is essential to maximize crop production. However, around half of the applied N is lost to the environment, causing water and air pollution and contributing to climate change. Understanding the natural genetic and metabolic basis underlying plants N use efficiency is of great interest to attain an agriculture with less N demand and thus more sustainable. The study of ammonium (NH4+) nutrition is of particular interest, because it mitigates N losses due to nitrate (NO3-) leaching or denitrification. In this work, we studied Brachypodium distachyon, the model plant for C3 grasses, grown with NH4+ or NO3- supply. We performed gene expression analysis in the root of the B. distachyon reference accession Bd21 and examined the phenotypic variation across 52 natural accessions through analyzing plant growth and a panel of 22 metabolic traits in leaf and root. We found that the adjustment of primary metabolism to NH4+ nutrition is essential for the natural variation of NH4+ tolerance, notably involving NH4+ assimilation and phosphoenolpyruvate carboxylase (PEPC) activity. Additionally, genome-wide association studies (GWAS) indicated several loci associated with B. distachyon growth and metabolic adaptation to NH4+ nutrition. We found that the GDH2 gene was associated with the induction of root glutamate dehydrogenase activity under NH4+ nutrition and that two genes encoding malic enzyme were associated with leaf PEPC activity. Altogether, our work underlines the value of natural variation and the key role of primary metabolism to improve NH4+ tolerance.

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.

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.

An atlas of Brachypodium distachyon lateral root development

The root system of plants is a vital part for successful development and adaptation to different soil types and environments. A major determinant of the shape of a plant root system is the formation of lateral roots, allowing for expansion of the root system. Arabidopsis thaliana, with its simple root anatomy, has been extensively studied to reveal the genetic program underlying root branching. However, to get a more general understanding of lateral root development, comparative studies in species with a more complex root anatomy are required. Here, by combining optimized clearing methods and histology, we describe an atlas of lateral root development in Brachypodium distachyon, a wild, temperate grass species. We show that lateral roots initiate from enlarged phloem pole pericycle cells and that the overlying endodermis reactivates its cell cycle and eventually forms the root cap. In addition, auxin signaling reported by the DR5 reporter was not detected in the phloem pole pericycle cells or young primordia. In contrast, auxin signaling was activated in the overlying cortical cell layers, including the exodermis. Thus, Brachypodium is a valuable model to investigate how signaling pathways and cellular responses have been repurposed to facilitate lateral root organogenesis.

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.

The First Observation of the Filamentous Fungus Neurospora crassa Growing in the Roots of the Grass Brachypodium distachyon

The model organism Neurospora crassa has been cultivated in laboratories since the 1920s and its saprotrophic lifestyle has been established for decades. However, beyond their role as saprotrophs, fungi engage in intricate relationships with plants, showcasing diverse connections ranging from mutualistic to pathogenic. Although N. crassa has been extensively investigated under laboratory conditions, its ecological characteristics remain largely unknown. In contrast, Brachypodium distachyon, a sweet grass closely related to significant crops, demonstrates remarkable ecological flexibility and participates in a variety of fungal interactions, encompassing both mutualistic and harmful associations. Through a comprehensive microscopic analysis using electron, fluorescence, and confocal laser scanning microscopy, we discovered a novel endophytic interaction between N. crassa and B. distachyon roots, where fungal hyphae not only thrive in the apoplastic space and vascular bundle but also may colonize plant root cells. This new and so far hidden trait of one of the most important fungal model organisms greatly enhances our view of N. crassa, opening new perspectives concerning the fungus' ecological role. In addition, we present a new tool for studying plant-fungus interspecies communication, combining two well-established model systems, which improves our possibilities of experimental design on the molecular level.

Expansions and contractions of repetitive DNA elements reveal contrasting evolutionary responses to the polyploid genome shock hypothesis in Brachypodium model grasses

Brachypodium grass species have been selected as model plants for functional genomics of grass crops, and to elucidate the origins of allopolyploidy and perenniality in monocots, due to their small genome sizes and feasibility of cultivation. However, genome sizes differ greatly between diploid or polyploid Brachypodium lineages. We have used genome skimming sequencing data to uncover the composition, abundance, and phylogenetic value of repetitive elements in 44 representatives of the major Brachypodium lineages and cytotypes. We also aimed to test the possible mechanisms and consequences of the "polyploid genome shock hypothesis" (PGSH) under three different evolutionary scenarios of variation in repeats and genome sizes of Brachypodium allopolyploids. Our data indicated that the proportion of the genome covered by the repeatome in the Brachypodium species showed a 3.3-fold difference between the highest content of B. mexicanum-4x (67.97%) and the lowest of B. stacei-2x (20.77%), and that changes in the sizes of their genomes were a consequence of gains or losses in their repeat elements. LTR-Retand and Tekay retrotransposons were the most frequent repeat elements in the Brachypodium genomes, while Ogre retrotransposons were found exclusively in B. mexicanum. The repeatome phylogenetic network showed a high topological congruence with plastome and nuclear rDNA and transcriptome trees, differentiating the ancestral outcore lineages from the recently evolved core-perennial lineages. The 5S rDNA graph topologies had a strong match with the ploidy levels and nature of the subgenomes of the Brachypodium polyploids. The core-perennial B. sylvaticum presents a large repeatome and characteristics of a potential post-polyploid diploidized origin. Our study evidenced that expansions and contractions in the repeatome were responsible for the three contrasting responses to the PGSH. The exacerbated genome expansion of the ancestral allotetraploid B. mexicanum was a consequence of chromosome-wide proliferation of TEs and not of WGD, the additive repeatome pattern of young allotetraploid B. hybridum of stabilized post-WGD genome evolution, and the genomecontraction of recent core-perennials polyploids (B. pinnatum, B. phoenicoides) of repeat losses through recombination of these highly hybridizing lineages. Our analyses have contributed to unraveling the evolution of the repeatome and the genome size variation in model Brachypodium grasses.

Building blocks toward sustainable biofertilizers: variation in arbuscular mycorrhizal spore germination when immobilized with diazotrophic bacteria in biodegradable hydrogel beads

AIM

We investigated whether there was interspecies and intraspecies variation in spore germination of 12 strains of arbuscular mycorrhizal fungi when co-entrapped with the diazotrophic plant growth-promoting bacteria, Azospirillum brasilense Sp7 in alginate hydrogel beads.

METHODS AND RESULTS

Twelve Rhizophagus irregularis, Rhizophagus intraradices, and Funneliformis mosseae strains were separately combined with a live culture of Azospirillum brasilense Sp7. Each fungal-bacterial consortia was supplemented with sodium alginate to a 2% concentration (v/v) and cross-linked in calcium chloride (2% w/v) to form biodegradable hydrogel beads. One hundred beads from each combination (total of 1200) were fixed in solidified modified Strullu and Romand media. Beads were observed for successful spore germination and bacterial growth over 14 days. In all cases, successful growth of A. brasilense was observed. For arbuscular mycorrhizal fungi, interspecies variation in spore germination was observed, with R. intraradices having the highest germination rate (64.3%), followed by R. irregularis (45.5%) and F. mosseae (40.3%). However, a difference in intraspecies germination was only observed among strains of R. irregularis and F. mosseae. Despite having varying levels of germination, even the strains with the lowest potential were still able to establish with the plant host Brachypodium distachyon in a model system.

CONCLUSIONS

Arbuscular mycorrhizal spore germination varied across strains when co-entrapped with a diazotrophic plant growth-promoting bacteria. This demonstrates that hydrogel beads containing a mixed consortium hold potential as a sustainable biofertilizer and that compatibility tests remain an important building block when aiming to create a hydrogel biofertilizer that encases a diversity of bacteria and fungi. Moving forward, further studies should be conducted to test the efficacy of these hydrogel biofertilizers on different crops across varying climatic conditions in order to optimize their potential.

The Brachypodium distachyon DREB transcription factor BdDREB-39 confers oxidative stress tolerance in transgenic tobacco

Key message BdDREB-39 is a DREB/CBF transcription factor, localized in the nucleus with transactivation activity, and BdDREB-39-overexpressing transgenic yeasts and tobacco enhanced the tolerance to oxidative stress.Abstract The DREB/CBF transcription factors are generally recognized to play an important factor in plant growth, development and response to various abiotic stresses. However, the mechanism of DREB/CBFs in oxidative stress response is largely unknown. This study isolated a DREB/CBF gene BdDREB-39 from Brachypodium distachyon (B. distachyon). Multiple sequence alignment and phylogenetic analysis showed that BdDREB-39 was closely related to the DREB proteins of oats, barley, wheat and rye and therefore its study can provide a reference for the excavation and genetic improvement of BdDREB-39 or its homologs in its closely related species. The transcript levels of BdDREB-39 were significantly up-regulated under H2O2 stress. BdDREB-39 was localised in the nucleus and functioned as a transcriptional activator. Overexpression of BdDREB-39 enhanced H2O2 tolerance in yeast. Transgenic tobaccos with BdDREB-39 had higher germination rates, longer root, better growth status, lesser reactive oxygen species (ROS) and malondialdehyde (MDA), and higher superoxide dismutase (SOD) and peroxidase (POD) activities than wild type (WT). The expression levels of ROS-related and stress-related genes were improved by BdDREB-39. In summary, these results revealed that BdDREB-39 can improve the viability of tobacco by regulating the expression of ROS and stress-related genes, allowing transgenic tobacco to accumulate lower levels of ROS and reducing the damage caused by ROS to cells. The BdDREB-39 gene has the potential for developing plant varieties tolerant to stress.

Genome-wide identification of the plant homeodomain-finger family in rye and ScPHD5 functions in cold tolerance and flowering time

KEY MESSAGE

111 PHD genes were newly identified in rye genome and ScPHD5's role in regulating cold tolerance and flowering time was suggested. Plant homeodomain (PHD)-finger proteins regulate the physical properties of chromatin and control plant development and stress tolerance. Although rye (Secale cereale L.) is a major winter crop, PHD-finger proteins in rye have not been studied. Here, we identified 111 PHD genes in the rye genome that exhibited diverse gene and protein sequence structures. Phylogenetic tree analysis revealed that PHDs were genetically close in monocots and diverged from those in dicots. Duplication and synteny analyses demonstrated that ScPHDs have undergone several duplications during evolution and that high synteny is conserved among the Triticeae species. Tissue-specific and abiotic stress-responsive gene expression analyses indicated that ScPHDs were highly expressed in spikelets and developing seeds and were responsive to cold and drought stress. One of these genes, ScPHD5, was selected for further functional characterization. ScPHD5 was highly expressed in the spike tissues and was localized in the nuclei of rye protoplasts and tobacco leaves. ScPHD5-overexpressing Brachypodium was more tolerant to freezing stress than wild-type (WT), with increased CBF and COR gene expression. Additionally, these transgenic plants displayed an extremely early flowering phenotype that flowered more than two weeks earlier than the WT, and vernalization genes, rather than photoperiod genes, were increased in the WT. RNA-seq analysis revealed that diverse stress response genes, including HSPs, HSFs, LEAs, and MADS-box genes, were also upregulated in transgenic plants. Our study will help elucidate the roles of PHD genes in plant development and abiotic stress tolerance in rye.

TaMAPK3 phosphorylates TaCBF and TaICE and plays a negative role in wheat freezing tolerance

Freezing temperature during overwintering often kills plants; plants have thus, developed a defense mechanism called 'cold acclimation', in which a number of genes are involved in increasing cell protection and gene expression. Mitogen-activated protein kinase (MAPK) controls proteins' activities by phosphorylation and is involved in numerous metabolic pathways. In this study, we identified the protein interaction between TaMAPK3 and the proteins in the cold response pathway, ICE41, ICE87, and CBFIVd-D9. The subcellular localization and bimolecular fluorescence complement (BiFC) assays revealed that these proteins interact in the nucleus or in the plasma membrane. Furthermore, MAPK3-mediated phosphorylation of ICE41, ICE87, and CBFIVd-D9 was verified using an in vitro phosphorylation assay. TaMAPK3-overexpressing transgenic Brachypodium showed a lower survival rate upon freezing stress and lower proline content during cold acclimation, compared to wild-type plants. Furthermore, cold response gene expression analysis revealed that the expression of these genes was suppressed in the transgenic lines under cold treatment. It was further elucidated that MAPK3 mediates the degradation of ICE and CBF proteins, which implies the negative impact of MAPK3 on the freezing tolerance of plants. This study will help to elucidate the molecular mechanisms of cold tolerance and the activity of MAPK3 in wheat.

Three near-complete genome assemblies reveal substantial centromere dynamics from diploid to tetraploid in Brachypodium genus

BACKGROUND

Centromeres are critical for maintaining genomic stability in eukaryotes, and their turnover shapes genome architectures and drives karyotype evolution. However, the co-evolution of centromeres from different species in allopolyploids over millions of years remains largely unknown.

RESULTS

Here, we generate three near-complete genome assemblies, a tetraploid Brachypodium hybridum and its two diploid ancestors, Brachypodium distachyon and Brachypodium stacei. We detect high degrees of sequence, structural, and epigenetic variations of centromeres at base-pair resolution between closely related Brachypodium genomes, indicating the appearance and accumulation of species-specific centromere repeats from a common origin during evolution. We also find that centromere homogenization is accompanied by local satellite repeats bursting and retrotransposon purging, and the frequency of retrotransposon invasions drives the degree of interspecies centromere diversification. We further investigate the dynamics of centromeres during alloploidization process, and find that dramatic genetics and epigenetics architecture variations are associated with the turnover of centromeres between homologous chromosomal pairs from diploid to tetraploid. Additionally, our pangenomes analysis reveals the ongoing variations of satellite repeats and stable evolutionary homeostasis within centromeres among individuals of each Brachypodium genome with different polyploidy levels.

CONCLUSIONS

Our results provide unprecedented information on the genomic, epigenomic, and functional diversity of highly repetitive DNA between closely related species and their allopolyploid genomes at both coarse and fine scale.

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.

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.

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.

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.

Plant virology: an RNA treasure trove

Key principles pertaining to RNA biology not infrequently have their origins in plant virology. Examples have arisen from studies on viral RNA-intrinsic properties and the infection process from gene expression, replication, movement, and defense evasion to biotechnological applications. Since RNA is at the core of the central dogma in molecular biology, how plant virology assisted in the reinforcement or adaptations of this concept, while at other instances shook up elements of the doctrine, is discussed. Moreover, despite the negative effects of viral diseases in agriculture worldwide, plant viruses can be considered a scientific treasure trove. Today they remain tools of discovery for biotechnology, studying evolution, cell biology, and host-microbe interactions.

MIF-like domain containing protein orchestrates cellular differentiation and virulence in the fungal pathogen Magnaporthe oryzae

Macrophage migration inhibitory factor (MIF) is a pleiotropic protein with chemotactic, pro-inflammatory, and growth-promoting activities first discovered in mammals. In parasites, MIF homologs are involved in immune evasion and pathogenesis. Here, we present the first comprehensive analysis of an MIF protein from the devastating plant pathogen Magnaporthe oryzae (Mo). The fungal genome encodes a single MIF protein (MoMIF1) that, unlike the human homolog, harbors multiple low-complexity regions (LCRs) and is unique to Ascomycota. Following infection, MoMIF1 is expressed in the biotrophic phase of the fungus, and is strongly down-regulated during subsequent necrotrophic growth in leaves and roots. We show that MoMIF1 is secreted during plant infection, affects the production of the mycotoxin tenuazonic acid and inhibits plant cell death. Our results suggest that MoMIF1 is a novel key regulator of fungal virulence that maintains the balance between biotrophy and necrotrophy during the different phases of fungal infection.

Impact of high atmospheric carbon dioxide on the biotic stress response of the model cereal species Brachypodium distachyon

Losses due to disease and climate change are among the most important issues currently facing crop production. It is therefore important to establish the impact of climate change, and particularly of high carbon dioxide (hCO2), on plant immunity in cereals, which provide 60% of human calories. The aim of this study was to determine if hCO2 impacts Brachypodium distachyon immunity, a model plant for temperate cereals. Plants were grown in air (430 ppm CO2) and at two high CO2 conditions, one that is relevant to projections within the coming century (1000 ppm) and a concentration sufficient to saturate photosynthesis (3000 ppm). The following measurements were performed: phenotyping and growth, salicylic acid contents, pathogen resistance tests, and RNAseq analysis of the transcriptome. Improved shoot development was observed at both 1000 and 3000 ppm. A transcriptomic analysis pointed to an increase in primary metabolism capacity under hCO2. Alongside this effect, up-regulation of genes associated with secondary metabolism was also observed. This effect was especially evident for the terpenoid and phenylpropanoid pathways, and was accompanied by enhanced expression of immunity-related genes and accumulation of salicylic acid. Pathogen tests using the fungus Magnaporthe oryzae revealed that hCO2 had a complex effect, with enhanced susceptibility to infection but no increase in fungal development. The study reveals that immunity in B. distachyon is modulated by growth at hCO2 and allows identification of pathways that might play a role in this effect.

Brachypodium BdCHS is a homolog of Arabidopsis AtCHS involved in the synthesis of flavonoids and lateral root development

Flavonoids are a kind of plant-specific secondary metabolites, which play an important role in regulating plant growth and development, stress response, and also have medicinal value. Chalcone synthase is the key enzyme in the synthesis of flavonoids. The function of chalcone synthase in Arabidopsis thaliana has been well studied, but its homologous protein in Brachypodium distachyon has not been reported. In this study, we identified a homolog of AtCHS in B. distachyon, named BdCHS, and described its function. Phylogenetic tree analysis showed that BdCHS was most closely related to CHS in Triticum aestivum. Transgene analysis revealed that BdCHS protein was localized in the cytoplasm of Arabidopsis root cells. BdCHS protein can complement the phenotype of AtCHS mutants with lighter seed coat color and increased lateral root density. The content of superoxide anion in the cortical cells above the lateral root primordium in AtCHS mutants was higher than that in the wild-type, and BdCHS protein could restore the content of superoxide anion in AtCHS mutant to the level of that in the wild-type. The results showed that BdCHS was a functional homolog of AtCHS, which laid a foundation for the subsequent application of BdCHS in genetic breeding and crop improvement.

Complexity and regulation of age-dependent alternative splicing in Brachypodium distachyon

Alternative splicing (AS) is a gene regulatory mechanism that generates multiple transcripts of the same gene precursor by the spliceosome complex, promoting messenger RNA complexity and proteome diversity. Although AS is extensively studied in response to environmental stresses, whether it mediates age-dependent development and how it is adjusted by growth transitions are largely unknown. Here, we comprehensively explored the AS landscape at different development stages in the grass model plant Brachypodium (Brachypodium distachyon). We identified abundant coding genes and non-coding transcripts subject to dynamic AS regulation during juvenile, adult, and reproductive transitions. Moreover, we revealed that SC35-LIKE SPLICING FACTOR 33 (SCL33), a serine/arginine-rich splicing factor in spliceosomes, plays a redundant and antagonistic role with its putative paralog, SCL33L, in regulating intron assembly across distinct developmental stages. In addition, we determined global AS variations in microRNA156 (miR156)-overproducing plants, in which growth transitions are delayed, and found that SPLs were regulated by miR156 in intron retention alteration in addition to mRNA clearance and translation inhibition manners. Finally, we demonstrated a complex regulatory process of age-dependent AS events in B. distachyon that were coincidently or separately regulated by miR156 and SCL33/SCL33L. These results illustrate a substantial -machinery of AS that mediates phase transitions in plants.

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.

Fast-track transformation and genome editing in Brachypodium distachyon

BACKGROUND

Even for easy-to-transform species or genotypes, the creation of transgenic or edited plant lines remains a significant bottleneck. Thus, any technical advance that accelerates the regeneration and transformation process is welcome. So far, methods to produce Brachypodium distachyon (Bd) transgenics span at least 14 weeks from the start of tissue culture to the recovery of regenerated plantlets.

RESULTS

We have previously shown that embryogenic somatic tissues grow in the scutellum of immature zygotic Bd embryos within 3 days of in vitro induction with exogenous auxin and that the development of secondary embryos can be initiated immediately thereafter. Here, we further demonstrate that such pluripotent reactive tissues can be genetically transformed with Agrobacterium tumefaciens right after the onset of somatic embryogenesis. In brief, immature zygotic embryos are induced for callogenesis for one week, co-cultured with Agrobacterium for three days, then incubated on callogenesis selective medium for three weeks, and finally transferred on selective regeneration medium for up to three weeks to obtain plantlets ready for rooting. This 7-to-8-week procedure requires only three subcultures. Its validation includes the molecular and phenotype characterization of Bd lines carrying transgenic cassettes and novel CRISPR/Cas9-generated mutations in two independent loci coding for nitrate reductase enzymes (BdNR1 and BdNR2).

CONCLUSIONS

With a short callogenesis stage and streamlined in vitro regeneration following co-cultivation with Agrobacterium, transgenic and edited T0 Bd plantlets can be produced in about 8 weeks, a gain of one to two months compared to previously published methods, with no reduction in transformation efficiency and at lower costs.

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 CC-NB-LRR protein BSR1 from Brachypodium confers resistance to Barley stripe mosaic virus in gramineous plants by recognising TGB1 movement protein

Although some nucleotide binding, leucine-rich repeat immune receptor (NLR) proteins conferring resistance to specific viruses have been identified in dicot plants, NLR proteins involved in viral resistance have not been described in monocots. We have used map-based cloning to isolate the CC-NB-LRR (CNL) Barley stripe mosaic virus (BSMV) resistance gene barley stripe resistance 1 (BSR1) from Brachypodium distachyon Bd3-1 inbred line. Stable BSR1 transgenic Brachypodium line Bd21-3, barley (Golden Promise) and wheat (Kenong 199) plants developed resistance against BSMV ND18 strain. Allelic variation analyses indicated that BSR1 is present in several Brachypodium accessions collected from countries in the Middle East. Protein domain swaps revealed that the intact LRR domain and the C-terminus of BSR1 are required for resistance. BSR1 interacts with the BSMV ND18 TGB1 protein in planta and shows temperature-sensitive antiviral resistance. The R³⁹⁰ and T³⁹² residues of TGB1_ND (ND18 strain) and the G¹⁹⁶ and K¹⁹⁷ residues within the BSR1 P-loop motif are key amino acids required for immune activation. BSR1 is the first cloned virus resistance gene encoding a typical CNL protein in monocots, highlighting the utility of the Brachypodium model for isolation and analysis of agronomically important genes for crop improvement.

Plasticity's role in adaptive evolution depends on environmental change components

To forecast extinction risks of natural populations under climate change and direct human impacts, an integrative understanding of both phenotypic plasticity and adaptive evolution is essential. To date, the evidence for whether, when, and how much plasticity facilitates adaptive responses in changing environments is contradictory. We argue that explicitly considering three key environmental change components - rate of change, variance, and temporal autocorrelation - affords a unifying framework of the impact of plasticity on adaptive evolution. These environmental components each distinctively effect evolutionary and ecological processes underpinning population viability. Using this framework, we develop expectations regarding the interplay between plasticity and adaptive evolution in natural populations. This framework has the potential to improve predictions of population viability in a changing world.

Brachypodium: 20 years as a grass biology model system; the way forward?

It has been 20 years since Brachypodium distachyon was suggested as a model grass species, but ongoing research now encompasses the entire genus. Extensive Brachypodium genome sequencing programmes have provided resources to explore the determinants and drivers of population diversity. This has been accompanied by cytomolecular studies to make Brachypodium a platform to investigate speciation, polyploidisation, perenniality, and various aspects of chromosome and interphase nucleus organisation. The value of Brachypodium as a functional genomic platform has been underscored by the identification of key genes for development, biotic and abiotic stress, and cell wall structure and function. While Brachypodium is relevant to the biofuel industry, its impact goes far beyond that as an intriguing model to study climate change and combinatorial stress.

Patterns of pan-genome occupancy and gene coexpression under water-deficit in Brachypodium distachyon

Natural populations are characterized by abundant genetic diversity driven by a range of different types of mutation. The tractability of sequencing complete genomes has allowed new insights into the variable composition of genomes, summarized as a species pan-genome. These analyses demonstrate that many genes are absent from the first reference genomes, whose analysis dominated the initial years of the genomic era. Our field now turns towards understanding the functional consequence of these highly variable genomes. Here, we analysed weighted gene coexpression networks from leaf transcriptome data for drought response in the purple false brome Brachypodium distachyon and the differential expression of genes putatively involved in adaptation to this stressor. We specifically asked whether genes with variable "occupancy" in the pan-genome - genes which are either present in all studied genotypes or missing in some genotypes - show different distributions among coexpression modules. Coexpression analysis united genes expressed in drought-stressed plants into nine modules covering 72 hub genes (87 hub isoforms), and genes expressed under controlled water conditions into 13 modules, covering 190 hub genes (251 hub isoforms). We find that low occupancy pan-genes are under-represented among several modules, while other modules are over-enriched for low-occupancy pan-genes. We also provide new insight into the regulation of drought response in B. distachyon, specifically identifying one module with an apparent role in primary metabolism that is strongly responsive to drought. Our work shows the power of integrating pan-genomic analysis with transcriptomic data using factorial experiments to understand the functional genomics of environmental response.

Improving bread wheat yield through modulating an unselected AP2/ERF gene

Crop breeding heavily relies on natural genetic variation. However, additional new variations are desired to meet the increasing human demand. Inflorescence architecture determines grain number per spike, a major determinant of bread wheat (Triticum aestivum L.) yield. Here, using Brachypodium distachyon as a wheat proxy, we identified DUO-B1, encoding an APETALA2/ethylene response factor (AP2/ERF) transcription factor, regulating spike inflorescence architecture in bread wheat. Mutations of DUO-B1 lead to mild supernumerary spikelets, increased grain number per spike and, importantly, increased yield under field conditions without affecting other major agronomic traits. DUO-B1 suppresses cell division and promotes the expression of BH^t/WFZP, whose mutations could lead to branched 'miracle-wheat'. Pan-genome analysis indicated that DUO-B1 has not been utilized in breeding, and holds promise to increase wheat yield further.

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.

Manipulation of Insect Vectors’ Host Selection Behavior by Barley Yellow Dwarf Virus Is Dependent on the Host Plant Species and Viral Co-Infection

Previous studies have shown that vector-borne viruses can manipulate the host selection behavior of insect vectors, yet the tripartite interactions of pathogens, host plants and insect vectors have been documented only in a limited number of pathosystems. Here, we report that the host selection behavior of the insect vector of barley yellow dwarf virus-PAV (BYDV-PAV) and cereal yellow dwarf virus-RPS (CYDV-RPS) is dependent on the host plant species and viral co-infection. This study shows that a model cereal plant, Brachypodium distachyon, is a suitable host plant for examining tripartite interactions with BYDV-PAV and CYDV-RPS. We reveal that BYDV-PAV has a different effect on the host selection behavior of its insect vector depending on the host plant species. Viruliferous aphids significantly prefer non-infected plants to virus-infected wheat plants, whereas viral infection on a novel host plant, B. distachyon, is not implicated in the attraction of either viruliferous or nonviruliferous aphids. Furthermore, our findings show that multiple virus infections of wheat with BYDV-PAV and CYDV-RPS alter the preference of their vector aphid. This result indicates that BYDV-PAV acquisition alters the insect vector’s host selection, thereby varying the spread of multiple viruses.

Transcriptional repression of TaNOX10 by TaWRKY19 compromises ROS generation and enhances wheat susceptibility to stripe rust

Reactive oxygen species (ROS) are vital for plant immunity and regulation of their production is crucial for plant health. While the mechanisms that elicit ROS production have been relatively well studied, those that repress ROS generation are less well understood. Here, via screening Brachypodium distachyon RNA interference mutants, we identified BdWRKY19 as a negative regulator of ROS generation whose knockdown confers elevated resistance to the rust fungus Puccinia brachypodii. The three wheat paralogous genes TaWRKY19 are induced during infection by virulent P. striiformis f. sp. tritici (Pst) and have partially redundant roles in resistance. The stable overexpression of TaWRKY19 in wheat increased susceptibility to an avirulent Pst race, while mutations in all three TaWRKY19 copies conferred strong resistance to Pst by enhancing host plant ROS accumulation. We show that TaWRKY19 is a transcriptional repressor that binds to a W-box element in the promoter of TaNOX10, which encodes an NADPH oxidase and is required for ROS generation and host resistance to Pst. Collectively, our findings reveal that TaWRKY19 compromises wheat resistance to the fungal pathogen and suggest TaWRKY19 as a potential target to improve wheat resistance to the commercially important wheat stripe rust fungus.

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

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

The wild grass Brachypodium distachyon as a developmental model system

The arrival of cheap and high-throughput sequencing paired with efficient gene editing technologies allows us to use non-traditional model systems and mechanistically approach biological phenomena beyond what was conceivable just a decade ago. Venturing into different model systems enables us to explore for example clade-specific environmental responses to changing climates or the genetics and development of clade-specific organs, tissues and cell types. We-both early career researchers working with the wild grass model Brachypodium distachyon-want to use this review to (1) highlight why we think B. distachyon is a fantastic grass developmental model system, (2) summarize the tools and resources that have enabled discoveries made in B. distachyon, and (3) discuss a handful of developmental biology vignettes made possible by using B. distachyon as a model system. Finally, we want to conclude by (4) relating our personal stories with this emerging model system and (5) share what we think is important to consider before starting work with an emerging model system.

Brachypodium distachyon UNICULME4 and LAXATUM-A are redundantly required for development

In cultivated grasses, tillering, leaf, and inflorescence architecture, as well as abscission ability, are major agronomical traits. In barley (Hordeum vulgare), maize (Zea mays), rice (Oryza sativa), and brachypodium (Brachypodium distachyon), NOOT-BOP-COCH-LIKE (NBCL) genes are essential regulators of vegetative and reproductive development. Grass species usually possess 2-4 NBCL copies and until now a single study in O. sativa showed that the disruption of all NBCL genes strongly altered O. sativa leaf development. To improve our understanding of the role of NBCL genes in grasses, we extended the study of the two NBCL paralogs BdUNICULME4 (CUL4) and BdLAXATUM-A (LAXA) in the nondomesticated grass B. distachyon. For this, we applied reversed genetics and generated original B. distachyon single and double nbcl mutants by clustered regularly interspaced short palindromic repeats - CRISPR associated protein 9 (CRISPR-Cas9) approaches and genetic crossing between nbcl targeting induced local lesions in genomes (TILLING) mutants. Through the study of original single laxa CRISPR-Cas9 null alleles, we validated functions previously proposed for LAXA in tillering, leaf patterning, inflorescence, and flower development and also unveiled roles for these genes in seed yield. Furthermore, the characterization of cul4laxa double mutants revealed essential functions for nbcl genes in B. distachyon development, especially in the regulation of tillering, stem cell elongation and secondary cell wall composition as well as for the transition toward the reproductive phase. Our results also highlight recurrent antagonist interactions between NBCLs occurring in multiple aspects of B. distachyon development.

Migration without interbreeding: Evolutionary history of a highly selfing Mediterranean grass inferred from whole genomes

Wild plant populations show extensive genetic subdivision and are far from the ideal of panmixia which permeates population genetic theory. Understanding the spatial and temporal scale of population structure is therefore fundamental for empirical population genetics - and of interest in itself, as it yields insights into the history and biology of a species. In this study we extend the genomic resources for the wild Mediterranean grass Brachypodium distachyon to investigate the scale of population structure and its underlying history at whole-genome resolution. A total of 86 accessions were sampled at local and regional scales in Italy and France, which closes a conspicuous gap in the collection for this model organism. The analysis of 196 accessions, spanning the Mediterranean from Spain to Iraq, suggests that the interplay of high selfing and seed dispersal rates has shaped genetic structure in B. distachyon. At the continental scale, the evolution in B. distachyon is characterized by the independent expansion of three lineages during the Upper Pleistocene. Today, these lineages may occur on the same meadow yet do not interbreed. At the regional scale, dispersal and selfing interact and maintain high genotypic diversity, thus challenging the textbook notion that selfing in finite populations implies reduced diversity. Our study extends the population genomic resources for B. distachyon and suggests that an important use of this wild plant model is to investigate how selfing and dispersal, two processes typically studied separately, interact in colonizing plant species.

Enriched root bacterial microbiome in invaded vs native ranges of the model grass allotetraploid Brachypodium hybridum

Invasive species can shift the composition of key soil microbial groups, thus creating novel soil microbial communities. To better understand the biological drivers of invasion, we studied plant-microbial interactions in species of the Brachypodium distachyon complex, a model system for functional genomic studies of temperate grasses and bioenergy crops. While Brachypodium hybridum invasion in California is in an incipient stage, threatening natural and agricultural systems, its diploid progenitor species B. distachyon is not invasive in California. We investigated the root, soil, and rhizosphere bacterial composition of Brachypodium hybridum in both its native and invaded range, and of B. distachyon in the native range. We used high-throughput, amplicon sequencing to evaluate if the bacteria associated with these plants differ, and whether biotic controls may be driving B. hybridum invasion. Bacterial community composition of B. hybridum differed based on provenance (native or invaded range) for root, rhizosphere, and bulk soils, as did the abundance of dominant bacterial taxa. Bacteroidetes, Cyanobacteria and Bacillus spp. (species) were significantly more abundant in B. hybridum roots from the invaded range, whereas Proteobacteria, Firmicutes, Erwinia and Pseudomonas were more abundant in the native range roots. Brachypodium hybridum forms novel biotic interactions with a diverse suite of rhizosphere microbes from the invaded range, which may not exert a similar influence within its native range, ostensibly contributing to B. hybridum’s invasiveness. These associated plant microbiomes could inform future management approaches for B. hybridum in its invaded range and could be key to understanding, predicting, and preventing future plant invasions.

Comparative Genomics, Evolution, and Drought-Induced Expression of Dehydrin Genes in Model Brachypodium Grasses

Dehydration proteins (dehydrins, DHNs) confer tolerance to water-stress deficit in plants. We performed a comparative genomics and evolutionary study of DHN genes in four model Brachypodium grass species. Due to limited knowledge on dehydrin expression under water deprivation stress in Brachypodium, we also performed a drought-induced gene expression analysis in 32 ecotypes of the genus' flagship species B. distachyon showing different hydric requirements. Genomic sequence analysis detected 10 types of dehydrin genes (Bdhn) across the Brachypodium species. Domain and conserved motif contents of peptides encoded by Bdhn genes revealed eight protein architectures. Bdhn genes were spread across several chromosomes. Selection analysis indicated that all the Bdhn genes were constrained by purifying selection. Three upstream cis-regulatory motifs (BES1, MYB124, ZAT) were detected in several Bdhn genes. Gene expression analysis demonstrated that only four Bdhn1-Bdhn2, Bdhn3, and Bdhn7 genes, orthologs of wheat, barley, rice, sorghum, and maize genes, were expressed in mature leaves of B. distachyon and that all of them were more highly expressed in plants under drought conditions. Brachypodium dehydrin expression was significantly correlated with drought-response phenotypic traits (plant biomass, leaf carbon and proline contents and water use efficiency increases, and leaf water and nitrogen content decreases) being more pronounced in drought-tolerant ecotypes. Our results indicate that dehydrin type and regulation could be a key factor determining the acquisition of water-stress tolerance in grasses.

The origin and evolution of a two-component system of paralogous genes encoding the centromeric histone CENH3 in cereals

BACKGROUND

The cereal family Poaceae is one of the largest and most diverse angiosperm families. The central component of centromere specification and function is the centromere-specific histone H3 (CENH3). Some cereal species (maize, rice) have one copy of the gene encoding this protein, while some (wheat, barley, rye) have two. We applied a homology-based approach to sequenced cereal genomes, in order to finally trace the mutual evolution of the structure of the CENH3 genes and the nearby regions in various tribes.

RESULTS

We have established that the syntenic group or the CENH3 locus with the CENH3 gene and the boundaries defined by the CDPK2 and bZIP genes first appeared around 50 Mya in a common ancestor of the subfamilies Bambusoideae, Oryzoideae and Pooideae. This locus came to Pooideae with one copy of CENH3 in the most ancient tribes Nardeae and Meliceae. The βCENH3 gene as a part of the locus appeared in the tribes Stipeae and Brachypodieae around 35-40 Mya. The duplication was accompanied by changes in the exon-intron structure. Purifying selection acts mostly on αCENH3s, while βCENH3s form more heterogeneous structures, in which clade-specific amino acid motifs are present. In barley species, the βCENH3 gene assumed an inverted orientation relative to αCENH3 and the CDPK2 gene was substituted with LHCB-l. As the evolution and domestication of plant species went on, the locus was growing in size due to an increasing distance between αCENH3 and βCENH3 because of a massive insertion of the main LTR-containing retrotransposon superfamilies, gypsy and copia, without any evolutionary preference on either of them. A comparison of the molecular structure of the locus in the A, B and D subgenomes of the hexaploid wheat T. aestivum showed that invasion by mobile elements and concomitant rearrangements took place in an independent way even in evolutionarily close species.

CONCLUSIONS

The CENH3 duplication in cereals was accompanied by changes in the exon-intron structure of the βCENH3 paralog. The observed general tendency towards the expansion of the CENH3 locus reveals an amazing diversity of ways in which different species implement the scenario described in this paper.

Differential Gene Expression of Brachypodium distachyon Roots Colonized by Gluconacetobacter diazotrophicus and the Role of BdCESA8 in the Colonization

Alternatives to synthetic nitrogen fertilizer are needed to reduce the costs of crop production and offset environmental damage. Nitrogen-fixing bacterium Gluconacetobacter diazotrophicus has been proposed as a possible biofertilizer for monocot crop production. However, the colonization of G. diazotrophicus in most monocot crops is limited and deep understanding of the response of host plants to G. diazotrophicus colonization is still lacking. In this study, the molecular response of the monocot plant model Brachypodium distachyon was studied during G. diazotrophicus root colonization. The gene expression profiles of B. distachyon root tissues colonized by G. diazotrophicus were generated via next-generation RNA sequencing, and investigated through gene ontology and metabolic pathway analysis. The RNA sequencing results indicated that Brachypodium is actively involved in G. diazotrophicus colonization via cell wall synthesis. Jasmonic acid, ethylene, gibberellin biosynthesis. nitrogen assimilation, and primary and secondary metabolite pathways are also modulated to accommodate and control the extent of G. diazotrophicus colonization. Cellulose synthesis is significantly downregulated during colonization. The loss of function mutant for Brachypodium cellulose synthase 8 (BdCESA8) showed decreased cellulose content in xylem and increased resistance to G. diazotrophicus colonization. This result suggested that the cellulose synthesis of the secondary cell wall is involved in G. diazotrophicus colonization. The results of this study provide insights for future research in regard to gene manipulation for efficient colonization of nitrogen-fixing bacteria in Brachypodium and monocot crops.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Mitigating tradeoffs in plant breeding

Tradeoffs among plant traits help maintain relative fitness under unpredictable conditions and maximize reproductive success. However, modifying tradeoffs is a breeding challenge since many genes of minor effect are involved. The intensive crosstalk and fine-tuning between growth and defense responsive phytohormones via transcription factors optimizes growth, reproduction, and stress tolerance. There are regulating genes in grain crops that deploy diverse functions to overcome tradeoffs, e.g., miR-156-IPA1 regulates crosstalk between growth and defense to achieve high disease resistance and yield, while OsALDH2B1 loss of function causes imbalance among defense, growth, and reproduction in rice. GNI-A1 regulates seed number and weight in wheat by suppressing distal florets and altering assimilate distribution of proximal seeds in spikelets. Knocking out ABA-induced transcription repressors (AITRs) enhances abiotic stress adaptation without fitness cost in Arabidopsis. Deploying AITRs homologs in grain crops may facilitate breeding. This knowledge suggests overcoming tradeoffs through breeding may expose new ones.

Analysis of Codon Usage Patterns of Six Sequenced Brachypodium distachyon Lines Reveals a Declining CG Skew of the CDSs from the 5'-ends to the 3'-ends

Brachypodium distachyon, a new monocotyledonous model plant, has received wide attention in biological research due to its small genome and numerous genetic resources. Codon usage bias is an important feature of genes and genomes, and it can be used in transgenic and evolutionary studies. In this study, the nucleotide compositions and patterns of codon usage bias were calculated using Codon W. Additionally, an ENC plot, Parity rule 2 and correspondence analyses were used to explore the major factors influencing codon usage bias patterns. The numbers of hydrogen bonds and skews were used to analyze the GC trend in the 5'-ends of the coding sequences. The results showed that minor differences in the codon usage bias patterns were revealed by the ENC plot, Parity rule 2 and correspondence analyses. The analyses of the CG-skew and the number of hydrogen bonds showed a declining trend in the number of cytosines at the 5'-ends of the CDSs (from the 5'-ends to the 3'-ends), indicating that GC may play a major role in codon usage bias. In addition, our results laid a foundation for the study of codon usage bias patterns in Brachypodium genus and suggested that the GC plays a major role in determining these patterns.

Conserved, divergent and heterochronic gene expression during Brachypodium and Arabidopsis embryo development

KEY MESSAGE

Developmental and transcriptomic analysis of Brachypodium embryogenesis and comparison with Arabidopsis identifies conserved and divergent phases of embryogenesis and reveals widespread heterochrony of developmental gene expression. Embryogenesis, transforming the zygote into the mature embryo, represents a fundamental process for all flowering plants. Current knowledge of cell specification and differentiation during plant embryogenesis is largely based on studies of the dicot model plant Arabidopsis thaliana. However, the major crops are monocots and the transcriptional programs associated with the differentiation processes during embryogenesis in this clade were largely unknown. Here, we combined analysis of cell division patterns with development of a temporal transcriptomic resource during embryogenesis of the monocot model plant Brachypodium distachyon. We found that early divisions of the Brachypodium embryo were highly regular, while later stages were marked by less stereotypic patterns. Comparative transcriptomic analysis between Brachypodium and Arabidopsis revealed that early and late embryogenesis shared a common transcriptional program, whereas mid-embryogenesis was divergent between species. Analysis of orthology groups revealed widespread heterochronic expression of potential developmental regulators between the species. Interestingly, Brachypodium genes tend to be expressed at earlier stages than Arabidopsis counterparts, which suggests that embryo patterning may occur early during Brachypodium embryogenesis. Detailed investigation of auxin-related genes shows that the capacity to synthesize, transport and respond to auxin is established early in the embryo. However, while early PIN1 polarity could be confirmed, it is unclear if an active response is mounted. This study presents a resource for studying Brachypodium and grass embryogenesis and shows that divergent angiosperms share a conserved genetic program that is marked by heterochronic gene expression.

TaWRKY13-A Serves as a Mediator of Jasmonic Acid-Related Leaf Senescence by Modulating Jasmonic Acid Biosynthesis

Leaf senescence is crucial for crop yield and quality. Transcriptional regulation is a key step for integrating various senescence-related signals into the nucleus. However, few regulators of senescence implicating transcriptional events have been functionally characterized in wheat. Based on our RNA-seq data, we identified a WRKY transcription factor, TaWRKY13-A, that predominately expresses at senescent stages. By using the virus-induced gene silencing (VIGS) method, we manifested impaired transcription of TaWRKY13-A leading to a delayed leaf senescence phenotype in wheat. Moreover, the overexpression (OE) of TaWRKY13-A accelerated the onset of leaf senescence under both natural growth condition and darkness in Brachypodium distachyon and Arabidopsis thaliana. Furthermore, by physiological and molecular investigations, we verified that TaWRKY13-A participates in the regulation of leaf senescence via jasmonic acid (JA) pathway. The expression of JA biosynthetic genes, including AtLOX6, was altered in TaWRKY13-A-overexpressing Arabidopsis. We also demonstrated that TaWRKY13-A can interact with the promoter of AtLOX6 and TaLOX6 by using the electrophoretic mobility shift assay (EMSA) and luciferase reporter system. Consistently, we detected a higher JA level in TaWRKY13-A-overexpressing lines than that in Col-0. Moreover, our data suggested that TaWRKY13-A is partially functional conserved with AtWRKY53 in age-dependent leaf senescence. Collectively, this study manifests TaWRKY13-A as a positive regulator of JA-related leaf senescence, which could be a new clue for molecular breeding in wheat.

Arabis alpina: A perennial model plant for ecological genomics and life-history evolution

Many model organisms were chosen and achieved prominence because of an advantageous combination of their life‐history characteristics, genetic properties and also practical considerations. Discoveries made in Arabidopsis thaliana, the most renowned noncrop plant model species, have markedly stimulated studies in other species with different biology. Within the family Brassicaceae, the arctic–alpine Arabis alpina has become a model complementary to Arabidopsis thaliana to study the evolution of life‐history traits, such as perenniality, and ecological genomics in harsh environments. In this review, we provide an overview of the properties that facilitated the rapid emergence of A. alpina as a plant model. We summarize the evolutionary history of A. alpina, including genomic aspects, the diversification of its mating system and demographic properties, and we discuss recent progress in the molecular dissection of developmental traits that are related to its perennial life history and environmental adaptation. From this published knowledge, we derive open questions that might inspire future research in A. alpina, other Brassicaceae species or more distantly related plant families.

YSL3-mediated copper distribution is required for fertility, seed size and protein accumulation in Brachypodium

Addressing the looming global food security crisis requires the development of high-yielding crops. In agricultural soils, deficiency in the micronutrient copper significantly decreases grain yield in wheat (Triticum aestivum), a globally important crop. In cereals, grain yield is determined by inflorescence architecture, flower fertility, grain size, and weight. Whether copper is involved in these processes, and how it is delivered to the reproductive organs is not well understood. We show that copper deficiency alters not only the grain set but also flower development in both wheat and its recognized model, Brachypodium distachyon. We then show that the Brachypodium yellow stripe-like 3 (YSL3) transporter localizes to the phloem, transports copper in frog (Xenopus laevis) oocytes, and facilitates copper delivery to reproductive organs and grains. Failure to deliver copper, but not iron, zinc, or manganese to these structures in the ysl3 CRISPR-Cas9 mutant results in delayed flowering, altered inflorescence architecture, reduced floret fertility, grain size, weight, and protein accumulation. These defects are rescued by copper supplementation and are complemented by YSL3 cDNA. This knowledge will help to devise sustainable approaches for improving grain yield in regions where soil quality is a major obstacle for crop production. Copper distribution by a phloem-localized transporter is essential for the transition to flowering, inflorescence architecture, floret fertility, size, weight, and protein accumulation in seeds.

Plant thermotropism: an underexplored thermal engagement and avoidance strategy

Various strategies evolved in plants to adjust the position of organs relative to the prevailing temperature condition, which allows optimal plant growth and performance. Such responses are classically separated into nastic and tropic responses. During plant thermotropic responses, organs move towards (engage) or away (avoid) from a directional temperature cue. Despite thermotropism being a classic botanical concept, the underlying ecological function and molecular and biophysical mechanisms remain poorly understood to this day. This contrasts to the relatively well-studied thermonastic movements (hyponasty) of e.g., rosette leaves. In this review, we provide an update on the current knowledge on plant thermotropisms and propose directions for future research and application.