Update on the genomics and basic biology of Brachypodium: International Brachypodium Initiative (IBI)

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.

Identification, expression, and functional analysis of Hsf and Hsp20 gene families in Brachypodium distachyon under heat stress

Background

The heat shock factor (Hsf) and small heat shock protein (sHsp, also called Hsp20) complex has been identified as a primary component in the protection of plant cells from ubiquitous stresses, particularly heat stress. Our study aimed to characterize and analyze the Hsf and Hsp genes in Brachypodium distachyon, an annual temperate grass and model plant in cereal and grass studies.

Results

We identified 24 Hsf and 18 Hsp20 genes in B. distachyon and explored their evolution in gene organization, sequence features, chromosomal localization, and gene duplication. Our phylogenetic analysis showed that BdHsfs could be divided into three categories and BdHsp20s into ten subfamilies. Further analysis showed that the 3’UTR length of BdHsp20 genes had a negative relationship with their expression under heat stress. Expression analyses indicated that BdHsp20s and BdHsfs were strongly and rapidly induced by high-temperature treatment. Additionally, we constructed a complex regulatory network based on their expression patterns under heat stress. Morphological analysis suggested that the overexpression of five BdHsp20 genes enhanced the seed germination rate and decreased cell death under high temperatures.

Conclusion

Ultimately, our study provided important evolutionary and functional characterizations for future research on the regulatory mechanisms of BdHsp20s and BdHsfs in herbaceous plants under environmental stress.

DNA Methylation-An Epigenetic Mark in Mutagen-Treated Brachypodium distachyon Cells

The chromatin structure is significantly influenced by some epigenetic modifications including DNA methylation. The nuclear organization plays an essential role in the cell response to external stresses including mutagens. We present an analysis of the correlation between epigenetic modifications and the instability of the Brachypodium distachyon genome, which are observed as micronuclei, following maleic hydrazide (MH) and nitroso-N-methylurea (MNU) treatments. We compared the level of DNA methylation in the control (untreated) and mutagen-treated B. distachyon nuclei. An immunostaining method using specific antibodies against modified DNA anti-5-methylcytosine was used for the evaluation of DNA methylation in a single nucleus and micronucleus. Interestingly, we showed an alteration of DNA methylation in cells after mutagenic treatments. The results indicate that DNA methylation might be involved in the response of the B. distachyon genome to mutagenic treatments. This demonstrates that analyses of the epigenetic modifications should be integrated into current plant genetic toxicology in order to explain the mechanisms of DNA damage and repair in plants.

The Metabolic Response of Brachypodium Roots to the Interaction with Beneficial Bacteria Is Affected by the Plant Nutritional Status

The potential of plant growth promoting (PGP) bacteria in improving the performance of plants in suboptimal environments is increasingly acknowledged, but little information is available on the mechanisms underlying this interaction, particularly when plants are subjected to a combination of stresses. In this study, we investigated the effects of the inoculation with the PGP bacteria Azospirillum brasilense (Azospirillum) on the metabolism of the model cereal Brachypodium distachyon (Brachypodium) grown at low temperatures and supplied with insufficient phosphorus. Investigating polar metabolite and lipid fluctuations during early plant development, we found that the bacteria initially elicited a defense response in Brachypodium roots, while at later stages Azospirillum reduced the stress caused by phosphorus deficiency and improved root development of inoculated plants, particularly by stimulating the growth of branch roots. We propose that the interaction of the plant with Azospirillum was influenced by its nutritional status: bacteria were sensed as pathogens while plants were still phosphorus sufficient, but the interaction became increasingly beneficial for the plants as their phosphorus levels decreased. Our results provide new insights on the dynamics of the cereal-PGP bacteria interaction, and contribute to our understanding of the role of beneficial microorganisms in the growth of cereal crops in suboptimal environments.

Comparative analysis of MAPK and MKK gene families reveals differential evolutionary patterns in Brachypodium distachyon inbred lines

BACKGROUND

Mitogen-activated protein kinase (MAPK) cascades are involved with signal transduction in almost every aspect of plant growth and development, as well as biotic and abiotic stress responses. The evolutionary analysis of MAPKs and MKKs in individual or entire plant species has been reported, but the evolutionary patterns in the diverse inbred lines of Brachypodium distachyon are still unclear.

RESULTS

We conducted the systematical molecular evolutionary analysis of B. distachyon. A total of 799 MAPKs and 618 MKKs were identified from 53 B. distachyon inbred lines. Remarkably, only three inbred lines had 16 MPKs and most of those inbred lines lacked MPK7-2 members, whereas 12 MKKs existed in almost all B. distachyon inbred lines. Phylogenetic analysis indicated that MAPKs and MKKs were divided into four groups as previously reported, grouping them in the same branch as corresponding members. MPK21-2 was the exception and fell into two groups, which may be due to their exon-intron patterns, especially the untranslated regions (UTRs). We also found that differential evolution patterns of MKK10 paralogues from ancient tandem duplicates may have undergone functional divergence. Expression analyses suggested that MAPKs and MKKs likely played different roles in different genetic contexts within various tissues and with abiotic stresses.

CONCLUSION

Our study revealed that UTRs affected the structure and evolution of MPK21-2 genes and the differential evolution of MKK10 paralogues with ancient tandem duplication might have functional divergences. Our findings provide new insights into the functional evolution of genes in closely inbred lines.

Perspective: 50 years of plant chromosome biology

The past 50 years has been the greatest era of plant science discovery, and most of the discoveries have emerged from or been facilitated by our knowledge of plant chromosomes. At last we have descriptive and mechanistic outlines of the information in chromosomes that programs plant life. We had almost no such information 50 years ago when few had isolated DNA from any plant species. The important features of genes have been revealed through whole genome comparative genomics and testing of variants using transgenesis. Progress has been enabled by the development of technologies that had to be invented and then become widely available. Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa) have played extraordinary roles as model species. Unexpected evolutionary dramas were uncovered when learning that chromosomes have to manage constantly the vast numbers of potentially mutagenic families of transposons and other repeated sequences. The chromatin-based transcriptional and epigenetic mechanisms that co-evolved to manage the evolutionary drama as well as gene expression and 3-D nuclear architecture have been elucidated these past 20 years. This perspective traces some of the major developments with which I have become particularly familiar while seeking ways to improve crop plants. I draw some conclusions from this look-back over 50 years during which the scientific community has (i) exposed how chromosomes guard, readout, control, recombine, and transmit information that programs plant species, large and small, weed and crop, and (ii) modified the information in chromosomes for the purposes of genetic, physiological, and developmental analyses and plant improvement.

Genome Wide Identification and Comparative Analysis of the Serpin Gene Family in Brachypodium and Barley

Serpins (serine protease inhibitors) constitute one of the largest and most widely distributed superfamilies of protease inhibitors and have been identified in nearly all organisms. To gain significant insights, a comprehensive in silico analysis of the serpin gene family was carried out in the model plant for temperate grasses Brachypodium distachyon and barley Hordeum vulgare using bioinformatic tools at the genome level for the first time. We identified a total of 27 BdSRPs and 25 HvSRP genes in Brachypodium and barley, respectively, showing an unexpectedly high gene number in these model plants. Gene structure, conserved motifs and phylogenetic comparisons of serpin genes supported the role of duplication events in the expansion and evolution of serpin gene family. Further, purifying selection pressure was found to be a main driving force in the evolution of serpin genes. Genome synteny analysis indicated that BdSRP genes were present in syntenic regions of barley, rice, sorghum and maize, suggesting that they evolved before the divergence of these species from common ancestor. The distinct expression pattern in specific tissues further suggested a specialization of functions during development and in plant defense. These results suggest that the LR serpins (serpins with Leu-Arg residues at P2-P1') identified here can be utilized as candidates for exploitation in disease resistance, pest control and preventing stress-induced cell death. Additionally, serpins were identified that could lead to further research aimed at validating and functionally characterizing the role of potential serpin genes from other plants.

Gradual polyploid genome evolution revealed by pan-genomic analysis of Brachypodium hybridum and its diploid progenitors

Our understanding of polyploid genome evolution is constrained because we cannot know the exact founders of a particular polyploid. To differentiate between founder effects and post polyploidization evolution, we use a pan-genomic approach to study the allotetraploid Brachypodium hybridum and its diploid progenitors. Comparative analysis suggests that most B. hybridum whole gene presence/absence variation is part of the standing variation in its diploid progenitors. Analysis of nuclear single nucleotide variants, plastomes and k-mers associated with retrotransposons reveals two independent origins for B. hybridum, ~1.4 and ~0.14 million years ago. Examination of gene expression in the younger B. hybridum lineage reveals no bias in overall subgenome expression. Our results are consistent with a gradual accumulation of genomic changes after polyploidization and a lack of subgenome expression dominance. Significantly, if we did not use a pan-genomic approach, we would grossly overestimate the number of genomic changes attributable to post polyploidization evolution.

Comparatively Barcoded Chromosomes of Brachypodium Perennials Tell the Story of Their Karyotype Structure and Evolution

The Brachypodium genus is an informative model system for studying grass karyotype organization. Previous studies of a limited number of species and reference chromosomes have not provided a comprehensive picture of the enigmatic phylogenetic relationships in the genus. Comparative chromosome barcoding, which enables the reconstruction of the evolutionary history of individual chromosomes and their segments, allowed us to infer the relationships between putative ancestral karyotypes of extinct species and extant karyotypes of current species. We used over 80 chromosome-specific BAC (bacterial artificial chromosome) clones derived from five reference chromosomes of B. distachyon as probes against the karyotypes of twelve accessions representing five diploid and polyploid Brachypodium perennials. The results showed that descending dysploidy is common in Brachypodium and occurs primarily via nested chromosome fusions. Brachypodiumdistachyon was rejected as a putative ancestor for allotetraploid perennials and B. stacei for B. mexicanum. We propose two alternative models of perennial polyploid evolution involving either the incorporation of a putative x = 5 ancestral karyotype with different descending dysploidy patterns compared to B. distachyon chromosomes or hybridization of two x = 9 ancestors followed by genome doubling and descending dysploidy. Details of the karyotype structure and evolution in several Brachypodium perennials are revealed for the first time.

Detecting Brachypodium distachyon Chromosomes Bd4 and Bd5 in MH- and X-Ray-Induced Micronuclei Using mcFISH

Micronuclei are biomarkers of genotoxic effects and chromosomal instability. They are formed when chromosome fragments or whole chromosomes fail to disjoin into daughter nuclei. We present qualitative and quantitative analyses of the involvement of specific chromosome regions of chromosomes Bd4 and Bd5 in the formation of micronuclei of Brachypodium distachyon root tip cells following maleic hydrazide (MH) treatment and X-radiation. This is visualised by cytomolecular approaches using bacterial artificial chromosome (BAC)-based multicolour fluorescence in situ hybridisation (mcFISH) in combination with 5S and 25S rDNA probes. The results showed that the long arm of submetacentric chromosome Bd4 forms micronuclei at twice the frequency of its short arm, suggesting that the former is more prone to double-strand breaks (DSBs). In contrast, no difference was observed in the frequency of micronuclei derived from the long and short arms of submetacentric chromosome Bd5. Interestingly, the proximal region of the short arm of Bd5 is more prone to DSBs than its distal part. This demonstrates that 5S rDNA and 35S rDNA loci are not "hot spots" for DNA breaks after the application of these mutagens.

Low Lignin Mutants and Reduction of Lignin Content in Grasses for Increased Utilisation of Lignocellulose

Biomass rich in lignocellulose from grasses is a major source for biofuel production and animal feed. However, the presence of lignin in cell walls limits its efficient utilisation such as in its bioconversion to biofuel. Reduction of the lignin content or alteration of its structure in crop plants have been pursued, either by regulating genes encoding enzymes in the lignin biosynthetic pathway using biotechnological techniques or by breeding naturally-occurring low lignin mutant lines. The aim of this review is to provide a summary of these studies, focusing on lignin (monolignol) biosynthesis and composition in grasses and, where possible, the impact on recalcitrance to bioconversion. An overview of transgenic crops of the grass family with regulated gene expression in lignin biosynthesis is presented, including the effect on lignin content and changes in the ratio of p-hydroxyphenyl (H), guaiacyl (G) and syringyl (S) units. Furthermore, a survey is provided of low-lignin mutants in grasses, including cereals in particular, summarising their origin and phenotypic traits together with genetics and the molecular function of the various genes identified.

Homoeolog-specific activation of genes for heat acclimation in the allopolyploid grass Brachypodium hybridum

Background

Allopolyploid plants often show wider environmental tolerances than their ancestors; this is expected to be due to the merger of multiple distinct genomes with a fixed heterozygosity. The complex homoeologous gene expression could have been evolutionarily advantageous for the adaptation of allopolyploid plants. Despite multiple previous studies reporting homoeolog-specific gene expression in allopolyploid species, there are no clear examples of homoeolog-specific function in acclimation to a long-term stress condition.

Results

We found that the allopolyploid grass Brachypodium hybridum and its ancestor Brachypodium stacei show long-term heat stress tolerance, unlike its other ancestor, Brachypodium distachyon. To understand the physiological traits of B. hybridum, we compared the transcriptome of the 3 Brachypodium species grown under normal and heat stress conditions. We found that the expression patterns of approximately 26% and approximately 38% of the homoeolog groups in B. hybridum changed toward nonadditive expression and nonancestral expression, respectively, under normal condition. Moreover, we found that B. distachyon showed similar expression patterns between normal and heat stress conditions, whereas B. hybridum and B. stacei significantly altered their transcriptome in response to heat after 3 days of stress exposure, and homoeologs that were inherited from B. stacei may have contributed to the transcriptional stress response to heat in B. hybridum. After 15 days of heat exposure, B. hybridum and B. stacei maintained transcriptional states similar to those under normal conditions. These results suggest that an earlier response to heat that was specific to homoeologs originating from B. stacei contributed to cellular homeostasis under long-term heat stress in B. hybridum.

Conclusions

Our results provide insights into different regulatory events of the homoeo-transcriptome that are associated with stress acclimation in allopolyploid plants.

Germination and the Early Stages of Seedling Development in Brachypodium distachyon

Successful germination and seedling development are crucial steps in the growth of a new plant. In this study, we investigated the course of the cell cycle during germination in relation to grain hydration in the model grass Brachypodium distachyon (Brachypodium) for the first time. Flow cytometry was performed to monitor the cell cycle progression during germination and to estimate DNA content in embryo tissues. The analyses of whole zygotic embryos revealed that the relative DNA content was 2C, 4C, 8C, and 16C. Endoreplicated nuclei were detected in the scutellum and coleorhiza cells, whereas the rest of the embryo tissues only had nuclei with a 2C and 4C DNA content. This study was accompanied by a spatiotemporal profile analysis of the DNA synthetic activity in the organs of Brachypodium embryos during germination using EdU labelling. Upon imbibition, nuclear DNA replication was initiated in the radicle within 11 h and subsequently spread towards the plumule. The first EdU-labelled prophases were observed after 14 h of imbibition. Analysis of selected genes that are involved in the regulation of the cell cycle, such as those encoding cyclin-dependent kinases and cyclins, demonstrated an increase in their expression profiles.

Chromosome identification and reconstruction of evolutionary rearrangements in Brachypodium distachyon, B. stacei and B. hybridum

BACKGROUND AND AIMS

The Brachypodium genus represents a useful model system to study grass genome organization. Palaeogenomic analyses (e.g. Murat F, Armero A, Pont C, Klopp C, Salse J. 2017. Reconstructing the genome of the most recent common ancestor of flowering plants. Nature Genetics49: 490-496) have identified polyploidization and dysploidy as the prime mechanisms driving the diversity of plant karyotypes and nested chromosome fusions (NCFs) crucial for shaping grass chromosomes. This study compares the karyotype structure and evolution in B. distachyon (genome Bd), B. stacei (genome Bs) and in their putative allotetraploid B. hybridum (genomes BdBs).

METHODS

Brachypodium chromosomes were measured and identified using multicolour fluorescence in situ hybridization (mcFISH). For higher resolution, comparative chromosome barcoding was developed using sets of low-repeat, physically mapped B. distachyon-derived bacterial artificial chromosome (BAC) clones.

KEY RESULTS

All species had rather small chromosomes, and essentially all in the Bs genome were morphometrically indistinguishable. Seven BACs combined with two rDNA-based probes provided unambiguous and reproducible chromosome discrimination. Comparative chromosome barcoding revealed NCFs that contributed to the reduction in the x = 12 chromosome number that has been suggested for the intermediate ancestral grass karyotype. Chromosome Bd3 derives from two NCFs of three ancestral chromosomes (Os2, Os8, Os10). Chromosome Bs6 shows an ancient Os8/Os10 NCF, whilst Bs4 represents Os2 only. Chromosome Bd4 originated from a descending dysploidy that involves two NCFs of Os12, Os9 and Os11. The specific distribution of BACs along Bs9 and Bs5, in both B. stacei and B. hybridum, suggests a Bs genome-specific Robertsonian rearrangement.

CONCLUSIONS

mcFISH-based karyotyping identifies all chromosomes in Brachypodium annuals. Comparative chromosome barcoding reveals rearrangements responsible for the diverse organization of Bd and Bs genomes and provides new data regarding karyotype evolution since the split of the two diploids. The fact that no chromosome rearrangements were observed in B. hybridum compared with the karyotypes of its phylogenetic ancestors suggests prolonged genome stasis after the formation of the allotetraploid.

Reconstructing the origins and the biogeography of species' genomes in the highly reticulate allopolyploid-rich model grass genus Brachypodium using minimum evolution, coalescence and maximum likelihood approaches

The identification of homeologous genomes and the biogeographical analyses of highly reticulate allopolyploid-rich groups face the challenge of incorrectly inferring the genomic origins and the biogeographical patterns of the polyploids from unreliable strictly bifurcating trees. Here we reconstruct a plausible evolutionary scenario of the diverging and merging genomes inherited by the diploid and allopolyploid species and cytotypes of the model grass genus Brachypodium. We have identified the ancestral Brachypodium genomes and inferred the paleogeographical ranges for potential hybridization events that originated its allopolyploid taxa. We also constructed a comprehensive phylogeny of Brachypodium from five nuclear and plastid genes using Species Tree Minimum Evolution allele grafting and Species Network analysis. The divergence ages of the lineages were estimated from a consensus maximum clade credibility tree using fossil calibrations, whereas ages of origin of the diploid and allopolyploid species were inferred from coalescence Bayesian methods. The biogeographical events of the genomes were reconstructed using a stratified Dispersal-Extinction-Colonization model with three temporal windows. Our combined Minimum Evolution-coalescence-Bayesian approach allowed us to infer the origins and the identities of the homeologous genomes of the Brachypodium allopolyploids, matching the expected ploidy levels of the hybrids. To date, the current extant progenitor genomes (species) are only known for B. hybridum. Putative ancestral homeologous genome have been inherited by B. mexicanum, ancestral and recent genomes by B. boissieri, and only recently evolved genomes by B. retusum and the core perennial clade allopolyploids (B. phoenicoides, B. pinnatum 4x, B. rupestre 4x). We dissected the complex spatio-temporal evolution of ancestral and recent genomes and have detected successive splitting, dispersal and merging events for dysploid homeologous genomes in diverse geographical scenarios that have led to the current extant taxa. Our data support Mid-Miocene splits of the Holarctic ancestral genomes that preceded the Late Miocene origins of Brachypodium ancestors of the modern diploid species. Successive divergences of the annual B. stacei and B. distachyon diploid genomes were implied to have occurred in the Mediterranean region during the Late Miocene-Pliocene. By contrast, a profusion of splits, range expansions and different genome mergings were inferred for the perennial diploid genomes in the Mediterranean and Eurasian regions, with sporadic colonizations and further mergings in other continents during the Quaternary. A reliable biogeographical scenario was obtained for the Brachypodium genomes and allopolyploids where homeologous genomes split from their respective diploid counterpart lineages in the same ancestral areas, showing similar or distinct dispersals. By contrast, the allopolyploid taxa remained in the same ancestral ranges after hybridization and genome doubling events. Our approach should have utility in deciphering the genomic composition and the historical biogeography of other allopolyploid-rich organismal groups, which are predominant in eukaryotes.

Methods for Cytogenetic Chromosome Barcoding and Chromosome Painting in Brachypodium distachyon and Its Relative Species

Brachypodium distachyon provides a particularly appealing object for molecular cytogenetic analysis due to its compact genome and low repetitive DNA content, as well as low (x = 5) basic number of chromosomes easily identifiable on the basis of their morphometric features. Some of these features, such as genome compactness, are shared by the other members of the genus, thus making them amenable for comparative cytogenetic mapping. Cytogenetic infrastructure established for B. distachyon was initially based on fluorescence in situ hybridization with various tandemly repeated sequences as probes. The molecular cytogenetic studies advanced greatly with the development of B. distachyon large DNA insert genomic libraries. These resources coupled with the access to the fully sequenced genome of B. distachyon enabled chromosome painting in monocots for the first time. This pioneering work was subsequently extended to other Brachypodium species, allowing insight into grass karyotype evolution. In this protocol we describe the methods of making somatic and meiotic chromosome preparations, probe labeling, FISH with BAC clones, a strategy for chromosome barcoding and chromosome painting in B. distachyon, and comparative chromosome painting in the other Brachypodium species.

Comparative plastome genomics and phylogenomics of Brachypodium: flowering time signatures, introgression and recombination in recently diverged ecotypes

Few pan-genomic studies have been conducted in plants, and none of them have focused on the intraspecific diversity and evolution of their plastid genomes. We address this issue in Brachypodium distachyon and its close relatives B. stacei and B. hybridum, for which a large genomic data set has been compiled. We analyze inter- and intraspecific plastid comparative genomics and phylogenomic relationships within a family-wide framework. Major indel differences were detected between Brachypodium plastomes. Within B. distachyon, we detected two main lineages, a mostly Extremely Delayed Flowering (EDF+) clade and a mostly Spanish (S+) - Turkish (T+) clade, plus nine chloroplast capture and two plastid DNA (ptDNA) introgression and micro-recombination events. Early Oligocene (30.9 million yr ago (Ma)) and Late Miocene (10.1 Ma) divergence times were inferred for the respective stem and crown nodes of Brachypodium and a very recent Mid-Pleistocene (0.9 Ma) time for the B. distachyon split. Flowering time variation is a main factor driving rapid intraspecific divergence in B. distachyon, although it is counterbalanced by repeated introgression between previously isolated lineages. Swapping of plastomes between the three different genomic groups, EDF+, T+, S+, probably resulted from random backcrossing followed by stabilization through selection pressure.

TILLING in Brachypodium distachyon

TILLING is a low-cost screening method that allows for identification of mutations in a gene-of-interest within a range of few base pairs. TILLING can be applied to mutant populations or to plant collections of cultivars, landraces or crop wild relatives (Eco-TILLING). The method is based on the Cel1 enzyme cleavage of mismatches in PCR products amplified with labeled primers. The cleavage can be detected due to the labeled primers by different methods including capillary electrophoresis. Here, we introduce the development of the mutant population BRACHYLIFE and present a Brachypodium TILLING protocol based on fluorescing primers for PCR, enzymatic cleavage, and detection with Applied Biosystems 3130xl Genetic Analyzer.

The interplay between aridity and competition determines colonization ability, exclusion and ecological segregation in the heteroploid Brachypodium distachyon species complex

A higher competitive advantage of polyploid plants compared with their parental diploids is frequently invoked to explain their establishment success, colonization of novel environments and cytotypic ecological segregation, yet there is scarce experimental evidence supporting such hypotheses. Here, we investigated whether differential competitive ability of species of the Brachypodium distachyon (Poaceae) species complex, a model system for genomic, ecological and evolutionary studies of temperate grasses, contributes to explaining their ecological segregation as well as their coexistence in diploid/allotetraploid contact zones. We conducted two field experiments in dry and humid localities to evaluate the tolerance to competition of diploids and allotetraploids in densely occupied environments, and to parameterize models of intra- and intercytotype competition as a mechanism for species exclusion/coexistence. We provide experimental evidence supporting the hypothesis that, under natural field conditions, allotetraploids have superior ecological success compared with one of their parental diploids in terms of both colonizing competitive habitats and intercytotypic competition, with the balance of intra/intercytotype competition favoring polyploid population establishment. These findings, together with previous data on ecogeographic segregation and adaptive response to water stress, suggest that the interplay between aridity and competitive outcome determines the ability to colonize competitive environments, the exclusion of diploids, especially in arid localities, and species geographic segregation.

Environmental isolation explains Iberian genetic diversity in the highly homozygous model grass Brachypodium distachyon

BACKGROUND

Brachypodium distachyon (Poaceae), an annual Mediterranean Aluminum (Al)-sensitive grass, is currently being used as a model species to provide new information on cereals and biofuel crops. The plant has a short life cycle and one of the smallest genomes in the grasses being well suited to experimental manipulation. Its genome has been fully sequenced and several genomic resources are being developed to elucidate key traits and gene functions. A reliable germplasm collection that reflects the natural diversity of this species is therefore needed for all these genomic resources. However, despite being a model plant, we still know very little about its genetic diversity. As a first step to overcome this gap, we used nuclear Simple Sequence Repeats (nSSR) to study the patterns of genetic diversity and population structure of B. distachyon in 14 populations sampled across the Iberian Peninsula (Spain), one of its best known areas.

RESULTS

We found very low levels of genetic diversity, allelic number and heterozygosity in B. distachyon, congruent with a highly selfing system. Our results indicate the existence of at least three genetic clusters providing additional evidence for the existence of a significant genetic structure in the Iberian Peninsula and supporting this geographical area as an important genetic reservoir. Several hotspots of genetic diversity were detected and populations growing on basic soils were significantly more diverse than those growing in acidic soils. A partial Mantel test confirmed a statistically significant Isolation-By-Distance (IBD) among all studied populations, as well as a statistically significant Isolation-By-Environment (IBE) revealing the presence of environmental-driven isolation as one explanation for the genetic patterns found in the Iberian Peninsula.

CONCLUSIONS

The finding of higher genetic diversity in eastern Iberian populations occurring in basic soils suggests that these populations can be better adapted than those occurring in western areas of the Iberian Peninsula where the soils are more acidic and accumulate toxic Al ions. This suggests that the western Iberian acidic soils might prevent the establishment of Al-sensitive B. distachyon populations, potentially causing the existence of more genetically depauperated individuals.

Diversity and association of phenotypic and metabolomic traits in the close model grasses Brachypodium distachyon, B. stacei and B. hybridum

BACKGROUND AND AIMS

Morphological traits in combination with metabolite fingerprinting were used to investigate inter- and intraspecies diversity within the model annual grasses Brachypodium distachyon, Brachypodium stacei and Brachypodium hybridum .

METHODS

Phenotypic variation of 15 morphological characters and 2219 nominal mass ( m / z ) signals generated using flow infusion electrospray ionization-mass spectrometry (FIE-MS) were evaluated in individuals from a total of 174 wild populations and six inbred lines, and 12 lines, of the three species, respectively. Basic statistics and multivariate principal component analysis and discriminant analysis were used to differentiate inter- and intraspecific variability of the two types of variable, and their association was assayed with the rcorr function.

KEY RESULTS

Basic statistics and analysis of variance detected eight phenotypic characters [(stomata) leaf guard cell length, pollen grain length, (plant) height, second leaf width, inflorescence length, number of spikelets per inflorescence, lemma length, awn length] and 434 tentatively annotated metabolite signals that significantly discriminated the three species. Three phenotypic traits (pollen grain length, spikelet length, number of flowers per inflorescence) might be genetically fixed. The three species showed different metabolomic profiles. Discriminant analysis significantly discriminated the three taxa with both morphometric and metabolome traits and the intraspecific phenotypic diversity within B. distachyon and B. stacei . The populations of B. hybridum were considerably less differentiated.

CONCLUSIONS

Highly explanatory metabolite signals together with morphological characters revealed concordant patterns of differentiation of the three taxa. Intraspecific phenotypic diversity was observed between northern and southern Iberian populations of B. distachyon and between eastern Mediterranean/south-western Asian and western Mediterranean populations of B. stacei . Significant association was found for pollen grain length and lemma length and ten and six metabolomic signals, respectively. These results would guide the selection of new germplasm lines of the three model grasses in ongoing genome-wide association studies.

Brachypodium distachyon - A Useful Model in the Qualification of Mutagen-Induced Micronuclei Using Multicolor FISH

Brachypodium distachyon (Brachypodium) is now intensively utilized as a model grass species in various biological studies. Its favorable cytological features create a unique foundation for a convenient system in mutagenesis, thereby potentially enabling the 'hot spots' and 'cold spots' of DNA damage in its genome to be analyzed. The aim of this study was to analyze the involvement of 5S rDNA, 25S rDNA, the Arabidopsis-type (TTTAGGG)n telomeric sequence and the Brachypodium-originated centromeric BAC clone CB33J12 in the micronuclei formation in Brachypodium root tip cells that were subjected to the chemical clastogenic agent maleic hydrazide (MH). To the best of our knowledge, this is the first use of a multicolor fluorescence in situ hybridization (mFISH) with four different DNA probes being used simultaneously to study plant mutagenesis. A quantitative analysis allowed ten types of micronuclei, which were characterized by the presence or absence of specific FISH signal(s), to be distinguished, thus enabling some specific rules governing the composition of the MH-induced micronuclei with the majority of them originating from the terminal regions of chromosomes, to be identified. The application of rDNA sequences as probes showed that 5S rDNA-bearing chromosomes are involved in micronuclei formation more frequently than the 25S rDNA-bearing chromosomes. These findings demonstrate the promising potential of Brachypodium to be a useful model organism to analyze the effects of various genotoxic agents on the plant nuclear genome stability, especially when the complex FISH-based and chromosome-specific approaches such as chromosome barcoding and chromosome painting will be applied in future studies.

Characterization of the LBD gene family in Brachypodium: a phylogenetic and transcriptional study

An unambiguous nomenclature is proposed for the twenty-eight-member LOB domain transcription factor family in Brachypodium . Expression analysis provides unique transcript patterns that are characteristic of a wide range of organs and plant parts. LOB (lateral organ boundaries)-domain proteins define a family of plant-specific transcription factors involved in developmental processes from embryogenesis to seed production. They play a crucial role in shaping the plant architecture through coordinating cell fate at meristem to organ boundaries. Despite their high potential importance, our knowledge of them is limited, especially in the case of monocots. In this study, we characterized LOB domain protein coding genes (LBDs) of Brachypodium distachyon, a model plant for grasses, and present their phylogenetic relationships and an overall spatial expression study. In the Brachypodium genome database, 28 LBDs were found and then classified based on the presence of highly conserved LOB domain motif. Their transcript amounts were measured via quantitative real-time RT-PCR in 37 different plant parts from root tip to generative organs. Comprehensive phylogenetic analysis suggests that there are neither Brachypodium- nor monocot-specific lineages among LBDs, but there are differences in terms of complexity of subclasses between monocots and dicots. Although LBDs in Brachypodium have wide variation of tissue-specific expression and relative transcript levels, overall expression patterns show similarity to their counterparts in other species. The varying transcript profiles we observed support the hypothesis that Brachypodium LBDs have diverse but conserved functions in plant organogenesis.

Recreating Stable Brachypodium hybridum Allotetraploids by Uniting the Divergent Genomes of B. distachyon and B. stacei

Brachypodium hybridum (2n = 30) is a natural allopolyploid with highly divergent sub-genomes derived from two extant diploid species, B. distachyon (2n = 10) and B. stacei (2n = 20) that differ in chromosome evolution and number. We created synthetic B. hybridum allotetraploids by hybridizing various lines of B. distachyon and B. stacei. The initial amphihaploid F1 interspecific hybrids were obtained at low frequencies when B. distachyon was used as the maternal parent (0.15% or 0.245% depending on the line used) and were sterile. No hybrids were obtained from reciprocal crosses or when autotetraploids of the parental species were crossed. Colchicine treatment was used to double the genome of the F1 amphihaploid lines leading to allotetraploids. The genome-doubled F1 plants produced a few S1 (first selfed generation) seeds after self-pollination. S1 plants from one parental combination (Bd3-1×Bsta5) were fertile and gave rise to further generations whereas those of another parental combination (Bd21×ABR114) were sterile, illustrating the importance of the parental lineages crossed. The synthetic allotetraploids were stable and resembled the natural B. hybridum at the phenotypic, cytogenetic and genomic levels. The successful creation of synthetic B. hybridum offers the possibility to study changes in genome structure and regulation at the earliest stages of allopolyploid formation in comparison with the parental species and natural B. hybridum.

Spatial Distribution of Selected Chemical Cell Wall Components in the Embryogenic Callus of Brachypodium distachyon

Brachypodium distachyon L. Beauv. (Brachypodium) is a species that has become an excellent model system for gaining a better understanding of various areas of grass biology and improving plant breeding. Although there are some studies of an in vitro Brachypodium culture including somatic embryogenesis, detailed knowledge of the composition of the main cell wall components in the embryogenic callus in this species is missing. Therefore, using the immunocytochemical approach, we targeted 17 different antigens of which five were against the arabinogalactan proteins (AGP), three were against extensins, six recognised pectic epitopes and two recognised hemicelluloses. These studies were complemented by histological and scanning electron microscopy (SEM) analyses. We revealed that the characteristic cell wall components of Brachypodium embryogenic calli are AGP epitopes that are recognised by the JIM16 and LM2 antibodies, an extensin epitope that is recognised by the JIM11 antibody and a pectic epitopes that is recognised by the LM6 antibody. Furthermore, we demonstrated that AGPs and pectins are the components of the extracellular matrix network in Brachypodium embryogenic culture. Additionally, SEM analysis demonstrated the presence of an extracellular matrix on the surface of the calli cells. In conclusion, the chemical compositions of the cell walls and ECMSN of Brachypodium callus show spatial differences that correlate with the embryogenic character of the cells. Thus, the distribution of pectins, AGPs and hemicelluloses can be used as molecular markers of embryogenic cells. The presented data extends the knowledge about the chemical composition of the embryogenic callus cells of Brachypodium.

Expansion and stress responses of AP2/EREBP superfamily in Brachypodium distachyon

APETALA2/ethylene-responsive element binding protein (AP2/EREBP) transcription factors constitute one of the largest and most conserved gene families in plant, and play essential roles in growth, development and stress response. Except a few members, the AP2/EREBP family has not been characterized in Brachypodium distachyon, a model plant of Poaceae. We performed a genome-wide study of this family in B. distachyon by phylogenetic analyses, transactivation assays and transcript profiling. A total of 149 AP2/EREBP genes were identified and divided into four subfamilies, i.e., ERF (ethylene responsive factor), DREB (dehydration responsive element binding gene), RAV (related to ABI3/VP) and AP2. Tandem duplication was a major force in expanding B. distachyon AP2/EREBP (BdAP2/EREBP) family. Despite a significant expansion, genomic organizations of BdAP2/EREBPs were monotonous as the majority of them, except those of AP2 subfamily, had no intron. An analysis of transcription activities of several closely related and duplicated BdDREB genes showed their functional divergence and redundancy in evolution. The expression of BdAP2/EREBPs in different tissues and the expression of DREB/ERF subfamilies in B. distachyon, wheat and rice under abiotic stresses were investigated by next-generation sequencing and microarray profiling. Our results are valuable for further function analysis of stress tolerant AP2/EREBP genes in B. distachyon.

Molecular Cloning and Functional Characterization of Two Brachypodium distachyon UBC13 Genes Whose Products Promote K63-Linked Polyubiquitination

Living organisms are constantly subject to DNA damage from environmental sources. Due to the sessile nature of plants, UV irradiation is a major genotoxic agent and imposes a significant threat on plant survival, genome stability and crop yield. In addition, other environmental chemicals can also influence the stability of the plant genome. Eukaryotic organisms have evolved a mechanism to cope with replication-blocking lesions and stabilize the genome. This mechanism is known as error-free DNA damage tolerance, and is mediated by K63-linked PCNA polyubiquitination. Genes related to K63-linked polyubiquitination have been isolated recently from model plants like Arabidopsis and rice, but we are unaware of such reports on the crop model Brachypodium distachyon. Here, we report the identification and functional characterization of two B. distachyon UBC13 genes. Both Ubc13s form heterodimers with Uevs from other species, which are capable of catalyzing K63 polyubiquitination in vitro. Both genes can functionally rescue the yeast ubc13 null mutant from killing by DNA-damaging agents. These results suggest that Ubc13-Uev-promoted K63-linked polyubiquitination is highly conserved in eukaryotes including B. distachyon. Consistent with recent findings that K63-linked polyubiquitination is involved in several developmental and stress-responsive pathways, the expression of BdUbc13s appears to be constitutive and is regulated by abnormal temperatures.

Three Brachypodium distachyon Uev1s Promote Ubc13-Mediated Lys63-Linked Polyubiquitination and Confer Different Functions

In this study, we report the identification and functional characterization of three Brachypodium distachyon UEV genes. All three BdUev1s form heterodimers with BdUbc13s, which are capable of catalyzing Lys63-linked polyubiquitination in vitro. The three BdUEV1 genes are also able to functionally complement the budding yeast mms2 mutant defective in DNA-damage tolerance. BdUev1A differs from the other two BdUev1s in that it contains an 18-amino acid tail, which appears to compromise its function in yeast, as deletion of this tail restores full function. BdUev1A is excluded from the nucleus, whereas BdUev1B, BdUev1C and the C-terminally truncated BdUev1A are mainly found in the nucleus. These and the BdUEV1 gene expression analysis allow us to speculate that although all three BdUev1s function by promoting Lys63-linked polyubiquitination, BdUev1B and BdUev1C are involved in DNA-damage response and possibly other nuclear functions, while BdUev1A is required for non-nuclear function(s).

Natural variation, differentiation, and genetic trade-offs of ecophysiological traits in response to water limitation in Brachypodium distachyon and its descendent allotetraploid B. hybridum (Poaceae)

Differences in tolerance to water stress may underlie ecological divergence of closely related ploidy lineages. However, the mechanistic basis of physiological variation governing ecogeographical cytotype segregation is not well understood. Here, using Brachypodium distachyon and its derived allotetraploid B. hybridum as model, we test the hypothesis that, for heteroploid annuals, ecological divergence of polyploids in drier environments is based on trait differentiation enabling drought escape. We demonstrate that under water limitation allotetraploids maintain higher photosynthesis and stomatal conductance and show earlier flowering than diploids, concordant with a drought-escape strategy to cope with water stress. Increased heterozygosity and greater genetic variability and plasticity of polyploids could confer a superior adaptive capability. Consistent with these predictions, we document (1) greater standing within-population genetic variation in water-use efficiency (WUE) and flowering time in allotetraploids, and (2) the existence of (nonlinear) environmental clines in physiology across allotetraploid populations. Increased gas exchange and diminished WUE occurred at the driest end of the gradient, consistent with a drought-escape strategy. Finally, we found that allotetraploids showed weaker genetic correlations than diploids congruous with the expectation of relaxed pleiotropic constraints in polyploids. Our results suggest evolutionary divergence of ecophysiological traits in each ploidy lineage.

Influence of rhizobacterial volatiles on the root system architecture and the production and allocation of biomass in the model grass Brachypodium distachyon (L.) P. Beauv

BACKGROUND

Plant growth-promoting rhizobacteria are increasingly being seen as a way of complementing conventional inputs in agricultural systems. The effects on their host plants are diverse and include volatile-mediated growth enhancement. This study sought to assess the effects of bacterial volatiles on the biomass production and root system architecture of the model grass Brachypodium distachyon (L.) Beauv.

RESULTS

An in vitro experiment allowing plant-bacteria interaction throughout the gaseous phase without any physical contact was used to screen 19 bacterial strains for their growth-promotion ability over a 10-day co-cultivation period. Five groups of bacteria were defined and characterised based on their combined influence on biomass production and root system architecture. The observed effects ranged from unchanged to greatly increased biomass production coupled with increased root length and branching. Primary root length was increased only by the volatile compounds emitted by Enterobacter cloacae JM22 and Bacillus pumilus T4. Overall, the most significant results were obtained with Bacillus subtilis GB03, which induced an 81 % increase in total biomass, as well as enhancing total root length, total secondary root length and total adventitious root length by 88.5, 201.5 and 474.5 %, respectively.

CONCLUSIONS

This study is the first report on bacterial volatile-mediated growth promotion of a grass plant. Contrasting modulations of biomass production coupled with changes in root system architecture were observed. Most of the strains that increased total plant biomass also modulated adventitious root growth. Under our screening conditions, total biomass production was strongly correlated with the length and branching of the root system components, except for primary root length. An analysis of the emission kinetics of the bacterial volatile compounds is being undertaken and should lead to the identification of the compounds responsible for the observed growth-promotion effects. Within the context of the inherent characteristics of our in vitro system, this paper identifies the next critical experimental steps and discusses them from both a fundamental and an applied perspective.

Genomic architecture and functional relationships of intronless, constitutively- and alternatively-spliced genes in Brachypodium distachyon

Splicing and alternative splicing (AS) are widespread co- and post-transcriptional regulatory processes in plants. Recently, we characterized genome-wide AS landscapes and virus-induced AS patterns in Brachypodium distachyon (Brachypodium), a C3 model grass. Brachypodium plants infected with Panicum mosaic virus (PMV) alone or in mixed infections with its satellite virus (SPMV) were used for high-throughput, paired-end RNA sequencing. Here, using gene attributes of ∼5,655 intronless genes, ∼13,302 constitutively spliced, and ∼7,564 alternatively spliced genes, we analyzed the influence of genomic features on splicing incidence and AS frequency. In Brachypodium, gene length, coding sequence length, and exon and intron number were positively correlated to splicing incidence and AS frequency. In contrast, exon length and the percentage composition of GC (%GC) content were inversely correlated with splicing incidence and AS frequency. Although gene expression status had little correlation with splicing occurrence per se, it negatively correlated to AS frequency: i.e., genes with ≥5 alternatively spliced transcripts were significantly less expressed compared to genes encoding <5 alternative transcripts. Further gene set enrichment analysis uncovered unique functional relationships among nonspliced, constitutively spliced and alternatively spliced genes.