# The Plant Short-Chain Dehydrogenase (SDR) superfamily:genome-wide inventory and diversification patterns

Hanane Moummou, Yvonne Kallberg, Libert Brice Tonfack, Bengt Persson, Benoît van Der Rest

<figure><img src="image_1.png" /></figure>

To cite this version:

Hanane Moummou, Yvonne Kallberg, Libert Brice Tonfack, Bengt Persson, Benoît van Der Rest. The Plant Short-Chain Dehydrogenase (SDR) superfamily:genome-wide inventory and diversification patterns. BMC Plant Biology, 2012, vol. 12, pp.1-17. 10.1186/1471-2229-12-219. http://hal-00793421

HAL Id: hal-00793421

https://hal.science/hal-00793421

Submitted on 22 Feb 2013

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
**To cite this version:** Moummou, Hanane and Kallberg, Yvonne and Tonfack, Libert Brice and Persson, Bengt and Van der Rest, Benoît *The Plant Short-Chain Dehydrogenase (SDR) superfamily: genome-wide inventory and diversification patterns*. (2012) BMC Plant Biology, vol. 12 (n° 1). pp. 219. ISSN 1471-2229
# The Plant Short-Chain Dehydrogenase (SDR) superfamily: genome-wide inventory and diversification patterns

doi:10.1186/1471-2229-12-219

Hanane Moummou (flowernuina@yahoo.fr)

Yvonne Kallberg (yvonne.kallberg@ki.se)

Libert Brice Tonfack (libricetonfack@yahoo.fr)

Bengt Persson (bpn@ifm.liu.se)

Benoît van der Rest (benoit.van-der-rest@ensat.fr)
# The Plant Short-Chain Dehydrogenase (SDR) superfamily: genome-wide inventory and diversification patterns

Hanane Moummou<fnref n="1" />,<fnref n="2" />

Email: flowernuina@yahoo.fr

Yvonne Kallberg<fnref n="3" />

Email: yvonne.kallberg@ki.se

Libert Brice Tonfack<fnref n="1" />,<fnref n="4" />

Email: libricetonfack@yahoo.fr

Bengt Persson<fnref n="5" />,<fnref n="6" />

Email: bpn@ifm.liu.se

Benoît van der Rest<fnref n="1" />,<fnref n="7" />,*

Email: benoit.van-der-rest@ensat.fr

<footnote n="1">¹ Université de Toulouse, INPT-ENSAT, UMR990 Génomique et Biotechnologie des Fruits, Avenue de l'Agrobiopole, BP 32607, Castanet-Tolosan F-31326, France</footnote>

<footnote n="2">² Laboratory of Food Science, Faculty of Science Semlalia, University CADI AYYAD, Marrakech, Morocco</footnote>

<footnote n="3">³ Bioinformatics Infrastructure for Life Sciences, Science for Life Laboratory, Centre for Molecular Medicine, Karolinska Institutet, SE-171 77 Stockholm, Sweden</footnote>

<footnote n="4">⁴ Laboratory of Biotechnology and Environment, Unit of Plant Physiology and Improvement, Department of Plant Biology, Faculty of Science, University of Yaoundé 1, PO BOX 812, Yaoundé, Cameroon</footnote>

<footnote n="5">⁵ Science for Life Laboratory, Department of Cell and Molecular Biology (CMB), Karolinska Institutet, SE-17177 Stockholm, Sweden</footnote>

<footnote n="6">⁶ IFM Bioinformatics and Swedish e-Science Research Centre (SeRC), Linköping University, SE-58183 Linköping, Sweden</footnote>

<footnote n="7">⁷ INRA, UMR990 Génomique et Biotechnologie des Fruits, 24 Chemin de Borde Rouge, Castanet-Tolosan F-31326, France</footnote>

<footnote>* Corresponding author. INRA, UMR990 Génomique et Biotechnologie des Fruits, 24 Chemin de Borde Rouge, Castanet-Tolosan F-31326, France</footnote>
**Abstract**

**Background**

Short-chain dehydrogenases/reductases (SDRs) form one of the largest and oldest NAD(P)(H) dependent oxidoreductase families. Despite a conserved ‘Rossmann-fold’ structure, members of the SDR superfamily exhibit low sequence similarities, which constituted a bottleneck in terms of identification. Recent classification methods, relying on hidden-Markov models (HMMs), improved identification and enabled the construction of a nomenclature. However, functional annotations of plant SDRs remain scarce.

**Results**

Wide-scale analyses were performed on ten plant genomes. The combination of hidden Markov model (HMM) based analyses and similarity searches led to the construction of an exhaustive inventory of plant SDR. With 68 to 315 members found in each analysed genome, the inventory confirmed the over-representation of SDRs in plants compared to animals, fungi and prokaryotes. The plant SDRs were first classified into three major types — ‘classical’, ‘extended’ and ‘divergent’ — but a minority (10 % of the predicted SDRs) could not be classified into these general types (‘unknown’ or ‘atypical’ types). In a second step, we could categorize the vast majority of land plant SDRs into a set of 49 families. Out of these 49 families, 35 appeared early during evolution since they are commonly found through all the Green Lineage. Yet, some SDR families — tropinone reductase-like proteins (SDR65C), ‘ABA2-like’-NAD dehydrogenase (SDR110C), ‘salutaridine/menthone-reductase-like’ proteins (SDR114C), ‘dihydroflavonol 4-reductase’-like proteins (SDR108E) and ‘isoflavone-reductase-like’ (SDR460A) proteins — have undergone significant functional diversification within vascular plants since they diverged from Bryophytes. Interestingly, these diversified families are either involved in the secondary metabolism routes (terpenoids, alkaloids, phenolics) or participate in developmental processes (hormone biosynthesis or catabolism, flower development), in opposition to SDR families involved in primary metabolism which are poorly diversified.

**Conclusion**

The application of HMMs to plant genomes enabled us to identify 49 families that encompass all Angiosperms (‘higher plants’) SDRs, each family being sufficiently conserved to enable simpler analyses based only on overall sequence similarity. The multiplicity of SDRs in plant kingdom is mainly explained by the diversification of large families involved in different secondary metabolism pathways, suggesting that the chemical diversification that accompanied the emergence of vascular plants acted as a driving force for SDR evolution.

**Keywords**

Short-chain dehydrogenase/reductase (SDRs), SDR Nomenclature Initiative, Hidden Markov Model, Multigenic family, Plant
# Background

Short-chain dehydrogenases/reductases (SDRs) constitute one of the largest and oldest protein superfamilies known to date. This ancient family, found in all domains of life (Archea, Eukaryotes, Prokaryotes and viruses), is characterized by large sequence divergences but several common properties: (i) a conserved 3D structure consisting of ‘Rossmann-fold’ β-sheet with α-helices on both sides, (ii) an N-terminal dinucleotide cofactor binding motif, (iii) an active site with a catalytical residue motif YxxxK [1,2]. With the release of genome sequences of numerous living organisms, the availability of around 300 crystal structures and the identification of many enzymatic functions, much attention has been given to classify the members of the SDR superfamily. A first discrimination was established between five types of SDR: the ‘classical’ type, consisting of approximately 250 amino acids, the ‘extended’ type that has an additional 100-residue domain in the C-terminal region, the ‘intermediate’ type that displays a specific G/AxxGxxG/A cofactor binding motif, the ‘divergent’ type that comprises enoyl-reductases from plant and bacteria and harbours modifications both in the cofactor binding site and active site motifs and the ‘complex’ SDR which are usually part of large multi-domain enzymes, such as mammalian fatty acid synthases or bacterial polyketide synthases [2-4]. Moreover, the discovery of new oxidoreductase structures harbouring the SDR ‘Rossmann-fold’ motif revealed the existence of uncommon types, often referred to as ‘unknown’ or ‘atypical’ types. More recently, the diversity of SDRs, either their amino acid sequences or their functions, led to the development of a second classification effort: the ‘SDR Nomenclature Initiative’ that aims at being more informative regarding SDRs functions and at establishing a sustainable and expandable nomenclature system based on the use of a large set of hidden Markov models (HMM) [5]. Nowadays, 449 families have been listed in this nomenclature [6].

Although mentioned by several authors [2,4], the diversity of SDRs in plants has never been investigated thoroughly. The recent advances in sequencing techniques and the still-increasing speed of genome releases now facilitate an exhaustive review of complex multigenic families. In the case of SDRs, a second challenge for plant scientific community is to unravel the functions of these oxidoreductases. Indeed, in the TAIR10 annotation of *Arabidopsis thaliana* genome, a large majority of ‘classical’ SDRs (two thirds) are merely annotated as NAD(P)-binding Rossmann-fold superfamily protein oxidoreductase [7]. This lack of information prompted us to adopt an exhaustive approach on plant SDRs. In a previous paper, we reviewed the involvement of different SDRs in primary and secondary metabolism [8]. In the present paper, we combined the use of HMMs and phylogenetic analyses on a set of genomes representative of plant diversity, in order to conduct a global inventory of plant SDRs coherent with the current SDR classification and nomenclature. This inventory was integrated into a functional classification of plant SDRs. Since this genome-wide inventory confirmed the high diversity of plant SDRs, the distribution and evolution of the different SDR families was examined, notably to investigate the link between SDR diversification and the emergence of secondary metabolism in vascular plants.

# Methods

## Analysed genomes

Genome analyses were performed on ten distinct genomes comprising four Dicots, three Monocots, the Pteridophytae *Selagniella moellendorffii*, the moss *Physcomitrella patens* and
the Alga *Chlamydomonas reinhardtii*. Sequences and most annotations were downloaded from the Joint Genome Institute website. The predicted proteomes analysed [7,9-17] were deduced from the annotations given in Table 1.

<table><caption>**Table 1 Reference and size of the analyzed genomes**</caption><thead><tr><th>Species</th><th>Taxa</th><th>Annotation used</th><th>Number of loci</th><th>Reference</th></tr></thead><tbody><tr><td>Chlamydomonas reinhardtii</td><td>Chlorophyte</td><td>Chlre4.1_Augustus9</td><td>15935</td><td>[9]</td></tr><tr><td>Physcomitrella patens</td><td>Moss</td><td>proteins. Phypa1_1.FilteredModels</td><td>35938</td><td>[10]</td></tr><tr><td>Selaginella moellendorffii</td><td>Lycophyte</td><td>Selmo1_GeneModels_FilteredModels3</td><td>22285</td><td>[11]</td></tr><tr><td>Arabidopsis thaliana</td><td>Eudicot</td><td>TAIR9</td><td>27379</td><td>[7]</td></tr><tr><td>Populus trichocarpa</td><td>Eudicot</td><td>Populus.trichocarpa.v2.0</td><td>41377</td><td>[12]</td></tr><tr><td>Vitis vinifera</td><td>Eudicot</td><td>12X March 2010 release</td><td>26346</td><td>[13]</td></tr><tr><td>Glycine max</td><td>Eudicot</td><td>Glyma1_pacId</td><td>46367</td><td>[14]</td></tr><tr><td>Oryza sativa</td><td>Monocot</td><td>MSU Rice Genome Annotation (Osa1) Release 6.1</td><td>40577</td><td>[15]</td></tr><tr><td>Zea mays</td><td>Monocot</td><td>ZmB73_4a.53 working translations</td><td>102202</td><td>[16]</td></tr><tr><td>Sorghum bicolor</td><td>Monocot</td><td>Sorbi1_GeneModels_Sbi1_4_aa</td><td>34496</td><td>[17]</td></tr></tbody></table>

The number of loci corresponds to the protein coding genes predicted by the annotation

## HMM-based analyses of plant genomes

Genomic sets of predicted proteins were challenged with three Pfams HMMs [18]: PF00106, PF01370 and PF01073 using HMMER3. SDR Nomenclature Initiative HMMs were defined and updated as described previously [5]. The five SDR types ('classical', 'extended', 'intermediate', 'divergent', and 'complex') each has an HMM trained to identify sequences of respective type. The HMMs were created using HMMER3, with manually adjusted alignments of representative sequences as seed. Cutoffs are used to decide if a hit is significant or not: 'classical' — 138, 'extended' — 108, 'intermediate' — 162, 'divergent' — 160, and 'complex' — 140. In addition to the five types, an 'unknown' label is used for sequences with scores lower than these cutoffs but still high enough to safely predict the sequence as an SDR: 'classical' — 29, 'extended' — 75, and 'divergent' — 100. Scores below the cutoffs are considered not positive.

For the PLR/IFR family, an HMM was created and incorporated to the 'SDR Nomenclature Initiative' set of family HMMs. The procedure for training the HMM was the same as previously developed with iterative refinement of the model until no new members were found [5].

## Decision rules for SDR inventory

For each sequence recognized by an HMM (hit), a score was assigned. Yet, several sequences were only recognized by one or two HMMs (either the Pfam derived HMM or the SDR-type HMM) and sometimes with a very low score. Thus, we defined a series of rules schematized in a decision tree (Figure 1). Sequences that were recognized by SDR nomenclature initiative were directly considered as positive. Sequences identified with both remaining sets of HMMs were also considered as positive. For the remaining sequences recognized only by one HMM, we first looked at the existence of strong homology with positive hits identified in the previous steps in order to include putative 'truncated' proteins (see Inventory refinement, below). Alternatively, we checked individually the existence of structural data in the
scientific literature, which allowed either including some hits in our inventory or discarding certain families of enzymes, notably the medium-chain dehydrogenases that display distinct structural motifs. In absence of structural data, the sequences recognized by a single HMM could not be classified and were included in a list of ambiguous sequences (Additional file 1: Table S1) that contains oxido-reductases that still await structural data before confirming or infirming their belonging to the SDR superfamily.

<caption>**Figure 1 Decision rules used to make an inventory of plant SDRs using three sets of HMM.** All the HMM sets were run independently on the 10 predicted proteomes. The complete inventory and the ambiguous predictions are included as supplementary material (Additional file 2: Table S2 and Additional file 1: Table S1)</caption>

## Inventory refinement

For the gene loci that are associated with several gene models and therefore with different protein predictions, a sole amino acid sequence was selected according two criteria: (1) the maximum HMM score and (2) the maximum alignment score deduced from a BlastP performed on other plant genomes. When the HMM and BlastP analyses led to contradictory predictions, a single protein prediction was manually selected after aligning the different gene models with its closest homologues. To include in the SDR classification the truncated proteins that failed to be recognized by the HMMs, a BlastP sequence search was performed on each genome using as query sequences the complete list of SDRs recognized in the first round of HMM searches. All sequences that displayed a segment of 60 amino acids with more than 50 % identity were classified in the same type or family as its closest homologue.

## Distance matrices and phylogenetic analyses

Phylogenetic analyses and distance matrices were built using the Mega5 package [19]. Full length amino acid sequences were aligned using the ClustalW algorithm. Distance matrices evaluating the percentage of sequence identity were calculated on the basis of p-distance with the pairwise deletion option. Unless stated differently, phylogenetic trees were built using the Neighbor-Joining method. The percentage of replicate trees in which the associated taxa clustered together was calculated in the bootstrap test (500 replicates). Trees were drawn to scale, with branch lengths in the same units as those of the evolutionary distances used to infer the phylogenetic tree. The evolutionary distances were computed using the Poisson correction method and were expressed as the units of the number of amino acid substitutions per site.

## Statistical analyses

Principal Component Analysis was performed on the distribution matrix given in Figure 2 using the R 2.14.1 software [20] with in-house developed scripts (Elie Maza, personal communication). The robustness of the conclusions was checked by carrying the same analysis after removal of the individual exhibiting extreme values (SDR108E).

<caption>**Figure 2 Distribution of the SDR families in the analyzed plant genomes represented as a heat map.** The heat map was built on a distribution matrix deduced from the inventory classification shown in Additional file 2: Table S2. The blue to red color gradient reflects the number of SDR listed in each family in the different genomes; the absence of family is</caption>
indicated with a white square. The names of the families were deduced from the 'SDR Nomenclature Initiative' HMMs or by a representative gene accession for orphan families not recognized by a specific HMM

# Results and discussion

## HMM-driven inventory of plant SDR

Initial HMM analyses were performed on ten complete genomes: 4 Eudicots (Arabidopsis thaliana, Populus trichocarpa, Vitis vinifera, Glycine max), 3 Monocots (Zea mays, Oryza sativa and Sorghum bicolor), the lycophyte Selaginella moellendorffii, the moss Physcomitrella patens and the unicellular green algae Chlamydomonas reinhardtii (Table 1). The predicted 'proteomes' deduced from the genome annotations were searched against three distinct sets of HMMs: the Pfam HMMs considered to encompass most SDR (PF00106, PF01370, PF01073), HMMs developed in the framework of the SDR nomenclature initiative [5,21] and a set of HMMs developed to predict the type ('classical', 'extended', 'intermediary', 'divergent' and 'complex') of SDR (see Methods).

This first analysis led to an exhaustive inventory of plant SDRs presented in supplemental data (Additional file 2: Table S2 and Additional file 1: Table S1). This inventory was divided into a main list (Additional file 2: Table S2), where the HMM scores or the high similarity with known SDRs were sufficient to establish a good prediction, and a complementary list of ambiguous SDR predictions (Additional file 1: Table S1), containing proteins with low HMM scores and absence of structural data (see decision tree in Figure 1 and Methods). Despite its very low HMM scores, we included in the main list a large family that comprises pinoresinol reductase (PLR), isoflavone reductase (IFR), vestitone reductase, phenylcoumaran benzylic ether reductase and eugenol synthase. Indeed, the structures of several members of this family were resolved by crystallography and the data revealed the presence of a SDR-typical Rossmann-fold [22-25]. Subsequently, an HMM was created and incorporated to the 'SDR Nomenclature Initiative' set of HMMs. The PLR/IFR family was named SDR460A, where the 'A' stands for 'atypical'.

## Distribution of plant SDRs

The number and the distribution of plant SDRs of different types are summarized in Table 2. As in other Eukaryotes, the major types consist of 'classical' and 'extended' SDRs. 'Divergent' SDRs are limited to one conserved family: an enoyl-ACP reductase (ENR) involved in lipid biosynthesis (AT2G05990.1 in *Arabidopsis*) [8,26]. While neither 'intermediate' nor 'complex' types are found in plants, we can notice a high number of 'unknown' types, meaning that the sequence patterns clearly differ from the other types. As previously noticed [2], the SDR family is highly represented in the plant kingdom: while 73 SDRs were numbered in the human genome [27] and 39 in the cyanobacteria *Synechocystis* sp. PCC 6803 [28], the number of SDRs in land plants vary from 126 in the moss *P. patens* to 315 in soybean (*G. max*). Even if we consider the variations due to the genome sizes (Table 1), SDRs are more represented in the Angiosperms than in the algae *C. reinhardtii* or in *P. patens*, suggesting a relationship between the emergence of vascular plant and the apparent multiplicity of plant SDRs.
<table><caption>Table 2 Distribution of SDRs in different plants</caption><thead><tr><th rowspan="2"></th><th rowspan="2">Total SDR</th><th colspan="5">Types of SDR</th></tr><tr><th>Classical</th><th>Divergent</th><th>Extended</th><th>Atypical (PLR-IFR)</th><th>Unknown</th></tr></thead><tbody><tr><td>Arabidopsis</td><td>178</td><td>90</td><td>1</td><td>72</td><td>8</td><td>7</td></tr><tr><td>Poplar</td><td>268</td><td>122</td><td>4</td><td>106</td><td>16</td><td>20</td></tr><tr><td>Grapevine</td><td>205</td><td>95</td><td>2</td><td>88</td><td>14</td><td>6</td></tr><tr><td>Soybean</td><td>315</td><td>145</td><td>4</td><td>138</td><td>15</td><td>13</td></tr><tr><td>Rice</td><td>227</td><td>110</td><td>2</td><td>95</td><td>10</td><td>10</td></tr><tr><td>Maize</td><td>230</td><td>97</td><td>3</td><td>113</td><td>7</td><td>10</td></tr><tr><td>Sorghum</td><td>237</td><td>106</td><td>2*</td><td>114</td><td>7</td><td>8</td></tr><tr><td>Selaginella</td><td>142</td><td>64</td><td>1</td><td>55</td><td>5</td><td>17</td></tr><tr><td>Physcomitrella</td><td>126</td><td>59</td><td>2</td><td>55</td><td>0</td><td>10</td></tr><tr><td>Chlamydomonas</td><td>68</td><td>41</td><td>1</td><td>21</td><td>1</td><td>4</td></tr></tbody></table>

Families with low scores and no structural data (listed in Table S2) were omitted. *The presence of divergent SDRs in *Sorghum bicolor* was deduced from the Sbicolor_79_peptide annotation

## Sub-classification of plant SDR

The HMMs developed by the SDR nomenclature initiative [5] aim at classifying the SDR superfamily into a large number of families, at a level where this classification would be informative regarding the functions of SDRs. In a first analysis, the HMMs defined in the frame of the SDR nomenclature initiative directly recognized 74 % of plant SDRs. After performing similarity searches and associating truncated proteins with its closest homologues (see Methods), the proportion of non-classified SDRs dropped markedly since 94.5 % of plant SDRs were categorized into 49 families. While the majority of these families are found in most Tracheophytes (Table 3), seven (SDR58C, 59C, 74C, 86C, 90C, 103U, 107C; Additional file 2: Table S2) are found only in *P. patens* or in *C. reinhardtii*. The occurrence of different family sequences in the analysed genomes was represented as a heat map (Figure 2).

<table><caption>Table 3 Classification of plant SDRs</caption><thead><tr><th>Representative gene</th><th>SDR nomenclature initiative</th><th>Known functions</th><th>Occurence</th><th>Average identity (%)</th></tr></thead><tbody><tr><td>AT4G23420</td><td>SDR7C</td><td><i>Pisum sativum</i> Tic32 (chloroplast protein import translocon)</td><td>ViridP</td><td>49,4</td></tr><tr><td>AT1G67730</td><td>SDR12C</td><td>β-ketoacyl reductase (fatty acids elongation)</td><td>LandP</td><td>48,4</td></tr><tr><td>AT3G12800</td><td>SDR17C</td><td>-</td><td>ViridP</td><td>64,1</td></tr><tr><td>AT4G05530</td><td>SDR25C</td><td>SDRA-IBR1 (indole-3-butyric acid response 1)</td><td>ViridP</td><td>67,7</td></tr><tr><td>AT3G03330</td><td>SDR34C</td><td>-</td><td>ViridP</td><td>56,1</td></tr><tr><td>AT3G06060</td><td>SDR35C</td><td>-</td><td>ViridP</td><td>47,9</td></tr><tr><td>AT4G09750</td><td>SDR40C</td><td>-</td><td>ViridP*</td><td>70,8</td></tr><tr><td>AT1G54870</td><td>SDR57C</td><td>-</td><td>ViridP</td><td>58,0</td></tr><tr><td>AT5G06060</td><td>SDR65C</td><td>Tropinone Reductase</td><td>ViridP</td><td>53,3</td></tr><tr><td>AT3G03980</td><td>SDR68C</td><td>-</td><td>TracheoP</td><td>57,0</td></tr></tbody></table>
<table>
  <tr>
    <td>AT5G54190</td>
    <td>SDR73C</td>
    <td>Protochlorophyllide Oxidoreductase</td>
    <td>ViridP</td>
    <td>74,5</td>
  </tr>
  <tr>
    <td>AT3G50560</td>
    <td>SDR84C</td>
    <td>-</td>
    <td>ViridP</td>
    <td>60,4</td>
  </tr>
  <tr>
    <td>AT1G52340</td>
    <td>SDR110C</td>
    <td>ABA2 (xanthoxin oxidase), Tasselseed2, Secoisolariciresinol dehydrogenase, Momilactone A synthase, Isopiperitenol dehydrogenase</td>
    <td>LandP</td>
    <td>47,1</td>
  </tr>
  <tr>
    <td>AT3G61220</td>
    <td>SDR114C</td>
    <td>Salutaridine reductase, Menthone reductase, Isopiperitenone reductase</td>
    <td>ViridP</td>
    <td>45,4</td>
  </tr>
  <tr>
    <td>AT5G50600</td>
    <td>SDR119C</td>
    <td>Hydroxysteroid Dehydrogenase</td>
    <td>LandP</td>
    <td>44,4</td>
  </tr>
  <tr>
    <td>AT3G55290</td>
    <td>SDR132C</td>
    <td>Solanum tuberosum TDF511</td>
    <td>ViridP</td>
    <td>62,4</td>
  </tr>
  <tr>
    <td>AT1G24360</td>
    <td>SDR152C</td>
    <td>FAS-II- β-ketoacyl reductase (FabG)</td>
    <td>ViridP</td>
    <td>68,3</td>
  </tr>
  <tr>
    <td>AT1G10310</td>
    <td>SDR357C</td>
    <td>Pterin aldehyde reductase (folate salvage)</td>
    <td>TracheoP</td>
    <td>70,0</td>
  </tr>
  <tr>
    <td>AT5G10050</td>
    <td>SDR368C</td>
    <td>-</td>
    <td>ViridP</td>
    <td>45,8</td>
  </tr>
  <tr>
    <td>AT4G27760</td>
    <td>SDR369C</td>
    <td>Arabidopsis thaliana Forever Young</td>
    <td>ViridP</td>
    <td>57,2</td>
  </tr>
  <tr>
    <td>AT2G05990</td>
    <td>SDR87D</td>
    <td>Enoyl-ACP reductase (ENR)</td>
    <td>ViridP</td>
    <td>75,0</td>
  </tr>
  <tr>
    <td>AT1G49670</td>
    <td>-</td>
    <td>-</td>
    <td>ViridP</td>
    <td>50,6</td>
  </tr>
  <tr>
    <td>AT3G01980</td>
    <td>-</td>
    <td>Cucumis melo ADH2</td>
    <td>LandP</td>
    <td>57,8</td>
  </tr>
  <tr>
    <td>AT4G13250</td>
    <td>-</td>
    <td>NYC1/NOL (chlorophyll b reductase)</td>
    <td>ViridP</td>
    <td>48,1</td>
  </tr>
  <tr>
    <td>AT4G20760</td>
    <td>-</td>
    <td>-</td>
    <td>ViridP</td>
    <td>61,8</td>
  </tr>
  <tr>
    <td>AT5G04070</td>
    <td>-</td>
    <td>-</td>
    <td>LandP</td>
    <td>52,7</td>
  </tr>
  <tr>
    <td>AT4G10960</td>
    <td>SDR1E</td>
    <td>UDP-D-glucose/UDP-D-galactose 4-epimerase, UDP-arabinose 4-epimerase</td>
    <td>Virid</td>
    <td>55,4</td>
  </tr>
  <tr>
    <td>AT1G78570</td>
    <td>SDR2E</td>
    <td>NDP-L-rhamnose synthase/epimerase</td>
    <td>ViridP</td>
    <td>74,7</td>
  </tr>
  <tr>
    <td>AT5G66280</td>
    <td>SDR3E</td>
    <td>GDP-mannose 4,6-dehydratase</td>
    <td>LandP</td>
    <td>72,3</td>
  </tr>
  <tr>
    <td>AT1G17890</td>
    <td>SDR4E</td>
    <td>GDP-4-keto-6-deoxymannose-3,5-epimerase-4-reductase</td>
    <td>LandP</td>
    <td>73,1</td>
  </tr>
  <tr>
    <td>AT2G28760</td>
    <td>SDR6E</td>
    <td>UDP-xylose synthase, UDP-glucuronic acid decarboxylase</td>
    <td>ViridP</td>
    <td>69,7</td>
  </tr>
  <tr>
    <td>AT2G20360</td>
    <td>SDR22E</td>
    <td>-</td>
    <td>ViridP</td>
    <td>60,1</td>
  </tr>
  <tr>
    <td>AT1G47290</td>
    <td>SDR31E</td>
    <td>3β-hydroxysteroid-dehydrogenase/decarboxylase</td>
    <td>ViridP</td>
    <td>48,2</td>
  </tr>
  <tr>
    <td>AT2G33630</td>
    <td>SDR42E</td>
    <td>-</td>
    <td>ViridP*</td>
    <td>66,2</td>
  </tr>
  <tr>
    <td>AT4G30440</td>
    <td>SDR50E</td>
    <td>UDP-D-glucuronate 4-epimerase</td>
    <td>ViridP</td>
    <td>61,3</td>
  </tr>
  <tr>
    <td>AT4G33030</td>
    <td>SDR52E</td>
    <td>UDP-sulfoquinovose synthase</td>
    <td>ViridP</td>
    <td>73,8</td>
  </tr>
  <tr>
    <td>AT1G08200</td>
    <td>SDR67E</td>
    <td>UDP-D-apiose/UDP-D-xylose synthase</td>
    <td>LandP</td>
    <td>81,9</td>
  </tr>
  <tr>
    <td>AT5G28840</td>
    <td>SDR93E</td>
    <td>GDP-D-mannose 3',5'-epimerase</td>
    <td>ViridP</td>
    <td>87,4</td>
  </tr>
  <tr>
    <td>AT5G42800</td>
    <td>SDR108E</td>
    <td>Dihydroflavonol 4-reductase, Anthocyanidin reductase, Cinnamoyl-CoA reductase, Phenylacetaldehyde reductase, Eutypine reductase</td>
    <td>ViridP</td>
    <td>36,6</td>
  </tr>
  <tr>
    <td>GRMZM2G086773</td>
    <td>SDR115E</td>
    <td>HC-toxin reductase</td>
    <td>FlowerP</td>
    <td>55,0</td>
  </tr>
  <tr>
    <td>AT5G22500</td>
    <td>SDR117E</td>
    <td>fatty-acyl-CoA reductase</td>
    <td>LandP</td>
    <td>46,8</td>
  </tr>
  <tr>
    <td>AT4G24220</td>
    <td>SDR75U</td>
    <td>VEIN PATTERNING 1 (VEP1), progesterone 5β-reductase</td>
    <td>LandP</td>
    <td>53,4</td>
  </tr>
</table>
<table><tr><td>AT4G35250</td><td>SDR81U</td><td>-</td><td>ViridP</td><td>76,4</td></tr><tr><td>AT1G09340</td><td>SDR83U</td><td>Chloroplast stem-loop binding protein</td><td>ViridP</td><td>50,7</td></tr><tr><td>AT5G18660</td><td>SDR98U</td><td>3,8-divinyl protochlorophyllide a 8-vinyl reductase</td><td>ViridP</td><td>62,5</td></tr><tr><td>AT5G02240</td><td>SDR358U</td><td>-</td><td>ViridP*</td><td>68,9</td></tr><tr><td>AT1G32100 (PLR-IFR)</td><td>SDR460A</td><td>Pinoresinol reductase, Isoflavone reductase, Vestitone reductase, Phenylcoumaran benzylic ether reductase, Eugenol synthase</td><td>TracheoP</td><td>45,3</td></tr><tr><td>AT4G33360</td><td>-</td><td>Farnesol NAD dehydrogenase</td><td>LandP</td><td>63,2</td></tr><tr><td>AT4G00560</td><td>-</td><td>-</td><td>ViridP</td><td>56,5</td></tr></table>

Each family was associated with a representative gene and, when possible, with a specific SDR nomenclature initiative HMM. Information on the occurrence of SDRs in different genomes are reported by the taxon name (ViridP: Viridiplantae; LandP: Embryophytae; TracheoP: Tracheophytae; FlowerP: Magnioliophyta). Average pairwise identities were calculated from the sequences of plant genomes. Ambiguously predicted SDRs and families absent in flowering plants were omitted. *: occurrence in Viriplantae was deduced from the presence of homologues in other Green Algae genomes

On the opposite, 5.5 % of plant SDRs (from 4 % in Angiosperms to 29 % in *C. reinhardtii*) remained unclassified. The existence of these orphan SDRs lays in the conception of the 'SDR nomenclature initiative' HMMs. In order to achieve robust HMMs, the authors considered only families with sufficient number of representative and non-redundant sequences [5], thus excluding SDR families with too few members. To circumvent this difficulty, we examined the possibility to define new families on the sole basis of amino-acid sequence conservation. Therefore, all the unclassified SDRs from the main inventory (Additional file 2: Table S2) were associated to its closest homologues using BlastP searches and sequence alignments. Interestingly, all the unclassified sequences from Angiosperms clearly matched with at least one *Arabidopsis* SDR, the e value obtained from a BlastP against *Arabidopsis* predicted proteome never exceeding 1 e-40. Thus, seven new clusters were defined on the basis of sequence conservation, four being common to all the Viridiplantae genomes while three were found only in land plants (Figure 2 and Table 3). Within these clusters, the average pairwise sequence identities ranged from 48 % to 62 %. These conservation rates are consistent with the average pairwise identities observed for the families defined by a 'SDR nomenclature initiative'-HMM, that ranges from 37 % to 82 % identity (Table 3). All these clusters were represented by a limited number of sequences in each genome, supporting the explanation that the lack of 'SDR-nomenclature-initiative'-HMMs is simply the consequence of an insufficient set of sequences and that these families might be defined in the future, with the release of new sequences in the UNIPROT database. To complete the plant SDR classification, each new cluster was assigned a representative gene, based on an *Arabidopsis thaliana* identifier. While all angiosperms SDRs could be categorized in a family, defined either by a specific HMM or by primary structure conservation, 15 sequences from *C. reinhardtii*, 4 sequences from *P. patens* and one sequence from *S. moellendorffii* were too distant to other SDR sequences and remained unclassified.

By extension, the ambiguous SDR sequences were also clustered on the basis of sequence homologies, allowing the definition of nine potential families (Additional file 1: Table S1). Yet, in absence of structural data confirming the existence of typical SDR structures, these sequences were not analysed further.
In a last step, plant SDR classification was combined with functional information. Taking
advantage of our previous bibliographic research [8] and of the annotations found for
*Arabidopsis* (TAIR10), we completed the classification by mentioning all the known
functions described in the scientific literature in Table 3. Also, to each family, a
representative gene was chosen according to three criteria: (1) favour *Arabidopsis* accessions
with respect to the quality of TAIR annotations and its pertinence as a model plant; (2) when
possible, opt for genes that have been functionally characterised; otherwise (3), priority was
given to the accession that displayed the lowest average distance with other members of its
family.

Evolution and diversification of plant SDR as a potential trait of land plant emergence?

The distribution of the different families in the different taxa was further examined to
understand the evolution of the plant SDR superfamily. We first addressed the question of
potential origins of the different SDR families. Out of the 49 families listed in Table 3, 32
were found both in the algae *C. reinhardtii* and in the majority of land plants, suggesting that
most plant SDRs families emerged prior to land plant radiation that started -460 Myear ago,
in the Ordovician period [29]. For three additional families (SDR40C, SDR42E and
SDR358U), the absence of a member in *C. reinhardtii* or even in *P. patens* predicted
proteomes masked the occurrence of these families in other genus of green algae (*Volvox*,
*Micromonas*, *Chlorella* and *Ostreococcus*), suggesting that the families were ancestral, but
that the genes might have been lost in some taxa. In addition, 10 families absent in green
algae are common to all land plants (Figure 2 and Table 3), indicating that 45 families are
shared among land plants (embryophytes). 48 families were common to vascular plants as 3
additional families were specific to *S. moellendorffii* and Angiosperms. At last, a sole family,
SDR115E, was found only in Angiosperms. The origins of some families may be very
ancient: SDR1E, 2E, 6E and 7C families are found in all domains of life (Archea, Eukaryote,
Prokaryote) while the SDR12C, 17C, 25C, 34C, 35C, 22E and 31E families are common to
the majority of Eukaryotes [5]. Besides, several ancestral SDR families are close to
Prokaryotic ‘homologues’. For example, the origin of the plastids is illustrated by the
presence of chloroplastic SDRs similar to its cyanobacterial homologue. In a recent paper,
Kramm et al. [28] listed 39 SDRs in the genome of the cyanobacteria *Synechocystis* sp. PCC
6803. 20 of these SDRs show clear homologies (>35 % identity) with plant SDRs (data not
shown). The SDRs clusters present both in cyanobacteria and plant genomes include the very
ancient families (SDR1E, 2E, 3E, 6E) and several plastidial proteins involved in primary
metabolism, such as sulfolipid biosynthesis protein (SDR52E), protochlorophyllide
oxidoreductase (SDR73C), 3,8-divinyl protochlorophyllide a 8-vinyl reductase (SDR98U) or
the members of the fatty acid synthase (FasII) complex (SDR152C and SDR87D).

The origin of these taxon-specificities probably results from three evolutionary mechanisms: horizontal gene transfers, differentiation of a novel family from a pre-existing SDR family and loss of genes. Indeed, Tarrio *et al.* [30] established that the Vein Patterning 1 (SDR75U) gene family had undergone five lateral gene transfer events, one occurring from bacteria to an ancestor of land plants. Conversely, extensive search of SDR homologues in the Genbank database revealed clear homologies between independent taxa, such as the similarities between the Tracheophyte SDR68C members and its Proteobacteria homologues or the close relationship between plant PLR-IFR family and Bacteria or Ascomycete isoflavone reductase-like proteins (data not shown), thus illustrating the possible importance of horizontal gene transfers. An original example of SDR differentiation is illustrated by the
emergence of the Angiosperm-specific HC-toxin reductase (SDR115E) family, involved in
the pathogen Helminthosporium carbonum (HC) toxin reduction [31]. Since previous
phylogenetic analyses [32] showed the existence of significant homologies between HC-toxin
reductase (SDR115E) and the large dihydroflavonol 4-reductase (4-DFR, SDR108E) family,
we integrated SDR108E and SDR115E amino acids sequences in the same alignment and
phylogenetic analysis (Figure 3 and Additional file 3: Figure S1A). The topology of the
deduced tree (Figure 3 and Additional file 3: Figure S1A) suggests that the SDR115E branch
belongs to a larger clade that includes 4-DFR (AT5G42800.1 cluster, [33]), anthocyanidin
reductase (AT1G61720.1 cluster, [34,35]) and the brassinosteroid related 4-DFR-like protein
BEN1 (AT2G45400.1 cluster, [36]). The robustness of this topology was further checked
using different phylogeny algorithms (Neighbour-Joining and Maximum Likelihood) or
rooting the tree with external sequences from other SDR families (SDR1E, SDR6E,
SDR31E). All trees displayed similar topologies, SDR115E members always clustering with
4-DFR, anthocyanidin reductase and BEN1 (data not shown), thus supporting the view that
the HC-toxin reductase (SDR115E) branch evolved from an ancestor belonging to the
SDR108E family. The divergences of sequences within the SDR108E-115E ‘clade’ were
sufficient to establish two distinct HMM profiles. At last, two distinct features may illustrate
the role of loss of genes in SDR evolution: (i) although found in Monocots, grapevine, poplar
and soybean genomes, the SDR115E family is absent in Arabidopsis genome or ESTs
database; (ii) some families found in *P. patens* or in *S. moellendorffii* genomes (SDR74C,
86C, 103U, Additional file 2: Table S2) are absent in all the Angiosperms genomes,
suggesting that genes might have been lost during before flowering plants radiation.

**Figure 3 Phylogenetic tree of the SDR108E and SDR115E families.** The blue arrow indicates the node at the origin of the ‘AnR, 4-DFR and SDR115E’ branch. Amino acid sequences recognized by the SDR108E and SDR115E HMMs were aligned with ClustalW algorithm. The evolutionary history was inferred using the Neighbor-Joining method. The percentages of replicate trees in which the associated taxa clustered together in the bootstrap test (500 replicates) are shown next to the branches. Full references of sequences compressed in different clusters are provided as supplemental data (Additional file 3: figure 1A). Consistent trees were obtained using the Maximum Likelihood method or rooting the tree with other SDR families (SDR1E, SDR6E, SDR31E) as outgroups

The second obvious feature, when observing the distribution of SDR families (Figure 2), is
the expansion pattern of the different families. A Principal Component Analysis (PCA) was
performed on the distribution matrix used to build the heat map presented in Figure 2. It
allowed the individualization of ten families displaying high values on the first axis (Figure
4A). All these families are characterized by a large number of members in contrast to the
majority of SDR families represented in plant genomes with a limited set of sequences.
Interestingly, the second axis is mainly driven by the vectors formed by *P. patens* and *C.
reinhardtii* genomes (Figure 4B) and it discriminates two patterns of diversification: families
expanded both in the moss *P. patens* and in vascular plants (SDR1E, SDR2E, SDR6E,
SDR7C, SDR50E) and families expanded in vascular and flowering plants (SDR65C,
SDR108E, SDR110C, SDR114C and SDR460A).

**Figure 4 Diversification patterns of plant SDR families deduced from Principal Component Analysis.** PCA was calculated on the distribution matrix shown in Figure 2. A) Scatter plot deduced from the two first components: the first and second axes respectively participate for 79 % and 9 % of the diversity. B) Contribution of different genomes (expressed as vectors) in the first and second axes values. Angiosperms genomes follow the
order (anticlockwise): *G. max*, *Z. mays*, *A. thaliana*, *S. bicolor*, *P. trichocarpa*, *O. sativa*, *V. vinifera*

Remarkably, all the five families expanded in vascular plants comprise enzymes involved in secondary metabolism (Table 3): tropinone reductases (SDR65C) are known for their involvement in alkaloids biosynthesis; SDR110C NAD-dehydrogenases oxidize various phenolic or terpenic compounds, including xanthoxin, a precursor of abscissic acid (ABA); SDR114C menthone and salutaridine reductase, are involved in monterpene and alkaloid metabolism respectively; the large SDR108E family members catalyze the reduction of several phenolic precursors (4-dihydroflavonol, anthocyanidin, cinnamoyl-CoA, phenylacetaldehyde or eutypine) and last, the atypical PLR/IFR family (SDR460A) is also involved in phenolic metabolism. On the opposite, several poorly diversified clusters (SDR52E, 73C, 152C, 87D, 357C) that contain highly conserved sequences participate in primary metabolism such as chlorophyll synthesis or degradation, lipid metabolism or vitamin synthesis.

## Identification of functional clusters within SDR families

For multigenic SDR families, the analyses can be conducted further with phylogenetic calculations. To illustrate the importance of this complementary approach, we focused on two large families involved in secondary metabolism: SDR110C (ABA2 xanthoxin dehydrogenase family) and SDR108E (4-DFR) family. For tropinone reductase (SDR65C) and menthone/salutaridine reductase (SDR114C) families, readers are referred respectively to Brock *et al.* [37] and Ziegler *et al.* [38] for complete phylogenetic analyses.

In our previous review [8], we listed six different functions described for SDR110C in the scientific literature: Arabidopsis ABA2 xanthoxin dehydrogenase (abscissic acid biosynthesis) [39,40], rice diterpenoid momilactone synthase A [41], mint (*Mintha sativa*) isopiperitenol dehydrogenase [42], *Forsythia intermedia* secoisolariciresinol dehydrogenase (lignan biosynthesis), the maize or rice feminization gene *TASSELSEED2* [43] and the Arabidopsis *AtATA1* gene, involved in pollen and anther tapetal cells development [44]. A phylogenetic tree was established from the analysed genome sequences and completed with mint isopiperitenol dehydrogenase and forsythia secoisolariciresinol dehydrogenase sequences (Figure 5 and Additional file 3: Figure S1B). Remarkably, the six functions described in the literature were distributed on five different clades, thus giving valuable hypotheses regarding the putative function of the orthologues or paralogues in other Angiosperm species. On the different clades, we also observe that highly homologous SDRs are often clustered in specific chromosomal regions, illustrating the importance of gene duplication events in the diversification process. At last, the different accessions of Selaginella are distributed in three distinct clades. However, the bootstrap values of the different nodes were too low to clearly establish a clear relationship between lycophytes and Angiosperms SDR110C sequences.

**Figure 5 Phylogenetic tree of the SDR110C family.** Amino acids sequences recognized by the SDR110C HMM were aligned with ClustalW algorithm. The evolutionary history and the bootstrap test (500 replicates) were computed as described for SDR108E (Figure 3). Full references of sequences compressed in different clusters are provided as supplemental data (Additional file 3: Figure S1B)
For the highly variable SDR108E family, we included the SDR115E family in the analysis as both family are closely related (see above). As reported for the SDR110C family, several branches can be associated with functions described in the literature [8]: 4-DFR [33], anthocyanidin reductase (AnR) [34,35], HC-toxin reductase [31], phenylacetaldehyde reductase [45], cinnamoyl-CoA reductase (CCR) [46] or eutypine reductase [47] (Figure 3). In contrast to SDR110C, the tree is also informative concerning the evolution of SDRs among land plant since distinct sequences from *S. moellendorffii* and *P. patens* are clearly associated with independent clades. These associations are of special interest for certain classes of enzymes such as the CCR catalysing the first irreversible oxidation step leading to monolignol synthesis. Indeed, several enzymes involved in the lignin biosynthesis pathway appeared early in land plant evolution and the moss are believed to accumulate uncondensed monolignols [48]. Thus, the association on the same branch of sequences from *P. patens* and *S. moellendorffii* with Angiosperms bona fide CCR suggests that the enzyme anciently acquired its specificity and diverged rapidly from other SDR108E members. Last, as observed for SDR110C, several highly similar genes are clustered in specific chromosomal regions. Hence, with numerous members and a low conservation rate of amino acid sequences, the SDR108E family and its daughter branch SDR115E constitute a good example of a gradual and fast evolution of a multigenic family. Since the majority of the described enzymes reduce phenolic compounds, we may hypothesize that the SDR108E evolution accompanied tightly the complexification of phenolic and phenylpropanoid metabolism during land plant radiation.

Although essential for the functional study of large SDR families, phylogenetic analyses may be very informative for smaller families as well. This is the case of the Non-Yellow-Coloring 1 (NYC1) chlorophyllase b family, where the phylogenetic analyses clearly divide the family in two distinct clades: NYC1 and NOL (Non-yellow-coloring-Like) that diverged from a common ancestor (Figure 6). It was suggested that during evolution, the divergence led to the emergence of a functional hetero-oligomer, since both genes are necessary for chlorophyll b degradation [49,50].

<caption>**Figure 6** Phylogenetic trees of the chlorophyllase b (NYC1/NOL) family. Amino acids sequences from the AT4G13250 family were aligned with ClustalW algorithm. The evolutionary history and the bootstrap test (500 replicates) were computed as described for SDR108E (Figure 3)</caption>

At last, we carried analyses on three SDR families involved in lipids primary metabolism: fatty-acid synthase (FAS-II)-β-ketoacyl synthase (SDR152C) [51], (FAS-II)-enoyl-ACP reductase (SDR87D) [52] and the UDP-sulfoquinovose synthase (SDR52E) involved in sulfolipids biosynthesis [53]. By contrast with the families involved in secondary metabolism discussed above, the average sequence identity is high (Table 3), ranging from 68 % to 75 %. When the N-terminal chloroplast peptide signals are removed from sequence alignments, the average identities reach the scores of 79 % (SDR152C) and 84 % (SDR87D and SDR 52E). Phylogenetic trees were deduced from sequence alignments (Figure 7). Despite the presence of some duplication events observed for SDR152C (Figure 7A) and SDR87D (Figure 7B), the tree topologies are in good agreement with plant taxonomy for all the three SDR families, thus suggesting that the primary structure has been conserved under a high pressure of selection.
**Figure 7 Phylogenetic trees of three families involved in lipid primary metabolism.** (A) SDR152C-FasII-β-keto-reductase (β-KR); (B): SDR87D-FasII-Enoyl-ACP-reductase (ENR); (C) SDR52E UDP-sulfoquinovose synthase (SQD1). Amino acids sequences recognized by the SDR152C, SDR87D and SDR52E HMMs were aligned with ClustalW algorithm. The evolutionary history and the bootstrap test (500 replicates) were computed as described for SDR108E (Figure 3)

Conclusions

Work presented in this paper aimed at providing a full picture of plant SDRs using the current classification, especially the recent SDR nomenclature initiative. The combination of HMM models and similarity searches enabled us to classify most of the plant SDRs into a core of 49 families. Of these 49 families, 42 could be associated to an HMM, while the other 7 families being only defined on the basis of amino acids sequence conservation. Remarkably, all predicted SDRs from Angiosperms or *S. Moellendorffii* (corresponding to the so-called ‘higher plants’) could be categorized within these families. As all families exhibit a high degree of primary structure conservation, the average amino acid identities ranging from 37 % to 87 % among plant genomes, all SDRs sequences from Angiosperms can be analysed easily on the sole basis of sequence alignment, using very classical software (Blast, Multialin, ClustalW). For moss *P. patens* and green algae *C. reinhardtii* sequences, the predictions are less accurate, 3 % and 20 % of predicted SDRs remain unclassified. This limitation probably results from the under-representation of bryophyte and chlorophyte sequences compared to Angiosperms. In addition, the development of genome sequencing on more distant taxa (for example charophytes, liverworts or hornworts) should increase the number of UNIPROT sequences with sufficient divergences, thus improving the quality of HMM and allowing, in a mid-term, the definition of HMMs for the orphan SDR families.

Strikingly, the number of families found in Angiosperms (49) does not differ much from the 47 SDR families listed in the human genome [5]. The large proportion of families (35 out of 49) found in all Viridiplantae, from Algae to Angiosperms, is consistent with the view that most SDR sub-branches diverged early during evolution [54]. Plants possess either SDRs common to all Eukaryotes or SDRs of bacterial origin, in particular SDRs deriving from the plastidial endosymbiosis. However, the major difference between plants and other eukaryotes, that explains the high number of SDRs in ‘higher plants’, lies in the existence of large multigenic families. These families expanded much later during evolution, as attested by their under-representation in moss and algae. Because of their involvement in secondary metabolism routes (including hormone biosynthesis), they can be considered as an adaptive character that emerged during land colonization and emergence of the vascular apparatus.

Abbreviations

AnR, Anthocyanidin reductase; A. thaliana, Arabidopsis thaliana; C. reinhardtii, Chlamydomonas reinhardtii; 4-DFR, Dihydroflavonol 4-reductase; G. max, Glycine max; HMM, Hidden Markov model; IFR, Isoflavone reductase; O. sativa, Oryza sativa; P. patens, Physcomitrella patens; P. trichocarpa, Populus trichocarpa; PLR, Pinoresinol reductase; S. bicolor, Sorghum bicolor; S. moellendorffii, Selaginella moellendorffii; SDR, Short-chain dehydrogenase/reductase; V. vinifera, Vitis vinifera; Z. mays, Zea mays.
**Competing interests**

The authors declare that they have no competing interests.

**Author' contributions**

HM participated in sequence alignments, phylogenetic analyses and collection of functional annotations. LB initiated the SDR inventory and participated in sequence alignments. BVDR conceived the study, participated in its design and coordination and drafted the manuscript. YK and BP led the HMM analyses and helped to draft the manuscript. All authors read and approved the final manuscript.

**Acknowledgements**

Hanane Moummou has received a doctoral fellowship from EGIDE within the frame of a French-Morocco Volubilis project (MA-06-155) and Libert Brice Tonfack from "Service de Coopération et d'Action Culturelle" of the French Embassy (Cameroon). The authors are grateful to Jean-Claude Pech (INP-ENSAT), Mohamed Benichou (CADI AYYAD-Marrakech University) and Emmanuel Youmbi (University of Yaoundé) for the coordination of the exchange programs and their remarks during the manuscript preparation and acknowledge Elie Maza (INP-ENSAT) and Christine Rousseau (INP-ENSAT) for their help in statistical and bioinformatic analyses. Bioinformatic analyses benefited from the Bioinfo-GenoToul facilities and those of Bioinformatics Infrastructure for Life Sciences (BILS).

**References**

1. Filling C, Berndt KD, Benach J, Knapp S, Prozorovski T, Nordling E, Ladenstein R, Jörnvall H, Oppermann U: *Critical residues for structure and catalysis in short-chain dehydrogenases/reductases*. *J Biol Chem* 2002, **277**:25677-25684.

2. Kavanagh KL, Jörnvall H, Persson B, Oppermann U: *Medium- and short-chain dehydrogenase/reductase gene and protein families: the SDR superfamily: functional and structural diversity within a family of metabolic and regulatory enzymes*. *Cell Mol Life Sci* 2008, **65**:3895-3906.

3. Kallberg Y, Oppermann U, Jörnvall H, Persson B: *Short-chain dehydrogenase/reductase (SDR) relationships: a large family with eight clusters common to human, animal, and plant genomes*. *Protein Sci* 2002, **11**:636-641.

4. Kallberg Y, Oppermann U, Jörnvall H, Persson B: *Short-chain dehydrogenases/reductases (SDRs)*. *Eur J Biochem* 2002, **269**:4409-4417.

5. Kallberg Y, Oppermann U, Persson B: *Classification of the short-chain dehydrogenase/reductase superfamily using hidden Markov models*. *FEBS J* 2010, **277**:2375-2386.

6. Persson B, Kallberg Y, Oppermann U: SDR (short-chain dehydrogenases/reductases). [http://www.sdr-enzymes.org](http://www.sdr-enzymes.org).
7. Swarbreck D, Wilks C, Lamesch P, Berardini TZ, Garcia-Hernandez M, Foerster H, Li D, Meyer T, Muller R, Ploetz L, Radenbaugh A, Singh S, Swing V, Tissier C, Zhang P, Huala E: The Arabidopsis Information Resource (TAIR): gene structure and function annotation. *Nucleic Acids Res* 2008, **36**:D1009–D1014.

8. Tonfack LB, Moummou H, Latché A, Youmbi E, Benichou M, Pech J-C, van der Rest B: The plant SDR superfamily: involvement in primary and secondary metabolism. *Curr Top Plant Biol*, **12**:41–53.

9. Merchant SS, Prochnik SE, Vallon O, Harris EH, Karpowicz SJ, Witman GB, Terry A, Salamov A, Fritz-Laylin LK, Maréchal-Drouard L, Marshall WF, Qu L-H, Nelson DR, Sanderfoot AA, Spalding MH, Kapitonov VV, Ren Q, Ferris P, Lindquist E, Shapiro H, Lucas SM, Grimwood J, Schmutz J, Cardol P, Cerutti H, Chanfreau G, Chen C-L, Cognat V, Croft MT, Dent R, Dutcher S, Fernández E, et al: The Chlamydomonas genome reveals the evolution of key animal and plant functions. *Science* 2007, **318**:245–250.

10. Rensing SA, Lang D, Zimmer AD, Terry A, Salamov A, Shapiro H, Nishiyama T, Perroud P-F, Lindquist EA, Kamisugi Y, Tanahashi T, Sakakibara K, Fujita T, Oishi K, Shin-I T, Kuroki Y, Toyoda A, Suzuki Y, Hashimoto S-I, Yamaguchi K, Sugano S, Kohara Y, Fujiyama A, Anterola A, Aoki S, Ashton N, Barbazuk WB, Barker E, Bennetzen JL, Blankenship R, Cho SH, Dutcher SK, et al: The Physcomitrella genome reveals evolutionary insights into the conquest of land by plants. *Science* 2008, **319**:64–69.

11. Banks JA, Nishiyama T, Hasebe M, Bowman JL, Gribskov M, dePamphilis C, Albert VA, Aono N, Aoyama T, Ambrose BA, Ashton NW, Axtell MJ, Barker E, Barker MS, Bennetzen JL, Bonawitz ND, Chapple C, Cheng C, Correa LGG, Dacre M, DeBarry J, Dreyer I, Elias M, Engstrom EM, Estelle M, Feng L, Finet C, Floyd SK, Frommer WB, Fujita T, Gramzow L, Gutensohn M, et al: The Selaginella genome identifies genetic changes associated with the evolution of vascular plants. *Science* 2011, **332**:960–963.

12. Tuskan GA, Difazio S, Jansson S, Bohlmann J, Grigoriev I, Hellsten U, Putnam N, Ralph S, Rombauts S, Salamov A, Schein J, Sterck L, Aerts A, Bhalerao RR, Bhalerao RP, Blaudez D, Boerjan W, Brun A, Brunner A, Busov V, Campbell M, Carlson J, Chalot M, Chapman J, Chen G-L, Cooper D, Coutinho PM, Couturier J, Covert S, Cronk Q, Cunningham R, Davis J, et al: The genome of black cottonwood, Populus trichocarpa (Torr. & Gray). *Science* 2006, **313**:1596–1604.

13. Jaillon O, Aury J-M, Noel B, Policriti A, Clepet C, Casagrande A, Choisne N, Aubourg S, Vitulo N, Jubin C, Vezzi A, Legeai F, Hugueney P, Dasilva C, Horner D, Mica E, Jublot D, Poulain J, Bruyère C, Billault A, Segurens B, Gouyvenoux M, Ugarte E, Cattonaro F, Anthouard V, Vico V, Del Fabbro C, Alaux M, Di Gaspero G, Dumas V, Felice N, Paillard S, et al: The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. *Nature* 2007, **449**:463–467.

14. Schmutz J, Cannon SB, Schlueter J, Ma J, Mitros T, Nelson W, Hyten DL, Song Q, Thelen JJ, Cheng J, Xu D, Hellsten U, May GD, Yu Y, Sakurai T, Umezawa T, Bhattacharyya MK, Sandhu D, Valliyodan B, Lindquist E, Peto M, Grant D, Shu S, Goodstein D, Barry K, Futrell-Griggs M, Abernathy B, Du J, Tian Z, Zhu L, Gill N, Joshi T, Libault M, Sethuraman A, Zhang X-C, Shinozaki K, Nguyen HT, Wing RA, Cregan P,
Specht J, Grimwood J, Rokhsar D, Stacey G, Shoemaker RC, Jackson SA: **Genome sequence of the palaeopolyploid soybean.** *Nature* 2010, **463**:178–183.

15. Ouyang S, Zhu W, Hamilton J, Lin H, Campbell M, Childs K, Thibaud-Nissen F, Malek RL, Lee Y, Zheng L, Orvis J, Haas B, Wortman J, Buell CR: **The TIGR Rice Genome Annotation Resource: improvements and new features.** *Nucleic Acids Res* 2007, **35**:D883–D887.

16. Schnable PS, Ware D, Fulton RS, Stein JC, Wei F, Pasternak S, Liang C, Zhang J, Fulton L, Graves TA, Minx P, Reily AD, Courtney L, Kruchowski SS, Tomlinson C, Strong C, Delehaunty K, Fronick C, Courtney B, Rock SM, Belter E, Du F, Kim K, Abbott RM, Cotton M, Levy A, Marchetto P, Ochoa K, Jackson SM, Gillam B, Chen W, Yan L, et al: **The B73 maize genome: complexity, diversity, and dynamics.** *Science* 2009, **326**:1112–1115.

17. Paterson AH, Bowers JE, Bruggmann R, Dubchak I, Grimwood J, Gundlach H, Haberer G, Hellsten U, Mitros T, Poliakov A, Schmutz J, Spannagl M, Tang H, Wang X, Wicker T, Bharti AK, Chapman J, Feltus FA, Gowik U, Grigoriev IV, Lyons E, Maher CA, Martis M, Narechania A, Otillar RP, Penning BW, Salamov AA, Wang Y, Zhang L, Carpita NC, Freeling M, Gingle AR, Hash CT, Keller B, Klein P, Kresovich S, McCann MC, Ming R, Peterson DG, Mehboob-ur-Rahman, Ware D, Westhoff P, Mayer KFX, Messing J, Rokhsar DS: **The Sorghum bicolor genome and the diversification of grasses.** *Nature* 2009, **457**:551–556.

18. Finn RD, Mistry J, Tate J, Coggill P, Heger A, Pollington JE, Gavin OL, Gunasekaran P, Ceric G, Forslund K, Holm L, Sonnhammer ELL, Eddy SR, Bateman A: **The Pfam protein families database.** *Nucleic Acids Res* 2010, **38**:D211–D222.

19. 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**:2731–2739.

20. R Core Team: *R: A language and environment for statistical computing*. [http://www.R-project.org/].

21. Persson B, Kallberg Y, Bray JE, Bruford E, Dellaporta SL, Favia AD, Duarte RG, Jörnvall H, Kavanagh KL, Kedishvili N, Kisiela M, Maser E, Mindnich R, Orchard S, Penning TM, Thornton JM, Adamski J, Oppermann U: **The SDR (short-chain dehydrogenase/reductase and related enzymes) nomenclature initiative.** *Chem Biol Interact* 2009, **178**:94–98.

22. Min T, Kasahara H, Bedgar DL, Youn B, Lawrence PK, Gang DR, Halls SC, Park H, Hilsenbeck JL, Davin LB, Lewis NG, Kang C: **Crystal structures of pinoresinol-lariciresinol and phenylcoumaran benzylic ether reductases and their relationship to isoflavone reductases.** *J Biol Chem* 2003, **278**:50714–50723.

23. Shao H, Dixon RA, Wang X: **Crystal structure of vestitone reductase from alfalfa (Medicago sativa L.).** *J Mol Biol* 2007, **369**:265–276.

24. Wang X, He X, Lin J, Shao H, Chang Z, Dixon RA: **Crystal structure of isoflavone reductase from alfalfa (Medicago sativa L.).** *J Mol Biol* 2006, **358**:1341–1352.
25. Louie GV, Baiga TJ, Bowman ME, Koeduka T, Taylor JH, Spassova SM, Pichersky E, Noel JP: **Structure and reaction mechanism of basil eugenol synthase.** *PLoS One* 2007, **2:e993**.

26. Rafferty JB, Simon JW, Baldock C, Artymiuk PJ, Baker PJ, Stuitje AR, Slabas AR, Rice DW: **Common themes in redox chemistry emerge from the X-ray structure of oilseed rape (Brassica napus) enoyl acyl carrier protein reductase.** *Structure* 1995, **3**:927–938.

27. Bray JE, Marsden BD, Oppermann U: **The human short-chain dehydrogenase/reductase (SDR) superfamily: a bioinformatics summary.** *Chem Biol Interact* 2009, **178**:99–109.

28. Kramm A, Kisiela M, Schulz R, Maser E: **Short-chain dehydrogenases/reductases in cyanobacteria.** *FEBS J* 2012, **279**:1030–1043.

29. Steemans P, Hérissé AL, Melvin J, Miller MA, Paris F, Verniers J, Wellman CH: **Origin and radiation of the earliest vascular land plants.** *Science* 2009, **324**:353.

30. Tarrío R, Ayala FJ, Rodríguez-Trelles F: **The Vein Patterning 1 (VEP1) gene family laterally spread through an ecological network.** *PLoS One* 2011, **6:e22279**.

31. Meeley RB, Johal GS, Briggs SP, Walton JD: **A biochemical phenotype for a disease resistance gene of maize.** *Plant Cell* 1992, **4**:71–77.

32. Sindhu A, Chintamanani S, Brandt AS, Zanis M, Scofield SR, Johal GS: **A guardian of grasses: specific origin and conservation of a unique disease-resistance gene in the grass lineage.** *Proc Natl Acad Sci U S A* 2008, **105**:1762–1767.

33. Shirley BW, Hanley S, Goodman HM: **Effects of ionizing radiation on a plant genome: analysis of two Arabidopsis transparent testa mutations.** *Plant Cell* 1992, **4**:333–347.

34. Devic M, Guilleminot J, Debeaujon I, Bechtold N, Bensaude E, Koornneef M, Pelletier G, Delseny M: **The BANYULS gene encodes a DFR-like protein and is a marker of early seed coat development.** *Plant J* 1999, **19**:387–398.

35. Xie D-Y, Sharma SB, Paiva NL, Ferreira D, Dixon RA: **Role of anthocyanidin reductase, encoded by BANYULS in plant flavonoid biosynthesis.** *Science* 2003, **299**:396–399.

36. Yuan T, Fujioka S, Takatsuto S, Matsumoto S, Gou X, He K, Russell SD, Li J: **BEN1, a gene encoding a dihydroflavonol 4-reductase (DFR)-like protein, regulates the levels of brassinosteroids in Arabidopsis thaliana.** *Plant J* 2007, **51**:220–233.

37. Brock A, Brandt W, Dräger B: **The functional divergence of short-chain dehydrogenases involved in tropinone reduction.** *Plant J* 2008, **54**:388–401.

38. Ziegler J, Facchini PJ, Geissler R, Schmidt J, Ammer C, Kramell R, Voigtländer S, Gesell A, Pienkny S, Brandt W: **Evolution of morphine biosynthesis in opium poppy.** *Phytochemistry* 2009, **70**:1696–1707.
39. González-Guzmán M, Apostolova N, Bellés JM, Barrero JM, Piqueras P, Ponce MR, Micol JL, Serrano R, Rodríguez PL: The short-chain alcohol dehydrogenase ABA2 catalyzes the conversion of xanthoxin to abscisic aldehyde. *Plant Cell* 2002, **14**:1833–1846.

40. Cheng W-H, Endo A, Zhou L, Penney J, Chen H-C, Arroyo A, Leon P, Nambara E, Asami T, Seo M, Koshiba T, Sheen J: A unique short-chain dehydrogenase/reductase in Arabidopsis glucose signaling and abscisic acid biosynthesis and functions. *Plant Cell* 2002, **14**:2723-2743.

41. Shimura K, Okada A, Okada K, Jikumaru Y, Ko K-W, Toyomasu T, Sassa T, Hasegawa M, Kodama O, Shibuya N, Koga J, Nojiri H, Yamane H: Identification of a biosynthetic gene cluster in rice for momilactones. *J Biol Chem* 2007, **282**:34013-34018.

42. Ringer KL, Davis EM, Croteau R: Monoterpene metabolism. *Cloning, expression, and characterization of (-)-isopiperitenol/(-)-carveol dehydrogenase of peppermint and spearmint. Plant Physiol* 2005, **137**:863-872.

43. DeLong A, Calderon-Urrea A, Dellaporta SL: Sex determination gene TASSELSEED2 of maize encodes a short-chain alcohol dehydrogenase required for stage-specific floral organ abortion. *Cell* 1993, **74**:757-768.

44. Lebel-Hardenack S, Ye D, Koutnikova H, Saedler H, Grant SR: Conserved expression of a TASSELSEED2 homolog in the tapetum of the dioecious Silene latifolia and Arabidopsis thaliana. *Plant J* 1997, **12**:515-526.

45. Tieman DM, Loucas HM, Kim JY, Clark DG, Klee HJ: Tomato phenylacetaldehyde reductases catalyze the last step in the synthesis of the aroma volatile 2-phenylethanol. *Phytochemistry* 2007, **68**:2660-2669.

46. Lacombe E, Hawkins S, Van Doorsselaere J, Piquemal J, Goffner D, Poeydomenge O, Boudet AM, Grima-Pettenati J: Cinnamoyl CoA reductase, the first committed enzyme of the lignin branch biosynthetic pathway: cloning, expression and phylogenetic relationships. *Plant J* 1997, **11**:429-441.

47. Guillén P, Guis M, Martínez-Reina G, Colrat S, Dalmayrac S, Deswarte C, Bouzayen M, Roustan JP, Fallot J, Pech JC, Latché A: A novel NADPH-dependent aldehyde reductase gene from Vigna radiata confers resistance to the grapevine fungal toxin eutypine. *Plant J* 1998, **16**:335-343.

48. Xu Z, Zhang D, Hu J, Zhou X, Ye X, Reichel KL, Stewart NR, Syrenne RD, Yang X, Gao P, Shi W, Doeppke C, Sykes RW, Burris JN, Bozell JJ, Cheng MZ-M, Hayes DG, Labbe N, Davis M, Stewart CN Jr, Yuan JS: Comparative genome analysis of lignin biosynthesis gene families across the plant kingdom. *BMC Bioinforma* 2009, **10**(Suppl 11):S3.

49. Kusaba M, Ito H, Morita R, Iida S, Sato Y, Fujimoto M, Kawasaki S, Tanaka R, Hirochika H, Nishimura M, Tanaka A: Rice NON-YELLOW COLORING1 is involved in light-harvesting complex II and grana degradation during leaf senescence. *Plant Cell* 2007, **19**:1362-1375.
50. Sato Y, Morita R, Katsuma S, Nishimura M, Tanaka A, Kusaba M: Two short-chain dehydrogenase/reductases, NON-YELLOW COLORING 1 and NYC1-LIKE, are required for chlorophyll b and light-harvesting complex II degradation during senescence in rice. *Plant J* 2009, **57**:120–131.

51. Slabas AR, Chase D, Nishida I, Murata N, Sidebottom C, Safford R, Sheldon PS, Kekwick RG, Hardie DG, Mackintosh RW: Molecular cloning of higher-plant 3-oxoacyl-(acyl carrier protein) reductase. Sequence identities with the nodG-gene product of the nitrogen-fixing soil bacterium Rhizobium meliloti. *Biochem J* 1992, **283**:321–326.

52. Kater MM, Koningstein GM, Nijkamp HJ, Stuitje AR: cDNA cloning and expression of Brassica napus enoyl-acyl carrier protein reductase in Escherichia coli. *Plant Mol Biol* 1991, **17**:895–909.

53. Sanda S, Leustek T, Theisen MJ, Garavito RM, Benning C: Recombinant Arabidopsis SQD1 converts udp-glucose and sulfite to the sulfolipid head group precursor UDP-sulfoquinovose in vitro. *J Biol Chem* 2001, **276**:3941–3946.

54. Jörnvall H, Hedlund J, Bergman T, Oppermann U, Persson B: Superfamilies SDR and MDR: from early ancestry to present forms. Emergence of three lines, a Zn-metalloenzyme, and distinct variabilities. *Biochem Biophys Res Commun* 2010, **396**:125–130.

## Additional files

### Additional_file_1 as XLS

Additional file 1: Table S1 List of ambiguous predictions. Proteins were only recognized with one set of HMMs with a low score.

### Additional_file_2 as XLS

Additional file 2: Table S2 Exhaustive inventory of SDR.

### Additional_file_3 as PPT

Additional file 3: Figure S1 Full phylogenetic trees of SDR108E (A) and SDR110C (B) families. The evolutionary history was inferred using the Neighbor-Joining method. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (500 replicates) are shown next to the branches.
<figure><img src="image_3.png" /><figcaption>Figure 1</figcaption></figure>
<table><thead><tr><th></th><th>C. reinhardtii</th><th>P. patens</th><th>S. moellendorfii</th><th>O. sativa</th><th>Z. mays</th><th>S. bicolor</th><th>V. vinifera</th><th>A. thaliana</th><th>P. trichocarpa</th><th>G. max</th></tr></thead><tbody><tr><td>SDR7C</td><td>3</td><td>9</td><td>12</td><td>11</td><td>9</td><td>8</td><td>7</td><td>8</td><td>11</td><td>23</td></tr><tr><td>SDR12C</td><td>1</td><td>1</td><td>2</td><td>10</td><td>5</td><td>9</td><td>4</td><td>2</td><td>5</td><td>5</td></tr><tr><td>SDR17C</td><td>1</td><td>2</td><td>2</td><td>1</td><td>2</td><td>3</td><td>2</td><td>2</td><td>2</td><td>2</td></tr><tr><td>SDR25C</td><td>1</td><td>1</td><td>1</td><td>1</td><td>2</td><td>2</td><td>2</td><td>1</td><td>2</td><td>3</td></tr><tr><td>SDR34C</td><td>1</td><td>2</td><td>0</td><td>1</td><td>1</td><td>1</td><td>1</td><td>2</td><td>2</td><td>2</td></tr><tr><td>SDR35C</td><td>1</td><td>4</td><td>2</td><td>1</td><td>2</td><td>2</td><td>1</td><td>2</td><td>1</td><td>2</td></tr><tr><td>SDR40C</td><td>0</td><td>0</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td></tr><tr><td>SDR57C</td><td>3</td><td>3</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>2</td><td>4</td><td>7</td></tr><tr><td>SDR65C</td><td>2</td><td>3</td><td>2</td><td>18</td><td>7</td><td>13</td><td>21</td><td>16</td><td>14</td><td>19</td></tr><tr><td>SDR68C</td><td>0</td><td>0</td><td>1</td><td>6</td><td>5</td><td>7</td><td>2</td><td>4</td><td>5</td><td>1</td></tr><tr><td>SDR73C</td><td>1</td><td>2</td><td>1</td><td>2</td><td>3</td><td>2</td><td>2</td><td>3</td><td>3</td><td>3</td></tr><tr><td>SDR84C</td><td>0</td><td>2</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>2</td><td>3</td></tr><tr><td>SDR110C</td><td>0</td><td>2</td><td>6</td><td>24</td><td>12</td><td>14</td><td>15</td><td>12</td><td>28</td><td>27</td></tr><tr><td>SDR114C</td><td>1</td><td>3</td><td>11</td><td>11</td><td>17</td><td>20</td><td>13</td><td>6</td><td>16</td><td>21</td></tr><tr><td>SDR119C</td><td>0</td><td>2</td><td>4</td><td>7</td><td>7</td><td>7</td><td>6</td><td>8</td><td>9</td><td>4</td></tr><tr><td>SDR132C</td><td>1</td><td>1</td><td>2</td><td>1</td><td>1</td><td>1</td><td>3</td><td>6</td><td>1</td><td>2</td></tr><tr><td>SDR152C</td><td>1</td><td>2</td><td>1</td><td>3</td><td>4</td><td>3</td><td>2</td><td>1</td><td>2</td><td>4</td></tr><tr><td>SDR357C</td><td>0</td><td>0</td><td>1</td><td>1</td><td>4</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td></tr><tr><td>SDR368C</td><td>1</td><td>2</td><td>2</td><td>2</td><td>3</td><td>2</td><td>3</td><td>2</td><td>4</td><td>4</td></tr><tr><td>SDR369C</td><td>1</td><td>2</td><td>1</td><td>1</td><td>2</td><td>1</td><td>1</td><td>3</td><td>2</td><td>2</td></tr><tr><td>AT1G49670.1</td><td>1</td><td>2</td><td>4</td><td>1</td><td>2</td><td>2</td><td>1</td><td>1</td><td>1</td><td>1</td></tr><tr><td>AT3G01980.1</td><td>0</td><td>1</td><td>1</td><td>1</td><td>2</td><td>1</td><td>1</td><td>1</td><td>2</td><td>1</td></tr><tr><td>AT4G13250.1</td><td>2</td><td>3</td><td>1</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>3</td></tr><tr><td>AT4G20760.1</td><td>1</td><td>2</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>2</td></tr><tr><td>AT5G04070.1</td><td>0</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>2</td><td>1</td><td>2</td></tr><tr><td>SDR87D</td><td>1</td><td>2</td><td>1</td><td>2</td><td>3</td><td>2</td><td>2</td><td>1</td><td>4</td><td>4</td></tr><tr><td>SDR1E</td><td>1</td><td>9</td><td>2</td><td>7</td><td>14</td><td>7</td><td>5</td><td>9</td><td>7</td><td>12</td></tr><tr><td>SDR2E</td><td>2</td><td>6</td><td>3</td><td>3</td><td>4</td><td>4</td><td>1</td><td>4</td><td>5</td><td>11</td></tr><tr><td>SDR3E</td><td>0</td><td>2</td><td>1</td><td>2</td><td>2</td><td>4</td><td>1</td><td>2</td><td>2</td><td>2</td></tr><tr><td>SDR4E</td><td>0</td><td>1</td><td>1</td><td>3</td><td>1</td><td>2</td><td>1</td><td>2</td><td>2</td><td>2</td></tr><tr><td>SDR6E</td><td>1</td><td>5</td><td>3</td><td>6</td><td>13</td><td>6</td><td>6</td><td>6</td><td>7</td><td>14</td></tr><tr><td>SDR22E</td><td>1</td><td>2</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>2</td><td>1</td></tr><tr><td>SDR31E</td><td>2</td><td>2</td><td>1</td><td>3</td><td>4</td><td>11</td><td>2</td><td>3</td><td>2</td><td>9</td></tr><tr><td>SDR42E</td><td>1</td><td>1</td><td>1</td><td>1</td><td>2</td><td>2</td><td>1</td><td>1</td><td>2</td><td>2</td></tr><tr><td>SDR50E</td><td>1</td><td>9</td><td>9</td><td>5</td><td>12</td><td>5</td><td>2</td><td>6</td><td>9</td><td>11</td></tr><tr><td>SDR52E</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>2</td></tr><tr><td>SDR67E</td><td>0</td><td>1</td><td>2</td><td>1</td><td>2</td><td>1</td><td>6</td><td>2</td><td>2</td><td>6</td></tr><tr><td>SDR75U</td><td>0</td><td>1</td><td>8</td><td>3</td><td>3</td><td>3</td><td>1</td><td>1</td><td>12</td><td>3</td></tr><tr><td>SDR81U</td><td>1</td><td>2</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>2</td><td>2</td></tr><tr><td>SDR83U</td><td>2</td><td>2</td><td>4</td><td>2</td><td>3</td><td>2</td><td>2</td><td>2</td><td>4</td><td>4</td></tr><tr><td>SDR93E</td><td>1</td><td>3</td><td>9</td><td>2</td><td>2</td><td>2</td><td>2</td><td>1</td><td>2</td><td>5</td></tr><tr><td>SDR98U</td><td>1</td><td>1</td><td>1</td><td>2</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>2</td></tr><tr><td>SDR108E</td><td>8</td><td>7</td><td>20</td><td>44</td><td>44</td><td>49</td><td>48</td><td>24</td><td>49</td><td>45</td></tr><tr><td>SDR115E</td><td>0</td><td>0</td><td>0</td><td>7</td><td>4</td><td>4</td><td>4</td><td>0</td><td>2</td><td>3</td></tr><tr><td>SDR117E</td><td>0</td><td>2</td><td>2</td><td>7</td><td>5</td><td>12</td><td>5</td><td>8</td><td>9</td><td>10</td></tr><tr><td>SDR358U</td><td>0</td><td>3</td><td>2</td><td>2</td><td>2</td><td>1</td><td>1</td><td>2</td><td>1</td><td>2</td></tr><tr><td>SDR460A (PLR-IFR)</td><td>0</td><td>0</td><td>5</td><td>10</td><td>7</td><td>7</td><td>14</td><td>8</td><td>16</td><td>15</td></tr><tr><td>AT4G00560.1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>2</td></tr><tr><td>AT4G03360.1</td><td>0</td><td>1</td><td>1</td><td>1</td><td>1</td><td>2</td><td>1</td><td>1</td><td>2</td><td>1</td></tr></tbody></table>

<figure><img src="image_4.png" /></figure>
<figure><img src="image_5.png" /><figcaption>Figure 3</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>Figure 4</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>Figure 5</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>Figure 6.5</figcaption></figure>
<figure><img src="image_9.png" /><figcaption>(A) SDR152C
(FasII-β-KR)</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>(B) SDR87D
(FasII-ENR)</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>(C) SDR52E
(SQD1)</figcaption></figure>
Additional files provided with this submission:

Additional file 1: 6226069457362749_add1.xls, 27K

http://www.biomedcentral.com/imedia/1710856356853177/supp1.xls

Additional file 2: 6226069457362749_add2.xls, 478K

http://www.biomedcentral.com/imedia/1953275481853177/supp2.xls

Additional file 3: 6226069457362749_add3.ppt, 170K

http://www.biomedcentral.com/imedia/1683065970853177/supp3.ppt