Comparative analyses reveal potential uses of Brachypodium distachyon as a model for cold stress responses in temperate grasses

doi:10.1186/1471-2229-12-65

BMC Plant Biol. 2012; 12: 65.

Published online 2012 May 8. doi: 10.1186/1471-2229-12-65

PMCID: PMC3487962

Comparative analyses reveal potential uses of Brachypodium distachyon as a model for cold stress responses in temperate grasses

Chuan Li,^1,² Heidi Rudi,² Eric J Stockinger,³ Hongmei Cheng,⁴ Moju Cao,¹ Samuel E Fox,⁵ Todd C Mockler,⁶ Bjørge Westereng,⁷ Siri Fjellheim,² Odd Arne Rognli,² and Simen R Sandve²

¹Maize Research Institute, Sichuan Agricultural University, Sichuan, China

²Department of Plant and Environmental Sciences, Norwegian University of Life Sciences, ÅS, Norway

³Department of Horticulture and Crop Science, The Ohio State University/OARDC, Wooster, OH, 44691, USA

⁴Biotechnology Research Institute, Chinese Academy of Agricultural Sciences, Beijing, 100081, China

⁵Department of Botany and Plant Pathology and Center for Genome Research and Biocomputing, Oregon State University, Corvallis, OR, USA

⁶Donald Danforth Plant Science Center, Saint Louis, MO, 63132, USA

⁷Department of Chemistry, Biotechnology and Food Science, Norwegian University of Life Sciences, Ås, Norway

Corresponding author.

Chuan Li: lichuan172/at/163.com; Heidi Rudi: heidi.rudi/at/umb.no; Eric J Stockinger: stockinger.4/at/osu.edu; Hongmei Cheng: chenghm/at/caas.net.cn; Moju Cao: mojupp/at/163.com; Samuel E Fox: foxsa/at/onid.orst.edu; Todd C Mockler: tmockler/at/danforthcenter.org; Bjørge Westereng: bjorge.westereng/at/umb.no; Siri Fjellheim: siri.fjellheim/at/umb.no; Odd Arne Rognli: odd-arne.rognli/at/umb.no; Simen R Sandve: simen.sandve/at/umb.no

Received December 23, 2011; Accepted April 27, 2012.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License ( http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

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.

Results

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.

Conclusions

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.

Keywords: Brachypodium distachyon, Cold climate adaptation, Ice recrystallization inhibition protein, Gene expression, Fructosyltransferase, C-repeat binding factor, Gene family evolution

Background

Brachypodium distachyon became the first Pooideae grass species to have its genome fully sequenced [1]. 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 [2]; 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) [1], and it can be genetically transformed via Agrobacterium-mediated transformation [3]. 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 [4], and disease resistance [5]. 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 [6] (Figure (Figure1).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 [9]. Hence, the adaptation of the Pooideae to cooler climates makes this group an ideal model system for studying adaptive evolution in plants [10]. 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 (Figure1),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 (Figure11).

Figure 1

Worldwide distribution of members of Pooideae species as percentage of the total grass flora. The native distributions of Brachypodium distachyon and a typical cool season Pooideae forage grass, Phleum pratense, are indicated in two different green shadings. (more ...)

Phylogenetic studies suggest that B. distachyon diverged from the core Pooideae approximately 35 million years ago [2] while key Pooideae-specific adaptations to cold climates evolved during the Eocene-Oligocene cooling period (34–26 Mya), after the B. distachyon-core Pooideae split [9]. If this is correct, shared ancestral molecular mechanisms involved in cold and freezing stress might differ between B. distachyon and agriculturally important species of the Triticeae and Poeae tribes (referred to hereafter as core Pooideae).

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 [17] by gaining an ice binding domain [16] through an expansion of a repeat motif (NxVxG/NxVxxG) [18]. 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 [30] 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.

Methods

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. [36] 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 distachyon IRIP-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 [37] to verify that B. distachyon IRIP genes contained the typical NxVxG/NxVxxG-repeat ice binding domain. IRIP-paralog specific primers were designed using primer3 [38] 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 [39]. 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; [40]). 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.

Phylogenetic analyses

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. [41] as queries. The CBF3/4 phylogeny were constructed using CBF3/4 genes from Triticeae [41], Poaea [42], and rice [30]. 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 [37] using a codon model at the Guidance web server ( http://guidance.tau.ac.il) [43]. 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 [44] on amino acid residues which were back translated to codons for phylogenetic analyses. All alignments were manually checked and edited in BioEdit [45].

Maximum Likelihood (ML) phylogenetic reconstruction using the best evolutionary model according to Akaike Information Criterion (AIC) was carried out in Treefinder [46]. The trees from ML phylogenetic reconstructions were checked for consistency with alternative Bayesian phylogenetic reconstruction estimated with MrBayes [47]. The model used in the MrBayes analyses were GTR + G + I (lset nst = 6 rates = invgamma), and simulations were run for 1 × 10⁵ or 1 × 10⁶ 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 [48].

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 [49]. 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).

Fructan measurement

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. [50]. For the colorimetric quantification we made slight modifications to the method described in Pollock [51]; 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

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

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

Seven Brachypodium distachyon IRIP-like genes from two gene clusters containing two and five genes were identified on chromosome 5 (Figure (Figure2a).2a). Bradi5g22870.1 and Bradi5g22880.1 in the proximal cluster have both truncated ice binding domains (Additional file 1). To test whether the B. distachyon IRIPs are induced by low temperature, paralog specific qRT-PCR primers were designed for four IRIP paralogs with non-truncated ice-binding domains (Bradi5g27300.1, Bradi5g27310.1, Bradi5g27330.1, Bradi5g27350.1) (Table (Table1).1). Strong cold induction of all genes except Bradi5g27300.1 was observed after one day of cold acclimation in both winter types and the 21–1 line (Figure (Figure3).3). Bradi5g27300.1 was cold induced in three genotypes, but the level of expression was generally lower than the other IRIP genes (Figure (Figure3,3, Additional file 2). Expression data from the microarray experiment confirmed the general patterns of IRIP gene cold induction observed in the qRT-PCR experiment; little or no cold induction of Bradi5g27300.1 and medium to strong (2.5-25 fold) cold induction of the other IRIP genes within 24 hours cold treatment (Table (Table22).

Figure 2

Localization and distribution of seven IRIP genes identified in Brachypodium distachyon on chromosome 5 and B. distachyon IRIP genes phylogeny. (a) The location of seven Brachypodium distachyon IRIP genes on chromosome 5. (b) 50% consensus ML phylogeny (more ...)

Table 1

Primer sequences used for qRT-PCR experiments

Figure 3

Quantitative reverse transcript PCR analyses of Brachypodium distachyon IRIP genes at three timepoints after cold acclimation. (a) Bradi5g27300.1 (b) Bradi5g27310.1 (c) Bradi5g27330.1 (d) Bradi5g27350.1 The expression values are given as 2^-ΔΔCt (more ...)

Table 2

Gene expression values from microarray study of cold induction of IRIP, CBF3and FST-like homologs inBrachypodium distachyon

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 (Figure3a-b,d),3a-b,d), while the fourth IRIP gene (Figure (Figure3c)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 (Table2).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 (Figure2b)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 (Figure3).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 [4]. 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 [54]. 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 [4].

Brachypodium distachyon has a large cold responsive CBF3 family but lack CBF4 genes

Fourteen out of the total 18 B. distachyon CBF3/4 gene homologs identified in the blast search could be classified as CBF3 members according to the phylogenetic analysis of all CBF homologs (Additional file 3). Thirteen of these genes belong to the CBF3c/d clades. Figure

4 shows the phylogeny of the CBF3c/d genes in B. distachyon and two core Pooideae species. One B. distachyon gene (Bradi4g35630.1) belongs to the CBF3c subgroup, while the other 12 belong to the CBF3d subgroup. Both the single CBF3c (Bradi4g35630.1) gene member and the majority of the CBF3d genes in B. distachyon belong to a gene cluster on chromosome 4. This chromosome is in large parts syntenic with Triticeae chromosome 5 [1] which contains clusters of cold induced wheat and barley CBF3c/d genes [55,56]. Microarray expression data shows that all but two CBF3d genes (Bradi1g57970.1 and Bradi3g57360.1) are cold induced (>2 fold) during 24 hours of cold exposure (Table (Table2).2). Together with the phylogenetic analyses this expression data suggests that both cold responsive CBF3c and CBF3d genes were present in a Pooideae ancestor prior to Brachypodium and core Pooideae divergence.

Figure 4

Brachypodium distachyon CBF3c/d and core Pooideae CBF3c/d phylogenetic analysis. Fifty % consensus ML phylogeny (J2

G model) with 100 bootstrap replicates showing the relationships between the B. distachyon CBF3c/d members and core Pooideae (more ...)

A monophyletic origin of all B. distachyon CBF3d paralogs is supported by both ML and Bayesian phylogenetic reconstructions (Figure (Figure4).4). Moreover, several CBF3d members are found on chromosome 2 and 3 (Figure (Figure4),4), which does not conform to the syntenic Triticeae 5 relationship [1]. 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 [57]. 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 distachyon lacks FST genes and differs from core pooideae species in fructan accumulation during cold stress

In total, fifteen genes with some homology to the core Pooideae FSTs were identified. Thirteen of these were more distantly related than the closest FST homolog in rice (an invertase-like gene), and were thus not considered in the analyses. Core Pooideae FST genes encode a diagnostic h(A/G)Y/F motif [58] and are induced by low temperatures [59,60]. Neither of the two FST-like homologs in B. distachyon encoded the FST-motif, but just the invertase motif (Figure (Figure5).5). The most distant homolog to the FSTs (Bradi3g00910.1) was cold induced in Bd21 (Table (Table2)2) and belongs to a monophyletic cluster of invertase-like genes containing gene members from core Pooideae species (Figure (Figure5).5). Bradi1g52210.1 was placed closest to the FSTs in the phylogeny, but was not cold induced in Bd21-1 within 24 hours of low temperature treatment (Table (Table2).2). Some FST-genes are known to be induced by cold only after several days of CA [60], but unfortunately we do not have experimental data for B. distachyon cold treatments longer than 24

hrs. Hence it is possible that also Bradi1g52210.1 is induced later in the CA process.

Figure 5

Phylogeny analysing Brachypodium distachyon FST and invertase genes compare to core Pooideae FST and invertase genes. Fifty % consensus ML phylogeny (GTR + G + I model) of FST genes from core Pooideae species, B. distachyon FST homologs and two rice invertase-like (more ...)

FST gain probably evolved through a mutation in the substrate binding site of a vacuolar invertase gene [61]. A vacuolar invertase in rice is induced by low temperatures [62] 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.

To test whether B. distachyon accumulated fructans in response to low temperatures, total carbohydrate and fructan levels were measured before and after cold acclimation. Both core Pooideae species and B. distachyon showed marked responses in carbohydrate accumulation during cold stress (Figure (Figure6)6) as has been shown for many plant species [63-67]. Interestingly, cold induced fructan content increases were much higher in core Pooideae (0.8-1.2 fold) than in B. distachyon (0.2-0.3 fold). Qualitative analysis by HPAEC confirmed that both B. distachyon and L. perenne induce a large increase in fructose, glucose and sucrose during low temperature. However, the profile of short oligosaccharides that accumulated in B. distachyon was different from core Pooideae (Additional file 4). For example, significant levels of kestose, nystose and other unidentified oligosaccharides are present after cold acclimation of L. perenne but not in cold acclimated B. distachyon (Additional file 4a). Moreover, sugars with higher degree of polymerization (DP) are also present to a much larger extent in cold treated core Pooideae compared to B. distachyon (Additional file 4b). The correlation of cold induced modifications in fructan content and the phylogenetic clustering of presence/absence of amino acid motifs conferring fructosyltransferase activity support that B. distachyon does not possess orthologs of the core Pooideae FST enzymes (Figure (Figure5).5).

Figure 6

Carbohydrate accumulation in Brachypodium distachyon and core Pooideae species in response to four days of cold acclimation at 2°C.(a) Fructan contents (mg

ml^-1) and (b) total carbohydrate content (mg

ml^-1) from colorimetric quantification. (more ...)

The distribution of fructan synthesising plants is skewed towards ecosystems characterized by intermittent drought and low temperature stress [68]. Furthermore, both correlative [69-72] and transgenic studies [73] 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 [74]. 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.

Conclusion

It is evident from our comparative analyses that B. distachyon and the core Pooideae differ in key cold stress pathways (Figure (Figure7).7). Even though this difference limits the use of B. distachyon as a holistic model for the molecular biology of low temperature stress in core Pooideae species, B. distachyon will be useful to study specific genes and pathways, such as CBF3 or IRIP genes. For example, IRIP RNAi knockout/knockdown lines can be used to test the importance of IRIP gene function for cold and freezing tolerance in Pooideae, and IRIP promoter deletion constructs will be able to shed light on IRIP transcriptional regulation. B. distachyon could also be useful to understand the functional divergence between different CBF-gene families, such as the CBF3d and CBF3c group genes, and thus increase our general understanding of transcriptional control of cold stress responses in Pooideae. Since the CBF3c gene is present in single copy in B. distachyon, this could facilitate functional analyses and the pathways this gene affect at a mechanistic level.

Figure 7

Summary of the key steps in the evolution CBF3/4, IRIP, and FST genes in Pooideae after the divergence from rice. The evolutionary steps fall into three categories, ‘function’ (i.e. change of gene function), ‘expansion’ (more ...)

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 (Figure3),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 [4]. 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.

Authors’ contributions

CL carried out the experimental work on qRT-PCR gene expression and analyses of total sugar/fructan levels, as well as participating in drafting the manuscript. HR helped design the study, carried out experimental work on gene expression. EJS and HC carried out the screening of B. distachyon genomic libraries. EJS also participated in drafting the manuscript. MC was involved in the data acquisition and manuscript preparation. SF collected and analyzed the species distributions and helped draft the manuscript. SEF and TCM carried out microarray based gene expression analysis and participated in drafting the manuscript. BW was responsible for all sugar content analyses and carried out the HPAEC analyses. OAR participated in the design of the study and helped draft the manuscript. SRS participated in the study design, carried out the phylogenetic analyses, and was responsible for the final manuscript preparation.

Supplementary Material

Additional file 1

Amino-acid sequences alignment analysis of Brachypodium distachyon IRI proteins with Lolium perenne IRI proteins sequences. Sites with with black shade are highly conserved (>70% of sequences). Bradi5g27350.1 and Bradi5g22870.1 have truncated ice binding-domains.

Click here for file^{(1.6M, pdf)}

Additional file 2

qRT-PCR expression levels in fold change and p-values of Brachypodium distachyon IRIP genes.

Click here for file^{(58K, pdf)}

Additional file 3

Minimum evolution CBF3/4 gene phylogeny including all Brachypodium distachyon homologs. The Tamura 3-parameter method with gamma distributed rate variation and pairwise deletion was used to calculate evolutionary relatedness. All non-CBF3c/d and CBF4 genes as classified by the phylogeny are in green and red colour, respectively, while other CBF-homologs are in black. Species abbreviations: Os, Oryza sativa; Hv, Hordeum vulgare; Lp, Lolium perenne.

Click here for file^{(188K, pdf)}

Additional file 4

HPAEC results for carbohydrate contents before and after low temperature treatment. Detector intensity given as nano coulomb (nC). Black curve and Blue curve are the cold treated and none treated Lolium perenne plants carbohydrate extraction elution separately. Purple curve and brown curve are the cold treated and none treated Brachypodium distachyon plant carbohydrate extraction elution. (a) Low degree of polymerization oligosaccharides (DP two to five) (glucose, fructose, sucrose, kestose, nystose and other unidentified oligosaccharides) are detected in L. perenne and B. distachyon using HPAEC. (b) Higher degree of polymerization oligosaccharides are detected in L. perenne than in B. distachyon using HPAEC.

Click here for file^{(222K, pdf)}

Acknowledgements

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.

References

The International Brachypodium Initiative. Genome sequencing and analysis of the model grassBrachypodium distachyon. Nature. 2010;463(7282):763–768. doi: 10.1038/nature08747. [PubMed] [Cross Ref]
Bevan MW, Garvin DF, Vogel JP. Brachypodium distachyon genomics for sustainable food and fuel production. Curr Opin Biotech. 2010;21(2):211–217. doi: 10.1016/j.copbio.2010.03.006. [PubMed] [Cross Ref]
Vogel J, Hill T. High-efficiency Agrobacterium-mediated transformation of Brachypodium distachyon inbred line Bd21-3. Plant Cell Rep. 2008;27(3):471–478. doi: 10.1007/s00299-007-0472-y. [PubMed] [Cross Ref]
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. doi: 10.1007/s12155-009-9069-3. [Cross Ref]
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. doi: 10.1104/pp.010196. [PubMed] [Cross Ref]
Hartley W. Studies on origin, evolution, and distribution of gramineae. V.the subfamily Festucoideae. Aust J Bot. 1973;21(2):201–234. doi: 10.1071/BT9730201. [Cross Ref]
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. doi: 10.2307/3298585. [Cross Ref]
Stromberg CAE. Evolution of grasses and grassland ecosystems. Annu Rev Earth Pl Sc. 2011;39:517–544. doi: 10.1146/annurev-earth-040809-152402. [Cross Ref]
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. doi: 10.1111/j.1365-294X.2010.04629.x. [PubMed] [Cross Ref]
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. doi: 10.1016/j.plantsci.2010.07.011. [PubMed] [Cross Ref]
Hulten E, Fries M. Atlas of North European vascular plants - North of the Tropic of Cancer. Federal Republic of Germany: Koeltz Scientific Books, Königstein; 1986.
Tutin TG. In: Flora Europaea. Tutin DHH TG, Burges NA, Moore DM, Valentine DH, Walters SM, Webb DA, editor. Cambridge: Cambridge University Press; 1980. CXCIII Graminea; pp. 118–267.
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. doi: 10.1139/G09-062. [PubMed] [Cross Ref]
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. doi: 10.1007/s00425-010-1163-4. [PubMed] [Cross Ref]
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. doi: 10.1111/j.1365-3040.2009.01925.x. [PubMed] [Cross Ref]
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. doi: 10.1038/35018639. [PubMed] [Cross Ref]
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. doi: 10.1093/pcp/pci093. [PubMed] [Cross Ref]
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. doi: 10.1186/1471-2148-8-245. [PMC free article] [PubMed] [Cross Ref]
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. doi: 10.1016/j.cryobiol.2008.09.005. [PubMed] [Cross Ref]
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. doi: 10.1073/pnas.92.25.11652. [PubMed] [Cross Ref]
Smeekens S, Vijn I. Fructan: more than a reserve carbohydrate? Plant Physiol. 1999;120(2):351–359. doi: 10.1104/pp.120.2.351. [PubMed] [Cross Ref]
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. doi: 10.1111/j.1469-8137.1994.tb04256.x. [Cross Ref]
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. doi: 10.1016/j.jplph.2010.07.019. [PubMed] [Cross Ref]
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. doi: 10.1271/bbb.66.2297. [PubMed] [Cross Ref]
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. doi: 10.1093/jxb/erm367. [PubMed] [Cross Ref]
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. doi: 10.1016/j.plantsci.2004.05.037. [Cross Ref]
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. doi: 10.1046/j.1365-313x.1998.00310.x. [PubMed] [Cross Ref]
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. doi: 10.1007/s11103-005-2498-2. [PubMed] [Cross Ref]
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. doi: 10.1007/s00425-004-1361-z. [PubMed] [Cross Ref]
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. doi: 10.1007/s00438-006-0206-9. [PMC free article] [PubMed] [Cross Ref]
Choi D-W, Rodriguez EM, Close TJ. Barley Cbf3 gene ientification, expression pattern, and map location. Plant Physiol. 2002;129(4):1781–1787. doi: 10.1104/pp.003046. [PubMed] [Cross Ref]
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. [PubMed]
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. doi: 10.1007/s00122-010-1288-7. [PubMed] [Cross Ref]
Vogel J, Bragg J. In: Genetics and Genomics of the Triticeae. vol. 7. Muehlbauer GJ, Feuillet C, editor. New York: Springer; 2009. Brachypodium distachyon, a new model for the Triticeae; pp. 427–449.
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. doi: 10.1371/journal.pone.0013461. [PMC free article] [PubMed] [Cross Ref]
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. doi: 10.1007/s00122-011-1590-z. [PubMed] [Cross Ref]
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. doi: 10.1093/nar/gki198. [PMC free article] [PubMed] [Cross Ref]
Rozen S, Skaletsky H. In: Bioinformatics methods and protocols: methods in molecular biology. 132. Misener S, Krawetz SA, editor. Humana Press; 1999. Primer3 on the WWW for general users and for biologist programmers; pp. 365–386.
Fox S, Filichkin S, Mockler TC. Applications of ultra-high-throughput sequencing. Methods Mol Biol. 2009;553:79–108. doi: 10.1007/978-1-60327-563-7_5. [PubMed] [Cross Ref]
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. doi: 10.1093/bioinformatics/19.2.185. [PubMed] [Cross Ref]
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. doi: 10.1007/s11103-005-2498-2. [PubMed] [Cross Ref]
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. [PubMed]
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. doi: 10.1093/molbev/msq066. [PMC free article] [PubMed] [Cross Ref]
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. doi: 10.1093/oxfordjournals.molbev.a003851. [PubMed] [Cross Ref]
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.
Jobb G, von Haeseler A, Strimmer K. TREEFINDER: a powerful graphical analysis environment for molecular phylogenetics. BMC Evol Biol. 2004;4(1):1–9. doi: 10.1186/1471-2148-4-1. [PMC free article] [PubMed] [Cross Ref]
Ronquist F, Huelsenbeck JP. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003;19(12):1572–1574. doi: 10.1093/bioinformatics/btg180. [PubMed] [Cross Ref]
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. doi: 10.1093/molbev/msr121. [PMC free article] [PubMed] [Cross Ref]
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. doi: 10.1046/j.1365-313X.2003.01661.x. [PubMed] [Cross Ref]
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. doi: 10.1078/0176-1617-00657. [Cross Ref]
Pollock CJ. Patterns of turnover of fructans in leaves of Dactylis glomerata L. New Phytol. 1982;90(4):645–650. doi: 10.1111/j.1469-8137.1982.tb03273.x. [Cross Ref]
Valluru R, Van den Ende W. Plant fructans in stress environments: emerging concepts and future prospects. J Exp Bot. 2008;59(11):2905–2916. doi: 10.1093/jxb/ern164. [PubMed] [Cross Ref]
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. doi: 10.1007/s00018-009-0002-x. [PMC free article] [PubMed] [Cross Ref]
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. doi: 10.1104/pp.110.159079. [PubMed] [Cross Ref]
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. doi: 10.1007/s00438-005-0076-6. [PubMed] [Cross Ref]
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. doi: 10.1007/s00122-005-0185-y. [PubMed] [Cross Ref]
Osada N, Innan H. Duplication and gene conversion in the Drosophila melanogaster genome. PLoS Genet. 2008;4(12):e1000305. doi: 10.1371/journal.pgen.1000305. [PMC free article] [PubMed] [Cross Ref]
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 6G-fructosyltransferase from perennial ryegrass into a sucrose:sucrose 1-fructosyltransferase. Plant Physiol. 2009;149(1):327–339. doi: 10.1104/pp.108.125559. [PubMed] [Cross Ref]
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. doi: 10.1093/pcp/pcp008. [PubMed] [Cross Ref]
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. doi: 10.1111/j.1469-8137.2008.02409.x. [PubMed] [Cross Ref]
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. doi: 10.1139/g06-066. [PubMed] [Cross Ref]
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. doi: 10.1111/j.1365-3040.2005.01390.x. [Cross Ref]
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. doi: 10.1016/S0168-9452(99)00058-8. [Cross Ref]
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. doi: 10.2135/cropsci2001.412443x. [Cross Ref]
Livingston DP. Nonstructural carbohydrate accumulation in winter oat crowns before and during cold hardening. Crop Sci. 1991;31(3):751–755. doi: 10.2135/cropsci1991.0011183X003100030043x. [Cross Ref]
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. doi: 10.1016/j.cryobiol.2005.11.001. [PubMed] [Cross Ref]
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. doi: 10.1111/j.1399-3054.1993.tb01390.x. [Cross Ref]
Hendry GAF. Evolutionary origins and natural functions of fructans - a climatological, biogeographic and mechanistic appraisal. New Phytol. 1993;123(1):3–14.
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.
Livingston DP III, 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. doi: 10.1104/pp.116.1.403. [Cross Ref]
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. doi: 10.2135/cropsci2008.12.0712. [Cross Ref]
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. doi: 10.1104/pp.103.022335. [PubMed] [Cross Ref]
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. doi: 10.1016/j.plantsci.2004.05.037. [Cross Ref]
Vijn I, Smeekens S. Fructan: more than a reserve carbohydrate? Plant Physiol. 1999;120(2):351–360. doi: 10.1104/pp.120.2.351. [PubMed] [Cross Ref]
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. [PubMed]
Hincha DK, Hellwege EM, Heyer AG, Crowe JH. Plant fructans stabilize phosphatidylcholine liposomes during freeze-drying. Eur J Biochem. 2000;267(2):535–540. doi: 10.1046/j.1432-1327.2000.01028.x. [PubMed] [Cross Ref]
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. [PubMed]

Articles from BMC Plant Biology are provided here courtesy of
BioMed Central