Evolution and taxonomic split of the model grass Brachypodium distachyon

Background and Aims

Brachypodium distachyon is being widely investigated across the world as a model plant for temperate cereals. This annual plant has three cytotypes (2n =  10, 20, 30) that are still regarded as part of a single species. Here, a multidisciplinary study has been conducted on a representative sampling of the three cytotypes to investigate their evolutionary relationships and origins, and to elucidate if they represent separate species.

Methods

Statistical analyses of 15 selected phenotypic traits were conducted in individuals from 36 lines or populations. Cytogenetic analyses were performed through flow cytometry, fluorescence in situ hybridization (FISH) with genomic (GISH) and multiple DNA sequences as probes, and comparative chromosome painting (CCP). Phylogenetic analyses were based on two plastid (ndhF, trnLF) and five nuclear (ITS, ETS, CAL, DGAT, GI) genes from different Brachypodium lineages, whose divergence times and evolutionary rates were estimated.

Key Results

The phenotypic analyses detected significant differences between the three cytotypes and demonstrated stability of characters in natural populations. Genome size estimations, GISH, FISH and CCP confirmed that the 2n = 10 and 2n = 20 cytotypes represent two different diploid taxa, whereas the 2n = 30 cytotype represents the allotetraploid derived from them. Phylogenetic analysis demonstrated that the 2n = 20 and 2n = 10 cytotypes emerged from two independent lineages that were, respectively, the maternal and paternal genome donors of the 2n = 30 cytotype. The 2n = 20 lineage was older and mutated significantly faster than the 2n = 10 lineage and all the core perennial Brachypodium species.

Conclusions

The substantial phenotypic, cytogenetic and molecular differences detected among the three B. distachyon sensu lato cytotypes are indicative of major speciation processes within this complex that allow their taxonomic separation into three distinct species. We have kept the name B. distachyon for the 2n = 10 cytotype and have described two novel species as B. stacei and B. hybridum for, respectively, the 2n = 20 and 2n = 30 cytotypes.     

Comparative cytogenetic analysis of the genomes of the model grass Brachypodium distachyon and its close relatives

Background and Aims

Brachypodium is a small genus of temperate grasses that comprises 12–15 species. Brachypodium distachyon is now well established as a model species for temperate cereals and forage grasses. In contrast to B. distachyon, other members of the genus have been poorly investigated at the chromosome level or not at all.

Methods

Twenty accessions comprising six species and two subspecies of Brachypodium were analysed cytogenetically. Measurements of nuclear genome size were made by flow cytometry. Chromosomal localization of 18–5·8–25S rDNA and 5S rDNA loci was performed by dual-colour fluorescence in situ hybridization (FISH) on enzymatically digested root-tip meristematic cells. For comparative phylogenetic analyses genomic in situ hybridization (GISH) applied to somatic chromosome preparations was used.

Key Results

All Brachypodium species examined have rather small genomes and chromosomes. Their chromosome numbers and genome sizes vary from 2n = 10 and 0·631 pg/2C in B. distachyon to 2n = 38 and 2·57 pg/2C in B. retusum, respectively. Genotypes with 18 and 28 chromosomes were found among B. pinnatum accessions. GISH analysis revealed that B. pinnatum with 28 chromosomes is most likely an interspecific hybrid between B. distachyon (2n = 10) and B. pinnatum (2n = 18). Two other species, B. phoenicoides and B. retusum, are also allopolyploids and B. distachyon or a close relative seems to be one of their putative ancestral species. In chromosomes of all species examined the 45S rDNA loci are distally distributed whereas loci for 5S rDNA are pericentromeric.

Conclusions

The increasing significance of B. distachyon as a model grass emphasizes the need to understand the evolutionary relationships in the genus Brachypodium and to ensure consistency in the biological nomenclature of its species. Modern molecular cytogenetic techniques such as FISH and GISH are suitable for comparative phylogenetic analyses and may provide informative chromosome- and/or genome-specific landmarks.     

Brachypodium distachyon as a new model system for understanding iron homeostasis in grasses: phylogenetic and expression analysis of Yellow Stripe-Like (YSL) transporters

Background and Aims

Brachypodium distachyon is a temperate grass with a small stature, rapid life cycle and completely sequenced genome that has great promise as a model system to study grass-specific traits for crop improvement. Under iron (Fe)-deficient conditions, grasses synthesize and secrete Fe(III)-chelating agents called phytosiderophores (PS). In Zea mays, Yellow Stripe1 (ZmYS1) is the transporter responsible for the uptake of Fe(III)–PS complexes from the soil. Some members of the family of related proteins called Yellow Stripe-Like (YSL) have roles in internal Fe translocation of plants, while the function of other members remains uninvestigated. The aim of this study is to establish brachypodium as a model system to study Fe homeostasis in grasses, identify YSL proteins in brachypodium and maize, and analyse their expression profiles in brachypodium in response to Fe deficiency.

Methods

The YSL family of proteins in brachypodium and maize were identified based on sequence similarity to ZmYS1. Expression patterns of the brachypodium YSL genes (BdYSL genes) were determined by quantitative RT–PCR under Fe-deficient and Fe-sufficient conditions. The types of PS secreted, and secretion pattern of PS in brachypodium were analysed by high-performance liquid chromatography.

Key Results

Eighteen YSL family members in maize and 19 members in brachypodium were identified. Phylogenetic analysis revealed that some YSLs group into a grass-specific clade. The Fe status of the plant can regulate expression of brachypodium YSL genes in both shoots and roots. 3-Hydroxy-2′-deoxymugineic acid (HDMA) is the dominant type of PS secreted by brachypodium, and its secretion is diurnally regulated.

Conclusions

PS secretion by brachypodium parallels that of related crop species such as barley and wheat. A single grass species-specific YSL clade is present, and expression of the BdYSL members of this clade could not be detected in shoots or roots, suggesting grass-specific functions in reproductive tissues. Finally, the Fe-responsive expression profiles of several YSLs suggest roles in Fe homeostasis.     

Genome-wide exploration of the molecular evolution and regulatory network of mitogen-activated protein kinase cascades upon multiple stresses in Brachypodium distachyon

Background

Brachypodium distachyon is emerging as a widely recognized model plant that has very close relations with several economically important Poaceae species. MAPK cascade is known to be an evolutionarily conserved signaling module involved in multiple stresses. Although the gene sequences of MAPK and MAPKK family have been fully identified in B. distachyon, the information related to the upstream MAPKKK gene family especially the regulatory network among MAPKs, MAPKKs and MAPKKKs upon multiple stresses remains to be understood.

Results

In this study, we have identified MAPKKKs which belong to the biggest gene family of MAPK cascade kinases. We have systematically investigated the evolution of whole MAPK cascade kinase gene family in terms of gene structures, protein structural organization, chromosomal localization, orthologs construction and gene duplication analysis. Our results showed that most BdMAPK cascade kinases were located at the low-CpG-density region, and the clustered members in each group shared similar structures of the genes and proteins. Synteny analysis showed that 62 or 21 pairs of duplicated orthologs were present between B. distachyon and Oryza sativa, or between B. distachyon and Arabidopsis thaliana respectively. Gene expression data revealed that BdMAPK cascade kinases were rapidly regulated by stresses and phytohormones. Importantly, we have constructed a regulation network based on co-expression patterns of the expression profiles upon multiple stresses performed in this study.

Conclusions

BdMAPK cascade kinases were involved in the signaling pathways of multiple stresses in B. distachyon. The network of co-expression regulation showed the most of duplicated BdMAPK cascade kinase gene orthologs demonstrated their convergent function, whereas few of them developed divergent function in the evolutionary process. The molecular evolution analysis of identified MAPK family genes and the constructed MAPK cascade regulation network under multiple stresses provide valuable information for further investigation of the functions of BdMAPK cascade kinase genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1452-1) contains supplementary material, which is available to authorized users.     

Vernalization response of Phleum pratense and its relationships to stem lignification and floral transition

Background

Timothy is a long-day grass species well adapted for cultivation in northern latitudes. It produces elongating tillers not only in spring growth but also later in summer. As the quantity and quality of harvested biomass is dictated by canopy architecture and the proportion of stem-forming flowering tillers, the regulation of flowering is of great interest in forage grass production.

Methods

Canopy architecture, stem morphology and freezing tolerance of vernalized timothy were investigated in greenhouse and field experiments. The molecular control of development was examined by analysing the relationship between apex development and expression of timothy homologues of the floral inducer VRN1 and repressor VRN2.

Key Results

True stem formation and lignification of the sclerenchyma ring occur in both vernalized and regrowing stems irrespective of the developmental stage of the apex. The stems had, however, divergent morphology. Vernalization enhanced flowering, and the expression of the VRN1 homologue was elevated when the apex had passed into the reproductive stage. High VRN1 homologue expression was not associated with reduction in freezing tolerance and the expression coincided with increased levels of the floral repressor VRN2 homologue. Field experiments supported the observed linkage between the upregulation of the VRN1 homologue and the transition to the reproductive stage in vernalized tillers. The upregulation of putative VRN1 or VRN2 genes was restricted to vernalized tillers in the spring yield and, thus, not detected in non-vernalized tillers of the second yield; so-called regrowth.

Conclusions

The formation of a lignified sclerenchyma ring that efficiently reduces the digestibility of the stem was not related to apex development but rather to a requirement for mechanical support. The observed good freezing tolerance of reproductive timothy tillers could be one important adaptation mechanism ensuring high yields in northern conditions. Both VRN1 and VRN2 homologues required a vernalization signal for expression so the development of yield-forming tillers in regrowth was regulated independently of the studied genes.     

Effects of acclimation and diapause on the thermal tolerance of the two-spotted spider mite Tetranychus urticae

The two-spotted spider mite, Tetranychus urticae, is a worldwide pest species that overwinters as diapausing females. Cold hardening is presumed to start during diapause development to ensure the successful overwintering of this species. To address this hypothesis, we compared cold tolerance between non-diapausing and diapausing females. We measured supercooling point (SCP) and survival to acute cold stress by exposing the mites at a range of sub-zero temperatures (from −4 to −28 °C for 2 h). The mean SCPs of non-diapausing and diapausing females were −19.6±0.5 and −24.7±0.3 °C respectively, and freezing killed the mites. Diapausing females were significantly more cold tolerant than non-diapausing ones, with LT₅₀ of −19.7 and −13.3 °C, respectively. Further, we also examined the effects of cold acclimation (10 d at 0 or 5 °C) in non-diapausing and diapausing females. Our findings indicated that diapause decreased SCP significantly, while cold acclimation had no effect on the SCP except for non-diapausing females that were acclimated at 5 °C. Acclimation at 5 °C enhanced survival to acute cold stress in diapausing and non-diapausing females, with LT₅₀ of −22.0 and −17.1 °C, respectively. Altogether, our results indicate that T. urticae is a chill tolerant species, and that diapause and cold acclimation elevate cold hardiness in this species.     

Differences in leaf proteome response to cold acclimation between Lolium perenne plants with distinct levels of frost tolerance

Perennial ryegrass (Lolium perenne) is a high quality forage and turf grass mainly due to its excellent nutritive values and rapid establishment rate. However, this species has limited ability to perform in harsh winter climates. Though winter hardiness is a complex trait, it is commonly agreed that frost tolerance (FT) is its main component. Species growing in temperate regions can acquire FT through exposure to low, non-lethal temperatures, a phenomenon known as cold acclimation (CA). The research on molecular basis of FT has been performed on the model plants, but they are not well adapted to extreme winter climates. Thus, the mechanisms of cell response to low temperature in winter crops and agronomically important perennial grasses have yet to be revealed. Here, two L. perenne plants with contrasting levels of FT, high frost tolerant (HFT) and low frost tolerant (LFT) plants, were selected for comparative proteomic research. The work focused on analyses of leaf protein accumulation before and after 2, 8, 26 h, and 3, 5, 7, 14 and 21 days of CA, using a high-throughput two-dimensional electrophoresis, and on the identification of proteins which were accumulated differentially between the selected plants by the application of mass spectrometry (MS). Analyses of 580 protein profiles revealed a total of 42 (7.2%) spots that showed at a minimum of 1.5-fold differences in protein abundance, at a minimum of at one time point of CA between HFT and LFT genotypes. It was shown that significant differences in profiles of protein accumulation between the analyzed plants appeared most often on the 5th (18 proteins) and the 7th (19 proteins) day of CA. The proteins derived from 35 (83.3%) spots were successfully identified by the use of MS and chloroplast proteins were shown to be the major group selected as differentially accumulated during CA. The functions of the identified proteins and their probable influence on the level of FT in L. perenne are discussed.