Comparative analyses reveal potential uses of Brachypodium distachyonas a model for cold stress responses in temperate grasses
© Li et al.; licensee BioMed Central Ltd. 2012
Received: 23 December 2011
Accepted: 27 April 2012
Published: 8 May 2012
Little is known about the potential of Brachypodium distachyon as a model for low temperature stress responses in Pooideae. The ice recrystallization inhibition protein (IRIP) genes, fructosyltransferase (FST) genes, and many C-repeat binding factor (CBF) genes are Pooideae specific and important in low temperature responses. Here we used comparative analyses to study conservation and evolution of these gene families in B. distachyon to better understand its potential as a model species for agriculturally important temperate grasses.
Brachypodium distachyon contains cold responsive IRIP genes which have evolved through Brachypodium specific gene family expansions. A large cold responsive CBF3 subfamily was identified in B. distachyon, while CBF4 homologs are absent from the genome. No B. distachyon FST gene homologs encode typical core Pooideae FST-motifs and low temperature induced fructan accumulation was dramatically different in B. distachyon compared to core Pooideae species.
We conclude that B. distachyon can serve as an interesting model for specific molecular mechanisms involved in low temperature responses in core Pooideae species. However, the evolutionary history of key genes involved in low temperature responses has been different in Brachypodium and core Pooideae species. These differences limit the use of B. distachyon as a model for holistic studies relevant for agricultural core Pooideae species.
KeywordsBrachypodium distachyon Cold climate adaptation Ice recrystallization inhibition protein Gene expression Fructosyltransferase C-repeat binding factor Gene family evolution
Brachypodium distachyon became the first Pooideae grass species to have its genome fully sequenced . The Brachypodium genus is a phylogenetic sister group to the Triticeae (cereals) and Poeae (forage grasses) tribes, which provided compelling rationale for sequencing the B. distachyon genome to develop a model more suitable for temperate grasses than rice (Oryza sativa). B. distachyon possesses features typical of a model plant ; it is of relatively short height (15–20 cm), there are inbred lines with an annual and rapid life cycle (eight to twelve weeks), it’s genome is one of the smallest among diploid grass genomes (about 300 Mbp) , and it can be genetically transformed via Agrobacterium-mediated transformation . Different ecotypes exhibit a range of adaptations to environments which also are important challenges faced in agricultural production systems, for example differences in flowering time, vernalization requirements , and disease resistance . Altogether, these features make B. distachyon a suitable model plant for studying agronomic traits in Pooideae grasses.
A characteristic feature of species in the Pooideae sub-family is their adaptation to temperate ecosystems, which is reflected in the global distribution of Pooideae grasses  (Figure 1). The most recent common ancestor of Pooideae grasses was adapted to tropical or sub-tropical climates [7, 8]. Subsequent radiation of Pooideae into cooler environments is thought to be associated with evolution of mechanisms involved in low temperature stress . Hence, the adaptation of the Pooideae to cooler climates makes this group an ideal model system for studying adaptive evolution in plants . Nonetheless, large intraspecific variation in tolerance to cold and freezing stress exists within Pooideae. Some Pooideae species (e.g. Phleum pratense) can tolerate extreme winter climates and has a species range which includes sub-arctic regions (Figure 1), while B. distachyon on the other hand is not adapted to extreme winter climates, which is reflected by the middle-eastern and Mediterranean geographical distribution (Figure 1).
Three Pooideae-specific genetic features involved in low temperature stress tolerance have been well described; (1) ice-recrystallization inhibition proteins (IRIPs), (2) fructosyltransferase genes (FSTs), and (3) expansions within the C-repeat binding factors (CBF) family. IRIPs are thought to minimize cell damage during ice formation by restricting ice crystal growth in the apoplastic space [14–16], a process known as ice crystal recrystallization. The ancestral IRIP gene in Pooideae is thought to have evolved from an LLR-protein kinase  by gaining an ice binding domain  through an expansion of a repeat motif (NxVxG/NxVxxG) . In core Pooideae species extensive gene duplications have given rise to large IRIP gene families [18, 19]. FSTs convert sucrose molecules into fructan sugars [20, 21] and low temperature stress induces FST gene expression and fructan accumulation in core Pooideae species [22–24]. Both the introduction of FSTs into plants lacking endogenous FST pathways and over-expression of endogenous FSTs in core Pooideae species has shown to increase freezing tolerance [25, 26]. CBFs are transcription factors that regulate suites of genes during drought and low temperature stress [27–29]. Two CBF subfamilies, CBF3 and CBF4, have undergone lineage specific duplications in the Pooideae  and the members of these CBF3/4-subfamilies are thought to play roles in Pooideae freezing stress [31–33].
Even though IRIP, FST, and Pooideae-specific CBF3/4 genes have been studied in great detail in core Pooideae species, a systematic study of homologs of these genes in B. distachyon is lacking. In this study we ask the question; to what extent are IRIP, FST, and Pooideae specific CBF genes conserved between B. distachyon and agriculturally important core Pooideae species? We answer this question by employing a suite of methods including comparative genomics, gene expression analyses, and characterization of carbohydrate metabolism. Our aims were to (1) assess the use of B. distachyon as a model to study mechanisms of low temperature stress responses in core Pooideae species, and (2) improve our understanding of the evolution of cold stress response in the Pooideae lineage.
Plant material, growth conditions and tissue sampling for gene expression studies
Four diploid inbred B. distachyon lines, Bd3-1, Bd21-1, Bd1-1, and Bd29-1 were used to characterize cold induced IRIP gene expression by quantitative real-time polymerase chain reaction (qRT-PCR). The seeds were kindly provided by Dr. David Garvin, University of Minnesota, USA. Bd3-1 and Bd21-1 originate from Iraq and are spring genotypes that do not require vernalization to induce flowering. Bd1-1 and Bd29-1 originate from northern Turkey and the Ukraine, respectively, and are winter genotypes which require long vernalization periods (six and 12 weeks respectively) to flower (http://www.ars.usda.gov/SP2UserFiles/person/1931/GarvinLabCoreBrachypodiumdistachyonLineSet(2).pdf) [34, 35]. For the microarray gene expression experiments only Bd21-1 was used.
In the qRT-PCR experiment 7 weeks old plants were used which had been established from seeds using the following growth conditions: 20/16°C day/night temperature and 16 h photoperiod with a photon flux density of 150 μmol m-2 s-1. Half of the plants were cold acclimated (CA) following the procedure outlined in Alm et al.  except that pre-acclimation was done at 12 h photoperiod, 12/6°C day/night temperature for 1 week. Plants were kept in CA conditions at 14 h photoperiod at a constant temperature of 1°C. Leaf tissues were collected from control plants (non-acclimated, NA) at the start of the pre-acclimation period and from CA plants at 4 h, 1 day and 10 days after the start of the cold acclimation period. To avoid experimental bias introduced by diurnal clock regulation of gene expression all samples were collected at the same hour of the day (in the morning) for each time point. Tissue for RNA extraction was sampled from different plants of each genotype at each time point. All sampled leaf tissues were frozen immediately in liquid nitrogen and stored at −80°C.
For the microarray experiment plants were grown in 16 hours photoperiod in a controlled growth room. The temperature was 23°C and the photon flux density was 200 μmol m-2 s-1. Cold experiments were conducted on three-week-old plants in a walk-in cold room at 4°C with a photon flux density of 200 μmol m-2 s-1. Control plants remained in the environmentally controlled growth room at 23°C. Experimental treatment began two hours post-dawn (10:00 am). Leaves and stems (total above ground tissues) from individual plants were collected at 1, 2, 5, 10, and 24 hours after experiments were initiated.
Identification of Brachypodium distachyonIRIP-homologs and design of paralog specific primers
Brachypodium distachyon IRIP homologs were identified through web-based blast search (http://www.Brachypodium.org) using IRIP genes from Lolium perenne (AY968588; EU680848; EU680850; EU680851) as queries. Multiple alignments of translated IRIP genes were made with default settings on the MAFFT web server  to verify that B. distachyon IRIP genes contained the typical NxVxG/NxVxxG-repeat ice binding domain. IRIP-paralog specific primers were designed using primer3  such that there were mismatches between IRIP-paralogs in the 5’ end of at least one primer for each primer-pair. IRIP paralog specificity was verified by cloning the PCR-product of the paralog specific primers pairs using a TOPO TA cloning kit (Invitrogen) and subsequent sequencing of five to ten clones. Final TaqMan MGB probes and primer sets for the quantitative reverse transcriptase PCR analyses (qRT-PCR) of IRIP genes were designed using Primer Express Software (Applied Biosystems).
Gene expression analysis by qRT-PCR
Total RNA was isolated with the RNeasy plant mini kit (Qiagen) and the RNA extraction was performed as described in the manufacturer`s protocol using 100 mg frozen tissue powder (ground with mortar in liquid nitrogen). DNase digestion was used to eliminate genomic DNA contamination. RNA quality was controlled with an Agilent 2100 Bioanalyzer (Aglilent Technologies) and RNA quantity measured on a Nanodrop ND-1000 UV–vis Spectrophotometer (Nanodrop Technologies). 2.5 μg of total RNA was reversed transcribed using SuperScriptVilo cDNA synthesis kit (Invitrogen). For qRT-PCR we used the EXPRESS two-step qRT-PCR universal kit (Invitrogen) with the superscript VILO cDNA synthesis kit (Invitrogen).
Two μl cDNA in a total reaction volume of 20 μl was used for each qRT-PCR reaction. Primers were used at a concentration of 0.5 μM and TaqMan probes at a concentration of 0.2 μM. Final ROX Reference Dye was 0.05 μM. Transcript levels were analysed using a ABI7500 real-time PCR machine (Applied Biosystems) with Fast Cycling Program; 95°C for 20 s and 40 cycles of 95°C for 3 s, and 60°C for 30 s. Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) was used as reference gene. Standard curves were made to control that primer and probe pairs had efficiency close to 100%. B. distachyon IRIP genes transcript levels were calculated relative to GADPH gene transcript levels using the comparative threshold cycle method (ΔCt method). Three biological replicates (leaf samples from three different plants) were used to estimate expression levels. Mean and standard deviation of 2-ΔΔCt was calculated for comparison of relative expression levels in CA compared to NA samples. A t-test was used to test if ΔCt values of CA samples were significantly different from NA samples (i.e. cold induced gene expression). Three no-template controls for each qRT-PCR plate per gene were performed to control for primer-dimer formation and DNA contamination.
Expression analysis by microarray
Leaf tissues were ground in liquid nitrogen and total cellular RNA was extracted using RNA Plant reagent (Invitrogen) and RNeasy kits (Qiagen) and treated with RNase-free DNase as previously described . RNA integrity was evaluated using an Agilent Bioanalyzer. Labeled target cDNA was prepared from 125 ng total RNA samples using the NuGen Applause WT-Amp PlusST RNA amplification system Kit protocols (Cat# 5510–24) and Encore Biotin module V2 (Cat# 4200–12). Approximately 4.55 μg fragmented cDNA from each sample was hybridized for 18 hours to an Affymetrix Brachypodium Genome Array (BradiAR1b520742). Hybridization was performed using GeneChip® Fluidics Station 450. Arrays were scanned using GeneChip® Scanner 3000 with autoloader at 570 nm and quality-controlled according to the standard Affymetrix protocols (Affymetrix GeneChip® Expression Analysis Technical Manual, 701021 Rev. 5; http://www.affymetrix.com) at the Oregon State University Center for Genome Research and Bioinformatics, Central Service Laboratory (detailed protocols are available at http://www.cgrb.oregonstate.edu/). Image processing and data extraction were performed using AGCC software version 3.0. The Affymetrix eukaryotic hybridization control kit and Poly-A RNA control kit were used to ensure efficiency of hybridization and cDNA amplification. All cDNAs from cold stress treatments and control samples were synthesized at the same time and microarray hybridizations were conducted simultaneously. Each array image was visually screened to discount for signal artifacts, scratches or debris.
Probe level normalization was done with Robust Multi-array Analysis (RMA) utilizing the Affymetrix Power Tools (APT) software package (http://www.affymetrix.com/partners_programs/programs/developer/tools/powertools.affx; ). Probe set summarization and expression estimates for each gene were conducted using the apt-probeset-summarize (1.14.3) program from Affymetrix. Data manipulations were performed using Perl scripts to calculate fold change between normalized treatment and control probe set values.
IRIP, FST, and CBF gene families contain recently derived paralogs (i.e. having few substitutions), hence all phylogenetic analyses were carried out on the nucleotide level. Sets of coding sequences (CDS) of IRIP, FST and CBF3/CBF4 genes from core Pooideae species were assembled from different sources. For the IRIP phylogeny, a representative collection of CDS from both Triticeae and Poeae tribe species were downloaded from NCBI and merged with the identified B. distachyon IRIP homologs (see text above). To identify CBF3/4 and FST homologs in B. distachyon, the CDS (v1.2) annotation was downloaded from (http://ftp.brachypodium.org/files/Annotation/) and local blast searches were performed using CBF3/4 genes identified in Skinner et al.  as queries. The CBF3/4 phylogeny were constructed using CBF3/4 genes from Triticeae , Poaea , and rice . Additional barley CDS sequences were collected from a large collection of full length (fl) cDNA and used in one analyses of the FST phylogeny. CDS were predicted from fl-cDNA sequences with orfpredictor (http://proteomics.ysu.edu/tools/OrfPredictor.html) using homology information from blastx searches against proteins from rice, maize, sorghum, and B. distachyon downloaded from ftp.plantbiology.msu.edu, ftp.brachypodium.org, ftp.maizesequence.org, and ftp://ftpmips.helmholtz-muenchen.de/plants/sorghum/, respectively. Only target sequences with a blast evalue <1e-10 were included in further analyses. Multiple sequence alignments of IRIPs and CBFs were made with MAFFT  using a codon model at the Guidance web server (http://guidance.tau.ac.il) . One hundred bootstrap replicates were calculated for each guidance alignment and sequences below 0.6 sequence score and 0.8 column score were automatically removed from the alignments. The FST alignment were produced with prankster  on amino acid residues which were back translated to codons for phylogenetic analyses. All alignments were manually checked and edited in BioEdit .
Maximum Likelihood (ML) phylogenetic reconstruction using the best evolutionary model according to Akaike Information Criterion (AIC) was carried out in Treefinder . The trees from ML phylogenetic reconstructions were checked for consistency with alternative Bayesian phylogenetic reconstruction estimated with MrBayes . The model used in the MrBayes analyses were GTR + G + I (lset nst = 6 rates = invgamma), and simulations were run for 1 × 105 or 1 × 106 generations with tree sampling every 100 or 1000 generations, respectively. Average deviations between the split frequencies were <0.01 in all analyses and a ‘burn in’ which excluded half the sampled trees was used to generate a consensus tree. Phylogenetic trees were visualized in MEGA v5 .
Bacteriophage library construction and screening
DNA of the inbred diploid line Bd3-1 was used to construct a bacteriophage lambda genomic library. Bd3-1 genomic DNA was a kind gift from Dr. David Garvin, University of Minnesota, USA. This library was screened with a rice cDNA probe encoding C-Repeat Binding Factor (CBF)/Dehydration Responsive Element Binding Protein (Os-DREB1A), and a probe encoding Os-DREB1B . DREB1A is an CBF3 subfamily CBF while DREB1B is an CBF4 subfamily CBF. The same clones cross-hybridizing to Os-DREB1A cross hybridized to Os-DREB1B, and no additional Os-DREB1B cross-hybridizing clones were identified. All clones fell into one of two classes based on restriction enzyme patterns. Two representative clones were sequenced and these sequences deposited in GenBank (accessions JQ180470 and JQ180471).
We measured total carbohydrate and fructan content in one B. distachyon spring (21–1) and winter type (29–1) and in the core Pooideae species Lolium perenne L. (perennial ryegrass) and Phleum pratense L. (timothy) before and after cold acclimation. Eight plants of each species were grown in the greenhouse under 16 h photoperiod. After eight weeks, half the plants were placed in a cold chamber at 2°C, while the other half was kept in the greenhouse as control. Leaf tissue from cold treated and control plants was harvested after four days and stored at −80°C. Extractions of total carbohydrate and fructan were carried out as described in Thorsteinsson et al. . For the colorimetric quantification we made slight modifications to the method described in Pollock ; because levan is the principal monocot fructan in grasses [52, 53] we chose levan as our fructan standard instead of inulin. Different concentrations (0.2, 0.4, 0.6, 0.8, 1, 2, 3, 4, and 5 mg ml-1) of levan and glucose were used to make separate standard curves.
Fructan samples were analyzed by high-performance anion-exchange chromatography (HPAEC) on a Dionex ICS3000 system (Dionex Corp.). Two μl of each filtrate was injected on a CarboPac PA1 column (2x250 mm analytical column with a 2x50 mm guard column) operated at 30°C, with 0.25 ml min-1 and analyses were detected with pulsed amperometric detection (PAD). Analyte separation was obtained by applying a gradient of eluent A (100 mM NaOH) and B (1.0 M NaOAc in 0.1 M NaOH) starting at 100% eluent A, followed by a two min linear gradient to 5% eluent B, then increasing to 25% eluent B in 23 min and a final increase to 50% eluent B reached at 45 min was kept for three min. Column reconditioning was obtained by returning to initial conditions in one min which was kept for 10 min. The following external standards were used for peak identification; levan (Sigma-Aldrich), fructose (Sigma-Aldrich), glucose (Sigma-Aldrich), sucrose (Sigma-Aldrich), and a set of fructooligosaccharides (1-Kestose, Nystose and 1-Fructofuranosylnystose) (Wako Chemicals).
Results and discussion
Cold responsive IRIP genes evolved early in the pooideae lineage
Primer sequences used for qRT-PCR experiments
forward primer 5′ > 3′
reverse primer 5′ > 3′
probe 5′ > 3′
product size (bp)
Gene expression values from microarray study of cold induction of IRIP, CBF 3 and FST-like homologs in Brachypodium distachyon
Fold change values (treatment/control) for each time point
The large standard deviations of the qRT-PCR expression values obscures detailed analyses of expression differences, nevertheless one interesting pattern in the IRIP expression data is apparent. In our limited dataset, the winter types express IRIPs more rapid than the spring types. For three of the four IRIP genes only the winter types have significant p-values at 4 hours (Figure 3a-b,d), while the fourth IRIP gene (Figure 3c) has marked higher expression in winter types at 4 hours, compared to Bd21-1. The pattern of low IRIP gene induction early in the CA treatment (<5 hrs) is also observed for Bd21-1 in the microarray experiments (Table 2). One hypothesis is therefore that winter types have a more rapid IRIP gene cold induction compared to spring types. It must be noted that expression data from many more lines of different flowering habits is needed to test this hypothesis.
Occurrence of bona fide IRIP genes in B. distachyon, containing the conserved ice-binding domain and being cold responsive, places the evolution of IRIP’s early in the Pooideae evolution prior to the Brachypodium-core Pooideae divergence. Moreover, the phylogenetic analysis supports a monophyletic origin of the B. distachyon IRIP genes (Figure 2b) which means that independent IRIP gene family expansions occurred in the Brachypodium lineage after the divergence from core Pooideae species.
Low IRIP induction in Bd3-1 could be explained by extreme spring type habits
The Bd3-1 spring type showed dramatically lower IRIP cold induction compared to the other lines, with only three qRT-PCR measurements being significantly higher than the NA samples (Figure 3). Since the spring type Bd21-1 has strong cold induced IRIP gene expression, the Bd3-1 IRIP expression phenotype cannot be related to the spring type life strategy per se. A recent vernalization response study showed that Bd3-1 is a rapid flowering spring type which expresses the flowering promoting genes VRN1 and VRN3 at very high levels early in the life cycle. In fact non-vernalized seedlings of Bd3-1 had approximately 4- and 6-fold higher expression levels of one of the two VRN1 paralogs and VRN3, respectively, compared to Bd21 . The physiological transition from vegetative to generative growth stage in cereals is associated with VRN1 induction, repression of cold induced gene expression, and loss of freezing tolerance . Thus, the strikingly low IRIP gene expression observed in Bd3-1 could be related to the very rapid transition from vegetative growth form to flowering observed in this line .
Brachypodium distachyonhas a large cold responsive CBF3 family but lack CBF4 genes
A monophyletic origin of all B. distachyon CBF3d paralogs is supported by both ML and Bayesian phylogenetic reconstructions (Figure 4). Moreover, several CBF3d members are found on chromosome 2 and 3 (Figure 4), which does not conform to the syntenic Triticeae 5 relationship . Taken together, this data indicate extensive B. distachyon specific duplications of CBF3d genes, both tandem and to other chromosomes, even though alternative scenarios could explain the observed CBF3d topology. First, what appears to be a B. distachyon-specific CBF3d clade could have evolved prior to Brachypodium divergence and later lost in core Pooideae. Second, gene conversion may homogenize gene sequences and create an ‘illusion’ of evolutionary relatedness . Lastly, all sequence orthologs of the B. distachyon CBF3d genes might not yet have been discovered in core Pooideae species.
Surprisingly, none of the 13 B. distachyon CBF homologs identified in the blast analyses belonged to the CBF4 group genes (Additional file 3). Screening a bacteriophage lambda genomic library for B. distachyon CBF gene content also failed to recover CBF4 homologs from B. distachyon. Two phage clones, λBd1C (JQ180470) and λBd5D (JQ180471), harbored four and three CBF3 genes, respectively, and comparison to the Bd21-1 genome showed that λBd1C and λBd5D correspond to two regions on Bd21-1 chromosome 4. The sequence encompassed by both clones was colinear with the Bd21-1sequence over their entire length. Because rice contains a single CBF4 gene (OsCBF4 AY785894), the most parsimonious model to explain the absence of CBF4 in B. distachyon is lineage specific gene loss.
Brachypodium distachyonlacks FST genes and differs from core pooideae species in fructan accumulation during cold stress
FST gain probably evolved through a mutation in the substrate binding site of a vacuolar invertase gene . A vacuolar invertase in rice is induced by low temperatures  and in this study the expression of a vacuolar invertase homolog in B. distachyon (Bradi3g00910.1) is shown to be induced by low temperatures. Hence, it is likely that the ancestral FST gene was responsive to cold prior to the evolution of FST function. The observation that the closest homolog outside the FST clade was a barley gene (NIASHV2001M19) supports an evolutionary model of gain of FST-function after the core Pooideae group divergence from Brachypodium. However, a less parsimonious scenario of a FST gene(s) loss in Brachypodium could also explain the phylogenetic results and thus cannot be ruled out. Compared to a FST gain in core Pooideae, a ‘loss in Brachypodium’ scenario would require a higher number of evolutionary changes (i.e. 2 changes) to be consistent with the topology in Figure 5; in addition to the FST loss, multiple FST gains must have occurred in the Pooideae lineage, or the NIASHV2001M19 barley gene must have reverted back from a FST gene to an invertase-like gene.
The distribution of fructan synthesising plants is skewed towards ecosystems characterized by intermittent drought and low temperature stress . Furthermore, both correlative [69–72] and transgenic studies  provide compelling evidence for an important role of fructans in drought and cold stress tolerance in core Pooideae grasses. Cold stress associated fructan accumulation was historically assumed to be linked to storage of easily accessible energy reserves as plants prepare for winter . However, results from functional studies have provided insights into a more direct role of fructans in abiotic stress protection, as part of stability enhancing complexes of the cell membrane lipid bi-layers during freezing stress [75–77]. It is thus possible that evolution of FST function, and subsequent increase and diversification of this enzyme family, was important for adaptation to environments with increased abiotic stress levels, such as colder ecosystems, in a core Pooideae ancestor.
Differences between B. distachyon and core Pooideae also reveal interesting, and potentially biological important, clues to understand the evolution and function of cold acclimation and freezing tolerance in the Poaceae. The absence of genes encoding enzymes in fructan biosynthesis provides a unique opportunity to carry out investigations, using transgenic techniques, to test hypotheses on the evolution of fructan metabolism in relation to adaptation to abiotic stress in grasses.
Understanding the underlying genetic factors controlling climate adaptation within a species is of great importance for breeding of abiotic stress resilience in crop plants. Our study revealed large differences in the transcriptional responses to cold stress among different B. distachyon lines. The Bd3-1 spring type had substantially lower levels of IRIP expression at all time-points during cold acclimation (Figure 3), compared to the other Bd-lines in the study. It is possible that this ‘cold non-responsiveness’ is related to the early expression of flowering pathway genes in Bd3-1 . Bd3-1 could therefore be an interesting model to study how mechanistically the transition to flowering (reproductive stage) interact with the CA pathways. Moreover, winter and spring types differed in the transcriptional response time for the IRIP genes. These differences could be related to local adaptations to climatic conditions, and thus provide an interesting model system to study general population differences and adaptation to cold stress responses in Pooideae. Because major transcriptional regulating pathways are conserved across highly divergent species, knowledge about mechanisms for local adaptation in B. distachyon populations could be transferable and valuable for agricultural important Pooideae crops.
Thanks to Øyvind Jørgensen for technical assistance and Dr. David Garvin for kindly providing Brachypodium distachyon germplasm resources. This work was supported by the Department of Plant and Environmental Sciences, Norwegian University of Life Sciences and the United States Department of Energy Plant Feedstock Genomics for Bioenergy grant DE-FG02-08ER64630 to TCM.
- The International Brachypodium Initiative: Genome sequencing and analysis of the model grassBrachypodium distachyon. Nature. 2010, 463 (7282): 763-768. 10.1038/nature08747.View ArticleGoogle Scholar
- Bevan MW, Garvin DF, Vogel JP: Brachypodium distachyon genomics for sustainable food and fuel production. Curr Opin Biotech. 2010, 21 (2): 211-217. 10.1016/j.copbio.2010.03.006.PubMedView ArticleGoogle Scholar
- Vogel J, Hill T: High-efficiency Agrobacterium-mediated transformation of Brachypodium distachyon inbred line Bd21-3. Plant Cell Rep. 2008, 27 (3): 471-478. 10.1007/s00299-007-0472-y.PubMedView ArticleGoogle Scholar
- Schwartz CJ, Doyle MR, Manzaneda AJ, Rey PJ, Mitchell-Olds T, Amasino RM: Natural variation of flowering time and vernalization responsiveness in Brachypodium distachyon. Bioenerg Res. 2010, 3 (1): 38-46. 10.1007/s12155-009-9069-3.View ArticleGoogle Scholar
- Draper J, Mur LAJ, Jenkins G, Ghosh-Biswas GC, Bablak P, Hasterok R, Routledge APM: Brachypodium distachyon: a new model system for functional genomics in grasses. Plant Physiol. 2001, 127 (4): 1539-1555. 10.1104/pp.010196.PubMedView ArticleGoogle Scholar
- Hartley W: Studies on origin, evolution, and distribution of gramineae. V.the subfamily Festucoideae. Aust J Bot. 1973, 21 (2): 201-234. 10.1071/BT9730201.View ArticleGoogle Scholar
- Barker NP, Clark LG, Davis JI, Duvall MR, Guala GF, Hsiao C, Kellogg EA, Linder HP, Grass Phylogeny Working Group G: Phylogeny and subfamilial classification of the grasses (Poaceae). Ann Mo Bot Gard. 2001, 88 (3): 373-457. 10.2307/3298585.View ArticleGoogle Scholar
- Stromberg CAE: Evolution of grasses and grassland ecosystems. Annu Rev Earth Pl Sc. 2011, 39: 517-544. 10.1146/annurev-earth-040809-152402.View ArticleGoogle Scholar
- Sandve SR, Fjellheim S: Did gene family expansions during the Eocene-Oligocene boundary climate cooling play a role in pooideae adaptation to cool climates?. Mol Ecol. 2010, 19 (10): 2075-2088. 10.1111/j.1365-294X.2010.04629.x.PubMedView ArticleGoogle Scholar
- Sandve SR, Kosmala A, Rudi H, Fjellheim S, Rapacz M, Yamada T, Rognli OA: Molecular mechanisms underlying frost tolerance in perennial grasses adapted to cold climates. Plant Sci. 2011, 180 (1): 69-77. 10.1016/j.plantsci.2010.07.011.PubMedView ArticleGoogle Scholar
- Hulten E, Fries M: Atlas of North European vascular plants - North of theTropic of Cancer, Volume 3. Königstein, Federal Republic of Germany: KoeltzScientific Books; 1986.Google Scholar
- Tutin TG: CXCIII Graminea. In Flora Europaea. Edited by Tutin DHH TG,Burges NA, Moore DM, Valentine DH, Walters SM, Webb DA. Cambridge:Cambridge University Press; 1980:118–267.Google Scholar
- Filiz E, Ozdemir BS, Budak F, Vogel JP, Tuna M, Budak H: Molecular, morphological, and cytological analysis of diverse Brachypodium distachyon inbred lines. Genome. 2009, 52 (10): 876-890. 10.1139/G09-062.PubMedView ArticleGoogle Scholar
- Zhang C, Fei S-z, Arora R, Hannapel D, Hannapel D: Ice recrystallization inhibition proteins of perennial ryegrass enhance freezing tolerance. Planta. 2010, 232 (1): 155-164. 10.1007/s00425-010-1163-4.PubMedView ArticleGoogle Scholar
- John UP, Polotnianka RM, Sivakumaran KA, Chew O, Mackin L, Kuiper MJ, Talbot JP, Nugent GD, Mautord J, Schrauf GE, et al: Ice recrystallization inhibition proteins (IRIPs) and freeze tolerance in the cryophilic antarctic hair grass Deschampsia antarctica E. Desv. Plant Cell Environ. 2009, 32 (4): 336-348. 10.1111/j.1365-3040.2009.01925.x.PubMedView ArticleGoogle Scholar
- Sidebottom C, Buckley S, Pudney P, Twigg S, Jarman C, Holt C, Telford J, McArthur A, Worrall D, Hubbard R, et al: Phytochemistry: heat-stable antifreeze protein from grass. Nature. 2000, 406 (6793): 256-256. 10.1038/35018639.PubMedView ArticleGoogle Scholar
- Tremblay K, Ouellet F, Fournier J, Danyluk J, Sarhan F: Molecular characterization and origin of novel bipartite cold-regulated ice recrystallization inhibition proteins from cereals. Plant Cell Physiol. 2005, 46 (6): 884-891. 10.1093/pcp/pci093.PubMedView ArticleGoogle Scholar
- Sandve SR, Rudi H, Asp T, Rognli OA: Tracking the evolution of a cold stress associated gene family in cold tolerant grasses. BMC Evol Biol. 2008, 8: 245-10.1186/1471-2148-8-245.PubMedView ArticleGoogle Scholar
- Kumble KD, Demmer J, Fish S, Hall C, Corrales S, DeAth A, Elton C, Prestidge R, Luxmanan S, Marshall CJ, et al: Characterization of a family of ice-active proteins from the Ryegrass, Lolium perenne. Cryobiology. 2008, 57 (3): 263-268. 10.1016/j.cryobiol.2008.09.005.PubMedView ArticleGoogle Scholar
- Sprenger N, Bortlik K, Brandt A, Boller T, Wiemken A: Purification, cloning, and functional expression of sucrose:fructan 6-fructosyltransferase, a key enzyme of fructan synthesis in barley. P Natl Acad Sci USA. 1995, 92 (25): 11652-11656. 10.1073/pnas.92.25.11652.View ArticleGoogle Scholar
- Smeekens S, Vijn I: Fructan: more than a reserve carbohydrate?. Plant Physiol. 1999, 120 (2): 351-359. 10.1104/pp.120.2.351.PubMedView ArticleGoogle Scholar
- Livingston DP, Knievel DP, Gildow FE: Fructan synthesis in oat :1. oligomer accumulation in stems during cold hardening and their in-vitro synthesis in a crude enzyme extract. New Phytol. 1994, 127 (1): 27-36. 10.1111/j.1469-8137.1994.tb04256.x.View ArticleGoogle Scholar
- Rao RSP, Andersen JR, Dionisio G, Boelt B: Fructan accumulation and transcription of candidate genes during cold acclimation in three varieties of Poa pratensis. J Plant Physiol. 2011, 168 (4): 344-351. 10.1016/j.jplph.2010.07.019.PubMedView ArticleGoogle Scholar
- Kawakami A, Yoshida M: Molecular characterization of sucrose:sucrose 1-fructosyltransferase and sucrose:fructan 6-fructosyltransferase associated with fructan accumulation in winter wheat during cold Hardening. Biosci Biotech Bioch. 2002, 66 (11): 2297-2305. 10.1271/bbb.66.2297.View ArticleGoogle Scholar
- Kawakami A, Sato Y, Yoshida M: Genetic engineering of rice capable of synthesizing fructans and enhancing chilling tolerance. J Exp Bot. 2008, 59 (4): 793-802. 10.1093/jxb/erm367.PubMedView ArticleGoogle Scholar
- Yamada T, Hisano H, Kanazawa A, Kawakami A, Yoshida M, Shimamoto Y: Transgenic perennial ryegrass plants expressing wheat fructosyltransferase genes accumulate increased amounts of fructan and acquire increased tolerance on a cellular level to freezing. Plant Sci. 2004, 167 (4): 861-868. 10.1016/j.plantsci.2004.05.037.View ArticleGoogle Scholar
- Thomashow MF, Gilmour SJ, Zarka DG, Stockinger EJ, Salazar MP, Houghton JM: Low temperature regulation of the Arabidopsis CBF family of AP2 transcriptional activators as an early step in cold-induced COR gene expression. Plant J. 1998, 16 (4): 433-442. 10.1046/j.1365-313x.1998.00310.x.PubMedView ArticleGoogle Scholar
- Hayes PM, Skinner JS, von Zitzewitz J, Szucs P, Marquez-Cedillo L, Filichkin T, Amundsen K, Stockinger EJ, Thomashow MF, Chen THH: Structural, functional, and phylogenetic characterization of a large CBF gene family in barley. Plant Mol Biol. 2005, 59 (4): 533-551. 10.1007/s11103-005-2498-2.PubMedView ArticleGoogle Scholar
- Minami A, Nagao M, Ikegami K, Koshiba T, Arakawa K, Fujikawa S, Takezawa D: Cold acclimation in bryophytes: low-temperature-induced freezing tolerance in Physcomitrella patens is associated with increases in expression levels of stress-related genes but not with increase in level of endogenous abscisic acid. Planta. 2005, 220 (3): 414-423. 10.1007/s00425-004-1361-z.PubMedView ArticleGoogle Scholar
- Badawi M, Danyluk J, Boucho B, Houde M, Sarhan F: The CBF gene family in hexaploid wheat and its relationship to the phylogenetic complexity of cereal CBFs. Mol Genet Genomics. 2007, 277 (5): 533-554. 10.1007/s00438-006-0206-9.PubMedView ArticleGoogle Scholar
- Choi D-W, Rodriguez EM, Close TJ: Barley Cbf3 gene ientification, expression pattern, and map location. Plant Physiol. 2002, 129 (4): 1781-1787. 10.1104/pp.003046.PubMedView ArticleGoogle Scholar
- Vágújfalvi A, Galiba G, Cattivelli L, Dubcovsky J: The cold-regulated transcriptional activator Cbf3 is linked to the frost-tolerance locus Fr-A2 on wheat chromosome 5A. Mol Genet Genomics. 2003, 269 (1): 60-67.PubMedGoogle Scholar
- Knox AK, Dhillon T, Cheng HM, Tondelli A, Pecchioni N, Stockinger EJ: CBF gene copy number variation at Frost Resistance-2 is associated with levels of freezing tolerance in temperate-climate cereals. Theor Appl Genet. 2010, 121 (1): 21-35. 10.1007/s00122-010-1288-7.PubMedView ArticleGoogle Scholar
- Vogel J, Bragg J: Brachypodium distachyon, a new model for the Triticeae.In Genetics and Genomics of the Triticeae. vol. 7th edition. Edited byMuehlbauer GJ, Feuillet C. New York: Springer; 2009:427–449.View ArticleGoogle Scholar
- Febrer M, Goicoechea JL, Wright J, McKenzie N, Song X, Lin J, Collura K, Wissotski M, Yu Y, Ammiraju JSS, et al: An Integrated physical, genetic and cytogenetic map of Brachypodium distachyon, a model system for grass research. PLoS ONE. 2010, 5 (10): e13461-10.1371/journal.pone.0013461.PubMedView ArticleGoogle Scholar
- Alm V, Busso CS, Ergon A, Rudi H, Larsen A, Humphreys MW, Rognli OA: QTL analyses and comparative genetic mapping of frost tolerance, winter survival and drought tolerance in meadow fescue (Festuca pratensis Huds.). Theor Appl Genet. 2011, 123 (3): 369-382. 10.1007/s00122-011-1590-z.PubMedView ArticleGoogle Scholar
- Katoh K, Kuma K-i, Toh H, Miyata T: MAFFT version 5: improvement in accuracy of multiple sequence alignment. Nucleic Acids Res. 2005, 33 (2): 511-518. 10.1093/nar/gki198.PubMedView ArticleGoogle Scholar
- Rozen S, Skaletsky H: Primer3 on the WWW for general users and forbiologist programmers. In Bioinformatics methods and protocols: methods inmolecular biology. 132nd edition. Edited by Misener S, Krawetz SA: HumanaPress; 1999:365–386.View ArticleGoogle Scholar
- Fox S, Filichkin S, Mockler TC: Applications of ultra-high-throughput sequencing. Methods Mol Biol. 2009, 553: 79-108. 10.1007/978-1-60327-563-7_5.PubMedView ArticleGoogle Scholar
- Bolstad BM, Irizarry RA, Astrand M, Speed TP: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003, 19 (2): 185-193. 10.1093/bioinformatics/19.2.185.PubMedView ArticleGoogle Scholar
- Skinner J, Zitzewitz J, Szűcs P, Marquez-Cedillo L, Filichkin T, Amundsen K, Stockinger E, Thomashow M, Chen T, Hayes P: Structural, functional, and phylogenetic characterization of a large CBF gene family in barley. Plant Mol Biol. 2005, 59 (4): 533-551. 10.1007/s11103-005-2498-2.PubMedView ArticleGoogle Scholar
- Tamura K, Yamada T: A perennial ryegrass CBF gene cluster is located in a region predicted by conserved synteny between Poaceae species. Theor Appl Genet. 2007, 114 (2): 273-283.PubMedView ArticleGoogle Scholar
- Penn O, Privman E, Landan G, Graur D, Pupko T: An alignment confidence score capturing robustness to guide tree uncertainty. Mol Biol Evol. 2010, 27 (8): 1759-1767. 10.1093/molbev/msq066.PubMedView ArticleGoogle Scholar
- Whelan S, Goldman N: A general empirical model of protein evolution derived from multiple protein families using a maximum-likelihood approach. Mol Biol Evol. 2001, 18 (5): 691-699. 10.1093/oxfordjournals.molbev.a003851.PubMedView ArticleGoogle Scholar
- Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser. 1999, 41: 95-98.Google Scholar
- Jobb G, von Haeseler A, Strimmer K: TREEFINDER: a powerful graphical analysis environment for molecular phylogenetics. BMC Evol Biol. 2004, 4 (1): 1-9. 10.1186/1471-2148-4-1.View ArticleGoogle Scholar
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19 (12): 1572-1574. 10.1093/bioinformatics/btg180.PubMedView ArticleGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28 (10): 2731-2739. 10.1093/molbev/msr121.PubMedView ArticleGoogle Scholar
- Dubouzet JG, Sakuma Y, Ito Y, Kasuga M, Dubouzet EG, Miura S, Seki M, Shinozaki K, Yamaguchi-Shinozaki K: OsDREB genes in rice, Oryza sativa L., encode transcription activators that function in drought-, high-salt- and cold-responsive gene expression. Plant J. 2003, 33 (4): 751-763. 10.1046/j.1365-313X.2003.01661.x.PubMedView ArticleGoogle Scholar
- Thorsteinsson B, Harrison PA, Chatterton NJ: Fructan and total carbohydrate accumulation in leaves of two cultivars of timothy (Phleum pratense Vega and Climax) as affected by temperature. J Plant Physiol. 2002, 159 (9): 999-1003. 10.1078/0176-1617-00657.View ArticleGoogle Scholar
- Pollock CJ: Patterns of turnover of fructans in leaves of Dactylis glomerata L. New Phytol. 1982, 90 (4): 645-650. 10.1111/j.1469-8137.1982.tb03273.x.View ArticleGoogle Scholar
- Valluru R, Van den Ende W: Plant fructans in stress environments: emerging concepts and future prospects. J Exp Bot. 2008, 59 (11): 2905-2916. 10.1093/jxb/ern164.PubMedView ArticleGoogle Scholar
- Livingston D, Hincha D, Heyer A: Fructan and its relationship to abiotic stress tolerance in plants. Cell Mol Life Sci. 2009, 66 (13): 2007-2023. 10.1007/s00018-009-0002-x.PubMedView ArticleGoogle Scholar
- Dhillon T, Pearce SP, Stockinger EJ, Distelfeld A, Li C, Knox AK, Vashegyi I, Vágújfalvi A, Galiba G, Dubcovsky J: Regulation of freezing tolerance and flowering in temperate cereals: The VRN-1 connection. Plant Physiol. 2010, 153 (4): 1846-1858. 10.1104/pp.110.159079.PubMedView ArticleGoogle Scholar
- Miller A, Galiba G, Dubcovsky J: A cluster of 11 CBF transcription factors is located at the frost tolerance locus Fr-Am2 in Triticum monococcum. Mol Genet Genomics. 2006, 275 (2): 193-203. 10.1007/s00438-005-0076-6.PubMedView ArticleGoogle Scholar
- Skinner J, Szűcs P, von Zitzewitz J, Marquez-Cedillo L, Filichkin T, Stockinger E, Thomashow M, Chen T, Hayes P: Mapping of barley homologs to genes that regulate low temperature tolerance in Arabidopsis. Theor Appl Genet. 2006, 112 (5): 832-842. 10.1007/s00122-005-0185-y.PubMedView ArticleGoogle Scholar
- Osada N, Innan H: Duplication and gene conversion in the Drosophila melanogaster genome. PLoS Genet. 2008, 4 (12): e1000305-10.1371/journal.pgen.1000305.PubMedView ArticleGoogle Scholar
- Lasseur B, Schroeven L, Lammens W, Le Roy K, Spangenberg G, Manduzio H, Vergauwen R, Lothier J, Prud'homme M-P, Van den Ende W: Transforming a fructan:fructan 6 G-fructosyltransferase from perennial ryegrass into a sucrose:sucrose 1-fructosyltransferase. Plant Physiol. 2009, 149 (1): 327-339. 10.1104/pp.108.125559.PubMedView ArticleGoogle Scholar
- del Viso F, Puebla AF, Fusari CM, Casabuono AC, Couto AS, Pontis HG, Hopp HE, Heinz RA: Molecular characterization of a putative sucrose:fructan 6-fructosyltransferase (6-SFT) of the cold-resistant Patagonian grass Bromus pictus associated with fructan accumulation under low temperatures. Plant Cell Physiol. 2009, 50 (3): 489-503. 10.1093/pcp/pcp008.PubMedView ArticleGoogle Scholar
- Hisano H, Kanazawa A, Yoshida M, Humphreys MO, Iizuka M, Kitamura K, Yamada T: Coordinated expression of functionally diverse fructosyltransferase genes is associated with fructan accumulation in response to low temperature in perennial ryegrass. New Phytol. 2008, 178: 766-780. 10.1111/j.1469-8137.2008.02409.x.PubMedView ArticleGoogle Scholar
- Francki MG, Walker E, Forster JW, Spangenberg G, Appels R: Fructosyltransferase and invertase genes evolved by gene duplication and rearrangements: rice, perennial ryegrass, and wheat gene families. Genome. 2006, 49 (9): 1081-1091. 10.1139/g06-066.PubMedView ArticleGoogle Scholar
- Oliver SN, Van Dongen JT, Alfred SC, Mamun EA, Zhao X, Saini HS, Fernandes SF, Blanchard CL, Sutton BG, Geigenberger P, et al: Cold-induced repression of the rice anther-specific cell wall invertase gene OSINV4 is correlated with sucrose accumulation and pollen sterility. Plant Cell Environ. 2005, 28 (12): 1534-1551. 10.1111/j.1365-3040.2005.01390.x.View ArticleGoogle Scholar
- Vágújfalvi A, Kerepesi I, Galiba G, Tischner T, Sutka J: Frost hardiness depending on carbohydrate changes during cold acclimation in wheat. Plant Sci. 1999, 144 (2): 85-92. 10.1016/S0168-9452(99)00058-8.View ArticleGoogle Scholar
- Dionne J, Castonguay Y, Nadeau P, Desjardins Y: Freezing tolerance and carbohydrate changes during cold acclimation of green-type annual bluegrass (Poa annua L.) ecotypes. Crop Sci. 2001, 41 (2): 443-451. 10.2135/cropsci2001.412443x.View ArticleGoogle Scholar
- Livingston DP: Nonstructural carbohydrate accumulation in winter oat crowns before and during cold hardening. Crop Sci. 1991, 31 (3): 751-755. 10.2135/cropsci1991.0011183X003100030043x.View ArticleGoogle Scholar
- Livingston DP, Premakumar R, Tallury SP: Carbohydrate partitioning between upper and lower regions of the crown in oat and rye during cold acclimation and freezing. Cryobiology. 2006, 52 (2): 200-208. 10.1016/j.cryobiol.2005.11.001.PubMedView ArticleGoogle Scholar
- Tronsmo AM, Kobro G, Morgenlie S, Flengsrud R: Carbohydrate content and glycosidase activities following cold hardening in two grass species. Physiol Plantarum. 1993, 88 (4): 689-695. 10.1111/j.1399-3054.1993.tb01390.x.View ArticleGoogle Scholar
- Hendry GAF: Evolutionary origins and natural functions of fructans - a climatological, biogeographic and mechanistic appraisal. New Phytol. 1993, 123 (1): 3-14.View ArticleGoogle Scholar
- Puebla AF, Salerno GL, Pontis HG: Fructan metabolism in two species of Bromus subjected to chilling and water stress. New Phytol. 1997, 136 (1): 123-129.View ArticleGoogle Scholar
- Livingston DP, Henson CA: Apoplastic sugars, fructans, fructan exohydrolase, and invertase in winter oat: responses to second-phase cold hardening. Plant Physiol. 1998, 116 (1): 403-408. 10.1104/pp.116.1.403.View ArticleGoogle Scholar
- Dionne J, Rochefort S, Huff DR, Desjardins Y, Bertrand A, Castonguay Y: Variability for freezing tolerance among 42 ecotypes of green-type annual Bluegrass. Crop Sci. 2010, 50 (1): 321-336. 10.2135/cropsci2008.12.0712.View ArticleGoogle Scholar
- Amiard V, Morvan-Bertrand A, Billard JP, Huault C, Keller F, Prud'homme MP: Fructans, but not the sucrosyl-galactosides, raffinose and loliose, are affected by drought stress in perennial ryegrass. Plant Physiol. 2003, 132 (4): 2218-2229. 10.1104/pp.103.022335.PubMedView ArticleGoogle Scholar
- Hisano H, Kanazawa A, Kawakami A, Yoshida M, Shimamoto Y, Yamada T: Transgenic perennial ryegrass plants expressing wheat fructosyltransferase genes accumulate increased amounts of fructan and acquire increased tolerance on a cellular level to freezing. Plant Sci. 2004, 167 (4): 861-868. 10.1016/j.plantsci.2004.05.037.View ArticleGoogle Scholar
- Vijn I, Smeekens S: Fructan: more than a reserve carbohydrate?. Plant Physiol. 1999, 120 (2): 351-360. 10.1104/pp.120.2.351.PubMedView ArticleGoogle Scholar
- Demel RA, Dorrepaal E, Ebskamp MJM, Smeekens JCM, de Kruijff B: Fructans interact strongly with model membranes. BBA-Biomembranes. 1998, 1375 (1–2): 36-42.PubMedView ArticleGoogle Scholar
- Hincha DK, Hellwege EM, Heyer AG, Crowe JH: Plant fructans stabilize phosphatidylcholine liposomes during freeze-drying. Eur J Biochem. 2000, 267 (2): 535-540. 10.1046/j.1432-1327.2000.01028.x.PubMedView ArticleGoogle Scholar
- Vereyken IJ, Chupin V, Demel RA, Smeekens SCM, De Kruijff B: Fructans insert between the headgroups of phospholipids. BBA-Biomembranes. 2001, 1510 (1–2): 307-320.PubMedView ArticleGoogle Scholar