<title>Background</title> <p>The Poaceae family (grasses) contains some of the most economically important and well studied plant species, e.g. maize, wheat, barley, and rice. Generally speaking the Pooideae subfamily, which includes wheat, barley and forage grasses, are adapted to cold seasons. Many species in this subfamily can withstand temperatures far below freezing and intercellular ice formation [##REF##11154332##1##,##REF##17024382##2##]. Rice and maize on the other hand belongs to the subfamilies Ehrhartoideae and Panicoideae, respectively, and are adapted to warm and tropical climates. Pooideae lineage (from now on referred to as cold tolerant grasses) adaptation to cold climates makes grasses an interesting model system for studying climatic adaptation at the physiological and molecular level.</p> <p>Frost tolerance adaptations are in many organisms associated with the evolution of antifreeze proteins (AFPs) [##REF##11222024##3##]. AFPs can affect freezing- and ice crystallisation related stress via different mechanisms. Thermal hysteresis (TH) depresses the freezing point at which ice crystallisation initiates, which render it possible for organisms to survive under freezing temperatures. Ice re-crystallisation inhibition (IRI) on the other hand does not hinder ice crystallisation but manipulates the growth of the ice crystals such that small ice crystals grow at the expense of larger ice crystals, and this has been suggested to prevent or minimize the cellular damage in plants [##REF##12171647##4##]. A third mode of AFP action is membrane stabilisation which has been reported for a fish AFP [##REF##11806929##5##]. Animal AFPs generally possess high thermal hysteresis (TH) characteristics and lower ice crystallisation initiation temperature by 1–5°C [##REF##15358271##6##,##UREF##0##7##]. Plant AFPs on the other hand have low TH-activity, but exhibits strong ice re-crystallisation inhibition (IRI) activity [##REF##15358271##6##].</p> <p>Genes encoding peptides with IRI capacity have evolved independently several times in different lineages of higher plants. These IRI peptides are homologous to diverse protein classes, e.g. thaumatin like proteins, endochitinases, endo-B-1,3-glucanase, and leucine rich repeat (LRR) containing proteins [##REF##15358271##6##,##REF##8552719##8##,##REF##10333479##9##]. Three LRR-domain containing IRI proteins (LRR-IRI) have been identified in plants, one in carrot (DcAFP; accession number <ext-link ext-link-type="gen" xlink:href="AAC6293">AAC6293</ext-link>) and two in wheat (TaIRI1 and TaIRI2 with accession numbers <ext-link ext-link-type="gen" xlink:href="AAX81542">AAX81542</ext-link> and <ext-link ext-link-type="gen" xlink:href="AAX81543">AAX81543</ext-link>) [##REF##15792959##10##,##REF##9756474##11##]. DcAFP has been classified as a polygalacturonase-inhibiting protein (PGIP) but does not exhibit PGIP activity [##REF##17112975##12##]. LRR motifs span across the entire processed DcAFP protein and form 10-loop beta-helix secondary structure with solvent exposed asparagine residues at putative ice binding sites [##REF##14531728##13##]. <italic>TaIRI1 </italic>and <italic>TaIRI2 </italic>genes (accession numbers <ext-link ext-link-type="gen" xlink:href="AY9968588">AY9968588</ext-link> and <ext-link ext-link-type="gen" xlink:href="AY968589">AY968589</ext-link>) have been identified as homologous to the LRR-domain coding region of a rice phytosulfokine LRR receptor kinase (<italic>OsLRR-PSR</italic>: <ext-link ext-link-type="gen" xlink:href="NP_001058711">NP_001058711</ext-link>) and an <italic>Arabidopsis </italic>trans-membrane protein kinase (<italic>AtLRR-TPK</italic>: <ext-link ext-link-type="gen" xlink:href="NP_200200">NP_200200</ext-link>). The wheat IRI peptides differ structurally from DcAFP in that the LRR-domain only comprises about half of the processed peptide [##REF##15792959##10##].</p> <p>In addition to the N-terminal LLR domain, wheat IRI proteins have a C-terminal repeat domain consisting of two similar A and B motifs, NxVxG and NxVxxG, respectively. This repeat domain has been reported to exhibit strong <italic>in vitro </italic>IRI capacity [##REF##10917518##14##]. Interestingly, blast search yields no sequences with homology to the IRI-domain outside the subfamily of cold tolerant grasses [##REF##15792959##10##]. Protein modelling has shown that the A and B repeated motifs of the IRI-domain folds into a B-roll with ice binding sites matching the prism face of ice [##REF##11721016##15##]. Expression studies have shown that increased expression levels in wheat [##REF##15792959##10##] and perennial ryegrass [Rudi et al, unpublished] are correlated to cold acclimation, but no <italic>in vivo </italic>studies to determine the localisation of these grass IRI peptides have been reported in the literature. However, <italic>TaIRI1 </italic>and <italic>TaIRI2 </italic>have been predicted to encode a N-terminal 20 amino acid signal peptide domain targeting the proteins to the secretory pathway, suggesting that the peptides could be located in the extracellular space [##REF##15792959##10##].</p> <p>While <italic>DcAFP </italic>is evolutionarily closely related to PGIPs; <italic>TaIRI1 </italic>and <italic>TaIRI2 </italic>genes are thought to have evolved from a LRR-PSR like ancestor gene. Furthermore, the evolutionary origin of the IRI-domain in grass IRI genes is much less obvious, because the IRI-domain is not homologous to any other sequences outside the cold tolerant grass lineage. Tremblay et al. [##REF##15792959##10##] proposed a "TE-hypothesis" to explain this apparent lack of homologous coding regions; that the IRI-domain had arisen by a transposable element (TE) insertion. However, no TE signature sequence could be identified surrounding the IRI-domain [##REF##15792959##10##], thus no empirical data supports the TE-hypothesis so far.</p> <p>Here we report the identification and characterisation of novel LRR-IRI homologous genes in cold tolerant grass species. We perform a detailed study of the evolutionary relationships between <italic>OsLRR-PSR </italic>and IRI-like genes by analysing sequence divergence at synonymous sites. We also use synonymous site divergence to trace the evolutionary history of the IRI-like gene family with respect to gene duplication events. The evolution of gene families <italic>per se </italic>is in itself a much debated topic, and gene family expansion and subsequent functional diversification is thought to have been a significant factor contributing to adaptations to new environments [##REF##18079367##16##,##REF##16242725##17##]. The evolution of the IRI-like gene family of cold tolerant grasses is discussed in the context of the Duplication-Degeneration-Complementation (DDC) model [##REF##10101175##18##]. Finally we address the unresolved matter of the evolutionary mechanism underlying the birth of the cold tolerant grass IRI-domain, and propose a novel hypothesis on the evolution of this IRI-domain.</p>
<title>Methods</title> <title><italic>In silico </italic>IRI-like sequence mining</title> <p>A blastn search in the NCBI database was performed using <italic>TaIRI1</italic>. All sequences with blast E-value &lt; 1*10<sup>-20 </sup>were downloaded from the EST and core nucleotide databases. Contigs were aligned with alignment parameters set to &gt; 97% identity and &gt; 40 nucleotides overlap using Sequencher (Gene Codes Corp., Ann Arbor, MI, USA). The 97% identity threshold was set to allow contig alignments to include different allelic forms and polymorphisms caused by EST sequencing errors. Non-coding nucleotides (i.e. promoter and 3'UTR) were removed after an initial prediction of open reading frame (ORF), and subsequently the sequences were realigned with identical parameters. All contigs were translated into their predicted amino acid sequence. Sequence contigs with lack of start and stop codon due to incomplete sequence coverage or putative sequence errors causing frame shift mutations were not included in the analysis. We validated the <italic>in silico </italic>mining method by aligning EST mined unigenes with four full length cDNA clones of grass IRI-like genes from the NCBI core nucleotide collection (barley; <ext-link ext-link-type="gen" xlink:href="AK252915">AK252915</ext-link>, <ext-link ext-link-type="gen" xlink:href="AK249041">AK249041</ext-link>/wheat; <ext-link ext-link-type="gen" xlink:href="AY968588">AY968588</ext-link>, <ext-link ext-link-type="gen" xlink:href="AY968589">AY968589</ext-link>).</p> <title>BAC identification and sequencing</title> <p>Two perennial ryegrass BAC libraries were used to identify novel IRI-like genes [##UREF##7##45##]. Primers for the initial identification of novel IRI-like genes were designed from coding sequences of <italic>LpAFP </italic>(<ext-link ext-link-type="gen" xlink:href="AJ277399">AJ277399</ext-link>) and a partial sequence of a <italic>Festuca pratensis </italic>IRI-like homolog (<ext-link ext-link-type="gen" xlink:href="EU684537">EU684537</ext-link>). The LpAFP primer pair had forward primer 5'GATGAACAGCCGAATACGATTTCT3' and reverse primer 5'GCTTCCAGATACAACGTGGTTGCT3', denaturing at 94°C for 4 minutes and then 35 cycles of 94°C 30 s, 60°C 45 s, 72°C 45 s, and 72°C 10 min. Primer pairs designed from the <italic>F. pratensis </italic>sequence were forward primer 5'TGTCATATCGGGGAACAACA3' and reverse primer 5'ACATGGTTTCGTCCGGATAC3' denaturing at 94°C for 4 minutes and then 40 cycles of 94°C 10 s, 60°C 45 s, and 72°C 1.30 min, and 70°C for 10 min. We also designed a third primer pair, referred to as LpIRIx primer pair, with forward primer 5'GAATGCCGTATCTGGGGACC3' and reverse primer 5'GTGGTTCCCGGATACGGTATT3', based on multiple sequence data acquired from sequencing of the above mentioned genes. This primer pair was used under the same conditions as the LpAFP primers. DNA maxi-preps of the BAC-clones were preformed using the NucloBond BAC 100 Kit (MACHEREY-NAGEL, Düren, Germany). For BAC-sequencing 500 ng BAC-DNA was combined with 20 μM of primer, 8 μl BigDye 3.1 ready mix and dH<sub>2</sub>0, to 20 μl total volume. Following 5 min of denaturising at 95°C, 50 cycles were performed with 30 s at 95°C, 10 s 50°C, and 4 min 60°C. Subsequently, the sequencing reactions were precipitated and sequenced on an ABI PRISM 3100 (Applied Biosystems, Foster City, CA, USA).</p> <title>Protein domain characterisation</title> <p>Predicted peptide LRR domains were characterised using Pfam [##REF##16381856##46##]. We verified the Pfam results by visual inspection of the sequences defining a LRR motif as LxxLxLxx, or variations of it where L is substituted with I, V, or A. To track the molecular evolution of the LRR-domain, <italic>OsLRR-PSR </italic>was used as a template for comparison to the predicted domains of IRI-like sequences. LRR motifs predicted by Pfam were considered significant if the Pfam E-value was lower than 0.05. IRI-like amino acid sequences were aligned with the LRR-domain of OsLRR-PSR and the LRR motifs in IRI-like sequences were named according to which of the LRR motif number in OsLRR-PSR they aligned to. IRI-domain characterisation was performed by visual scoring of the total number of repeat motifs (NxVxG/NxVxxG). IRI-domain repeat motifs were considered as "present", and counted, when they contained no more than one amino acid substitution compared to the consensus motifs. Signal peptide domains were predicted by TargetP [##REF##9051728##47##].</p> <title>Estimation of substitution rates</title> <p>To estimate divergence times between putatively paralogous and orthologous sequences we used the average dS obtained from three different methods, Nei &amp; Gojobori [##REF##3444411##48##], Kumar [##UREF##8##49##], and the Li-Wu-Luo [##REF##3916709##50##] method, in MEGA (4.0) [##REF##17488738##51##]. As a control for evolutionary rates of IRI-like genes we calculated the average dS values of ten randomly selected orthologs from wheat, perennial ryegrass, and rice. Maximum likelihood estimation of non-synonymous to synonymous substitution ratios (w) was performed using Codeml in the PAML software package (v 3.15) [##UREF##9##52##]. The 3 × 4 codon substitution model was chosen for Codeml w estimations. PAL2NAL [##REF##16845082##53##] was used to make codon based nucleotide alignments for the use in MEGA and PAML. The absolute time of divergence between orthologs and paralogs was estimated using a rate of 6.5*10<sup>-9 </sup>substitutions/synonymous site/year for grasses [##UREF##2##23##]. For estimation of mutation rates and absolute divergence times we used the relationship k = dS/2T, where k is the absolute rate of synonymous substitutions per year, dS is the synonymous substitution rate, and T is the absolute time since divergence. To identify putative alleles not grouped in the same contig due to methodological errors, a cut-off threshold of dS = 0.03 was used. This threshold was set on the basis of average inter-allelic dS for LRR-domains of 27 disease resistance like genes in <italic>A. thaliana </italic>[##REF##16798885##54##], and inter-allelic dS<sub>max </sub>of <italic>LpIRI1 </italic>calculated from twelve European perennial ryegrass genotypes (dS<sub>max </sub>= 0.015, data not published).</p> <title>Molecular and phylogenetic analysis</title> <p>All amino acid and nucleotide alignments were made by MAFFT [##REF##15661851##55##] and manually edited in BioEdit [##UREF##10##56##], and the phylogenetic trees were constructed in Treefinder [##REF##15222900##57##]. An AIC criteria test [##UREF##11##58##], implemented in the Modeltest option in Treefinder, was used to choose substitution model for the phylogenetic analysis. ML trees were bootstrapped with 1000 replicates. Synonymous distance based trees were inferred by UPGMA from a pairwise dS distance matrix in MEGA. Alignment figures were prepared by BoxShade <ext-link ext-link-type="uri" xlink:href="http://www.ch.embnet.org/software/BOX_faq.html">http://www.ch.embnet.org/software/BOX_faq.html</ext-link>.</p>
<title>Results</title> <title>Screening of perennial ryegrass BAC libraries</title> <p>Initial screening of the perennial ryegrass (<italic>Lolium perenne</italic>) LTS18 BAC library with the <italic>LpAFP </italic>primer pair produced two hits, from which <italic>LpIRI1 </italic>(<ext-link ext-link-type="gen" xlink:href="EU680848">EU680848</ext-link>) and <italic>LpIRI2 </italic>(<ext-link ext-link-type="gen" xlink:href="EU680849">EU680849</ext-link>) were isolated. The NV#20F1-30 BAC library produced four hits with the LpIRIx primer pair, and three hits with the LpAFP primer pair. <italic>LpIRI4 </italic>(<ext-link ext-link-type="gen" xlink:href="EU680851">EU680851</ext-link>) was subsequently isolated from one of the four positive LpIRIx hits and <italic>LpIRI3 </italic>(<ext-link ext-link-type="gen" xlink:href="EU680850">EU680850</ext-link>) was isolated from one of three positive LpAFP hits. All genes isolated from perennial ryegrass were intronless and encoded putative peptides with high identity to the wheat <italic>TaIRI1 </italic>and <italic>TaIRI2 </italic>genes (blastp &lt; 4e-10). <italic>LpIRI1</italic>, <italic>LpIRI4</italic>, and <italic>LpIRI3 </italic>were similar in size and encoded peptides of 285, 242, and 254 amino acids, respectively. <italic>LpIRI2 </italic>encoded a shorter ORF of 150 amino acids that was 94% identical to the LpIRI4 IRI-domain. The IRI-domain of LpIRI3 is identical to an earlier identified partial IRI peptide encoded by <italic>LpAFP </italic>(<ext-link ext-link-type="gen" xlink:href="AJ277399">AJ277399</ext-link>).</p> <p>Nucleotide alignments of the LpIRI-like genes showed that <italic>LpIRI2 </italic>has undergone a deletion of almost the entire LRR-domain coding region, the only remains of it being a 102 base pair (bp) fragment upstream of the <italic>LpIRI2 </italic>putative start codon. This could indicate that <italic>LpIRI2 </italic>is a pseudogene or a non-functional allele. Non-functional sequences are expected to evolve under neutral expectation, which means that the rate of non-synonymous to synonymous substitutions (w) is expected to be 1. Average w between <italic>LpIRI2 </italic>and the other perennial ryegrass sequences was estimated to be 0.56 which suggests that <italic>LpIRI2 </italic>is under selective constraints despite the major deletion in the LRR-domain.</p> <title><italic>In silico </italic>identification of IRI-like sequences</title> <p>The blastn EST search resulted in 189 wheat, 100 barley, 21 tall fescue (<italic>F. arundinacea</italic>), 5 Italian ryegrass (<italic>L. multiflorum</italic>) and 2 darnel ryegrass (<italic>L. temulentum</italic>) sequences. <italic>In silico </italic>mining produced from zero to eleven full length IRI-like sequences (i.e. with a start and stop codon) per species (Table ##TAB##0##1##), and sequences were annotated as follows; an initial two letters indicating latin species name, "C" indicating a contig of more than two ESTs, and lastly an identifier number. The number of ESTs in the contigs ranged from 2 to 40 (Table ##TAB##1##2##), with an average of 17 ESTs per contig. In addition we identified several partial IRI-like sequences; one partial wheat IRI-like contig of 10 ESTs, and three partial barley sequences of 2, 3, and 6 EST sequences. The partial sequences were not included in further analysis. A barley mRNA, <ext-link ext-link-type="gen" xlink:href="AK249041">AK249041</ext-link>, did not align to any barley contig but coded for a full length IRI-like ORF, hence we included this mRNA in our dataset. <italic>In silico </italic>mining with tall fescue, Italian ryegrass, and darnel ryegrass ESTs did not produce any full length IRI-like sequences. <italic>In silico </italic>EST mined sequences were considered validated if a core nucleotide sequence from NCBI had &gt; 99% identity over &gt; 200 bp. TaC3 was validated as being identical to <italic>TaIRI1</italic>, and HvC1 was validated by being identical to a full length mRNA sequence (<ext-link ext-link-type="gen" xlink:href="AK252915">AK252915</ext-link>). Accession numbers for ESTs belonging to full length IRI-like contigs acquired through <italic>in silico </italic>mining are included in additional material [see Additional file ##SUPPL##0##1##].</p> <title>Predicted protein structure characterisation</title> <p>The <italic>OsLRR-PSR </italic>peptide sequence was used as a template for comparisons of the LRR-domains. The <italic>OsLRR-PSR </italic>is suitable for this purpose because it is a putative homologue to the LRR-domain of the IRI-like sequences. In the following peptide structure characterization we assume that the most recent common ancestor (MRCA) of all IRI-like grass genes encoded the same domain architecture as <italic>OsLRR-PSR</italic>. <italic>OsLRR-PSR </italic>was predicted to encode one LRR N-terminal domain (LRR-NT), 21 internal LRR homologous motifs (1–21), and a protein kinase (PK) C-terminal domain. Only the LRR-NT and six internal LRR motifs were predicted with significant E-values by Pfam. Compared to <italic>OsLRR-PSR </italic>all predicted IRI-like peptides have reduced number of LRR motifs and lacks the PK-domain. Three internal LRR motifs did not align to any OsLRR-PSR LRR motifs. These have probably arisen through deletions in the LRR-domain causing two partial LRR motifs to merge into a novel LRR motif. Comparative sequence analysis suggests that deletions between LRR motifs 1–16, 2–16, and 12–16 of a <italic>OsLRR-PSR</italic>-like grass ancestor gene have resulted in three novel LRR-domains; LRR1b, 2b, and 16b (Fig. ##FIG##0##1##).</p> <p>The IRI-domain also varies in size by number of repeated motifs (Fig. ##FIG##0##1##). About 60% of all sequences with an IRI-domain have 15 repeat motifs or more. Six sequences were detected to have a reduced number of repeat motifs, or had completely lost the IRI-domain. Analysis of codon based nucleotide alignments revealed that frameshift (FS) mutations could be identified in four of the IRI-like sequences (TaC3, TaC10, TaC11, and <ext-link ext-link-type="gen" xlink:href="AK249041">AK249041</ext-link>) that showed reduced IRI-domain size (data not shown). HvC3 is the only IRI-like sequence with a completely reduced IRI-domain. For all sequences in the HvC3 contig additional information on abiotic conditions under which the plants had been grown were included in the EST files. Without exception all ESTs originated from tissue sampled from etiolated barley seedlings, and not from cold acclimated tissue. This is congruent with HvC3 lacking the entire ice binding IRI-domain, suggesting that IRI-like paralogs are involved in several different stress responses.</p> <p>Prediction of the subcellular location of the IRI-like peptides (see methods) predicted a signal peptide that targets the peptides to the secretory pathway present in all IRI-like peptides, except from LpIRI2. The lack of an LpIRI2 signal peptide, and the fact that LpIRI2 has undergone a deletion of almost the entire LRR-domain could suggest that <italic>LpIRI2 </italic>is in fact a non-functional allele or pseudogene. However the results from the w estimates contradict the non-functionality hypothesis. Alternatively the lack of a signal peptide can be interpreted as that LpIRI2 simply has evolved a different function than the IRI-like peptides with a conserved signal peptide.</p> <title>Phylogenetic analysis of IRI-like paralogs</title> <p><italic>OsLRR-PSR </italic>has highest homology to IRI-like genes outside the Pooideae subfamily (blastp E-value of 2e<sup>-46</sup>) and was therefore used as an out-group in the phylogenetic analysis. The paralog phylogenies of perennial ryegrass, barley, and wheat IRI-like peptide sequences are given in Figure ##FIG##1##2A, 2B##, and ##FIG##1##2C##, respectively. In the wheat phylogeny TaIRI2, the amino acid sequence with the highest number of conserved LRR-domains, diverge as a monophyletic branch, while all other wheat sequences with further reduced LRR-domain sizes form a second monophyletic group. This large monophyletic group can be subdivided into two smaller monophyletic clades (Fig. ##FIG##1##2C##). In group I a FS mutation in the IRI-domain of TaC3 shortens the ORF and separates TaC3 structurally from the other sequences. One might be tempted to speculate if this FS is due to alignment errors but the FS in the TaC3 contig have been validated independently by the full length mRNA <italic>TaIRI1</italic>. Group II contain IRI-like genes with high degree of divergent predicted peptide structure, i.e. putative peptides with strongly reduced and completely lost IRI-domain (TaC9 and TaC11) as well as predicted peptides with large IRI-domains. Alignments of predicted amino acid sequences with all FS mutations removed, that were used for the phylogenetic analysis, are presented in additional material [see Additional file ##SUPPL##1##2##, ##SUPPL##2##3##, and ##SUPPL##3##4##].</p> <title>Estimation of synonymous divergence of IRI-like sequences</title> <p>If we assume a molecular clock, synonymous substitution rates (dS) between two DNA sequences can be interpreted as a relative measurement of time since MRCA, thus for two paralogous genes dS can be interpreted as the time since gene duplication [##REF##12566392##19##]. Without being able to account for all the IRI-like paralogs existing in a genome we cannot infer if two paralogs descend from a single duplication event (i.e. being true paralogs) or if they are products of two separate duplication events. We therefore restricted our initial analysis of IRI-like gene duplication events to only comprise the dS<sub>max </sub>and dS<sub>min </sub>for all pairwise comparisons. The dS<sub>max</sub>-dS<sub>min </sub>range can be interpreted as the evolutionary time span in which all duplications of IRI-like genes in our dataset have occurred.</p> <p>The average dS between <italic>OsLRR-PSR </italic>and LRR-domains of all IRI-like sequences from cold tolerant grasses was estimated to 0.97 (standard deviation (SD) = 0.11). Hence if <italic>OsLRR-PSR </italic>is the true ortholog of IRI-like genes, the radiation of cold tolerant grasses from rice predates the initiation of IRI-like gene family expansion in our dataset (Fig. ##FIG##2##3A##). All paralog dS ranges for wheat, barley and perennial ryegrass overlapped, however there are large range differences. These range differences is caused by the low dS<sub>min </sub>range of wheat and perennial ryegrass IRI-like sequences (Fig. ##FIG##2##3A##). For perennial ryegrass the dS between <italic>LpIRI2 </italic>and <italic>LpIRI4 </italic>lowered the dS<sub>min </sub>estimate from 0.40 to 0.03. Dissection of wheat paralog dS estimates showed that the low wheat dS<sub>min </sub>range could be traced back to low pairwise dS within the monophyletic clades I and II (Fig. ##FIG##1##2C##). If dS estimates from sequence pairs within clade I and II were removed the wheat dS<sub>min </sub>shifted from 0.01 to 0.39.</p> <p>Very low dS between paralog pairs might reflect the inclusion of highly diverged alleles in our paralog dataset. Sequencing errors in ESTs and inclusion of highly diverged genotypes in our dataset could potentially give an inflated polymorphism level producing artificial contigs that are alleles rather than paralogs. To identify putative false paralogs we set an allelic dS threshold of dS &lt; 0.03 (see methods section). Based on this definition we identified two putative allelic sequence pairs TaC2-TaC11 (dS = 0.01) and <italic>LpIRI2</italic>-<italic>LpIRI4 </italic>(dS = 0.03).</p> <title>Estimation of synonymous divergence between control genes</title> <p>The ten control orthologs chosen for evolutionary rate control and their pairwise dS estimates are listed in Table ##TAB##2##3##. Based on these genes the average dS between rice and cold tolerant grasses were estimated to be 0.54 with an average SD between dS estimation methods of 0.08. The estimated divergence of rice and cold tolerant grasses based on the control genes (dS = 0.54) is substantially lower than the minimum equivalent estimate found by using <italic>OsLRR-PSR </italic>and IRI-like genes (dS = 0.72) (Fig. ##FIG##2##3A## and ##FIG##2##3B##). The average dS between perennial ryegrass and wheat genes were 0.32, with an average SD of 0.04 between the estimation methods.</p> <title>Molecular analysis of LRR and IRI-domains</title> <p>The alignment of the LRR-domain of <italic>OsLRR-PSR </italic>and the IRI-domains of a subset of the IRI-like peptide sequences shows two blocks of strikingly high conservation between the OsLRR-PSR and two IRI-domain motifs (Fig. ##FIG##3##4##). The two motifs have 2/3 of the motif sequence residues conserved between rice and cold tolerant grasses. The first motif has an Ile-Val substitution while the second motif has a Leu-Val substitution, neither of which are substitutions with large effect on hydrophilic properties.</p> <p>Relatively low level of conservation between DcAFP and grass IRI-like peptide sequences was found (blastp 1e-4). All sequences in the LRR-domain alignment (Fig. ##FIG##4##5##) have a blast E-value of &lt; 1e-18 to at least one other sequence in the alignment, but no larger blocks of conservation between DcAFP and any monocot LRR-domain can be identified. However we observe that several blocks of 2–4 conserved residues exist throughout the alignment of the LRR-domain between DcAFP and grass IRI-sequences.</p>
<title>Discussion</title> <p>Until now only three cold tolerant grass IRI protein coding genes have been reported; a partial coding sequence of an IRI-domain from perennial ryegrass [##REF##10917518##14##], and two highly identical full length mRNA paralogs from wheat [##REF##15792959##10##]. Through <italic>in silico </italic>mining and BAC sequencing we have identified 15 full length and 8 partial novel IRI-like genes in cold tolerant grasses. In addition, we have obtained the complete sequence of <italic>LpAFP</italic>. The data accumulated leaves no doubt: cold tolerant grasses of the Pooideae subfamily have evolved a lineage specific family of IRI-like genes.</p> <title>IRI-gene family radiation happened after the cold tolerant grass divergence</title> <p>The prevailing hypothesis on the evolution of LRR-IRI-like genes belonging to cold tolerant grasses is that they are lineage specific and that an <italic>OsLRR-PSR</italic>-like gene is the MRCA [##REF##15792959##10##]. This hypothesis was proposed based on sequence homology data only and we therefore re-examined this idea using more rigorous statistical methods by estimating synonymous divergence. When employing a commonly used mutation rate for grasses of 6.5*10<sup>-9 </sup>[##UREF##1##20##, ####REF##16381819##21##, ##REF##15208399##22####15208399##22##], estimated by Gaut et al. [##UREF##2##23##], the synonymous divergence level between <italic>OsLRR-PSR </italic>and IRI-like sequences suggested a MRCA about 75 Mya. This is slightly higher than upper thresholds of some published rice-Pooideae divergence estimates [##REF##15208399##22##]. However, our estimate of divergence time between rice and cold tolerant grasses based on the control genes suggest a rice-Pooideae divergence only 42 Mya. This is similar to divergence estimates published by Patterson et al. [##UREF##3##24##] and Salse et al. [##UREF##1##20##], dating back 41–47 and ~46 Mya, respectively. Our wheat and perennial ryegrass divergence estimates is dated ~10 million years prior to a previously published estimate of ~35 Mya [##UREF##4##25##].</p> <p>The observed discrepancy between the two divergence estimates of rice and cold tolerant grasses in our study (Fig. ##FIG##2##3A## and ##FIG##2##3B##) can be interpreted in two different ways. The <italic>OsLRR-PSR </italic>gene is the true ortholog of IRI-like genes and the incongruent divergence time estimates are caused by differences in molecular clock rates. Or alternatively, if the molecular clock rate is similar, it follows that <italic>OsLRR-PSR </italic>diverged from IRI-like genes long before rice and cold tolerant grasses diverged (Fig. ##FIG##2##3##). The burst of IRI-like sequence duplications must then have occurred in the ancestor genome of rice, and this implies that the rice genome subsequently must have lost all genes belonging to the IRI-like gene family. Even though loss of genes and whole gene families is not an uncommon feature of plant genome evolution [##REF##18079367##16##,##UREF##3##24##], elevated clock rate differences is a more parsimonious explanation to the divergence estimate differences seen in Figure ##FIG##2##3A## and ##FIG##2##3B##. Evolutionary rate differences are highly common among closely related species, different lineages of a species, and also within a genome [##REF##18056075##26##,##REF##17451608##27##]. When using a single gene family to estimate divergence between species, as with the IRI-like genes, deviation from the average genome clock rate would be expected. As an example the clock rate of the ten control genes varied from 4.1–7.1*10<sup>-9</sup>, with an average rate of 5.4*10<sup>-9</sup>, when using a divergence time between rice-cold tolerant grasses of 50 My.</p> <p>Assuming true orthologous relationship between <italic>OsLRR-PSR </italic>and IRI-like genes we can calibrate an average molecular clock rate for IRI-like genes using the dS = 0.97 and assuming an absolute divergence time of 50 Mya (see methods). This gives us an estimate of an IRI-like gene family specific clock rate of 9.7*10<sup>-9</sup>. Employing this adjusted clock rate pushes the estimates for the initiation of IRI-like gene duplications forward to 36, 27, and 39 Mya for wheat, perennial ryegrass and barley, respectively. This is approximately 3–14 My after our estimate for divergence between rice and cold tolerant grasses based on the control genes.</p> <title>Species specific differences in IRI-like sequence numbers</title> <p>Twice as many IRI-like sequences, partial and full length, were identified in wheat compared to barley (Table ##TAB##0##1##). <italic>In silico </italic>mining is vulnerable to methodologically introduced uncertainties. For example, the fact that the wheat EST database at NCBI is more than twice as large (1.2 M ESTs) than for barley (500 K ESTs) could be a contributing factor to the differences in numbers of IRI-like mined sequences because we expect that the EST database size is positively correlated to transcriptome coverage of an organism. A separate effect of a larger EST database will be the inclusion of ESTs from an increased number of genotypes, which could be a source of introduction of allelic polymorphisms.</p> <p>Even though methodological properties might elevate the number of wheat sequences identified to some extent, we believe that much of the difference in wheat and barley IRI-like sequence numbers are related to genomic ploidy level differences. Wheat (<italic>Triticum aestivum</italic>) is an allo-hexaploid originating about 8.000 years ago. It has three homoeologous genomes A, B, and D, which are estimated to have diverged 4.5–2.5 Mya [##UREF##4##25##]. Our results from the phylogenetic analyses, supported by the dS estimates, suggests that IRI-like sequences within monophyletic clade I and II (Fig. ##FIG##1##2##) could be homoeologous rather than paralogous. But low pairwise dS can alternatively reflect recent gene duplications. Consequently our inferences on the evolutionary relationship between putative wheat homoeologous sequences must be viewed in a critical manner.</p> <title>A model of wheat IRI-like sequence evolution</title> <p>As mentioned, wheat differs from the diploid grasses in our study in that it has many more and younger IRI-like paralogs, some of which probably can be accounted for by wheat polyploidisation. Using the IRI gene specific clock rate of 9.7*10<sup>-9 </sup>we built a model of the evolutionary history of wheat IRI-sequences using a phylogenetic approach (Fig. ##FIG##5##6A##). According to literature and the data acquired in this study wheat and barley specific duplications would be younger than ~35-25 My and wheat specific duplications would be younger than ~11 My [##UREF##4##25##]. The internal nodes of the phylogeny in Fig. ##FIG##5##6## can therefore be interpreted as putative ancient duplication events happening before divergence within Pooideae (D<sub>A</sub>), putative duplication events specific to wheat and barley (D<sub>WB</sub>), putative speciation event between wheat and barley (S<sub>WB</sub>), putative gene duplication events specific to wheat (D<sub>W</sub>), speciation events of the A, B, and D genomes of wheat (S<sub>ABD</sub>), and lastly allelic divergence (A).</p> <p>The evolutionary model with adjusted clock rates strengthens our hypothesis on the homoeologous relationship of the wheat clade I and II sequences (Fig. ##FIG##1##2##). The divergence between TaC6 and TaC2/TaC11 are predicted slightly earlier than the polyploidisation event, thus this internal node could not be classified unambiguously. Sub-trees of clade I and II, in which we have included the barley sequence with lowest synonymous distance to each clade, are presented in Figure ##FIG##5##6B## and ##FIG##5##6C##, respectively. The sub-trees further support the hypothesis that clade I and II represents genes that mainly have arisen through polyploidisation events. In both sub-trees the closest related barley sequence is estimated to have diverged from the wheat clade I and II about 15 Mya, coinciding with wheat-barley divergence [##UREF##4##25##]. However without a complete knowledge of the orthologous relationships of the IRI-like sequences the inferences on evolutionary relationships are somewhat speculative.</p> <title>Structural and functional diversification of the IRI-like gene family</title> <p>A striking feature of the IRI-like gene family is the structural differentiation between paralogs (Fig. ##FIG##0##1##). Structural diversification of IRI-like genes, as seen in our sequence collection, would be expected to affect the spectrum of IRI-like peptide function, because both LRR and IRI-domains are known to be involved in substrate binding [##REF##11721016##15##,##REF##10973092##28##]. One interpretation of this pattern is that IRI-like sequences with complementary combinations of LRR motifs and IRI-domain sizes are selected for and retained in the genome, which is what we expect from the duplication-degeneration-complementation (DDC) model of paralog evolution [##REF##10101175##18##]. DDC predicts that mutations in regulatory elements increase the probability of paralog retention because it leads to partitioning of ancestral functions (subfunctionalisation), and the model has proven to be an important contribution to understanding evolution of paralogous genes [##REF##15122255##29##,##REF##15785853##30##]. The DDC model has later been expanded to coding sequences [##UREF##5##31##,##REF##17418444##32##], and recently a combination of regulatory and structural DDC has been demonstrated [##REF##17956721##33##,##REF##17928853##34##].</p> <p>Regulatory subfunctionalisation in gene expression and tissue localisation has been demonstrated between <italic>TaIRI1 </italic>and <italic>TaIRI2 </italic>(TaC3) [##REF##15792959##10##], two genes coding for highly divergent LRR-domain structure and length. In our study we have also found evidence that peptide structure divergence has led to sub- or neofunctionalisation. A barley IRI-like sequence contig (HvC3) with no IRI-domain still seems to play a functional role under etiolation. This suggests that the LRR-domain of IRI-like genes may play a functional role in multiple stress responses. Other LRR-domain containing genes in plants have also been shown to be involved in stress responses under drought stress and as a key membrane-bound regulator of absiscic acid signalling [##REF##15772289##35##,##REF##15144382##36##].</p> <p>One interesting aspect of the structural divergence pattern is that all genes except <italic>LpIRI2 </italic>are predicted to encode a conserved N-terminus signalling domain targeting the proteins to the secretory pathway. Secretion to the apoplast is expected for proteins with ice interacting functions. In the light of these data, an interpretation of the structural variability of LRR-domains, combined with the apparent conservation of the N-terminal signalling domain, is that IRI-like genes might be under selective pressure for a continuous ORF from the signalling domain across the LRR-domain and into the IRI-domain, conserving the crucial function of apoplast export of IRI peptides. The LRR-domain itself might not be under functional conservation. As an example: the full length sequenced mRNA <ext-link ext-link-type="gen" xlink:href="AK249041">AK249041</ext-link> from barley has a N-terminal conserved predicted signal peptide motif, a completely reduced LRR-domain with no predicted LRR motifs, and an IRI-domain.</p> <p>Less dramatic polymorphisms between paralogs, such as single amino acid substitutions or small motif number differences could potentially have a large effect on the functionality. Single amino acid substitutions have been shown to radically change AFP functionality in both plant and animal AFPs [##REF##14531728##13##,##REF##9688560##37##]. Chakrabartty and co-workers [##REF##2738068##38##] showed that only small deletions in an AFP with repetitive structure from flounder altered the ice interacting properties dramatically. Thus, all the observed polymorphisms between IRI paralogs, even down to single amino acid substitutions, could potentially be of functional significance.</p> <title>Birth of an IRI repeat domain</title> <p>The molecular mechanisms underlying the metamorphosis from an <italic>OsLRR-PSR</italic>-like ancestor gene into the first bipartite IRI-like gene have been addressed by Trembley et al. [##REF##15792959##10##]. They proposed the "transposable element hypothesis" (TE) suggesting that the IRI-domain is a TE insertion that has resulted in a FS mutation and caused the loss of the PK-domain. However no TE signatures were found flanking the IRI-domain [##REF##15792959##10##]. Based on results from our sequence analysis (Fig. ##FIG##3##4##) we propose a competing hypothesis on the evolution of the IRI-domain, namely the repeated motif expansion (RME) hypothesis. It has been shown that expansions of domains by duplication of repeated motifs are common in genes of repetitive structure [##REF##17573804##39##]. We suggest that IRI motifs have increased in copy number by a yet unknown mechanism, possibly illegitimate recombination, slippage, or uneven crossing over. Contrary to the TE-hypothesis the RME hypothesis can explain the evolution of the IRI-domain and at the same time account for the existence of two IRI motif-like blocks in <italic>OsLRR-PSR </italic>(Fig ##FIG##3##4##). Lastly, if the entire IRI-domain is a TE-insertion we would expect this TE sequence to be found at other loci in grasses. However no such reports of TE-like sequences homologous to the IRI-domain are known to our knowledge.</p> <title>Convergent evolution of LRR containing AFPs</title> <p>LRR-domain containing proteins are extremely abundant in plants. The largest LRR containing plant peptide group is LRR receptor kinases (LRR-RK), having more than 200 members in the <italic>A. thaliana </italic>genome [##UREF##6##40##]. Plant disease resistance associated genes (NBS-LRR) comprise another large LRR containing functional group [##REF##15231261##41##]. Common for the function of LRR domains in any peptide is that they are associated with peptide-peptide recognition and binding interactions [##REF##18202283##42##, ####REF##9437864##43##, ##REF##9811792##44####9811792##44##].</p> <p>Through comparative protein domain analysis we have shown that LRR-domains of IRI-like genes are much less conserved compared to the predicted signal peptide motif flanking the N-terminus of the LRR-domains. We believe that this could be due to lack of selective constrains on the LRR-domain function itself, or perhaps selection for divergent LRR-domain functions as predicted by DDC. Whatever functional role today's IRI-like sequence LRR-domains might play; there is little doubt that the LRR-domain of IRI-like genes in cold tolerant grasses shares an ancient common ancestor with the LRR-domain of DcAFP (Fig. ##FIG##4##5##). However, while cold tolerant grass IRI-like proteins have evolved ice binding capacity through the evolution of an IRI-domain [##REF##10917518##14##,##REF##11721016##15##], DcAFP have evolved ice binding capacity through changes in the LRR-domain itself [##REF##14531728##13##]. <italic>DcAFP </italic>and grass LRR-IRI genes are therefore intriguing examples of parallel evolution of function by two completely different molecular mechanisms; evolutionary alterations of a pre-existing LRR-domain and evolution of a novel repeat domain with ice binding properties.</p>
<title>Conclusion</title> <p>The IRI-like genes identified by Sidebottom et al. [##REF##10917518##14##] and Tremblay et al. [##REF##15792959##10##], and in this study tell a tale of a complex evolutionary history that includes birth of an ice binding domain, a burst of gene duplication events after cold tolerant grasses radiated from rice, domain structure differentiation between paralogs, and sub- and/or neofunctionalisation of IRI-like proteins. Given more detailed functional studies, the IRI-like gene family can provide a valuable example of how duplicated genes evolve novel functional spectres. The hypothesis that evolution of IRI-like genes has been important for Pooideae grass adaptation to cold climate [##REF##15792959##10##] is strengthened by this study as we show that the evolution of the IRI-like gene family probably happened after the divergence from rice, and furthermore that the numbers of IRI-like genes are higher than earlier known.</p>
<title>Background</title> <p>Grasses are adapted to a wide range of climatic conditions. Species of the subfamily Pooideae, which includes wheat, barley and important forage grasses, have evolved extreme frost tolerance. A class of ice binding proteins that inhibit ice re-crystallisation, specific to the Pooideae subfamily lineage, have been identified in perennial ryegrass and wheat, and these proteins are thought to have evolved from a leucine-rich repeat phytosulfokine receptor kinase (<italic>LRR-PSR</italic>)-like ancestor gene. Even though the ice re-crystallisation inhibition function of these proteins has been studied extensively <italic>in vitro</italic>, little is known about the evolution of these genes on the molecular level.</p> <title>Results</title> <p>We identified 15 putative novel ice re-crystallisation inhibition (IRI)-like protein coding genes in perennial ryegrass, barley, and wheat. Using synonymous divergence estimates we reconstructed the evolution of the IRI-like gene family. We also explored the hypothesis that the IRI-domain has evolved through repeated motif expansion and investigated the evolutionary relationship between a LRR-domain containing IRI coding gene in carrot and the Pooideae IRI-like genes. Our analysis showed that the main expansion of the IRI-gene family happened ~36 million years ago (Mya). In addition to IRI-like paralogs, wheat contained several sequences that likely were products of polyploidisation events (homoeologs). Through sequence analysis we identified two short motifs in the rice <italic>LRR-PSR </italic>gene highly similar to the repeat motifs of the IRI-domain in cold tolerant grasses. Finally we show that the LRR-domain of carrot and grass IRI proteins both share homology to an <italic>Arabidopsis thaliana </italic>LRR-trans membrane protein kinase (<italic>LRR-TPK</italic>).</p> <title>Conclusion</title> <p>The diverse IRI-like genes identified in this study tell a tale of a complex evolutionary history including birth of an ice binding domain, a burst of gene duplication events after cold tolerant grasses radiated from rice, protein domain structure differentiation between paralogs, and sub- and/or neofunctionalisation of IRI-like proteins. From our sequence analysis we provide evidence for IRI-domain evolution probably occurring through increased copy number of a repeated motif. Finally, we discuss the possibility of parallel evolution of LRR domain containing IRI proteins in carrot and grasses through two completely different molecular adaptations.</p>
<title>Authors' contributions</title> <p>SRS conceived the study, carried out molecular genetics work and molecular analysis, and drafted the manuscript. HR carried out molecular genetics work, helped with data analysis, and participated in manuscript drafting. TA carried out molecular work and helped draft the manuscript. OAR participated in drafting the manuscript.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>We thank Kjetil Fosnes for assistance with BAC sequencing and Dr. Magnus Dehli Vigeland, Dr. Gordon Allison, and anonymous reviewers for valuable comments on the manuscript. This study was funded through the KMB project <italic>Festulolium </italic>with Improved Forage Quality and Winter Survival for Norwegian Farming, project number 173319/I10, funded by the Research Council of Norway and Graminor AS.</p>
<fig id="F1" position="float"><label>Figure 1</label><caption><p><bold>Predicted peptide structure of <italic>OsLRR-PSR </italic>and IRI-like sequences</bold>. Peptide structures of the LRR and IRI-domains of IRI-like sequences as predicted by Pfam and by visual characterisation, respectively. <italic>OsLRR-PSR </italic>was used as a template for peptide structure comparisons. Grey LRR-domains were predicted with significant P-values in Pfam (P &lt; 0.05). Numbers of IRI motifs in the IRI-domains and whether frame shift mutations were identified are given in parentheses.</p></caption></fig> <fig id="F2" position="float"><label>Figure 2</label><caption><p><bold>Maximum Likelihood paralog phylogenies of all full length IRI-like amino acid sequences</bold>. Bootstrap values from 1000 replicates are in percent at internal nodes. Bootstrap values &lt; 50 are not shown. A) Phylogeny of barley IRI-like full length sequences. B) Phylogeny of perennial ryegrass IRI-like full length sequences. C) Phylogeny of wheat IRI-like full length sequences. Two monophyletic clades of very low within-clade synonymous divergence (I, II) are boxed.</p></caption></fig> <fig id="F3" position="float"><label>Figure 3</label><caption><p><bold>Synonymous divergence between paralogous and orthologous sequences</bold>. A) Maximum and minimum dS range between all IRI-like paralogous sequences of wheat, barley and perennial ryegrass and the average dS between <italic>OsLRR-PSR </italic>and all IRI-like sequences. Dotted lines for the divergence estimate between OsLRR-PSR and IRI-like sequences represent SD of the ortholog dS. Dotted lines for wheat and perennial ryegrass indicate changes in dS range if putative alleles are included. B) A tree indicating the average synonymous divergence, given in dS, estimated from ten control genes from rice, wheat, and perennial ryegrass.</p></caption></fig> <fig id="F4" position="float"><label>Figure 4</label><caption><p><bold>Amino acid alignment between IRI-like sequences and <italic>OsLRR-PSR </italic>at the initiation of the IRI-domain</bold>. Shaded residues are identical in &gt; 90% of sequences. Black bars underlines IRI-domain A- and B motifs (NxVxG/NxVxxG). Boxed motifs are IRI-domain motifs shared between OsLRR-PSR and IRI-like sequences.</p></caption></fig> <fig id="F5" position="float"><label>Figure 5</label><caption><p><bold>Amino acid alignment of homologous LRR-domains from monocot and dicot species</bold>. The alignment shows homologous LRR-domains between Arabidopsis (AtLRR-TPK), carrot (DcAFP), rice (OsLRR-PSR), wheat (TaIRI1) and perennial ryegrass (LpIRI1). All sequences in the alignment has a blastp E-value of &gt; 1*10<sup>-18 </sup>to at least one other sequence in the alignment. Shaded residues are identical in &gt; 60% of sequences.</p></caption></fig> <fig id="F6" position="float"><label>Figure 6</label><caption><p><bold>Evolutionary relationships between the wheat IRI-like sequences inferred with an IRI-like gene specific molecular clock rate</bold>. Phylogenetic trees based on pairwise synonymous distances made with Unweighted Pair Group Method with Arithmetic Mean (UPGMA). Distance scale shows synonymous site divergence and estimated absolute divergence time in million years. Internal nodes represent duplication, speciation, or allelic divergence events. D<sub>A </sub>= ancient duplication events shared by wheat, barley, and perennial ryegrass; D<sub>WB </sub>= duplication events shared by wheat and barley; D<sub>W </sub>= duplication events exclusive for the wheat lineage; S<sub>WB </sub>= speciation event between wheat and barley; S<sub>ABD </sub>= speciation event for the A, B, and D genomes of wheat; A = allelic divergence event. A) Tree based on all wheat IRI-like sequences. B) Sub-tree of monophyletic clade II including the two most recently radiated IRI-like sequences in our dataset (HvC1 and TaC4). C) Sub-tree of monophyletic clade I including the single most recently radiated barley IRI-like sequence.</p></caption></fig>
<table-wrap id="T1" position="float"><label>Table 1</label><caption><p>All IRI-like sequences identified through EST <italic>in silico </italic>mining.</p></caption><table frame="hsides" rules="groups"><thead><tr><th align="left">Species</th><th align="center">No. EST sequences</th><th align="center" colspan="2">IRI-like sequences</th></tr><tr><th/><th/><th colspan="2"><hr/></th></tr><tr><th/><th/><th align="center">Full length</th><th align="center">Partial</th></tr></thead><tbody><tr><td align="left">Wheat</td><td align="center">189</td><td align="center">11</td><td align="center">1</td></tr><tr><td align="left">Barley</td><td align="center">100</td><td align="center">3</td><td align="center">3</td></tr><tr><td align="left">Tall fescue</td><td align="center">21</td><td align="center">0</td><td align="center">2</td></tr><tr><td align="left">Italian ryegrass</td><td align="center">5</td><td align="center">0</td><td align="center">1</td></tr><tr><td align="left">Darnel ryegrass</td><td align="center">2</td><td align="center">0</td><td align="center">1</td></tr></tbody></table></table-wrap> <table-wrap id="T2" position="float"><label>Table 2</label><caption><p>Full length IRI-like sequences identified through EST <italic>in silico </italic>mining and the number of ESTs per contig.</p></caption><table frame="hsides" rules="groups"><thead><tr><th align="center">Species</th><th align="center">Contig name</th><th align="center">ESTs in contig</th></tr></thead><tbody><tr><td align="center">Wheat</td><td align="center">TaC1</td><td align="center">16</td></tr><tr><td/><td align="center">TaC2</td><td align="center">39</td></tr><tr><td/><td align="center">TaC3</td><td align="center">19</td></tr><tr><td/><td align="center">TaC4</td><td align="center">21</td></tr><tr><td/><td align="center">TaC5</td><td align="center">16</td></tr><tr><td/><td align="center">TaC6</td><td align="center">12</td></tr><tr><td/><td align="center">TaC7</td><td align="center">16</td></tr><tr><td/><td align="center">TaC8</td><td align="center">17</td></tr><tr><td/><td align="center">TaC9</td><td align="center">3</td></tr><tr><td/><td align="center">TaC10</td><td align="center">4</td></tr><tr><td/><td align="center">TaC11</td><td align="center">2</td></tr><tr><td align="center">Barley</td><td align="center">HvC1</td><td align="center">23</td></tr><tr><td/><td align="center">HvC2</td><td align="center">40</td></tr><tr><td/><td align="center">HvC3</td><td align="center">15</td></tr></tbody></table></table-wrap> <table-wrap id="T3" position="float"><label>Table 3</label><caption><p>Evolutionary rate control genes and their pairwise synonymous distances.</p></caption><table frame="hsides" rules="groups"><thead><tr><th align="center">Gene</th><th align="center" colspan="3">Accession number</th><th align="center" colspan="2">Synonymous distance (dS)</th></tr><tr><th/><th colspan="5"><hr/></th></tr><tr><th/><th/><th/><th/><th align="center">Ta vs Lp</th><th align="center">TaLp vs Os</th></tr><tr><th/><th align="left">Ta</th><th align="left">Lp</th><th align="left">Os</th><th align="center">dS</th><th align="center">dS*</th></tr></thead><tbody><tr><td align="left">Cytosolic glyceraldehyde-3-phosphate dehydrogenase</td><td align="left"><ext-link ext-link-type="gen" xlink:href="EF592180">EF592180</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="EF463063">EF463063</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001059674">NM_001059674</ext-link></td><td align="center">0.33 (0.04)</td><td align="center">0.47 (0.06)</td></tr><tr><td align="left">Actin</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB181991">AB181991</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AY014279">AY014279</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001062196">NM_001062196</ext-link></td><td align="center">0.30 (0.03)</td><td align="center">0.49 (0.07)</td></tr><tr><td align="left">Gibberellin 20-oxidase</td><td align="left"><ext-link ext-link-type="gen" xlink:href="Y14008">Y14008</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AY014281">AY014281</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001058486">NM_001058486</ext-link></td><td align="center">0.37 (0.02)</td><td align="center">0.70 (0.13)</td></tr><tr><td align="left">Phytochrome B</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AF137331">AF137331</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AF137308">AF137308</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001056445">NM_001056445</ext-link></td><td align="center">0.39 (0.04)</td><td align="center">0.57 (0.07)</td></tr><tr><td align="left">Casein protein kinase 2 alpha subunit</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB052133">AB052133</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB213317">AB213317</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001065287">NM_001065287</ext-link></td><td align="center">0.30 (0.05)</td><td align="center">0.48 (0.09)</td></tr><tr><td align="left">Na+/H+ antiporter precursor</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AY461512">AY461512</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AY987047">AY987047</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001074903">NM_001074903</ext-link></td><td align="center">0.33 (0.04)</td><td align="center">0.58 (0.09)</td></tr><tr><td align="left">Putative plasma membrane Na+/H+ antiporter</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AY326952">AY326952</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AY987046">AY987046</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="CB634542">CB634542</ext-link></td><td align="center">0.14 (0.02)</td><td align="center">0.44 (0.06)</td></tr><tr><td align="left">Myo-inositol phosphate synthase</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AF542968">AF542968</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AY154382">AY154382</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001055777">NM_001055777</ext-link></td><td align="center">0.32 (0.04)</td><td align="center">0.57 (0.07)</td></tr><tr><td align="left">Cinnamoyl CoA reductase</td><td align="left"><ext-link ext-link-type="gen" xlink:href="CK161291">CK161291</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AF010290">AF010290</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001052667">NM_001052667</ext-link></td><td align="center">0.28 (0.01)</td><td align="center">0.41 (0.03)</td></tr><tr><td align="left">Fructan beta-(2,1) fructosidase</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AJ564996">AJ564996</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="DQ016297">DQ016297</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001052039">NM_001052039</ext-link></td><td align="center">0.39 (0.06)</td><td align="center">0.72 (0.11)</td></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p>EST accession numbers. Table giving all EST accession numbers included in the full length IRI-like <italic>in silico </italic>mined sequences.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional file 2</title><p>Amino acid alignment of perennial ryegrass IRI-like sequences. Amino acid alignment of perennial ryegrass IRI-like sequences used for phylogenetic analysis.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S3"><caption><title>Additional file 3</title><p>Amino acid alignment of wheat IRI-like sequences. Amino acid alignment of wheat IRI-like sequences used for phylogenetic analysis.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S4"><caption><title>Additional file 4</title><p>Amino acid alignment of barley IRI-like sequences. Amino acid alignment of barley IRI-like sequences used for phylogenetic analysis.</p></caption></supplementary-material>
<table-wrap-foot><p>IRI-like gene homologous ESTs were identified by a blastn search using <italic>TaIRI1</italic>.</p></table-wrap-foot> <table-wrap-foot><p>Accession numbers are given for putative orthologous genes from wheat (Ta), perennial ryegrass (Lp) and rice (Os) used for evolutionary rate control. Their synonymous pairwise divergence was calculated as a mean of three estimates (see methods section), and the standard deviation of all three estimates is given in parenthesis. *denotes a mean synonymous distance of two pairwise comparisons between perennial ryegrass-rice and wheat-rice (referred to as TaLp-Os).</p></table-wrap-foot>
<graphic xlink:href="1471-2148-8-245-1"/> <graphic xlink:href="1471-2148-8-245-2"/> <graphic xlink:href="1471-2148-8-245-3"/> <graphic xlink:href="1471-2148-8-245-4"/> <graphic xlink:href="1471-2148-8-245-5"/> <graphic xlink:href="1471-2148-8-245-6"/>
<media xlink:href="1471-2148-8-245-S1.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2148-8-245-S2.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2148-8-245-S3.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2148-8-245-S4.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Barrett"], "given-names": ["J"], "article-title": ["Thermal hysteresis proteins"], "source": ["The International Journal of Biochemistry & Cell Biology"], "year": ["2001"], "volume": ["33"], "issue": ["2"], "fpage": ["105"], "lpage": ["117"], "pub-id": ["10.1016/S1357-2725(00)00083-2"]}, {"surname": ["Salse", "Bolot", "Throude", "Jouffe", "Piegu", "Masood Quraishi", "Calcagno", "Cooke", "Delseny", "Feuillet"], "given-names": ["J", "S", "M", "V", "B", "U", "T", "R", "M", "C"], "article-title": ["Identification and Characterization of Shared Duplications between Rice and Wheat Provide New Insight into Grass Genome Evolution"], "source": ["Plant Cell"], "year": ["2008"], "comment": ["tpc.107.056309"]}, {"surname": ["Gaut", "Morton", "McCaig", "Clegg"], "given-names": ["BS", "BR", "BC", "MT"], "article-title": ["Substitution rate comparisons between grasses and palms: Synonymous rate differences at the nuclear gene Adh parallel rate differences at the plastid gene rbcL"], "source": ["Proceedings of the National Academy of Sciences"], "year": ["1996"], "volume": ["93"], "issue": ["19"], "fpage": ["10274"], "lpage": ["10279"], "pub-id": ["10.1073/pnas.93.19.10274"]}, {"surname": ["Paterson", "Bowers", "Chapman"], "given-names": ["AH", "JE", "BA"], "article-title": ["Ancient polyploidization predating divergence of the cereals, and its consequences for comparative genomics"], "source": ["Proceedings of the National Academy of Sciences"], "year": ["2004"], "volume": ["101"], "issue": ["26"], "fpage": ["9903"], "lpage": ["9908"], "pub-id": ["10.1073/pnas.0307901101"]}, {"surname": ["Huang", "Sirikhachornkit", "Su", "Faris", "Gill", "Haselkorn", "Gornicki"], "given-names": ["S", "A", "X", "J", "B", "R", "P"], "article-title": ["Genes encoding plastid acetyl-CoA carboxylase and 3-phosphoglycerate kinase of the Triticum/Aegilops complex and the evolutionary history of polyploid wheat"], "source": ["Proceedings of the National Academy of Sciences"], "year": ["2002"], "volume": ["99"], "issue": ["12"], "fpage": ["8133"], "lpage": ["8138"], "pub-id": ["10.1073/pnas.072223799"]}, {"surname": ["Tocchini-Valentini", "Fruscoloni", "Tocchini-Valentini"], "given-names": ["GD", "P", "GP"], "article-title": ["From the Cover: Structure, function, and evolution of the tRNA endonucleases of Archaea: An example of subfunctionalization"], "source": ["Proceedings of the National Academy of Sciences"], "year": ["2005"], "volume": ["102"], "issue": ["25"], "fpage": ["8933"], "lpage": ["8938"], "pub-id": ["10.1073/pnas.0502350102"]}, {"surname": ["Shiu", "Bleecker"], "given-names": ["S-H", "AB"], "article-title": ["Receptor-like kinases from Arabidopsis form a monophyletic gene family related to animal receptor kinases"], "source": ["Proceedings of the National Academy of Sciences"], "year": ["2001"], "volume": ["98"], "issue": ["19"], "fpage": ["10763"], "lpage": ["10768"], "pub-id": ["10.1073/pnas.181141598"]}, {"surname": ["Farrar", "Asp", "Lubberstedt", "Xu", "Thomas", "Christiansen", "Humphreys", "Donnison"], "given-names": ["K", "T", "T", "ML", "AM", "C", "MO", "IS"], "article-title": ["Construction of two Lolium perenne BAC libraries and identification of BACs containing candidate genes for disease resistance and forage quality"], "source": ["Molecular Breeding"], "year": ["2007"], "volume": ["19"], "issue": ["1"], "fpage": ["15"], "lpage": ["23"], "pub-id": ["10.1007/s11032-006-9036-z"]}, {"surname": ["Nei", "Kumar"], "given-names": ["M", "S"], "source": ["Molecular Evolution and Phylogenetics"], "year": ["2000"], "publisher-name": ["Oxford University Press"]}, {"surname": ["Yang"], "given-names": ["Z"], "article-title": ["PAML: a program package for phylogenetic analysis by maximum likelihood"], "source": ["Computational Applied Bioscience"], "year": ["1997"], "volume": ["13"], "issue": ["5"], "fpage": ["555"], "lpage": ["556"]}, {"surname": ["Hall"], "given-names": ["TA"], "article-title": ["BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT"], "source": ["Nucleic Acids Symposium Series"], "year": ["1999"], "volume": ["41"], "fpage": ["95"], "lpage": ["98"]}, {"surname": ["Aikake"], "given-names": ["H"], "article-title": ["A new look at the statistical model identification"], "source": ["IEEE Transactions on Automatic Control"], "year": ["1974"], "volume": ["19"], "issue": ["3"], "fpage": ["716"], "lpage": ["723"], "pub-id": ["10.1109/TAC.1974.1100705"]}]
BMC Evol Biol. 2008 Sep 5; 8:245