Natural Variation in Brachypodium Links Vernalization and Flowering Time Loci as Major Flowering Determinants

Comprehensive Genome-Wide Natural Variation and Expression Analysis of Tubby-like Proteins Gene Family in Brachypodium distachyon

The Tubby-like proteins (TLPs) gene family is a group of transcription factors found in both animals and plants. In this study, we identified twelve B. distachyon TLPs, divided into six groups based on conserved domains and evolutionary relationships. We predicted cis-regulatory elements involved in light, hormone, and biotic and abiotic stresses. The expression patterns in response to light and hormones revealed that BdTLP3, 4, 7, and 14 are involved in light responses, and BdTLP1 is involved in ABA responses. Furthermore, BdTLP2, 7, 9, and 13 are expressed throughout vegetative and reproductive stages, whereas BdTLP1, 3, 5, and 14 are expressed at germinating grains and early vegetative development, and BdTLP4, 6, 8, and 10 are expressed at the early reproduction stage. The natural variation in the eleven most diverged B. distachyon lines revealed high conservation levels of BdTLP1-6 to high variation in BdTLP7-14 proteins. Based on diversifying selection, we identified amino acids in BdTLP1, 3, 8, and 13, potentially substantially affecting protein functions. This analysis provided valuable information for further functional studies to understand the regulation, pathways involved, and mechanism of BdTLPs.

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.

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.

Mutations in the predicted DNA polymerase subunit POLD3 result in more rapid flowering of Brachypodium distachyon

The timing of reproduction is a critical developmental decision in the life cycle of many plant species. Fine mapping of a rapid-flowering mutant was done using whole-genome sequence data from bulked DNA from a segregating F2 mapping populations. The causative mutation maps to a gene orthologous with the third subunit of DNA polymerase δ (POLD3), a previously uncharacterized gene in plants. Expression analyses of POLD3 were conducted via real time qPCR to determine when and in what tissues the gene is expressed. To better understand the molecular basis of the rapid-flowering phenotype, transcriptomic analyses were conducted in the mutant vs wild-type. Consistent with the rapid-flowering mutant phenotype, a range of genes involved in floral induction and flower development are upregulated in the mutant. Our results provide the first characterization of the developmental and gene expression phenotypes that result from a lesion in POLD3 in plants.

Mechanical stimulation in Brachypodium distachyon: Implications for fitness, productivity, and cell wall properties

Mechanical stimulation, including exposure to wind, is a common environmental variable for plants. However, knowledge about the morphogenetic response of the grasses (Poaceae) to mechanical stimulation and impact on relevant agronomic traits is very limited. Two natural accessions of Brachypodium distachyon were exposed to wind and mechanical treatments. We surveyed a wide range of stem-related traits to determine the effect of the two treatments on plant growth, development, and stem biomass properties. Both treatments induced significant quantitative changes across multiple scales, from the whole plant down to cellular level. The two treatments resulted in shorter stems, reduced biomass, increased tissue rigidity, delayed flowering, and reduced seed yield in both accessions. Among changes in cell wall-related features, a substantial increase in lignin content and pectin methylesterase activity was most notable. Mechanical stimulation also reduced the enzymatic sugar release from the cell wall, thus increasing biomass recalcitrance. Notably, treatments had a distinct and opposite effect on vascular bundle area in the two accessions, suggesting genetic variation in modulating these responses to mechanical stimulation. Our findings highlight that exposure of grasses to mechanical stimulation is a relevant environmental factor affecting multiple traits important for their utilization in food, feed, and bioenergy applications.

Molecular Regulation of Plant Responses to Environmental Temperatures

Temperature is a key factor governing the growth and development, distribution, and seasonal behavior of plants. The entire plant life cycle is affected by environmental temperatures. Plants grow rapidly and exhibit specific changes in morphology under mild average temperature conditions, a response termed thermomorphogenesis. When exposed to chilling or moist chilling low temperatures, flowering or seed germination is accelerated in some plant species; these processes are known as vernalization and cold stratification, respectively. Interestingly, once many temperate plants are exposed to chilling temperatures for some time, they can acquire the ability to resist freezing stress, a process termed cold acclimation. In the face of global climate change, heat stress has emerged as a frequent challenge, which adversely affects plant growth and development. In this review, we summarize and discuss recent progress in dissecting the molecular mechanisms regulating plant thermomorphogenesis, vernalization, and responses to extreme temperatures. We also discuss the remaining issues that are crucial for understanding the interactions between plants and temperature.

The florigen interactor BdES43 represses flowering in the model temperate grass Brachypodium distachyon

FLOWERING LOCUS T (FT) protein, physiologically florigen, has been identified as a system integrator of numerous flowering time pathways in many studies, and its homologs are found throughout the plant lineage. It is important to uncover how precisely florigenic homologs contribute to flowering initiation and how these factors interact genetically. Here we dissected the function of Brachypodium FT orthologs BdFTL1 and BdFTL2 using overexpression and gene-editing experiments. Transgenic assays showed that both BdFTL1 and BdFTL2 could promote flowering, whereas BdFTL2 was essential for flowering initiation. Notably, BdFTL1 is subject to alternative splicing (AS), and its transcriptional level and AS are significantly affected by BdFTL2. Additionally, BdFTL2 could bind with the PHD-containing protein BdES43, an H3K4me3 reader. Furthermore, BdES43 was antagonistic to BdFTL2 in flowering initiation in a transcription-dependent manner and significantly affected BdFTL1 expression. BdFTL2, BdES43 and H3K4me3 also had highly similar distribution patterns within the BdFTL1 locus, indicating their interplay in regulating target genes. Taken together, florigen BdFTL2 functions as a potential epigenetic effector of BdFTL1 by interacting with a BdES43-H3K4me3 complex. This finding provides an additional insight for the regulatory mechanism underlying the multifaceted roles of florigen.

The Role of FLOWERING LOCUS C Relatives in Cereals

FLOWERING LOCUS C (FLC) is one of the best characterized genes in plant research and is integral to vernalization-dependent flowering time regulation. Yet, despite the abundance of information on this gene and its relatives in Arabidopsis thaliana, the role FLC genes play in other species, in particular cereal crops and temperate grasses, remains elusive. This has been due in part to the comparative reduced availability of bioinformatic and mutant resources in cereals but also on the dominant effect in cereals of the VERNALIZATION (VRN) genes on the developmental process most associated with FLC in Arabidopsis. The strong effect of the VRN genes has led researchers to believe that the entire process of vernalization must have evolved separately in Arabidopsis and cereals. Yet, since the confirmation of the existence of FLC-like genes in monocots, new light has been shed on the roles these genes play in both vernalization and other mechanisms to fine tune development in response to specific environmental conditions. Comparisons of FLC gene function and their genetic and epigenetic regulation can now be made between Arabidopsis and cereals and how they overlap and diversify is coming into focus. With the advancement of genome editing techniques, further study on these genes is becoming increasingly easier, enabling us to investigate just how essential FLC-like genes are to modulating flowering time behavior in cereals.

Experiencing winter for spring flowering: A molecular epigenetic perspective on vernalization

Many over-wintering plants, through vernalization, overcome a block to flowering and thus acquire competence to flower in the following spring after experiencing prolonged cold exposure or winter cold. The vernalization pathways in different angiosperm lineages appear to have convergently evolved to adapt to temperate climates. Molecular and epigenetic mechanisms for vernalization regulation have been well studied in the crucifer model plant Arabidopsis thaliana. Here, we review recent progresses on the vernalization pathway in Arabidopsis. In addition, we summarize current molecular and genetic understandings of vernalization regulation in temperate grasses including wheat and Brachypodium, two monocots from Pooideae, followed by a brief discussion on divergence of the vernalization pathways between Brassicaceae and Pooideae.

Population Genomic Approaches for Weed Science

Genomic approaches are opening avenues for understanding all aspects of biological life, especially as they begin to be applied to multiple individuals and populations. However, these approaches typically depend on the availability of a sequenced genome for the species of interest. While the number of genomes being sequenced is exploding, one group that has lagged behind are weeds. Although the power of genomic approaches for weed science has been recognized, what is needed to implement these approaches is unfamiliar to many weed scientists. In this review we attempt to address this problem by providing a primer on genome sequencing and provide examples of how genomics can help answer key questions in weed science such as: (1) Where do agricultural weeds come from; (2) what genes underlie herbicide resistance; and, more speculatively, (3) can we alter weed populations to make them easier to control? This review is intended as an introduction to orient weed scientists who are thinking about initiating genome sequencing projects to better understand weed populations, to highlight recent publications that illustrate the potential for these methods, and to provide direction to key tools and literature that will facilitate the development and execution of weed genomic projects.

Genetic Variation for Seed Metabolite Levels in Brachypodium distachyon

Metabolite composition and concentrations in seed grains are important traits of cereals. To identify the variation in the seed metabolotypes of a model grass, namely Brachypodium distachyon, we applied a widely targeted metabolome analysis to forty inbred lines of B. distachyon and examined the accumulation patterns of 183 compounds in the seeds. By comparing the metabolotypes with the population structure of these lines, we found signature metabolites that represent different accumulation patterns for each of the three B. distachyon subpopulations. Moreover, we found that thirty-seven metabolites exhibited significant differences in their accumulation between the lines Bd21 and Bd3-1. Using a recombinant inbred line (RIL) population from a cross between Bd3-1 and Bd21, we identified the quantitative trait loci (QTLs) linked with this variation in the accumulation of thirteen metabolites. Our metabolite QTL analysis illustrated that different genetic factors may presumably regulate the accumulation of 4-pyridoxate and pyridoxamine in vitamin B6 metabolism. Moreover, we found two QTLs on chromosomes 1 and 4 that affect the accumulation of an anthocyanin, chrysanthemin. These QTLs genetically interacted to regulate the accumulation of this compound. This study demonstrates the potential for metabolite QTL mapping in B. distachyon and provides new insights into the genetic dissection of metabolomic traits in temperate grasses.

Divergent roles of FT-like 9 in flowering transition under different day lengths in Brachypodium distachyon

Timing of reproductive transition is precisely modulated by environmental cues in flowering plants. Facultative long-day plants, including Arabidopsis and temperate grasses, trigger rapid flowering in long-day conditions (LDs) and delay flowering under short-day conditions (SDs). Here, we characterize a SD-induced FLOWERING LOCUS T ortholog, FT-like 9 (FTL9), that promotes flowering in SDs but inhibits flowering in LDs in Brachypodium distachyon. Mechanistically, like photoperiod-inductive FT1, FTL9 can interact with FD1 to form a flowering activation complex (FAC), but the floral initiation efficiency of FTL9-FAC is much lower than that of FT1-FAC, thereby resulting in a positive role for FTL9 in promoting floral transition when FT1 is not expressed, but a dominant-negative role when FT1 accumulates significantly. We also find that CONSTANS 1 (CO1) can suppress FTL9 in addition to stimulate FT1 to enhance accelerated flowering under LDs. Our findings on the antagonistic functions of FTL9 under different day-length environments will contribute to understanding the multifaceted roles of FT in fine-tune modulation of photoperiodic flowering in plants.

Plant flowering time is modified by day length. Here the authors show that the model grass Brachypodium distachyon expresses different homologs of FT in short and long days to produce floral activator complexes with altered activities contributing to photoperiod-dependence of flowering time.

Global Diversity of the Brachypodium Species Complex as a Resource for Genome-Wide Association Studies Demonstrated for Agronomic Traits in Response to Climate

The development of model systems requires a detailed assessment of standing genetic variation across natural populations. The Brachypodium species complex has been promoted as a plant model for grass genomics with translation to small grain and biomass crops. To capture the genetic diversity within this species complex, thousands of Brachypodium accessions from around the globe were collected and genotyped by sequencing. Overall, 1897 samples were classified into two diploid or allopolyploid species, and then further grouped into distinct inbred genotypes. A core set of diverse B. distachyon diploid lines was selected for whole genome sequencing and high resolution phenotyping. Genome-wide association studies across simulated seasonal environments was used to identify candidate genes and pathways tied to key life history and agronomic traits under current and future climatic conditions. A total of 8, 22, and 47 QTL were identified for flowering time, early vigor, and energy traits, respectively. The results highlight the genomic structure of the Brachypodium species complex, and the diploid lines provided a resource that allows complex trait dissection within this grass model species.

Detection of Race-Specific Resistance Against Puccinia coronata f. sp. avenae in Brachypodium Species

Oat crown rust caused by Puccinia coronata f. sp. avenae is the most destructive foliar disease of cultivated oat. Characterization of genetic factors controlling resistance responses to Puccinia coronata f. sp. avenae in nonhost species could provide new resources for developing disease protection strategies in oat. We examined symptom development and fungal colonization levels of a collection of Brachypodium distachyon and B. hybridum accessions infected with three North American P. coronata f. sp. avenae isolates. Our results demonstrated that colonization phenotypes are dependent on both host and pathogen genotypes, indicating a role for race-specific responses in these interactions. These responses were independent of the accumulation of reactive oxygen species. Expression analysis of several defense-related genes suggested that salicylic acid and ethylene-mediated signaling but not jasmonic acid are components of resistance reaction to P. coronata f. sp. avenae. Our findings provide the basis to conduct a genetic inheritance study to examine whether effector-triggered immunity contributes to nonhost resistance to P. coronata f. sp. avenae in Brachypodium spp.

The genetic architecture of colonization resistance in Brachypodium distachyon to non-adapted stripe rust (Puccinia striiformis) isolates

Multilayered defense responses ensure that plants are hosts to only a few adapted pathogens in the environment. The host range of a plant pathogen depends on its ability to fully overcome plant defense barriers, with failure at any single step sufficient to prevent life cycle completion of the pathogen. Puccinia striiformis, the causal agent of stripe rust (=yellow rust), is an agronomically important obligate biotrophic fungal pathogen of wheat and barley. It is generally unable to complete its life cycle on the non-adapted wild grass species Brachypodium distachyon, but natural variation exists for the degree of hyphal colonization by Puccinia striiformis. Using three B. distachyon mapping populations, we identified genetic loci conferring colonization resistance to wheat-adapted and barley-adapted isolates of P. striiformis. We observed a genetic architecture composed of two major effect QTLs (Yrr1 and Yrr3) restricting the colonization of P. striiformis. Isolate specificity was observed for Yrr1, whereas Yrr3 was effective against all tested P. striiformis isolates. Plant immune receptors of the nucleotide binding, leucine-rich repeat (NB-LRR) encoding gene family are present at the Yrr3 locus, whereas genes of this family were not identified at the Yrr1 locus. While it has been proposed that resistance to adapted and non-adapted pathogens are inherently different, the observation of (1) a simple genetic architecture of colonization resistance, (2) isolate specificity of major and minor effect QTLs, and (3) NB-LRR encoding genes at the Yrr3 locus suggest that factors associated with resistance to adapted pathogens are also critical for non-adapted pathogens.

An ortholog of CURLY LEAF/ENHANCER OF ZESTE like-1 is required for proper flowering in Brachypodium distachyon

Many plants require prolonged exposure to cold to acquire the competence to flower. The process by which cold exposure results in competence is known as vernalization. In Arabidopsis thaliana, vernalization leads to the stable repression of the floral repressor FLOWERING LOCUS C via chromatin modification, including an increase of trimethylation on lysine 27 of histone H3 (H3K27me3) by Polycomb Repressive Complex 2 (PRC2). Vernalization in pooids is associated with the stable induction of a floral promoter, VERNALIZATION 1 (VRN1). From a screen for mutants with a reduced vernalization requirement in the model grass Brachypodium distachyon, we identified two recessive alleles of ENHANCER OF ZESTE-LIKE 1 (EZL1). EZL1 is orthologous to A. thaliana CURLY LEAF 1, a gene that encodes the catalytic subunit of PRC2. B. distachyon ezl1 mutants flower rapidly without vernalization in long-day (LD) photoperiods; thus, EZL1 is required for the proper maintenance of the vegetative state prior to vernalization. Transcriptomic studies in ezl1 revealed mis-regulation of thousands of genes, including ectopic expression of several floral homeotic genes in leaves. Loss of EZL1 results in the global reduction of H3K27me3 and H3K27me2, consistent with this gene making a major contribution to PRC2 activity in B. distachyon. Furthermore, in ezl1 mutants, the flowering genes VRN1 and AGAMOUS (AG) are ectopically expressed and have reduced H3K27me3. Artificial microRNA knock-down of either VRN1 or AG in ezl1-1 mutants partially restores wild-type flowering behavior in non-vernalized plants, suggesting that ectopic expression in ezl1 mutants may contribute to the rapid-flowering phenotype.

Establishment of a vernalization requirement in Brachypodium distachyon requires REPRESSOR OF VERNALIZATION1

A requirement for vernalization, the process by which prolonged cold exposure provides competence to flower, is an important adaptation to temperate climates that ensures flowering does not occur before the onset of winter. In temperate grasses, vernalization results in the up-regulation of VERNALIZATION1 (VRN1) to establish competence to flower; however, little is known about the mechanism underlying repression of VRN1 in the fall season, which is necessary to establish a vernalization requirement. Here, we report that a plant-specific gene containing a bromo-adjacent homology and transcriptional elongation factor S-II domain, which we named REPRESSOR OF VERNALIZATION1 (RVR1), represses VRN1 before vernalization in Brachypodium distachyon That RVR1 is upstream of VRN1 is supported by the observations that VRN1 is precociously elevated in an rvr1 mutant, resulting in rapid flowering without cold exposure, and the rapid-flowering rvr1 phenotype is dependent on VRN1 The precocious VRN1 expression in rvr1 is associated with reduced levels of the repressive chromatin modification H3K27me3 at VRN1, which is similar to the reduced VRN1 H3K27me3 in vernalized plants. Furthermore, the transcriptome of vernalized wild-type plants overlaps with that of nonvernalized rvr1 plants, indicating loss of rvr1 is similar to the vernalized state at a molecular level. However, loss of rvr1 results in more differentially expressed genes than does vernalization, indicating that RVR1 may be involved in processes other than vernalization despite a lack of any obvious pleiotropy in the rvr1 mutant. This study provides an example of a role for this class of plant-specific genes.

Winter Memory throughout the Plant Kingdom: Different Paths to Flowering

Molecular mechanisms contribute to the memory of winter in different plant groups.

Genetic Architecture of Flowering-Time Variation in Brachypodium distachyon

The transition to reproductive development is a crucial step in the plant life cycle, and the timing of this transition is an important factor in crop yields. Here, we report new insights into the genetic control of natural variation in flowering time in Brachypodium distachyon, a nondomesticated pooid grass closely related to cereals such as wheat (Triticum spp.) and barley (Hordeum vulgare L.). A recombinant inbred line population derived from a cross between the rapid-flowering accession Bd21 and the delayed-flowering accession Bd1-1 were grown in a variety of environmental conditions to enable exploration of the genetic architecture of flowering time. A genotyping-by-sequencing approach was used to develop SNP markers for genetic map construction, and quantitative trait loci (QTLs) that control differences in flowering time were identified. Many of the flowering-time QTLs are detected across a range of photoperiod and vernalization conditions, suggesting that the genetic control of flowering within this population is robust. The two major QTLs identified in undomesticated B. distachyon colocalize with VERNALIZATION1/PHYTOCHROME C and VERNALIZATION2, loci identified as flowering regulators in the domesticated crops wheat and barley. This suggests that variation in flowering time is controlled in part by a set of genes broadly conserved within pooid grasses.

Linking Dynamic Phenotyping with Metabolite Analysis to Study Natural Variation in Drought Responses of Brachypodium distachyon

Drought is an important environmental stress limiting the productivity of major crops worldwide. Understanding drought tolerance and possible mechanisms for improving drought resistance is therefore a prerequisite to develop drought-tolerant crops that produce significant yields with reduced amounts of water. Brachypodium distachyon (Brachypodium) is a key model species for cereals, forage grasses, and energy grasses. In this study, initial screening of a Brachypodium germplasm collection consisting of 138 different ecotypes exposed to progressive drought, highlighted the natural variation in morphology, biomass accumulation, and responses to drought stress. A core set of ten ecotypes, classified as being either tolerant, susceptible or intermediate, in response to drought stress, were exposed to mild or severe (respectively, 15 and 0% soil water content) drought stress and phenomic parameters linked to growth and color changes were assessed. When exposed to severe drought stress, phenotypic data and metabolite profiling combined with multivariate analysis revealed a remarkable consistency in separating the selected ecotypes into their different pre-defined drought tolerance groups. Increases in several metabolites, including for the phytohormones jasmonic acid and salicylic acid, and TCA-cycle intermediates, were positively correlated with biomass yield and with reduced yellow pixel counts; suggestive of delayed senescence, both key target traits for crop improvement to drought stress. While metabolite analysis also separated ecotypes into the distinct tolerance groupings after exposure to mild drought stress, similar analysis of the phenotypic data failed to do so, confirming the value of metabolomics to investigate early responses to drought stress. The results highlight the potential of combining the analyses of phenotypic and metabolic responses to identify key mechanisms and markers associated with drought tolerance in both the Brachypodium model plant as well as agronomically important crops.