- Research article
- Open Access
Genome-wide analysis of the rice and arabidopsis non-specific lipid transfer protein (nsLtp) gene families and identification of wheat nsLtp genes by EST data mining
BMC Genomics volume 9, Article number: 86 (2008)
Plant non-specific lipid transfer proteins (nsLTPs) are encoded by multigene families and possess physiological functions that remain unclear. Our objective was to characterize the complete nsLtp gene family in rice and arabidopsis and to perform wheat EST database mining for nsLtp gene discovery.
In this study, we carried out a genome-wide analysis of nsLtp gene families in Oryza sativa and Arabidopsis thaliana and identified 52 rice nsLtp genes and 49 arabidopsis nsLtp genes. Here we present a complete overview of the genes and deduced protein features. Tandem duplication repeats, which represent 26 out of the 52 rice nsLtp genes and 18 out of the 49 arabidopsis nsLtp genes identified, support the complexity of the nsLtp gene families in these species. Phylogenetic analysis revealed that rice and arabidopsis nsLTPs are clustered in nine different clades. In addition, we performed comparative analysis of rice nsLtp genes and wheat (Triticum aestivum) EST sequences indexed in the UniGene database. We identified 156 putative wheat nsLtp genes, among which 91 were found in the 'Chinese Spring' cultivar. The 122 wheat non-redundant nsLTPs were organized in eight types and 33 subfamilies. Based on the observation that seven of these clades were present in arabidopsis, rice and wheat, we conclude that the major functional diversification within the nsLTP family predated the monocot/dicot divergence. In contrast, there is no type VII nsLTPs in arabidopsis and type IX nsLTPs were only identified in arabidopsis. The reason for the larger number of nsLtp genes in wheat may simply be due to the hexaploid state of wheat but may also reflect extensive duplication of gene clusters as observed on rice chromosomes 11 and 12 and arabidopsis chromosome 5.
Our current study provides fundamental information on the organization of the rice, arabidopsis and wheat nsLtp gene families. The multiplicity of nsLTP types provide new insights on arabidopsis, rice and wheat nsLtp gene families and will strongly support further transcript profiling or functional analyses of nsLtp genes. Until such time as specific physiological functions are defined, it seems relevant to categorize plant nsLTPs on the basis of sequence similarity and/or phylogenetic clustering.
Plant non-specific lipid transfer proteins (nsLTPs) were first isolated from spinach leaves and named for their ability to mediate the in vitro transfer of phospholipids between membranes . NsLTPs are widely distributed in the plant kingdom and form multigenic families of related proteins. However, in vitro lipid transfer or binding has been demonstrated only for a limited number of proteins and most nsLTPs have been identified on the basis of sequence homology, sequences deduced from cDNA clones or genes. All known plant nsLTPs are synthesized as precursors with a N-terminal signal peptide. Plant nsLTPs are small (usually 6.5 to 10.5 kDa) and basic (isoelectric point (pI) ranging usually from 8.5 to 12) proteins characterized by an eight cysteine motif (8 CM) backbone as follows: C-Xn-C-Xn-CC-Xn-CXC-Xn-C-Xn-C . The cysteine residues are engaged in four disulfide bonds that stabilize a hydrophobic cavity, which allows the binding of different lipids and hydrophobic compounds in vitro . Based on their molecular masses, plant nsLTPs were first separated into two types: type I (9 kDa) and type II (7 kDa) that are distinct both in terms of primary sequence identity (less than 30%) and lipid transfer efficiency . Although they have different cysteine pairing patterns, type I and type II nsLTPs constitute a structurally related family of proteins. Type I nsLTPs are characterized by a long tunnel-like cavity [4, 5] while a wheat type II nsLTP has two adjacent hydrophobic cavities . Several anther-specific proteins that display considerable homology with plant nsLTPs  have been proposed as a third type that differs from the two others by the number of amino acid residues interleaved in the 8 CM structure . To date, no structural data exists on the lipid transfer ability of type III nsLTPs.
Because they have been shown to transfer lipid molecules between membranes in vitro, plant nsLTPs were first suggested to be involved in membrane biogenesis . However, as they are synthesized with a N-terminal signal peptide , nsLTPs could not fulfill this function and were thought to be involved in secretion of extracellular lipophillic material, including cutin monomers . NsLTPs are possibly involved in a range of other biological processes, but their physiological function is not clearly understood. Like many other families of low molecular mass cysteine-rich proteins, nsLTPs display intrinsic antimicrobial properties and are thought to participate in plant defense mechanisms [11, 12]. This hypothetical function is also supported by the induction of the expression of many nsLtp genes in response to biotic infections or application of fungal elicitors [13–17] and by the enhanced tolerance to bacterial pathogens by overexpression of a barley nsLtp gene in transgenic arabidopsis . Due to their possible involvement in plant defense mechanisms, nsLTPs are recognized to be pathogenesis-related proteins and constitute the PR-14 family . Roles in plant defense signaling pathways have also been proposed since the disruption of the arabidopsis DIR1 gene, which encodes a nsLTP with an 8 CM distinct from those of types I, II or III, impairs the systemic acquired resistance signaling pathway . Similarly a wheat nsLTP competes with the fungal cryptogein for a same binding site in tobacco plasma membranes . A role in the mobilization of lipid reserves has also been suggested for germination-specific nsLTPs [22–24]. Finally, nsLTPs are thought to possess a function in male reproductive tissues . This role appears to be mainly related to type III nsLTPs whose genes display anther-specific expression , and to a few type I nsLtp genes including the rape E2 gene , the arabidopsis AtLtp12 gene (At3g51590)  and the rice t42 gene (Os01g12020)  that are also predominantly expressed at the early stage of anther development. It has been suggested that nsLTPs are involved in the deposition of material in the developing pollen wall ; however their precise function in pollen remains to be elucidated.
Plant nsLTPs are encoded by small multigene families but to date none has been extensively characterized. Six members have been identified in pepper , 11 in cotton , 14 in loblolly pine , 15 in arabidopsis , and 23 in wheat . The availability of the complete sequence of the arabidopsis , rice for both indica  and japonica subspecies , poplar  and grapevine  genomes has greatly enhanced our ability to characterize complex multigene families [38–40]. In polyploid genomes such as the allohexaploid wheat Triticum aestivum, the presence of multiple putative copies of each gene increases the complexity of the multigene families and the number of closely related sequences. With around 16,000 Mb , the genome of the hexaploid wheat is 128 times the size of the genome of the dicotyledonous model plant Arabidopsis thaliana and 38 times that of the monocotyledonous model plant Oryza sativa and has not been sequenced yet. Nevertheless, efforts made to generate wheat cDNA libraries [42–45] mean EST database mining can also be a successful strategy for the identification of multigene family members in complex genomes [46, 47]. In wheat, novel genes encoding polyphenol oxidases , storage proteins  and nsLTPs  were identified by EST database mining.
In the present study, we took advantage of the completion of the rice (japonica subspecies) and arabidopsis genome sequences to perform a genome-wide analysis of the nsLtp gene family in both species. In an effort to identify new members of the wheat nsLtp gene family, we searched the large public-domain collection of wheat ESTs for sequences displaying homologies with characterized rice nsLtp genes. In order to compare rice, arabidopsis and wheat nsLTP evolution, we performed phylogenetic analysis of the nsLTPs from these three plant species.
The Oryza sativa nsLtpgene family is composed of 52 members
Based on a conserved 8 CM, nsLTPs remain a structurally-related family of proteins. However, as a structural scaffold, this motif is also found in several plant protein families that are clustered in a single family (protease inhibitor/seed storage/LTP family) in the Pfam collection of protein families and domains . In order to identify the complete and non-redundant set of nsLtp genes in rice, we conducted an in silico analysis of the Oryza sativa subsp. japonica 'Nipponbare' genome. At the time of this study (November 2006), the Gramene database contained 101 genomic sequences annotated putative rice nsLtp genes. Each of the deduced protein sequences was manually assessed through the analysis of the cysteine residue patterns. The diversity of the retrieved 8 CM proteins enabled several cell wall glycoproteins to be distinguished including 23 glycosylphosphatidylinositol-anchored proteins characterized by a specific C-terminal sorting sequence , 21 proline-rich proteins and hybrid proline-rich proteins characterized by a high proportion of proline, histidine and glycine residues in the sequence comprised between the signal peptide and the 8 CM , and one glycine-rich protein  (Additional file 1). All these sequences displayed a supplementary motif (described above) not present in nsLTPs and were thus discarded. Other proteins were also discarded; they consist of three alpha-amylase/trypsin inhibitors which contain 10 cysteine residues engaged in five disulfide bonds , three prolamin storage proteins which lack the CXC motif and two 2S albumin storage proteins which present a molecular mass (MM) of about 20 kDa. Additionally, we eliminated two probable pseudogenes that have no corresponding transcripts indexed in the GenBank database and display mutation accumulations that result in the absence of the CC motif (Os04g09520) or a truncated 5' exon that curtails the signal peptide sequence (Os02g24720). As a result, only 46 out of the 101 genomic sequences initially annotated as putative nsLtp genes were found to encode proteins displaying the features of plant nsLTPs (Table 1). In addition to the presence of a signal peptide and the 8 CM (C-Xn-C-Xn-CC-Xn-CXC-Xn-C-Xn-C), the major feature we observed was a generally small MM (6.5 to 10.5 kDa), criteria that were those of type I and II nsLTPs described as having a lipid transfer activity [1, 56].
Next, a search for misannotated putative nsLtp genes was performed by blastn and tblastn searches of the TIGR Rice Pseudomolecules  using as query sequences the 46 rice genes and the 35 previously identified wheat nsLTPs and nsLtp genes . This approach resulted in the identification of six additional putative nsLtp genes leading to a total of 52 rice nsLtp genes (Table 1). These new genes were originally not annotated as putative nsLtp genes (Os01g58660, Os03g44000, Os09g35700, Os11g02424) or the presence of a frame shift in the coding region failed to identify the deduced proteins as putative nsLTPs (Os11g02330, Os11g02379.1).
The Arabidopsis thaliana nsLtpgene family is composed of 49 members
The same approach was used for arabidopsis. Locus annotations and protein domain descriptions allowed the identification of 112 loci that potentially encode nsLTPs. Analysis of protein primary sequences indicated that 31 of them encode glycosylphosphatidylinositol-anchored proteins, 25 encode hybrid proline-rich proteins and five encode 2S albumin storage proteins that were eliminated (Additional file 1). Three other loci were also discarded since the corresponding deduced protein failed to present an 8 CM (At1g21360, At2g33470, At3g21260). As a result, only 48 out of the 112 loci were found to encode putative nsLTPs (Table 2). Finally, blastn and tblastn searches allowed us to identify one new locus (At1g52415) that encodes an 8 CM protein with no homology with known Pfam domains.
Organization and structure of the rice and arabidopsis nsLtpgenes
Analysis of the physical chromosomal loci revealed that 26 out of the 52 rice nsLtp genes and 18 out of the 49 arabidopsis nsLtp genes are arranged in tandem duplication repeats (Figure 1). To cover nomenclature in different species, we named rice and arabidopsis nsLtp genes encoding nsLTPs OsLtp and AtLtp, respectively. Genes encoding mature proteins sharing more than 30% identity were grouped in the same type . Genes encoding rice and arabidopsis type I nsLTPs were named OsLtpI and AtLtpI respectively, and consecutive roman numbers were assigned for the other types.
In rice, two significant clusters of six type I nsLtp genes are found on chromosomes 11 and 12. A dot plot alignment of these two clusters clearly showed a co-linear segment that reveals high nucleotide sequence conservation, and indicated homologies between all nsLtp genes mainly limited to the ORFs (data not shown). Type II nsLtp genes are present as a cluster of six copies repeated in tandem on chromosome 10. Three direct repeat tandems were also identified on chromosome 1 (OsLtpII.1 and OsLtpII.2; OsLtpIV.1 and OsLtpIV.2; OsLtpVI.1 and OsLtpVI.2) and one on chromosome 4 (OsLtpV.2 and OsLtpV.3). Due to these duplications,nsLtp genes are over-represented on rice chromosomes 1, 10, 11 and 12, which carry 33 out of the 52 identified genes. On the contrary, no nsLtp genes were identified on chromosome 2.
In arabidopsis, 18 nsLtp genes were found organized in seven direct repeat tandems. Whereas one tandem of three repeats is present on chromosome 1 (AtLtpII.1, AtLtpII.2, and AtLtpII.3) and one tandem of two repeats is present on both chromosome 2 (AtLtpI.4 and AtLtpI.5) and 3 (AtLtpI.7 and AtLtpI.8), four direct repeat tandems are found on chromosome 5. With two to four repeats, these four tandems lead to the over-representation of nsLtp genes on arabidopsis chromosome 5.
With the exception of the AtLtpIV.3 and AtLtpIV.5 genes, no introns were identified in the coding regions of type II and IV rice and arabidopsis nsLtp genes and type IX arabidopsis nsLtp genes. On the contrary, all the type I, III, V and VI rice and arabidopsis nsLtp genes (except the AtLtpI.5 and AtLtpIII.2 genes) were predicted to be interrupted by a single intron positioned 2 to 73 bp upstream of the stop codon.
Identification of T. aestivum nsLtpgenes by EST database mining
Because the genome of T. aestivum has not yet been sequenced, we aimed to identify new members of the wheat nsLtp gene family by EST database mining. Since we observed strong homologies between many of the 52 rice nsLtp genes, the mismatches consented during the assembly of wheat ESTs in tentative consensus sequences or UniGene clusters (indexed in the TIGR Wheat Gene Index Database and in the NCBI UniGene database, respectively) make these last not appropriate for the identification of novel wheat nsLtp genes. Consequently, blast searches were performed against the wheat ESTs indexed in the GenBank database and collected from 239 T. aestivum cDNA libraries. To this end, we used the coding sequence of each of the 52 rice nsLtp genes listed in Table 1 and each of the 32 wheat genomic and cDNA sequences identified by Boutrot et al. 2007 .
ClustalW multiple-sequence alignments were performed for each blastn search. For each new putative wheat nsLtp gene identified, additional reiterative blastn searches were performed against the wheat EST database to identify additional related sequences. In total, this survey led to the identification of 156 putative wheat nsLtp genes (Table 3 and Additional file 2).
We applied to wheat nsLtp genes and proteins the nomenclature used for rice and arabidopsis (see above) and the eight types were named TaLtpI to TaLtpVIII. However, to consider the hexaploid status of the wheat genome we grouped wheat genes into subfamilies of putative homoeologous genes. This was based on the identity matrix (data not shown) calculated from the multiple sequence alignments and the nomenclature criteria that group mature proteins sharing more than 30% identity in a type and more than 75% identity in a subfamily . The 12 type I subfamilies were named TaLtpIa to TaLtpIl. Finally, the different members of each subfamily were differentiated by consecutive numbers, i.e. TaLtpIb.1 to TaLtpIb.39 for the 39 members of the type Ib subfamily. The correspondence between the previous nomenclature of wheat nsLtp genes  and the one used in this paper is shown in Additional file 2.
Since different T. aestivum cultivars were used to construct the cDNA libraries, the existence of probable variants of one gene may have resulted in overestimation of nsLtp gene diversity. Nevertheless, ESTs corresponding to at least 91 out of the 156 nsLtp genes were identified in the T. aestivum 'Chinese Spring' ('CS') cultivar. The identification of complete subfamily sets in single cultivars, such as the eight members of the TaLtpVa subfamily in the 'CS' cultivar, suggests that all the closely related genes of a subfamily reflect recent evolution of paralogous genes. We failed to identify any members of the TaLtpIe, TaLtpIf, TaLtpIi, TaLtpIk, TaLtpIl, TaLtpIVd, TaLtpVb, TaLtpVc, TaLtpVIIa and TaLtpVIIIa subfamilies in the 'CS' cultivar. However, most members of these subfamilies were identified in cDNA libraries prepared from specific plant material that were not used to construct 'CS' cDNA libraries.
Rice, arabidopsis and wheat nsLTP characteristics
The characteristics of the 52 rice and 49 arabidopsis putative nsLTPs are presented in Table 1 and Table 2, respectively. The MM and the theoretical pI of the 122 non-redundant wheat mature nsLTPs are summarized in Table 3 (details in Additional file 2).
Wheat, rice and arabidopsis nsLTPs are synthesized as pre-proteins that contain a putative signal peptide of 16 to 38 amino acids. The putative subcellular targeting of the 257 rice, arabidopsis and wheat nsLTP pre-protein sequences was analyzed using the TargetP 1.1 program and 255 of them present an N-terminal signal sequence that is thought to lead the mature protein through the secretory pathway. TaLTPIVb.3 and TaLTPIl.2 sequences have been predicted to contain a mitochondrial targeting peptide and a signal peptide. But, no conclusion could be drawn about the subcellular localization of these two mature proteins since the reliability of prediction was very weak.
At the pre-protein level, the OsLTPI.9 and OsLTPI.16 deduced proteins are identical. After cleavage of their signal peptide (predicted by the SignalP program), the OsLTPI.8 and OsLTPI.15 mature proteins are identical, as are the OsLTPI.12 and OsLTPI.19 mature proteins and the OsLTPI.13 and OsLTPI.20 mature proteins (Table 1). Therefore, before potential post-translational modifications, the 52 rice nsLtp genes encode 48 different mature nsLTPs. The 49 arabidopsis nsLtp genes encode proteins that are distinct in both their pre-protein and mature forms (Table 2). Thirty-four wheat proteins are redundant after cleavage of their signal peptide, 15 of them being redundant at the pre-protein level. Therefore, before potential post-translational modifications the 156 wheat putative nsLtp genes encode 122 different mature TaLTPs (Additional file 2). The TaLTPIf subfamily displays the strongest conservation since the four members have identical mature protein sequences. A high level of redundancy was also observed in genes of the TaLtpIg subfamily since five out of the eight members encode the same TaLTPIg.2 mature protein.
Since it allows all the cysteine residues to be maintained in a conserved position, the HMMalign program was preferred to ClustalW and was thus used to perform the multiple alignments of rice (Figure 2), arabidopsis (Figure 3) and wheat (Figure 4) nsLTPs. Based on the identity matrix (data not shown) calculated from the multiple sequence alignments and the nomenclature criteria that group mature proteins sharing more than 30% identity in a type , 49 out of the 52 rice nsLTPs, 45 out of the 49 arabidopsis nsLTPs and the 122 wheat nsLTPs were found to be clustered in nine types. The majority (147 out of 223) of the rice, arabidopsis and wheat nsLtp genes encode proteins that belong to the type I and type II nsLTPs. Fourteen rice, 15 arabidopsis and 34 wheat proteins described six new nsLTP types named types IV to IX. Three rice proteins and four arabidopsis proteins display less than 30% identity between themselves or with other nsLTPs to either make a type by themselves or be integrated in an already identified type. Therefore, these proteins were named OsLTPY.1 to OsLTPY.3 and AtLTPY.1 to AtLTPY.4.
Rice, wheat and arabidopsis nsLTPs are small proteins since their MMs usually range from 6636 Da to 10909 Da. However the OsLTPI.6 protein and the three members of the type VII wheat nsLTPs display unusual high MMs (13–15 kDa) due to the presence of supernumerary amino acid residues located at the C-terminal or N-terminal extremity of the deduced mature proteins. While the MM of nsLTPs previously allowed discrimination of the 9 kDa type I and the 7 kDa type II, type III nsLTPs were also found to present a MM of about 7 kDa. With nine nsLTP types identified, the relationship between MM and nsLTP type becomes more complex and is not anymore a good criterion to classify nsLTPs. The majority (199 out of 223) rice, wheat and arabidopsis non-redundant nsLTPs display a basic pI that is another characteristic of nsLTPs. In no case did nsLTPs with an acidic pI (3.92–5.50) form a specific type.
One characteristic of plant nsLTPs types I and II is the absence of tryptophane residues. Although this is usually the case, we found two type I (AtLTPI.2, AtLTPI.10), three type II (OsLTPII.1, AtLTPII.3, AtLTPII.11), four type IV (OsLTPIV.3, AtLTPIV.1, AtLTPIV.2, TaLTPIVb.1) and three nsLTPY proteins (OsLTPY.2, AtLTPY.1, AtLTPY.3) that contain one or two tryptophane residues.
The main characteristic of plant nsLTPs is the presence of eight cysteine residues in a strongly conserved position Cys1-Xn-Cys2-Xn-Cys3Cys4-Xn-Cys5XCys6-Xn-Cys7-Xn-Cys8. All the rice nsLTPs display this feature whereas two arabidopsis and two wheat nsLTPs present a different pattern. The Cys8 is missing in AtLTPI.1 and the Cys6 in AtLTPII.10. The TaLTPIVd.1 lacks Cys5 and Cys6 in the CXC motif and the TaLTPVIa.5 lacks the Cys7. Conversely, the members of the TaLTPIVa subfamilies, TaLTPIVc.1, OsLTPIV.1 and OsLTPIV.2 harbor an additional cysteine residue between Cys2 and Cys3, the TaLTPVIa subfamily members, OsLTPVI.1, OsLTPVI.2 OsLTPVI.4 and AtLTPII.10 between Cys6 and Cys7, AtLTPII.6 after Cys7, and the TaLTPVIIa subfamily members and OsLTPVII.1 after the Cys8 of the 8 CM.
The multiple alignment of the cysteine motifs of rice, arabidopsis and wheat nsLTPs also revealed a variable number of inter-cysteine amino acid residues (summarized in Figure 5). The AtLTPII.8 which is phylogenetically distant from all other type II nsLtp genes (see the phylogenetic analysis below) was not taken into consideration. In this way, seven nsLTP types can be identified through typical spacings for this motif. For example, type I nsLTPs contain 19 residues between the conserved Cys4 and Cys5 residues while types III, VII and VIII contain respectively 12, 27 and 25 residues between the conserved Cys6 and Cys7 residues. Similarly, types II, V and IX can be described with respectively 7, 14 and 13 residues between the conserved Cys1 and Cys2 residues. Only types IV and VI can not be distinguished based on this simple feature. A closer analysis of the sequences indicates that type VI nsLTPs are always characterized by a methionine and a valine residue present 10 and 4 aa before Cys7, respectively (Figures 2, 3, 4). At these positions, these two aa are always different in type IV nsLTPs and allow the direct distinction of type IV and VI nsLTPs.
Phylogenetic analysis of rice, arabidopsis and wheat nsLTPs
In order to analyze the phylogenetic organization of the nsLTP families, we constructed a phylogenetic tree from the alignment of respectively 45, 49 and 122 sequences of arabidopsis, rice and wheat nsLTPs, using the maximum-likelihood inference. Redundant mature wheat nsLTPs were eliminated but the arabidopsis and rice complete families were included. The solidity of the nodes was assessed by 100 bootstrap resampling repetitions. The seven arabidopsis and rice nsLTPY proteins were first included but due to the fact that their position was not well supported (nodes with weak bootstrap values) and consequently risked muddling the phylogenetic signal, they were excluded from the alignment. In the first attempt, several cysteine-rich protein sequences (metallothioneins, thionins and defensins from arabidopsis and rice) were tested as potential roots, but their position was different and none were supported by significant bootstrap values. Moreover, the phylogenetic relationships between types were not reliable whatever the root chosen. Consequently, we chose to present the complete condensed unrooted tree (Figure 6) where each of the subtrees (detailed in Figure 7) is rooted by all the other sequences.
The general organisation of the tree is coherent with the classification of nsLTPs in nine types. All the sequences belonging to the same type are grouped and constitute monophyletic groups (i.e. clades) except for type II nsLTPs. The bootstrap values supporting the clades corresponding to types III, V, VI, VII, VIII and IX are high, respectively 77, 100, 78, 95, 72 and 100. Types I and IV have lower bootstrap values, respectively 50 and 39. Based on the criteria that group mature proteins sharing more than 30% identity in a type, AtLTPIX.1 and AtLTPIX.2 were first included in type IV although their identity with other type IV nsLTPs was very low (12.6% to 30.1%). However, according to their position in the phylogenetic tree these sequences probably do not share the same common ancestor as other type IV nsLTPS and were classed in a new type named type IX. Type II nsLTPs are close in the tree but do not constitute a clade. This is mainly due to several A. thaliana nsLTPs (AtLTPII.1, AtLTPII.2, AtLTPII.3, AtLTPII.7, AtLTPII.8, AtLTPII.10, AtLTPII.11, AtLTPII.12, AtLTPII.13, AtLTPII.14 and AtLTPII.15), which appear to be more distantly related to other type II sequences. When the tree is built only with wheat and rice sequences, type II nsLTPs appear to be monophyletic and highly supported (bootstrap value 95; data not shown).
The distribution of nsLTPs in the tree is not either quantitatively or qualitatively homogeneous. As can be seen in Figure 7, there are significant differences in the number of sequences, with as few as two sequences for type IX nsLTPs and 90 for type I nsLTPs. Moreover, nsLTPs of each species are not homogeneously distributed within each type. Surprisingly, arabidopsis does not posses any type VII nsLTPs and no type IX nsLTPs were identified in rice and wheat.
Only type VIII nsLTPs displayed the simple organization that one would expect to be the most frequent between arabidopsis, rice and wheat, i.e. one sequence of each species (or three for the hexaploid wheat) with wheat and rice closer to each other and more distantly related to arabidopsis. Two other groups of sequences are organized in a similar way. The first group is composed of TaLTPVb.1, OsLTPV.1 and AtLTPV.1, however rice and arabidopsis are more closely related than wheat and rice. The second group is composed of AtLTPIV.1, AtLTPIV.2, OsLTPIV.3, TaLTPIVd.1 and TaLTPIVb.1. Even if a probably recent duplication in arabidopsis genome led to the presence of two copies, both are closely related to one copy of rice and two copies of wheat. In all the other cases, the arabidopsis sequences are either grouped and constitute a separated clade within a given type or branched close to the root of the type subtree. This is particularly true for AtLTPI.1, AtLTPI.4, AtLTPI.5, AtLTPI.6, AtLTPI.7, AtLTPI.8, AtLTPI.10, AtLTPI.11 and AtLTPI.12 or AtLTPIV.3, AtLTPIV.4 and AtLTPIV.5 or AtLTPVI.1, AtLTPVI.3 and AtLTPVI.4 or type II nsLTPs. In these cases, no obvious correspondence between arabidopsis and wheat/rice sequences exist and it is not possible to identify orthology relationships between nsLTP gene members of each species. A likely explanation may be that functions of nsLTPs are mostly due to a few conserved features indicating that functional domains or specific positions will be more conserved than others. Once these features are identified, it will become more relevant to perform fine phylogenetic analyses domain by domain.
The classification of the wheat nsLTP members in subfamilies when they share at least 75% amino acid identity appeared to agree with their phylogenetic relationships. Indeed, almost all the subfamilies appeared to be monophyletic (solid brackets in Figure 7) and are supported by high bootstrap values. Only two subfamilies present a more complex organization and are paraphyletic (i.e. they do not include all the members deriving from their common ancestor; dotted brackets in Figure 7). The TaLTPIIb subfamily clearly appears to be derived from TaLTPIIa. These two subfamilies share a common ancestor (node highly supported: 93), but TaLTPIIb members appear to have diverged from the others as the branch grouping them is longer and the node highly supported (98). Another subfamily, TaLTPIj, harbors surprising characteristics since the three wheat sequences are identical to three nsLTP rice copies (OsLTPI.10, OsLTPI.11 and OsLTPI.18). In contrast, we observed wheat nsLTP subfamilies (TaLTPIa, TaLTPIb, TaLTPIi, TaLTPIIa, TaLTPIIb, TaLTPIVa and TaLTPIVd indicated with green brackets in Figure 7) in which the closest related rice nsLTP is already closer to another wheat nsLTP subfamily. These wheat nsLTPs correspond either to groups in which a closer copy existed in rice and was subsequently deleted, or to wheat copies that are undergoing an evolution process specific to wheat. Because this concerns a large number of genes and the largest TaLtpIb subfamily, the second hypothesis is more likely.
Encoded by multigene families, plant nsLTPs were clustered in three clades based on their primary structure . Here we report the genome-wide analysis of the nsLtp gene family in O. sativa 'Nipponbare' and A. thaliana, which enabled us to identify six additional clades.
Gene structures and chromosomal locations indicate that the complexity of the arabidopsis and rice nsLtp gene families is mainly due to tandem duplication repeats representing 16 of the 49 arabidopsis nsLtp genes and 26 of the 52 rice nsLtp genes. The arabidopsis genome has undergone several rounds of genome-wide duplication events, including polyploidy  which likely support this nsLtp gene complexity. The rice genome is also the result of an ancient whole-genome duplication, a recent segmental duplication and massive ongoing individual gene duplications . Characterized by Wang et al. 2005 , a large-scale segmental duplication is observed in rice chromosomes 11 and 12 and consists of blocks of 5.44 Mb and 4.27 Mb, respectively. Due to this genomic segmental duplication, a cluster of six tandem duplicated copies is present in both chromosomes.
Based on sequence identity, 35 rice nsLTPs and 30 arabidopsis nsLTPs are clustered in the previously described type I, type II and type III clades. Fourteen rice nsLTPs and 15 arabidopsis nsLTPs are clustered in the six new types identified in this work. In wheat, 58 out of the 122 non-redundant nsLTPs are type I nsLTPs, 29 belong to the type II and three are type III nsLTPs. Finally, 32 wheat nsLTPs were clustered in five of the new types.
The wheat EST survey failed to identify transcripts corresponding to seven genes or protein previously identified. In the case of the TaLtpIa.2, TaLtpIb.1,TaLtpIg.1, TaLtpIg.5 and TaLtpIh.1 genes, effective transcription is supported by isolation of cDNAs or protein. However, without cDNA or protein identified, the TaLtpIIa.5 and TaLtpIb.2 genomic sequences could be pseudogenes. In both cases, these seven haplotypes are possibly not detected in the EST databases analyzed because of inter-varietal polymorphism, or because of restricted or specific-tissue expression.
The phylogenetic tree revealed that the classification of nsLTP family members in types and subfamilies according to respectively 30% and 75% of amino acid identity enables a good representation of the organization of the family. All the types (except type II) and most of the resulting subfamilies are monophyletic and supported by convincing bootstrap values. The three species have members in all the types except arabidopsis in type VII and rice and wheat in type IX. Either type VII appeared specifically in rice/wheat lineage or has disappeared in arabidopsis. It would be interesting to trace its evolution at the monocot/dicot scale. The absence of type IX nsLTPs in rice and wheat suggests that type IX could be specific to dicot species. Search for type IX nsLTPs in other species whose whole genome was sequenced should allow confirmation of this point.
The distribution of the sequences of the three species is not homogenous. First, arabidopsis nsLTPs are grouped within types or isolated and branched close to the root of the type subtree (type II). The main conclusion we can draw from these observations is that the ancestral nsLTP gene family already included eight (or nine) types before separation between the lineage leading to arabidopsis and the lineage leading to wheat and rice, but that each type was probably represented by only one or two ancestral members. Subsequently, the family evolved specifically in each lineage in terms of copy number and speed of duplication or mutation accumulation. The alternative to this scenario would be that several copies of each type pre-existed in the ancestral nsLTP gene family before monocots and dicots diverged but that a large number of copies was lost. It would be interesting to test these hypotheses by adding nsLTPs from other species to the analysis when their complete genomic sequences become available.
Our phylogenetic approach turned out to be more informative about the evolutionary relationships of certain subfamilies, especially when based on probabilistic methods instead of computed distances. Indeed, two subfamilies (TaLTPIIa and TaLTPIj) appear to be paraphyletic, i.e. they do not include all the members derived from the same common ancestor. In the case of the TaLTPIIa subfamily, this is due to the fact that some members underwent a process of divergence which resulted in them being grouped in a different subfamily (TaLTPIIb). The TaLTPIj subfamily members appear to be grouped because they evolved not far from their closest common ancestor. Their surprisingly high level of conservation with rice nsLTPs reinforces this assumption. This subfamily groups members with common characteristics (high amino acid identity, slow evolution rates) but does not include all the descendants of the same ancestor and consequently does not represent a phylogenetic group. In conclusion, although grouping according to percentage identity may make sense, it is nevertheless important to perform a precise phylogenetic analysis to understand the relationships between the gene members. Within this context, the identification of conserved domains or residues will allow to use these specific regions to perform functional phylogenetic analysis.
Within the wheat nsLtp gene subfamilies for which we did not identify a closely related rice gene, it is amazing to find the largest wheat TaLtpIb, TaLtpIIa/TaLtpIIb gene subfamilies. The larger number of genes in these subfamilies may be the evolutionary consequence of adaptation to wheat-specific functions or various environmental changes.
Since synteny between homoeologous chromosomes was shown to be widely conserved in the hexaploid wheat T. aestivum , each gene identified should be related to two other homoeologous copies. However we report that, in single cultivars, nine nsLtp gene subfamilies had more than three members. In spite of the relaxed selective constraint often exerted on duplicated genes, the members of the subfamily share more than 75% identity, suggesting that recent duplications of nsLtp genes also occurred in the wheat genome. Diverged from a common ancestor 46 millions years ago, Oryza and Triticum species display remarkably similar genomic organization . However, with more than three wheat homoeologous copies identified for most of the related rice genes, the nsLtp genes family appears to be much bigger in Triticum than in the Oryza genome. It has often been suggested that polyploidy offers genome plasticity, which, in turn, increases the potential ability of newly formed species to adapt to new environmental conditions . When a family already presenting a high copy number at the diploid level is duplicated twice, the complexity of the redundancy and the possibilities of evolution it offers are vast. To understand the evolutionary pattern of the wheat nsLtp gene family, correct identification of homoeologous genes and classification of paralogous sequences is essential. To this end, gene-specific PCR primers will be designed allowing to amplify the different members of a subfamily and to determine their chromosomal locations using Chinese Spring aneuploid and deletion lines.
The high number of nsLtp genes in the hexaploid wheat T. aestivum is probably mainly due to gene duplication by polyploidization. Whether this leads to retention of function of duplicated genes or to functional diversification either at the level of gene expression or protein function remains to be determined. Depending on the species or on the gene family, both phenomena have been observed following polyploid-induced gene duplication .
By analyzing the complete nsLtp gene family in both rice and arabidopsis genome we identified six new types leading to a total of nine types of nsLTPs. The type VII was found only in rice and wheat whereas the type IX was only identified in arabidopsis. Wheat EST data mining emphasized the higher number of nsLtp genes and complexity of certain subfamilies. The diversity of rice, arabidopsis and wheat nsLTPs suggests that nsLTPs support different functions in plants. However, until such time as specific biological functions or functional domains are defined, it seems relevant to categorize plant nsLTPs on the basis of sequence similarity and/or phylogenetic clustering.
In silico identification of rice and arabidopsis nsLtpgenes
The Gramene rice genome database (TIGR pseudomolecule assembly release 4 of IRGSP finished sequence)  was searched for nsLtp gene sequences using the gene annotations. The TAIR arabidopsis genome database (TAIR release 6.0)  was searched for nsLtp genes annotated as encoding lipid transfer proteins and the entire arabidopsis proteome was searched for proteins displaying a HMMPfam domain PF00234 (Plant lipid transfer/seed storage/trypsin-alpha amylase inhibitor). Blastn and tblastn searches were further performed against both databases using the retrieved annotated gene sequences, the wheat nsLtp gene sequences and previously identified nsLTPs , and the wheat nsLtp gene sequences identified in this work. The putative rice and arabidopsis nsLtp gene sequences retrieved were then curated for intron-exon junction positions using the NetGen2 program , and from comparison with related EST sequences in the Gramene rice genome database. The amino acid sequences deduced from the newly identified rice and arabidopsis nsLtp genes were finally assessed through the analysis of the cysteine residue patterns.
Wheat EST database searches
The search for Triticum aestivum ESTs was performed by comparing the coding sequences of wheat and rice nsLtp genes against EST sequences available at NCBI  in blastn searches. Sequence hits with E-values of less than 10-4 and a bit score of 100 or more were identified as putative nsLtp homologues and extracted. EST multiple alignments were performed using the ClustalW program . When their ORF alignment overlapped, multiple ESTs were considered as derived from a single gene and resolved to a single representative EST. An ORF was considered as a new gene if at least one mutation was observed and if it was represented by at least two ESTs covering the complete ORF. Then the EST displaying the most widely represented sequence in the 3'- and 5'-UTR regions was chosen as representative of the new wheat nsLtp gene. Singleton ESTs and ESTs presenting incomplete ORF were not considered except when several of them support a novel ORF. For a limited number of genes (11), single EST sequences displaying full ORF were nevertheless taken into account when they were supported by multiple and overlapping incomplete EST sequences.
Amino-acid sequence analysis
Pre-proteins translated from the ORF of all nsLTP sequences were analyzed for presence of potential signal peptide cleavage sites using the SignalP 3.0 program . The subcellular localization of the mature protein was predicted using the TargetP 1.1 program . Following signal peptide removal, theoretical pI and MM were computed using the program provided at . Amino acid sequences were efficiently aligned to the Pfam profile HMM (glocal model) defined from the protease inhibitor/seed storage/LTP family  using HMMalign from the HMMER package . A sequence identity matrix of the mature nsLTP sequences was computed using BioEdit v188.8.131.52  enabling us to determine the gene subfamily assignment and their nomenclature following the guidelines proposed by Boutrot et al. .
Rice, arabidopsis and wheat amino-acid sequences were aligned to the Pfam glocal model using HMMalign. Because they were not informative and created aberrant multi alignments during the re-samplings procedure, a total of 47 sites were removed from the alignment (12 of them were represented by only one sequence and the 35 others were non or few-informative sites, among them 29 were only represented by the three type VII wheat nsLTPs). Phylogenetic trees were built from the protein alignment with the maximum-likelihood method using the PHYML program . Maximum-likelihood inference analyses were conducted under the Jones Taylor Thornton substitution model  with estimation of the proportion of invariant sites and estimation of variation rate among the remaining sites according to a gamma distribution. The confidence level of each node was estimated by the bootstrap procedure using 100 resampling repetitions of the data. The unrooted phylogenetic trees were visualized using the Treeview 1.6.6 program .
Kader JC, Julienne M, Vergnolle C: Purification and characterization of a spinach-leaf protein capable of transferring phospholipids from liposomes to mitochondria or chloroplasts. Eur J Biochem. 1984, 139 (2): 411-416. 10.1111/j.1432-1033.1984.tb08020.x.
José-Estanyol M, Gomis-Rüth FX, Puigdomènech P: The eight-cysteine motif, a versatile structure in plant proteins. Plant Physiol Biochem. 2004, 42 (5): 355-365. 10.1016/j.plaphy.2004.03.009.
Douliez JP, Michon T, Elmorjani K, Marion D: Structure, biological and technological functions of lipid transfer proteins and indolines, the major lipid binding proteins from cereal kernels. J Cereal Sci. 2000, 32 (1): 1-20. 10.1006/jcrs.2000.0315.
Gincel E, Simorre JP, Caille A, Marion D, Ptak M, Vovelle F: Three-dimensional structure in solution of a wheat lipid-transfer protein from multidimensional 1H-NMR data. A new folding for lipid carriers. Eur J Biochem. 1994, 226 (2): 413-422. 10.1111/j.1432-1033.1994.tb20066.x.
Lerche MH, Poulsen FM: Solution structure of barley lipid transfer protein complexed with palmitate. Two different binding modes of palmitate in the homologous maize and barley nonspecific lipid transfer proteins. Protein Sci. 1998, 7 (12): 2490-2498.
Hoh F, Pons JL, Gautier MF, de Lamotte F, Dumas C: Structure of a liganded type 2 non-specific lipid-transfer protein from wheat and the molecular basis of lipid binding. Acta crystallogr, D Biol Crystallogr. 2005, 61: 397-406. 10.1107/S0907444905000417.
Lauga B, Charbonnel-Campaa L, Combes D: Characterization of MZm3-3, a Zea mays tapetum-specific transcript. Plant Sci. 2000, 157 (1): 65-75. 10.1016/S0168-9452(00)00267-3.
Boutrot F, Guirao A, Alary R, Joudrier P, Gautier MF: Wheat non-specific lipid transfer protein genes display a complex pattern of expression in developing seeds. Biochim Biophys Acta, Gene Struct Exp. 2005, 1730 (2): 114-125.
Kader JC: Lipid-transfer proteins in plants. Annu Rev Plant Physiol Plant Mol Biol. 1996, 47: 627-654. 10.1146/annurev.arplant.47.1.627.
Sterk P, Booij H, Schellekens GA, van Kammen A, de Vries SC: Cell-specific expression of the carrot EP2 lipid transfer protein gene. Plant Cell. 1991, 3 (9): 907-921. 10.1105/tpc.3.9.907.
Broekaert WF, Cammue BPA, de Bolle MFC, Thevissen K, de Samblanx GW, Osborn RW: Antimicrobial peptides from plants. Crit Rev Plant Sci. 1997, 16 (3): 297-323. 10.1080/713608148.
García-Olmedo F, Molina A, Alamillo JM, Rodríguez-Palenzuéla P: Plant defense peptides. Biopolymers, Pept Sci. 1998, 47 (6): 479-491. 10.1002/(SICI)1097-0282(1998)47:6<479::AID-BIP6>3.0.CO;2-K.
Molina A, García-Olmedo F: Developmental and pathogen-induced expression of three barley genes encoding lipid transfer proteins. Plant J. 1993, 4 (6): 983-991. 10.1046/j.1365-313X.1993.04060983.x.
Guiderdoni E, Cordero MJ, Vignols F, García-Garrido JM, Lescot M, Tharreau D, Meynard D, Ferrière N, Notteghem JL, Delseny M: Inducibility by pathogen attack and developmental regulation of the rice Ltp1 gene. Plant Mol Biol. 2002, 49 (6): 683-699. 10.1023/A:1015595100145.
Gomès E, Sagot E, Gaillard C, Laquitaine L, Poinssot B, Sanejouand YH, Delrot S, Coutos-Thévenot P: Nonspecific lipid-transfer protein genes expression in grape (Vitis sp.) cells in response to fungal elicitor treatments. Mol Plant Microbe Interact. 2003, 16 (5): 456-464. 10.1094/MPMI.2003.16.5.456.
Jung HW, Kim W, Hwang BK: Three pathogen-inducible genes encoding lipid transfer protein from pepper are differentially activated by pathogens, abiotic, and environmental stresses. Plant Cell Environ. 2003, 26 (6): 915-928. 10.1046/j.1365-3040.2003.01024.x.
Lu ZX, Gaudet DA, Frick M, Puchalski B, Genswein B, Laroche A: Identification and characterization of genes differentially expressed in the resistance reaction in wheat infected with Tilletia tritici, the common bunt pathogen. J Biochem Mol Biol. 2005, 38 (4): 420-431.
Molina A, García-Olmedo F: Enhanced tolerance to bacterial pathogens caused by the transgenic expression of barley lipid transfer protein LTP2. Plant J. 1997, 12 (3): 669-675. 10.1046/j.1365-313X.1997.00669.x.
van Loon LC, van Strien EA: The families of pathogenesis-related proteins, their activities, and comparative analysis of PR-1 type proteins. Physiol Mol Plant Pathol. 1999, 55 (2): 85-97. 10.1006/pmpp.1999.0213.
Maldonado AM, Doerner P, Dixon RA, Lamb CJ, Cameron RK: A putative lipid transfer protein involved in systemic resistance signalling in Arabidopsis. Nature. 2002, 419: 399-403. 10.1038/nature00962.
Buhot N, Douliez JP, Jacquemard A, Marion D, Tran V, Maume B, Milat ML, Ponchet M, Mikes V, Kader JC, Blein JP: A lipid transfer protein binds to a receptor involved in the control of plant defence responses. FEBS Lett. 2001, 509 (1): 27-30. 10.1016/S0014-5793(01)03116-7.
Edqvist J, Farbos I: Characterization of germination-specific lipid transfer proteins from Euphorbia lagascae. Planta. 2002, 215 (1): 41-50. 10.1007/s00425-001-0717-x.
Gonorazky AG, Regente MC, de la Canal L: Stress induction and antimicrobial properties of a lipid transfer protein in germinating sunflower seeds. J Plant Physiol. 2005, 162: 618-624. 10.1016/j.jplph.2004.10.006.
Soufleri IA, Vergnolle C, Miginiac E, Kader JC: Germination-specific lipid transfer protein cDNAs in Brassica napus L. Planta. 1996, 199 (2): 229-237. 10.1007/BF00196563.
Foster GD, Robinson SW, Blundell RP, Roberts MR, Hodge R, Draper J, Scott RJ: A Brassica napus mRNA encoding a protein homologous to phospholipid transfer proteins, is expressed specifically in the tapetum and developing microspores. Plant Sci. 1992, 84 (2): 187-192. 10.1016/0168-9452(92)90133-7.
Ariizumi T, Amagai M, Shibata D, Hatakeyama K, Watanabe M, Toriyama K: Comparative study of promoter activity of three anther-specific genes encoding lipid transfer protein, xyloglucan endotransglucosylase/hydrolase and polygalacturonase in transgenic Arabidopsis thaliana. Plant Cell Rep. 2002, 21 (1): 90-96. 10.1007/s00299-002-0487-3.
Imin N, Kerim T, Weinman JJ, Rolfe BG: Low temperature treatment at the young microspore stage induces protein changes in rice anthers. Mol Cell Proteomics. 2006, 5 (2): 274-292.
Liu K, Jiang H, Moore S, Watkins C, Jahn M: Isolation and characterization of a lipid transfer protein expressed in ripening fruit of Capsicum chinense. Planta. 2006, 223 (4): 672-683. 10.1007/s00425-005-0120-0.
Feng JX, Ji SJ, Shi YH, Wei G, Zhu YX: Analysis of five differentially expressed gene families in fast elongating cotton fiber. Acta Biochim Biophys Sin. 2004, 36 (1): 51-57.
Kinlaw CS, Gerttula SM, Carter MC: Lipid transfer protein genes of loblolly pine are members of a complex gene family. Plant Mol Biol. 1994, 26 (4): 1213-1216. 10.1007/BF00040702.
Arondel V, Vergnolle C, Cantrel C, Kader JC: Lipid transfer proteins are encoded by a small multigene family in Arabidopsis thaliana. Plant Sci. 2000, 157 (1): 1-12. 10.1016/S0168-9452(00)00232-6.
Boutrot F, Meynard D, Guiderdoni E, Joudrier P, Gautier MF: The Triticum aestivum non-specific lipid transfer protein (TaLtp) gene family: comparative promoter activity of six TaLtp genes in transgenic rice. Planta. 2007, 225 (4): 843-862. 10.1007/s00425-006-0397-7.
The Arabidopsis Genome Initiative: Analysis of the genome sequence of the flowering plant Arabidopsis thaliana. Nature. 2000, 408 (6814): 796-815. 10.1038/35048692.
Yu J, Hu S, Wang J, Wong GKS, Li S, Liu B, Deng Y, Dai L, Zhou Y, Zhang et : A draft sequence of the rice genome (Oryza sativa L. ssp. indica). Science. 2002, 296 (5565): 79-92. 10.1126/science.1068037.
International Rice Genome Sequencing Project: The map-based sequence of the rice genome. Nature. 2005, 436 (7052): 793-800. 10.1038/nature03895.
Populus trichocarpa genome assembly 1.0. [http://genome.jgi-psf.org/Poptr1/Poptr1.home.html]
The French-Italian Public Consortium for Grapevine Genome Characterization: The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature. 2007, 449 (7161): 463-467. 10.1038/nature06148.
Chen F, Li Q, Sun L, He Z: The rice 14-3-3 gene family and its involvement in responses to biotic and abiotic stress. DNA Res. 2006, 13 (2): 53-63. 10.1093/dnares/dsl001.
Englbrecht C, Schoof H, Bohm S: Conservation, diversification and expansion of C2H2 zinc finger proteins in the Arabidopsis thaliana genome. BMC Genomics. 2004, 5 (1): 39-10.1186/1471-2164-5-39.
Yuan J, Yang X, Lai J, Lin H, Cheng ZM, Nonogaki H, Chen F: The endo-beta-mannanase gene families in Arabidopsis, rice, and poplar. Funct Integr Genomics. 2007, 7 (1): 1-16. 10.1007/s10142-006-0034-3.
Arumuganathan K, Earle ED: Nuclear DNA content of some important plant species. Plant Mol Biol Rep. 1991, 211-215.
Ogihara Y, Mochida K, Nemoto Y, Murai K, Yamazaki Y, Shin-I T, Kohara Y: Correlated clustering and virtual display of gene expression patterns in the wheat life cycle by large-scale statistical analyses of expressed sequence tags. Plant J. 2003, 33 (6): 1001-1011. 10.1046/j.1365-313X.2003.01687.x.
Wilson ID, Barker GLA, Beswick RW, Shepherd SK, Lu C, Coghill JA, Edwards D, Owen P, Lyons R, Parker JS, Lenton JR, Holdsworth MJ, Shewry PR, Edwards KJ: A transcriptomics resource for wheat functional genomics. Plant Biotechnol J. 2004, 2 (6): 495-506. 10.1111/j.1467-7652.2004.00096.x.
Zhang D, Choi DW, Wanamaker S, Fenton RD, Chin A, Malatrasi M, Turuspekov Y, Walia H, Akhunov ED, Kianian P, Otto C, Simons K, Deal KR, Echenique V, Stamova B, Ross K, Butler GE, Strader L, Verhey SD, Johnson R, Altenbach S, Kothari K, Tanaka C, Shah MM, Laudencia-Chingcuanco D, Han P, Miller RE, Crossman CC, Chao S, Lazo GR, Klueva N, Gustafson JP, Kianian SF, Dubcovsky J, Walker-Simmons MK, Gill KS, Dvorak J, Anderson OD, Sorrells ME, McGuire PE, Qualset CO, Nguyen HT, Close TJ: Construction and evaluation of cDNA libraries for large-scale expressed sequence tag sequencing in wheat (Triticum aestivum L.). Genetics. 2004, 168 (2): 595-608. 10.1534/genetics.104.034785.
Mochida K, Kawaura K, Shimosaka E, Kawakami N, Shin-I T, Kohara Y, Yamazaki Y, Ogihara Y: Tissue expression map of a large number of expressed sequence tags and its application to in silico screening of stress response genes in common wheat. Mol Genet Genomics. 2006, 276 (3): 304-312. 10.1007/s00438-006-0120-1.
Adams MD, Kelley JM, Gocayne JD, Dubnick M, Polymeropoulos MH, Xiao H, Merril CR, Wu A, Olde B, Moreno RF, Kerlavage AR, McCombie WR, Venter JC: Complementary DNA sequencing: expressed sequence tags and human genome project. Science. 1991, 252 (5013): 1651-1656. 10.1126/science.2047873.
Boguski MS, Tolstoshev CM, Bassett DEJ: Gene discovery in dbEST. Science. 1994, 265 (5181): 1993-1994. 10.1126/science.8091218.
Jukanti AK, Bruckner PL, Fischer AM: Evaluation of wheat polyphenol oxidase genes. Cereal Chem. 2004, 81 (4): 481-485. 10.1094/CCHEM.2004.81.4.481.
Kawaura K, Mochida K, Ogihara Y: Expression profile of two storage-protein gene families in hexaploid wheat revealed by large-scale analysis of expressed sequence tags. Plant Physiol. 2005, 139 (4): 1870-1880. 10.1104/pp.105.070722.
Kruger WM, Pritsch C, Chao SM, Muehlbauer GJ: Functional and comparative bioinformatic analysis of expressed genes from wheat spikes infected with Fusarium graminearum. Mol Plant Microbe Interact. 2002, 15 (5): 445-455. 10.1094/MPMI.2002.15.5.445.
Pfam collection of protein families and domains. [http://www.sanger.ac.uk/Software/Pfam]
Borner GHH, Lilley KS, Stevens TJ, Dupree P: Identification of glycosylphosphatidylinositol-anchored proteins in Arabidopsis. A proteomic and genomic analysis. Plant Physiol. 2003, 132 (2): 568-577. 10.1104/pp.103.021170.
Jose-Estanyol M, Puigdomènech P: Plant cell wall glycoproteins and their genes. Plant Physiol Biochem. 2000, 38 (1-2): 97-108. 10.1016/S0981-9428(00)00165-0.
Sachetto-Martins G, Franco LO, de Oliveira DE: Plant glycine-rich proteins: a family or just proteins with a common motif?. Biochim Biophys Acta, Gene Struct Exp. 2000, 1492 (1): 1-14.
Franco OL, Rigden DJ, Melo FR, Grossi-de-Sá MF: Plant alpha-amylase inhibitors and their interaction with insect alpha-amylases. Structure, function and potential for crop protection. Eur J Biochem. 2002, 269 (2): 397-412. 10.1046/j.0014-2956.2001.02656.x.
Monnet FP, Dieryck W, Boutrot F, Joudrier P, Gautier MF: Purification, characterisation and cDNA cloning of a type 2 (7 kDa) lipid transfer protein from Triticum durum. Plant Sci. 2001, 161 (4): 747-755. 10.1016/S0168-9452(01)00459-9.
Ware D, Jaiswal P, Ni J, Pan X, Chang K, Clark K, Teytelman L, Schmidt S, Zhao W, Cartinhour S, McCouch S, Stein L: Gramene: a resource for comparative grass genomics. Nucleic Acids Res. 2002, 30 (1): 103-105. 10.1093/nar/30.1.103.
Blanc G, Hokamp K, Wolfe KH: A recent polyploidy superimposed on older large-scale duplications in the Arabidopsis genome. Genome Res. 2003, 13 (2): 137-144. 10.1101/gr.751803.
Yu J, Wang J, Lin W, Li S, Li H, Zhou J, Ni P, Dong W, Hu S, Zeng C, Zhang J, Zhang Y, Li R, Xu Z, Li S, Li X, Zheng H, Cong L, Lin L, Yin J, Geng J, Li G, Shi J, Liu J, Lv H, Li J, Wang J, Deng Y, Ran L, Shi X, Wang X, Wu Q, Li C, Ren X, Wang J, Wang X, Li D, Liu D, Zhang X, Ji Z, Zhao W, Sun Y, Zhang Z, Bao J, Han Y, Dong L, Ji J, Chen P, Wu S, Liu J, Xiao Y, Bu D, Tan J, Yang L, Ye C, Zhang J, Xu J, Zhou Y, Yu Y, Zhang B, Zhuang S, Wei H, Liu B, Lei M, Yu H, Li Y, Xu H, Wei S, He X, Fang L, Zhang Z, Zhang Y, Huang X, Su Z, Tong W, Li J, Tong Z, Li S, Ye J, Wang L, Fang L, Lei T, Chen C, Chen H, Xu Z, Li H, Huang H, Zhang F, Xu H, Li N, Zhao C, Li S, Dong L, Huang Y, Li L, Xi Y, Qi Q, Li W, Zhang B, Hu W, Zhang Y, Tian X, Jiao Y, Liang X, Jin J, Gao L, Zheng W, Hao B, Liu S, Wang W, Yuan L, Cao M, McDermott J, Samudrala R, Wang J, Wong GKS, Yang H: The genomes of Oryza sativa: A history of duplications. PLoS Biology. 2005, 3 (2): e38-10.1371/journal.pbio.0030038.
Wang X, Shi X, Hao B, Ge S, Luo J: Duplication and DNA segmental loss in the rice genome: implications for diploidization. New Phytol. 2005, 165 (3): 937-946. 10.1111/j.1469-8137.2004.01293.x.
Akhunov ED, Akhunova AR, Linkiewicz AM, Dubcovsky J, Hummel D, Lazo GR, Chao S, Anderson OD, David J, Qi L, Echalier B, Gill BS, Miftahudin, Gustafson JP, La Rota M, Sorrells ME, Zhang D, Nguyen HT, Kalavacharla V, Hossain K, Kianian SF, Peng J, Lapitan NLV, Wennerlind EJ, Nduati V, Anderson JA, Sidhu D, Gill KS, McGuire PE, Qualset CO, Dvorak J: Synteny perturbations between wheat homoeologous chromosomes caused by locus duplications and deletions correlate with recombination rates. Proc Natl Acad Sci USA. 2003, 100 (19): 10836-10841. 10.1073/pnas.1934431100.
Gaut BS: Evolutionary dynamics of grass genomes. New Phytol. 2002, 154 (1): 15-28. 10.1046/j.1469-8137.2002.00352.x.
Moore RC, Purugganan MD: The evolutionary dynamics of plant duplicate genes . Curr Opin Plant Biol. 2005, 8 (2): 122-128. 10.1016/j.pbi.2004.12.001.
Wendel JF: Genome evolution in polyploids. Plant Mol Biol. 2000, 42 (1): 225-249. 10.1023/A:1006392424384.
Gramene Rice Genome Database. [http://www.gramene.org]
The Arabidopsis Information Resource (TAIR). [http://www.arabidopsis.org]
Hebsgaard SM, Korning PG, Tolstrup N, Engelbrecht J, Rouze P, Brunak S: Splice site prediction in Arabidopsis thaliana pre-mRNA by combining local and global sequence information. Nucleic Acids Res. 1996, 24 (17): 3439-3452. 10.1093/nar/24.17.3439.
NCBI Expressed Sequence Tags database. [http://www.ncbi.nlm.nih.gov/dbEST/index.html]
Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22 (22): 4673-4680. 10.1093/nar/22.22.4673.
Bendtsen JD, Nielsen H, von Heijne G, Brunak S: Improved prediction of signal peptides: SignalP 3.0. J Mol Biol. 2004, 340 (4): 783-795. 10.1016/j.jmb.2004.05.028.
Emanuelsson O, Nielsen H, Brunak S, von Heijne G: Predicting subcellular localization of proteins based on their N-terminal aAmino acid sequence. J Mol Biol. 2000, 300 (4): 1005-1016. 10.1006/jmbi.2000.3903.
Masse moléculaire, pI, composition, courbe de titrage. [http://www.iut-arles.up.univ-mrs.fr/w3bb/d_abim/compo-p.html]
Eddy SR: Profile hidden Markov models. Bioinformatics. 1998, 14 (9): 755-763. 10.1093/bioinformatics/14.9.755.
Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser. 1999, 41: 95-98.
Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52 (5): 696-704. 10.1080/10635150390235520.
Jones DT, Taylor WR, Thornton JM: The rapid generation of mutation data matrices from protein sequences. Comput Appl Biosci. 1992, 8 (3): 275-282.
Page RDM: TREEVIEW: An application to display phylogenetic trees on personal computers. Comput Appl Biosci. 1996, 12: 357-358.
The authors wish to thank Jean-Pascal Sirven for his help in collecting the wheat EST sequences. FB was the recipient of a fellowship from the French Ministère de l'Education Nationale, de l'Enseignement Supérieur et de la Recherche. The authors also thank Alberto Cenci and Stéphane De Mita for helpful discussions.
FB carried out rice and wheat database searches, comparative genome analysis, gene structure prediction and nomenclature, and drafted the manuscript. NC carried out the phylogenetic analysis, contributed to the collection of the wheat EST sequences and to the writing of the manuscript. MFG coordinated the study and contributed to the writing of the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 3: Alignment of the rice, arabidopsis and wheat nsLTP sequences. The mature sequences of the 122 non-redundant wheat nsLTPs, the 49 rice nsLTPs, and the 45 arabidopsis nsLTPs were aligned using HMMalign and then manually refined. The phylogenetic tree was built from this protein alignment (fasta format). (DOC 96 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Boutrot, F., Chantret, N. & Gautier, MF. Genome-wide analysis of the rice and arabidopsis non-specific lipid transfer protein (nsLtp) gene families and identification of wheat nsLtp genes by EST data mining. BMC Genomics 9, 86 (2008). https://doi.org/10.1186/1471-2164-9-86
- Hexaploid Wheat
- Mature Protein
- Segmental Duplication
- nsLtp Gene
- Homoeologous Copy