This Article Abstract Full Text (PDF) Supplemental Data All Versions of this Article: 19/6/1947    most recent tpc.107.051391v1 Alert me when this article is cited Alert me if a correction is posted Services Email this article to a friend Similar articles in this journal Similar articles in ISI Web of Science Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via CrossRef Citing Articles via ISI Web of Science (13) Citing Articles via Google Scholar Google Scholar Articles by Baumann, M. J. Articles by Brumer, H. Search for Related Content PubMed PubMed Citation Articles by Baumann, M. J. Articles by Brumer, H., III Agricola Articles by Baumann, M. J. Articles by Brumer, H.

Structural Evidence for the Evolution of Xyloglucanase Activity from Xyloglucan Endo-Transglycosylases: Biological Implications for Cell Wall Metabolism^[W]

^a School of Biotechnology, Royal Institute of Technology, AlbaNova University Center, S-10691 Stockholm, Sweden
^b Equipe Glycobiologie Marine, Unité Mixte de Recherche 7139, Végétaux Marins et Biomolécules, Centre National de la Recherche Scientifique/Université Pierre et Marie Curie-Paris 6, Station Biologique, 29682 Roscoff, Bretagne, France

¹ To whom correspondence should be addressed. E-mail czjzek{at}sb-roscoff.fr or harry{at}biotech.kth.se; fax 33-298-292-324 or 46-8-5537-8367.

	ABSTRACT

TOP ABSTRACT INTRODUCTION RESULTS DISCUSSION METHODS REFERENCES

 
High-resolution, three-dimensional structures of the archetypal glycoside hydrolase family 16 (GH16) endo-xyloglucanases Tm-NXG1 and Tm-NXG2 from nasturtium (Tropaeolum majus) have been solved by x-ray crystallography. Key structural features that modulate the relative rates of substrate hydrolysis to transglycosylation in the GH16 xyloglucan-active enzymes were identified by structure–function studies of the recombinantly expressed enzymes in comparison with data for the strict xyloglucan endo-transglycosylase Ptt-XET16-34 from hybrid aspen (Populus tremula x Populus tremuloides). Production of the loop deletion variant Tm-NXG1-YNIIG yielded an enzyme that was structurally similar to Ptt-XET16-34 and had a greatly increased transglycosylation:hydrolysis ratio. Comprehensive bioinformatic analyses of XTH gene products, together with detailed kinetic data, strongly suggest that xyloglucanase activity has evolved as a gain of function in an ancestral GH16 XET to meet specific biological requirements during seed germination, fruit ripening, and rapid wall expansion.

	INTRODUCTION

TOP ABSTRACT INTRODUCTION RESULTS DISCUSSION METHODS REFERENCES

 
The xyloglucan (XG) family of polysaccharides continues to receive considerable attention due to the key structural role these matrix polymers play in the cell walls of higher plants. Xyloglucans are major components of the primary wall of all dicots and some monocots (Carpita and McCann, 2000), where they bridge paracrystalline cellulose microfibrils by surface adsorption and chain intercalation (Pauly et al., 1999a), thereby modulating wall mechanical properties and affecting morphology (Rose and Bennett, 1999; Chanliaud et al., 2004; Cosgrove, 2005; Thompson, 2005; Brummell, 2006, and references therein). Indeed, based on the broad distribution of xyloglucans in vascular plants, it has been suggested that the evolution of xyloglucan as a cellulose binding polysaccharide may have conferred a particular advantage in the colonization of drier habitats (Popper and Fry, 2003, 2004). In addition to performing structural roles, xyloglucans are also found as the principal seed-storage polysaccharides in a number of land plants (Kooiman, 1960; Rao and Srivastava, 1973; Edwards et al., 1985; Wang et al., 1996; Buckeridge et al., 2000; Ren et al., 2004).

The primary structures of xyloglucans possess a common ß(14)-D-glucan backbone that is regularly substituted with (16)-D-xylopyranosyl residues. Two predominant xylosylation patterns have been observed in primary wall and storage xyloglucans (Vincken et al., 1997), both based on Glc₄ backbone repeating structures: XXXG and XXGG (where X is [-D-Xylp(16)]-Glcpß(14) and G is Glcpß(14), according to the nomenclature of Fry et al. [1993]). More rarely, Glc₅- and Glc₆-based backbone repeats have also been observed (Sims et al., 1996; Buckeridge et al., 1997; Tine et al., 2006). These basic motifs are further decorated with combinations of galactopyranose, fucopyranose, arabinofuranose, and O-acetyl residues, dependent on the tissue and species of origin, to produce a diversity of xyloglucan structures (see Hoffman et al., 2005, for a recent overview of xyloglucan phylogeny). The archetypal seed xyloglucan from Tamarindus indica is composed of XXXG, XXLG, XLXG, and XLLG oligosaccharides (where L is [ß-D-Galp-(12)--D-Xylp(16)]-Glcpß(14)) in the molar ratio 1.4:3:1:5.4 (York et al., 1990), while the majority of primary wall xyloglucans also contain XXFG (where F is L bearing an -L-Fucp(12) substituent) (Carpita and McCann, 2000; Hoffman et al., 2005).

Considerable interest in the enzymes responsible for the metabolism of xyloglucans is sustained as a result of the physiological importance these polysaccharides have in cell wall morphogenesis. The nucleotide diphospho-sugar glycosyl transferases responsible for the anabolism of xyloglucan are beginning to be revealed (Reiter, 2002; Lerouxel et al., 2006), while recent three-dimensional structural analysis is highlighting the molecular details of xyloglucan depolymerization by microbial enzymes in glycoside hydrolase families (Henrissat and Davies, 1997) GH5 (Gloster et al., 2007), GH12 (Pauly et al., 1999b; Gloster et al., 2007), and GH74 (Yaoi et al., 2004; Martinez-Fleites et al., 2006). In plants, specific endolytic activity toward xyloglucan has only been demonstrated for members of glycoside hydrolase family GH16, encoded by the XTH (for xyloglucan endo-transglycosylase/hydrolase) gene subfamily (Rose et al., 2002). Like other GH16 enzymes, which include a range of endo-glucan and galactanases (Michel et al., 2001), the gene products of the XTH subfamily utilize a double displacement/retaining mechanism of glycosyl transfer (Figure 1 ).

View larger version (10K): [in this window] [in a new window] [as a PowerPoint slide]   Figure 1. The Canonical Retaining Mechanism of Glycosyl Transfer Showing Product Partitioning between Hydrolysis and Transglycosylation. (A) The chemical mechanism. (B) Full kinetic scheme for the minimal xylogluco-oligosaccharide donor substrate XXXGXXXG under steady state conditions in the absence of products. Substrates (A = XXXGXXXG, B = H2O), enzyme intermediates (E, EA, F, FA, and FB), and products (P1 to P3) are denoted with uppercase letters, and individual rate constants (kn) are shown.

The majority of the XTH gene products for which full kinetic data exist have been shown to have very little or no detectable hydrolytic activity (see Supplemental Table 1 online) and are thus predominant or strict xyloglucan endo-transglycosylases (XET; EC 2.4.1.207) (Fry et al., 1992; Nishitani and Tominaga, 1992). Ptt-XET16-34 (previously Ptt-XET16A) from the hybrid aspen Populus tremula x Populus tremuloides is one such example that has undetectable xyloglucanase activity (EC 3.2.1.151; alt. xyloglucan endo-ß-1,4-glucanase [XEG] or xyloglucan endo-hydrolase [XEH]) (Kallas et al., 2005). The three-dimensional structure of this enzyme has been solved (Johansson et al., 2004), and the catalytic mechanism has been probed using a range of synthetic xylogluco-oligosaccharides to elucidate details of enzyme–substrate interactions (Fauré et al., 2006). Intriguingly, xyloglucanase activity and gene sequence have been conclusively linked for only one member of GH16, Tm-NXG1 from Tropaeolum majus (nasturtium), which is a predominant endo-xyloglucanase that can also perform xyloglucan endo-transglycosylation at elevated concentrations of acceptor substrates (Edwards et al., 1986; de Silva et al., 1993; Fanutti et al., 1993, 1996; Chanliaud et al., 2004).

As a next step toward deciphering the determinants of xyloglucan transglycosylation versus hydrolysis in the plant XTH gene subfamily of GH16, we undertook a detailed structure–function study of Tm-NXG1 and the homologous Tm-NXG2. We report here the three-dimensional structures of Tm-NXG1 and Tm-NXG2 as examples of plant GH16 xyloglucanases and highlight key active site features that differ from the previously determined XET structure (Johansson et al., 2004). Detailed bioinformatic analysis was used to guide the construction of a hybrid enzyme that demonstrated intermediate kinetic activity, which effectively indicated that XET activity could be increased at the expense of xyloglucanase activity. These results are discussed in the context of the molecular evolution of GH16 xyloglucan-active enzymes to meet specific biological requirements for cell wall remodeling and storage polysaccharide mobilization.

	RESULTS

TOP ABSTRACT INTRODUCTION RESULTS DISCUSSION METHODS REFERENCES

 
Cloning of Tm-NXG2 cDNA and Mutagenesis to Produce Tm-NXG1 and the Deletion Mutant Tm-NXG1-YNIIG
Tm-NXG2 (GenBank accession number X68255) represents a C-terminal sequence fragment that displays eight nucleotide polymorphisms relative to the full-length Tm-NXG1 sequence (GenBank accession number X68254), only one of which gives rise to variation in the protein primary structure: Met-129 in Tm-NXG1 is replaced by Leu in Tm-NXG2. In this study, a unique full-length cDNA was cloned from a pool obtained from the epicotyl of nasturtium seedlings that was 99.3% identical to Tm-NXG1 but that encoded the Tm-NXG2–like Leu-129 variant. To study both proteins, site-directed mutagenesis was used to introduce a L129M point mutation to produce a cDNA encoding Tm-NXG1 from the novel Tm-NXG2 sequence. The historical abbreviations NXG1 and NXG2 are retained throughout for direct comparison with the extensive literature record of these enzymes (Edwards et al., 1986; de Silva et al., 1993; Fanutti et al., 1993, 1996; Chanliaud et al., 2004). NXG1 has also been referred to as nXGase (for nasturtium xyloglucanase) (Chanliaud et al., 2004).

Alignment of 130 full-length putative XTH gene sequences clearly indicated that the archetypal GH16 xyloglucanase, Tm-NXG1, contained unique loop extensions adjacent to the active site that are not present in the strict XETs from cauliflower (Brassica oleracea var botrytis) and hybrid aspen (Populus tremula x Populus tremuloides) (see Supplemental Figures 4 and 5 online). Neither Bob-XET16A nor Ptt-XET16-34 have detectable rates of xyloglucan hydrolysis using a sensitive reducing sugar assay (Henriksson et al., 2003; Kallas et al., 2005), whereas Tm-NXG1 has been demonstrated to be a predominant endo-hydrolase (EC 3.2.1.151) (Chanliaud et al., 2004) that nonetheless has some capacity to catalyze xyloglucan endo-transglycosylation (EC 2.4.1.207) (Fanutti et al., 1996). Therefore, deletion mutants of Tm-NXG1 were produced to test the hypothesis that the loop between ß-strands ß6 and ß7, as well as the loop between ß8 and ß9, might be primary determinants affecting the hydrolysis-to-transglycosylation ratio in this enzyme family. Although recombinant expression of constructs containing deletion variants of the loop sequence joining ß6 and ß7 has been unsuccessful to date, the ß8-ß9 loop variant Tm-NXG1-YNIIG was produced at levels similar to those of Tm-NXG1 and Tm-NXG2 (50 mg/L crude protein) in the methylotrophic yeast Pichia pastoris.

Catalytic Properties of the Recombinant Tm-NXG1, Tm-NXG2, and Tm-NXG1-YNIIG
To resolve any potential ambiguity resulting from previous characterization of native purified enzymes (Edwards et al., 1986; Fanutti et al., 1993) in the context of the single amino acid polymorphism observed between Tm-NXG1 and Tm-NXG2 (de Silva et al., 1993), a detailed analysis of the activity of these enzymes was performed. These data, in turn, provided a basis for an analysis of the effects of loop truncation on the catalytic function of Tm-NXG1-YNIIG.

Substrate Specificity
The activities of recombinantly produced and purified Tm-NXG1 and Tm-NXG2 were tested with a range of substituted and nonsubstituted ß(1-3)- and ß(1-4)-linked polysaccharides. Consistent with data reported >20 years ago (Edwards et al., 1986), these enzymes were specific xyloglucanases that exhibited no detectable activity on konjac glucomannan, Icelandic moss lichenan, barley (Hordeum vulgare) ß-glucan, Avicel, and birchwood (Betula spp) xylan. A trace amount of activity was observed with hydroxyethyl cellulose, as was shown previously for an extracted and purified enzyme preparation (Edwards et al., 1986). The pH rate profile for the hydrolysis of xyloglucan by Tm-NXG1 was classically bell-shaped, with apparent kinetic pK_a values of 4.3 ± 0.1 and 5.3 ± 0.1 (pH optimum, 4.8; see Supplemental Figure 1 online), similar to previous observations (Edwards et al., 1986). Tm-NXG2 exhibited an essentially identical pH rate profile (data not shown). Analysis of the recombinantly produced and purified enzymes yielded specific hydrolytic activities toward xyloglucan of 5.5 ± 0.5 mol XGO·min^–1·mol^–1 protein for Tm-NXG1 and 4.0 ± 0.3 mol XGO·min^–1·mol^–1 protein for Tm-NXG2 at the pH optimum (1 g/L xyloglucan). The hydrolytic activity of the loop deletion mutant Tm-NXG1-YNIIG was significantly lower: 0.23 ± 0.058 mol XGO·min^–1·mol^–1 protein. Neither Tm-NXG1 nor Tm-NXG2 showed detectable transfer of xyloglucan to radiolabeled xylogluco-oligosaccharides under standard assay conditions (Fry et al., 1992; Henriksson et al., 2003; Kallas et al., 2005). Further comparative kinetic analyses of the three enzymes were performed using size-exclusion chromatography (SEC) and high-performance anion-exchange chromatography–pulsed-amperometric detection (HPAEC-PAD) assays.

Mode of Action
The depolymerization of xyloglucan by Tm-NXG1 was followed by SEC to confirm that this enzyme cleaved xyloglucan endolytically and to examine the effect of the addition of Glc₄-based xylogluco-oligosaccharides (XGO_Glc4) on the reaction. The M_r of XG rapidly decreased, with no initial production of xylogluco-oligosaccharides, upon incubation with Tm-NXG1 (Figure 2A ), clearly indicating endo-xyloglucanase activity (EC 3.2.1.151). The addition of XGO_Glc4 resulted in a reduction of depolymerization velocity (Figure 2B), indicating product inhibition of Tm-NXG1 and limited XET activity (EC 2.4.1.207). For comparison, XG was incubated with a true XET, Ptt-XET16-34 from the hybrid aspen, which has undetectable hydrolytic activity (Kallas et al., 2005; this work). In the absence of XGO_Glc4, no change in the M_r distribution of XG was observed after overnight incubation (Figure 2E). Upon addition of XGO_Glc4, an effective depolymerization of XG can be observed due to Ptt-XET16-34 activity (Figure 2F). SEC analysis indicated that Tm-NXG1-YNIIG exhibited a qualitatively slower rate of XG hydrolysis than Tm-NXG1; after overnight incubation, a large amount of XG was not broken down to xylogluco-oligosaccharides (Figure 2C). However, addition of XGO_Glc4 increased the effective rate of XG depolymerization by Tm-NXG1-YNIIG, thus demonstrating improved XET activity relative to Tm-NXG1 (Figure 2D).

View larger version (26K): [in this window] [in a new window] [as a PowerPoint slide]   Figure 2. Size-Exclusion Chromatography of the Products of Tm-NXG1, Tm-NXG1-YNIIG, and Ptt-XET16-34 Operating on Xyloglucan Polysaccharide (Initial Mr 105 to 107) in the Absence and Presence of Xyloglucan Oligosaccharides (Mr 103). Incubation times are indicated above each chromatogram. The vertical dotted line in all panels indicates the elution of Glc4-based XGOs. Protein concentrations were as follows: Tm-NXG1 (9 mg/L), Tm-NXG1-YNIIG (10 mg/L), and Ptt-XET16-34 (10 mg/L). Other experimental conditions are described in Methods.

Initial rate kinetic analysis indicates that Tm-NXG1 is a predominant xyloglucan hydrolase with some capacity to catalyze xyloglucan endo-transglycosylation at increased substrate concentration. A small apparent transglycosylation rate leading to the production of XGO_Glc12 and XGO_Glc4 was observed at XGO_Glc8 concentrations > 0.4 mM, which reached a maximum at 1 mM, as did the hydrolytic rate (Figure 3A ). Notably, at these higher substrate concentrations, the observed ratio of transglycosylation to hydrolysis is nearly constant at 1:12. In stark contrast, Ptt-XET16-34 was shown to be an overwhelmingly predominant transglycosylase over the full concentration range from 10 µM to 3 mM XGO_Glc8, in keeping with previous analyses (Kallas et al., 2005); under initial rate conditions, the stoichiometry of XGO_Glc12 to XGO_Glc4 production was 1:1 (Figure 3B). Remarkably, the deletion mutant Tm-NXG1-YNIIG, which was predicted to produce a more XET-like enzyme relative to the wild type, exhibited increased transglycosylation ability and decreased hydrolytic ability, which become equivalent at high substrate concentrations (Figure 3C).

View larger version (8K): [in this window] [in a new window] [as a PowerPoint slide]   Figure 3. Initial Rate Kinetics as a Function of XGOGlc8 Concentration. Kinetics are shown for Tm-NXG1 (A), Tm-NXG1-YNIIG (B), and Ptt-XET16-34 (C). Closed circles, rate of XGOGlc12 production due to transglycosylation (2 XGOGlc8 XGOGlc12 + XGOGlc4); open squares, total rate of XGOGlc4 production; closed squares, corrected rate of XGOGlc4 production obtained by subtracting the contribution of XGOGlc4 release due to substrate transglycosylation. This observed rate is twice the actual catalytic rate, according to the stoichiometry of the hydrolysis reaction (XGOGlc8 2 XGOGlc4). Error bars indicate SD for triplicate (A), duplicate (B), and duplicate (C) measurements.

The steady state kinetic equations describing the initial rate of formation of products from the two possible reaction pathways have been solved (see Supplemental Text online, Equations A and B); both equations are complex, containing four macroscopic rate constants and only a single independent variable, [XGO_Glc8]. As a consequence, the equations are underdetermined, and meaningful fits to the data could not be obtained. Instead, the hydrolytic and transglycosylation rate data shown in Figure 3 were fit by the classical Michaelis-Menten expression, with or without a substrate inhibition parameter as appropriate (Brumer et al., 1999), to produce the apparent macroscopic kinetic constants shown in Table 1 . Using this simplified treatment, the derived constants provide only a qualitative measure of the relative rates of enzymic transglycosylation and hydrolysis.

View this table: [in this window] [in a new window]   Table 1. Apparent Kinetic Constants of Tm-NXG1, Ptt-XET16-34, and Tm-NXG1-YNIIG on XGOGlc8 Obtained by Nonlinear Least-Squares Curve Fitting

(1)

Therefore, a partitioning constant, K_P, can be defined as a concentration-independent parameter that reflects the ratio of the initial transglycosylation and hydrolysis rates [_{o(transglycosylation)}/_{o(hydrolysis)}]. K_P is infinite for true XETs such as Ptt-XET16-34, implying that either k₃ or k₄ is effectively zero; water either does not bind productively or the complex is catalytically incompetent. By contrast, Tm-NXG1 exhibits a K_P value of 11,500 ± 900, while the loop variant Tm-NXG1-YNIIG has a K_P value of 434,000 ± 2000 (see Supplemental Figures 2 and online). These values clearly highlight the significantly increased rate of transglycosylation to hydrolysis of Tm-NXG1-YNIIG, which results from truncation of the loop connecting ß-strands 8 and 9. Nevertheless, this effect is not simply due to a reduction of the hydrolytic activity; Tm-NXG1-YNIIG shows an apparent 5.7-fold reduction of the hydrolytic rate at saturation compared with Tm-NXG1, as well as a doubling of the transglycosylation rate to one-half that observed for Ptt-XET16-34 (cf. data at [XGO_Glc8] > 1 mM or apparent k_{cat(transglycosylation)} values; Table 1).

Transfer to Cello-Oligosaccharides
Crude nasturtium XET enzyme extracts were recently demonstrated to use cellulose derivatives as alternative donor and acceptor substrates, together with xyloglucan oligosaccharides and polysaccharides (Mohand and Farkas, 2006). This study has confirmed that Tm-NXG1 possesses limited activity toward hydroxyethyl cellulose (Edwards et al., 1986), while Ptt-XET16-34 has been shown to use GGGGXXXG as a glycosyl donor with a fluorescent XXXG-based acceptor (Fauré et al., 2006). The ability to catalyze hetero-transglycosylation reactions, therefore, may be a feature of some GH16 XETs and xyloglucanases. To test the catalytic promiscuity of Tm-NXG1 and Tm-NXG1-YNIIG, these enzymes were incubated in the presence of xyloglucan polysaccharide and a series of cello-oligosaccharide acceptors. The masses of the transfer products for cellobiose (XXXGGG, XLXGGG or XXLGGG, XLLGGG), cellotriose (XXXGGGG, XLXGGGG or XXLGGGG, XLLGGGG), and cellotetrose (XXXGGGGG, XLXGGGGG or XXLGGGGG, XLLGGGGG) were readily observed by electrospray ionization–time of flight–mass spectrometry (ESI-TOF MS; see Supplemental Table 3 online). Higher order xylogluco-oligosaccharide transglycosylation products (based on XGO_Glc8) containing cello-oligosaccharides were also detected (data not shown). No transglycosylation products were detected when lactose (8 mM) was used as an acceptor substrate, thus indicating that glycosyl transfer most likely involved the equitorial C-4 hydroxy group on the nonreducing end of the cello-oligosaccharides. In comparison, Ptt-XET16-34 did not produce any hetero-transglycosylation products with cello-oligosaccharides or lactose, which indicated a much more strict substrate tolerance. The increased ability of Tm-NXG1-YNIIG to catalyze transglycosylation relative to the wild-type enzyme (data not shown) was also reflected in an increased amount of obtained hetero-transglycosylation products produced by this mutant. Under preparative-scale conditions, hetero-transglycosylation was indeed significant; a reaction containing 25 mg of XG and 15 mg of cellobiose yielded 11 mg of pure XGO_Glc4-cellobiose conjugates (XXXGGG, XLXGGG or XXLGGG, XLLGGG) and 6 mg of pure XGO_Glc8-cellobiose conjugates after gel filtration through Bio-Gel P6.

Crystal Structures of Tm-NXG1 and Tm-NXG2
Overall Structure
The crystal structures of Tm-NXG1 and Tm-NXG2 are illustrated in Figure 4A , and full crystallographic data are given in Table 2 . The asymmetric unit of Tm-NXG1 crystals contained 443 water molecules, 4 glycerol molecules, and 3 copies of the protein molecule comprising 274 residues. This includes a three-residue extension (Ala-Tyr-Val), clearly defined by electron density at the N terminus, from the expression vector. Two hundred ninety-one water molecules and two copies of the enzyme were present in the asymmetric unit of Tm-NXG2, both including the 271 residues reported for the primary sequence of the mature protein (UniProt accession number Q07524). One loop was disordered in both molecules of the Tm-NXG2 crystal structure (residues Ala-191 to Tyr-197; Figure 4B). In accordance with the sequence similarity (39% identity; see Supplemental Figure 4 online), the overall structures are strikingly similar to that of their xyloglucan endo-transglycosylase relative, Ptt-XET16-34 (Protein Data Bank [PDB] identifier 1UN1). Both display the ß-jellyroll fold typical for family GH16 enzymes and superimpose onto the structure of Ptt-XET16-34 with a root mean square deviation of 1.35 Å for Tm-NXG1 (241 C atoms) and 1.18 Å for Tm-NXG2 (221 C atoms). Like Ptt-XET16-34, the Tm-NXG structures contain the helical extension at the C terminus stabilized by two disulfide bonds, although the extremity of the C terminus clearly displays a different conformation in the two isoforms of Tm-NXG. Both Tm-NXG enzyme structures lack the glycosylation site reported for most of the known GH16 XTH gene family members (Campbell and Braam, 1999a). Glu-94 was identified as the catalytic nucleophile and Glu-98 as the catalytic acid/base in Tm-NXG1 by three-dimensional structural alignment with Ptt-XET16-34 (Johansson et al., 2004) and by mechanistic analogy with other GH16 enzymes (Planas, 2000).

View larger version (62K): [in this window] [in a new window] [as a PowerPoint slide]   Figure 4. Three-Dimensional Structures of Tm-NXG1 and Tm-NXG2. (A) Ribbon representation of the crystallographic structure of Tm-NXG1. The polypeptide chain is colored from blue (N terminus) to red (C terminus). The three strictly conserved amino acids forming the catalytic machinery are labeled. (B) Superimposition of the structures of Tm-NXG1 (yellow) and Tm-NXG2 (blue) highlighting the different conformations of the loops that surround the active site groove. The catalytic residues, the three Trp residues that line the substrate binding cleft, and Tyr-122 and Ile-124 are labeled. Due to high disorder, the loop residues between Thr-192 and Lys-196 could not be modeled in the crystal structure of Tm-NXG2. (C) Wall-eyed stereo view of the final 2Fo-Fc electron density map displayed at a 1 level around the loop containing Trp-185 and Trp-190 in the crystal structure of Tm-NXG1.

View larger version (42K): [in this window] [in a new window] [as a PowerPoint slide]   Figure 5. Structural Comparison of Ptt-XET16-34 and Tm-NXG1. (A) Ribbon representation of the superimposition of Ptt-XET16-34 (green) onto Tm-NXG1 (yellow) illustrating the structural differences that appear mainly in loops 1 through 3 and in the C-terminal region. The catalytic machinery is labeled. (B) Close-up view into the positive subsites of Ptt-XET16-34. The superimposition was performed using the complex structure of Ptt-XET16-34 with bound XLLG (PDB identifier 1UMZ). Asp-178 in Ptt-XET16-34 forms a hydrogen bond to a xylosyl branch of the bound oligosaccharide. Ser-189 in Tm-NXG1 (or Tm-NXG2) is too distant to be able to form an equivalent bond. Gly-183 is replaced by Asn-194 in Tm-NXG1, which would collide with bound sugar in this loop conformation. Ile-124, which is part of the critical loop insertion in Tm-NXG1, would collide with the glucose unit bound in subsite +1 of the acceptor binding cleft. (C) Ribbon representation of the superimposition of Tm-NXG1-YNIIG (blue) onto Ptt-XET16-34 (green) and Tm-NXG1 (yellow) illustrating that the C trace of the truncated loop 2 in Tm-NXG1-YNIIG is now closer to that of Ptt-XET16-34. The catalytic machinery is labeled.

Crystal Structure of the Loop 2 Truncation Mutant, Tm-NXG1-YNIIG
The asymmetric unit of Tm-NXG1-YNIIG crystals contained 235 water molecules and two copies of the protein molecule comprising 272 residues. This includes a one-residue extension (Val), clearly defined by electron density at the N terminus, arising from the expression construct. The root mean square deviation value for the superimposition of Tm-NXG1 with Tm-NXG1-YNIIG is 0.44 Å for 264 C atoms and 1.0 Å for 225 C atoms with Ptt-XET16-34. As expected, Tm-NXG1-YNIIG closely resembles Tm-NXG1, except for the region around the loop 2 insertion, where the C trace of Tm-NXG1-YNIIG intimately follows that of Ptt-XET16-34 (Figure 5C).

Phylogenetic Analysis of XTH Gene Products
XTH gene products from GH16 have historically been classified into three major phylogenetic subgroups (groups I to III) on the basis of sequence similarity (reviewed in Campbell and Braam, 1999a). The obvious question thus arises whether the current phylogenetic classification can be used to predict biochemical function in the XTH subfamily. A burgeoning amount of genomic data, while providing unprecedented insights into the sequence diversity of XTH gene products, is at the same time blurring the distinction between the traditional subgroups. In particular, comparison of XTH genes in the Arabidopsis thaliana and rice (Oryza sativa) genomes has indicated that there is no longer a clear division between the previously distinct groups I and II (Yokoyama et al., 2004). XTH gene products are small (33 kD) and display a high degree of sequence similarity, so it is likely that as the number of available XTH gene sequences continues to rise, the current subgroup definitions will require revision. Therefore, a new bioinformatic analysis of this GH16 subfamily was performed to provide an updated phylogenetic framework in which to discuss the present three-dimensional structure-function analysis of Tm-NXG1, the archetypal mixed-function xyloglucanase/XET from the historical group III (Edwards et al., 1986; de Silva et al., 1993; Fanutti et al., 1993, 1996; Chanliaud et al., 2004). Approximately 130 full-length protein sequences were used to derive a phylogenetic tree from a structure-based sequence alignment and reconstruction using a maximum likelihood method (Guindon and Gascuel, 2003), including all Arabidopsis, rice, and black cottonwood genomic sequences, all tomato (Solanum lycopersicum) and hybrid aspen EST sequences, and selected individual sequences from ripening fruits and germinating seeds (specified in Methods).

Using a structurally characterized bacterial lichenase as an outlier, the tree forms two predominant branches with significant bootstrap values (Figure 6 ). The largest cluster affirmed previous observations about the merging of groups I and II (Yokoyama et al., 2004). However, this large cluster can, in fact, be subdivided into a number of smaller, statistically significant clades of protein sequences. To avoid further obfuscation, we will nonetheless refer to this ensemble broadly as group I/II. The second main branch encompasses fewer sequences and coincides with the historical group III. This analysis indicates that group III is not monolithic, as commonly represented (Campbell and Braam, 1999a; Yokoyama et al., 2004; Saladié et al., 2006), but can instead be subdivided into two predominant clades. For the purposes of this work, the first clade, which includes Tm-NXG1, will be designated group III-A. Group III-B includes, notably, the Sl-XTH5 and At-XTH27 gene products, which are XETs (Campbell and Braam, 1999b; Saladié et al., 2006). A small outlying group, close to the root and comprising five sequences (At-XTH1 to At-XTH3, At-XTH11, and Pt-XTH15), was also observed. Bootstrap values indicate that these sequences form an intermediate, possibly ancestral, clade between group I/II and group III. At-XTH11 is a notably curious sequence that is likely to encode a catalytically inactive enzyme; the catalytic acid/base and nucleophile-helper residues are mutated to noncarboxylic acid residues (Cys and Gln, respectively), and loop 3 is truncated by nine amino acid residues. Removal of this sequence significantly increases the bootstrap value for this clade.

View larger version (26K): [in this window] [in a new window] [as a PowerPoint slide]   Figure 6. Unrooted Phylogenetic Tree of 130 Full-Length XTH Gene Products and Bacillus licheniformis Lichenase (PDB identifier 1GBG; GenPept Accession Number CAA40547). Bootstrap values from 100 maximum likelihood resamplings are indicated.

	DISCUSSION

TOP ABSTRACT INTRODUCTION RESULTS DISCUSSION METHODS REFERENCES

 
What Determines the Fate of the Glycosyl Enzyme Intermediate in the XTH Subfamily of GH16?
The potential to catalyze xyloglucan endo-transglycosylation (EC 2.4.1.207) or endo-hydrolysis (EC 3.2.1.151) is a key feature of the retaining catalytic mechanism used by plant xyloglucan-active enzymes from GH16 (Sinnott, 1990). Our present structural and biochemical analysis indicates that the extension of loop 2 in Tm-NXG1 and Tm-NXG2 is a determinant of endo-xyloglucanase activity in XTH gene products. The strict transglycosylase Ptt-XET16-34 from hybrid aspen lacks this extension, which the three-dimensional structure of Tm-NXG1 shows is capable of interacting with substrate bound in the positive subsites of the enzyme, potentially modulating the binding of xylogluco-oligosaccharide acceptor substrates. Truncation of this loop to produce a hybrid, Ptt-XET16-34–like enzyme does in fact boost the rate of transglycosylation slightly, while hydrolytic activity is diminished (Figure 3). The variation in the length of loop 2 is clearly an important difference between the XETs and the xyloglucanases in GH16, but our inability to generate a strict transglycosylase from a predominant hydrolase indicates that this is not the only determinant.

Another obvious difference that distinguishes Tm-NXG1 and Tm-NXG2 from Ptt-XET16-34 is the presence of a three-residue extension of loop 1. We have been unsuccessful in expressing Tm-NXG1 variants that possess a truncation of this loop, either alone or in combination with loop 2 truncation, but can nonetheless infer from recent literature data that loop 1 alone is not likely to affect the transglycosylation/hydrolysis ratio. The tomato Sl-XTH5 gene is abundantly expressed during fruit ripening (Saladié et al., 2006) and, like Tm-NXG1, possesses the longer loop 1 variant (see Supplemental Figure 5 online). Viscometric and HPLC-based assays recently indicated that Sl-XTH5 is a predominant XET; under identical assay conditions, Tm-NXG1 exhibited predominant or exclusive endo-xyloglucanase activity (Saladié et al., 2006). Similarly, early work by Campbell and Braam (1999b) demonstrated that recombinant Arabidopsis At-XTH27 (EXGT-A3) was likewise a XET, devoid of hydrolytic activity. Unfortunately, the assays employed in these studies do not allow quantitative comparison of the transglycosylation and hydrolysis rates of these enzymes in terms of turnover numbers; nonetheless, they show qualitatively that the longer loop 1 observed in Sl-XTH5 and At-XTH27 (see Supplemental Figure 5 online) does not confer significant hydrolytic activities to these enzymes. Instead, our evidence suggests that the extension of loop 2 in Tm-NXG1–like enzymes is a major structural change responsible for the XET-to-xyloglucanase switch in GH16, as indicated by our mutagenesis and quantitative kinetic analyses. It is apparent, however, that a number of smaller, currently unexplored sequence variations (see Supplemental Figures 4 and 5 online) also conspire to dictate the fate of the glycosyl enzyme intermediate (Figure 1). Thus, further structure–function studies are clearly needed to unlock all of the secrets necessary to fully convert a xyloglucanase into a XET or vice versa.

The Evolution of Divergent Xyloglucan Activities in GH16
In view of the entire GH16 family, these phylogenetic analyses raise a fundamental evolutionary question. Because the ancestral enzyme of GH16 is likely to be a hydrolytic laminarinase (Barbeyron et al., 1998; Michel et al., 2001), is the ancestor of the XTH gene family a hydrolytic or transglycosylating enzyme? Our phylogenetic tree, based on maximum likelihood reconstructions, indicates that both main branches contain predominant XETs, while predominant xyloglucan hydrolases only appear in group III-A. Closer examination indicates that the putative ancestral group, which contains only sequences with XET-like loops, is likely to have emerged earlier than the Tm-NXG1–containing group III-A. Furthermore, the observation that group III-A is not very divergent (between 60 and 70% pairwise sequence identity) provides additional support for the hypothesis that the emergence of Tm-NXG–like enzymes is a relatively recent event. Moreover, Tm-NXG–like enzymes are rare, and their presence appears to be limited to germinating seeds, ripening fruits, and some fast-growing tissues (vide infra). This rareness is also exemplified, for instance, by a census of the Arabidopsis XTH family: from >30 homologous enzymes, only 2 (At-XTH31 and At-XTH32) are predicted to have xyloglucanase activity by our analysis (Figure 6).

Perhaps the strongest argument for the evolution of GH16 xyloglucanases from XETs involves a common structural characteristic that distinguishes the XTH gene family from all other GH16 enzymes. As illuminated by the first XET three-dimensional structure (Johansson et al., 2004), these enzymes have diverged from their hydrolytic cousins, the lichenases and laminarinases, by extension of the C terminus, thereby increasing the protein–substrate interactions in the positive subsites. Like the true XETs, Tm-NXG1 also has the C-terminal loop extension, which provides an extra ß-strand extending the length of the substrate binding cleft. Therefore, it is highly unlikely that Tm-NXG1–like enzymes have emerged directly from a hydrolytic ancestor that lacks this C-terminal extension. Indeed, here, we show that in Tm-NXG1, the additional hydrolytic activity is principally due to subtle changes in a surface loop that interacts with the positive subsites, probably affecting acceptor substrate binding. Altogether, these arguments favor the conclusion that xyloglucan endo-transglycosylation is the ancestral activity, which predates xyloglucanase activity arising from a secondary divergence in (historical) group III.

Biological Implications
Analysis of the tissue-specific expression of a range of GH16 XTH genes has been performed for a variety of species, including Arabidopsis, Populus, and rice, whose complete genomes have been sequenced (Yokoyama and Nishitani, 2001; Yokoyama et al., 2004; Aspeborg et al., 2005; Geisler-Lee et al., 2006). Recent comprehensive work by Braam and coworkers using XTH:GUS gene fusions has further illuminated the temporal and spatial expression of these genes in Arabidopsis (Becnel et al., 2006). With regard to the biological role of group III-A enzymes, which may represent a specific subgroup of predominant xyloglucanases in GH16, some interesting patterns emerge. The tissue-specific expression of two nasturtium XTH genes, Tm-NXG1 (group III-A) and Tm-XET1 (group I/II; GenBank accession number L43094), has been studied previously (Rose et al., 1996). XET1, which is 88% identical and 94% similar to the nonhydrolytic Ptt-XET16-34 (Kallas et al., 2005) at the protein level, exhibits a systemic expression pattern, with the exception of germinating seed cotyledons. Tm-NXG1, on the other hand, is expressed exclusively in this tissue, from which the active enzyme has been isolated (Edwards et al., 1986). Furthermore, Tm-NXG1 homologs are expressed specifically in a range of ripening fruits, including Litchi (Lu et al., 2006), papaya (Carica papaya) (GenBank accession number AY032600; R. Haireen, Z.M. Ali, R. Othman, and H. Lazan, unpublished data), grape (Vitis vinifera) (Ishimaru and Kobayashi, 2002), and apricot (Prunus armeniaca) (GenBank accession number AY578086; F. Geuna, R. Banfi, and D. Bassi, unpublished data). These observations are consistent with a general role of Tm-NXG–like enzymes in the degradation of xyloglucan during seed germination, on the one hand, and fruit ripening, on the other.

Two Arabidopsis XTH gene products are also found in group III-A: At-XTH31 (65% identical and 85% similar protein sequence to Tm-NXG1) is predicted to have a role in seed germination (Aubert and Herzog, 1996) and root morphogenesis (Aubert and Herzog, 1996; Becnel et al., 2006), whereas At-XTH32 (68% identical and 81% similar to Tm-NXG1) is expressed primarily in the shoot apex (Becnel et al., 2006). The observation of the expression of these genes in expanding tissues is not inconsistent with a potential role in xyloglucan degradation allowing rapid wall extension. Xyloglucanase activity is not likely to be a ubiquitous requirement for wall degradation in processes such as fruit ripening, however (Brummell, 2006); the Tm-NXG1 homolog Sl-XTH6 (group III-A) was not observed in EST libraries from ripening tomato, whereas the group III-B XET encoded by Sl-XTH5 was implicated in fruit development (Saladié et al., 2006).

As a notable side activity, it is currently unclear whether the ability of Tm-NXG1 to perform hetero-transglycosylation to cello-oligosaccharide acceptors has biological relevance. Indeed, it is likely that the in vitro cross-conjugation of xyloglucan and linear ß-glucans observed in crude nasturtium extracts (Mohand and Farkas, 2006) is due, at least in part, to Tm-NXG1, based on our results. However, in the germinating nasturtium seed, the effective physiological concentrations of internal XG bonds and water are likely to be much higher than the concentration of nonreducing ends of other ß-glucans, thus implying that the effective rate of hetero-transglycosylation would be low. Similarly, a barley XET isozyme was recently shown to have very low hetero-transglycosylation activity using alternating XG:ß-glucan and ß-glucan:XG donor/acceptor pairs in vitro (Hrmova et al., 2007). Although there is currently no evidence for the existence of analogous block copolysaccharide structures in muro (Hrmova et al., 2007), the predominance of (1,3:1,4)-ß-glucan over XG in barley primary cell wall raises the possibility that hetero-transglycosylation is kinetically feasible in that system. However, it should be borne in mind that promiscuous glycosyl transfer by retaining glycosidases to nonnatural acceptors is well documented, including for endo-(xylo)glucanases (York and Hawkins, 2000).

In this context, it is interesting that Tm-NXG1 is a comparatively feeble xyloglucanase, with a v_o/[E]_t value of 5.5 ± 0.5 min^–1 at the pH optimum and 1 g/L tamarind xyloglucan. Under similar conditions and an identical substrate concentration, a strict xyloglucanase from Paenibacillus pabuli (GH5, retaining) has a v_o/[E]_t value of 8700 ± 220 min^–1, while a nonspecific endo-glucanase from Bacillus licheniformis (GH12, retaining) has a v_o/[E]_t value of 54 ± 3 min^–1 (Gloster et al., 2007). Similarly, the inverting microbial xyloglucanases from GH74 have catalytic rates of at least 5 to 10 times greater than that of Tm-NXG1 (Irwin et al., 2003; Martinez-Fleites et al., 2006). Therefore, the poor hydrolytic activity of Tm-NXG1 may reflect the different physiological requirements for xyloglucan degradation by plants and microorganisms. During germination, it may be anticipated that the developing embryo requires metered, continuous release of sugar in the confines of the seed, whereas in rapidly expanding tissues, xyloglucan degradation must not be so rapid as to lead to irreparable wall loosening. On the other hand, the ability to rapidly degrade available carbohydrate sources would confer a selection advantage to microbes in more competitive environments. Tm-NXG1 nonetheless exhibits a rate of transglycosylation that is only fivefold lower than that of the strict transglycosylase Ptt-XET16-34 over a wide concentration range (Figure 3, Table 2). This supports bioinformatic analyses that indicate that mixed-function xyloglucanases/XETs may have evolved from true XETs; xyloglucanase activity has resulted from a gain of function in GH16, where glycosyl transfer to water has been favored by extension of specific active site loops. It is apparent that additional comprehensive biochemical studies are required to further illuminate the catalytic diversity hidden within this enzyme family.

	METHODS

TOP ABSTRACT INTRODUCTION RESULTS DISCUSSION METHODS REFERENCES

 
Material Sources
Ultrapure water was used in all experiments and refers to water purified on a Milli-Q system (Millipore) with a resistivity of > 18.2 M·cm. Polysaccharides were obtained from the following sources: konjac glucomannan, Icelandic moss lichenan, tamarind xyloglucan, and barley (Hordeum vulgare) ß-glucan (all from Megazyme); Avicel and hydroxyethyl cellulose (both from Fluka); and birchwood xylan (from Sigma-Aldrich). Xylogluco-oligosaccharides based on Glc₄, Glc₈, and Glc₁₂ backbones (XGO_Glc4, XGO_Glc8, and XGO_Glc12, respectively) were produced as reported by Martinez-Fleites et al. (2006). These samples represent mixtures of variably galactosylated oligosaccharides arising from the natural distribution of XXXG, XLXG, XXLG, and XLLG in tamarind kernel xyloglucan. The XGO_Glc8 preparation was contaminated with 0.075% XGO and 0.5% XGO_Glc12 (HPAEC-PAD analysis; vide infra); in all calculations of enzyme activities, these background values were subtracted. The average molecular weight of XGO_Glc8 (M_r 2530) was calculated from the composition of oligosaccharides XXXG, XLXG, XXLG, and XLLG (15:7:32:46) in tamarind xyloglucan.

Protein Analytical Methods
Proteins were analyzed by SDS-PAGE using an Xcell II Mini Cell (Novex). Precast NuPAGE gradient gels (4 to 12%) were purchased from Invitrogen. Molecular mass standards (LMW calibration kit) were obtained from Amersham Biosciences. Protein gels were stained by Coomassie Brilliant Blue R 250. Protein concentrations were determined spectrophotometrically by the Bradford method using BSA as a standard. DNA sequencing was performed by the Sanger method using an ABI Prism 377 DNA sequencer (Perkin-Elmer). The pH rate profiles of Tm-NXG1 and Tm-NXG2 were determined in duplicate using the method described by Nelson (1944) with XG (1 g/L) as substrate. The following 100 mM buffer systems were used: sodium acetate, pH 4.0 to 5.75; sodium phosphate, pH 6.0 and pH 6.5. The unfolding temperature of proteins was determined by the fluorescence-based thermal shift assay of Pantoliano et al. (2001), as described by Kallas et al. (2005). ESI-MS analysis of proteins was performed as described by Sundqvist et al. (2007).

Carbohydrate Analytical Methods
XET and xyloglucanase activities were measured with the colorimetric assay according to Sulova et al. (1995) and the radioactive assay devised by Fry et al. (1992). Activities on soluble polysaccharides other than xyloglucan were determined at concentrations of 0.5 g/L and an enzyme concentration of 0.02 g/L at 30°C. Activity on Avicel was measured on a 0.5 g/L suspension. Aliquots were taken at 10, 30, 60, and 130 min, and reducing sugars were quantified using the method of Nelson (1944). A six-point standard curve was made with glucose in the range 20 to 150 mg/L.

HPAEC-PAD analysis of xylogluco-oligosaccharides was performed using a Dionex Carbopac PA100 column on a Dionex ICS-3000 HPLC system controlled by Chromelion software, version 6.80 (Dionex). An optimized gradient program was used for simultaneous analysis of all xylogluco-oligosaccharides (XGO_Glc4, Glc₄-based xylogluco-oligosaccharides; XGO_Glc8, Glc₈-based xylogluco-oligosaccharides; and XGO_Glc12, Glc₁₂-based xylogluco-oligosaccharides). Solvent A was water; solvent B was 1 M sodium hydroxide; solvent C was 1 M sodium acetate. The gradient program was as follows: 0 to 3 min, 100 mM sodium hydroxide and 40 mM sodium acetate; 3 to 18 min, 100 mM sodium hydroxide and a linear gradient of 40 to 300 mM sodium acetate; 18 to 19 min, 500 mM sodium hydroxide and 500 mM sodium acetate. The system was then reequilibrated for 4 min with the initial conditions prior to the next injection. For kinetic experiments, standard curves were made for XGO (0 to 50 µM, 13 points) and the transglycosylation product XGO_Glc12 (Glc₁₂-based xylogluco-oligosaccharides; 0 to 50 µM, 13 points). Chromelion (version 6.80; Dionex) HPLC software was used to fit a cubic function to data from samples analyzed in triplicate.

Time Dependence of the Tm-NXG–Catalyzed Depolymerization of Xyloglucan
Samples (250 µL total volume) containing XG (1 g/L) and Tm-NXG1 (9 mg/L), Ptt-XET16-34 (10 mg/L), or Tm-NXG1-YNIIG (10 mg/L) in sodium acetate buffer (25 mM, pH 5.0) were incubated without or with XGO (0.35 g/L, 0.0035 g/L in the case of Tm-NXG1-YNIIG) at 30°C for 30, 60, or 150 min and overnight. The reactions were immediately frozen, lyophilized, and redissolved in 250 µL of DMSO prior to SEC analysis as described by Brumer et al. (2004). Control samples were treated under identical conditions but without enzyme.

Initial Rate Kinetics with XGO_Glc8 Xylogluco-Oligosaccharides
For each substrate concentration (10 to 3000 µM, 19 points), 20 µL of XGO_Glc8 dilutions was prewarmed for 5 min at 30°C. Enzyme solution in 25 mM sodium acetate buffer, pH 5.0 (5 µL), was added and the reaction was incubated at 30°C; after 7 min, 10 µL of sodium hydroxide (2.5 M) was added to stop the reaction. Enzyme activity was adjusted so that the substrate (XGO_Glc8) turnover was in all cases <10% of the initial concentration and typically <1%. All samples were made in duplicate. Samples were kept at 4°C for not longer than 6 h until analysis by HPAEC-PAD.

Transfer of Xyloglucan to Cello-Oligosaccharide Acceptor Substrates
For analytical reactions, Tm-NXG1 (0.02 g/L) was incubated for 18 h with XG (1 g/L) and lactose (8 mM), cellobiose (8 mM), cellotriose (8 mM), or cellotetraose (8 mM) in 16 mM sodium acetate buffer, pH 5.0. Samples (25 µL total volume) were analyzed by HPAEC-PAD and ESI-TOF-MS (Brumer et al., 2004). For preparative-scale synthesis of transfer products, Tm-NXG1-YNIIG (0.01 g/L) was incubated with cellobiose (3.6 M) and XG (2 g/L) in 16 mM sodium acetate buffer, pH 5.0 (12 mL), for 24 h. Products were purified by the gel-filtration method described by Martinez-Fleites et al. (2006) and identified by ESI-TOF-MS.

Bioinformatic Analysis
Bioinformatic analysis was performed on 130 full-length protein sequences that reflected the diversity of XETs and xyloglucanases in GH16. These included all genome-derived sequences from Arabidopsis thaliana (33 sequences) and Oryza sativa (29 sequences) and individual sequences from Carica papaya (GenPept accession number AAK51119; R. Haireen, Z.M. Ali, R. Othman, and H. Lazan, unpublished data), Tropaeolum majus (Tm-NXG1) (de Silva et al., 1993), Litchi chinensis (Lc-XET1) (Lu et al., 2006), and Vitis labrusca x Vitis vinifera (Vv-XET1) (Ishimaru and Kobayashi, 2002), selected from curated Carbohydrate-Active Enzyme database entries (http://www.cazy.org). Full-length Solanum lycopersicum ESTs (16 sequences) were selected from those reported by Saladié et al. (2006). Thirteen full-length Populus tremula x Populus tremuloides ESTs (Nishiku