Identification and functional characterization of a flax UDP-glycosyltransferase glucosylating secoisolariciresinol (SECO) into secoisolariciresinol monoglucoside (SMG) and diglucoside (SDG)
© Ghose et al.; licensee BioMed Central Ltd. 2014
Received: 4 August 2013
Accepted: 19 March 2014
Published: 28 March 2014
Lignans are a class of diphenolic nonsteroidal phytoestrogens often found glycosylated in planta. Flax seeds are a rich source of secoisolariciresinol diglucoside (SDG) lignans. Glycosylation is a process by which a glycosyl group is covalently attached to an aglycone substrate and is catalyzed by uridine diphosphate glycosyltransferases (UGTs). Until now, very little information was available on UGT genes that may play a role in flax SDG biosynthesis. Here we report on the identification, structural and functional characterization of 5 putative UGTs potentially involved in secoisolariciresinol (SECO) glucosylation in flax.
Five UGT genes belonging to the glycosyltransferases’ family 1 (EC 2.4.x.y) were cloned and characterized. They fall under four UGT families corresponding to five sub-families referred to as UGT74S1, UGT74T1, UGT89B3, UGT94H1, UGT712B1 that all display the characteristic plant secondary product glycosyltransferase (PSPG) conserved motif. However, diversity was observed within this 44 amino acid sequence, especially in the two peptide sequences WAPQV and HCGWNS known to play a key role in the recognition and binding of diverse aglycone substrates and in the sugar donor specificity. In developing flax seeds, UGT74S1 and UGT94H1 showed a coordinated gene expression with that of pinoresinol-lariciresinol reductase (PLR) and their gene expression patterns correlated with SDG biosynthesis. Enzyme assays of the five heterologously expressed UGTs identified UGT74S1 as the only one using SECO as substrate, forming SECO monoglucoside (SMG) and then SDG in a sequential manner.
We have cloned and characterized five flax UGTs and provided evidence that UGT74S1 uses SECO as substrate to form SDG in vitro. This study allowed us to propose a model for the missing step in SDG lignan biosynthesis.
KeywordsFlax Lignan UGTs SDG Secoisolariciresinol Glucosylation Glycosyltranferases
Lignans are a class of diphenolic nonsteroidal phytoestrogens with a wide variety of purported health benefits [1–4]. Different types of lignans have been reported in various plant species and include secoisolariciresinol diglucoside (SDG) found mainly in flax (Linum usitatissimum L.) [5–10]. Flax seeds are a rich source of SDG lignans that have been associated with positive roles in the prevention of chronic metabolic diseases in human [11–14].
In planta, lignans are usually found glycosylated in oligomeric chains . Glycosylation is a key mechanism that determines the chemical complexity and diversity of plant natural products [16, 17], ensures their chemical stability and water solubility while reducing chemical reactivity or toxicity , and facilitates their sorting, intercellular transport, storage and accumulation in plant cells [19, 20]. Glycosylation is one of the key modifications in the secondary metabolite biosynthesis and is catalyzed by carbohydrate active enzymes (CAZymes) which include the superfamily of glycosyltransferases (GTs) . The specific glycosylation position in biologically active compounds may serve to modulate their pharmacological activity and/or to alter and optimize their potential use as drugs . Members of the GT superfamily have been classified into 94 families where family 1 refers to the uridine glycosyl transferases (UGTs) [22, 23]. Plant UGTs are often characterized by a 44 amino acid consensus signature motif, the plant secondary product glycosyltransferase (PSPG) box [15, 23, 24]. UGTs transfer UDP-activated sugar moieties, including UDP-glucose, to specific acceptor molecules . Based on sequence homology, more than 120 UGTs have been reported in Arabidopsis and were grouped into 30 sub-families classified as UGT71 to UGT100 . In the course of this study, the flax draft genome was released . Barvkar et al.  probed this flax draft genome and reported 137 flax UGTs but did not assign functions to any of these UGTs.
To gain insights into SDG lignan glucosylation with potential applications in lignan metabolism engineering, we attempted to identify and characterize flax UGTs responsible for SECO glucosylation. Using database mining, molecular cloning, heterologous expression and enzyme assays, we isolated five putative UDP-glycosyltransferases from flax seeds and demonstrated that UGT74S1 glucosylated SECO, forming sequentially SECO monoglucoside (SMG) and then SDG. The findings, not only reported the first functional characterization of a SECO specific UGT in flax, but also pave the way for engineered SDG lignan metabolite species in vitro and in planta.
Library mining and UGT cloning
Using 19 NAPGEN EST library-derived gene-specific UGT primers and one degenerate (UGT-F2) primer, a total of 16 combinations produced unique PCR products of the expected sizes. The partial cDNA sequences were analyzed using BLASTx which confirmed the identity of each sequence as belonging to the UGT family. A ClustalW multiple sequence alignment showed that some of them were the same and a consensus phylogenetic tree revealed that eight were unique (Additional file 1). Subsequently, one representative sequence from each of the eight UGTs was selected for the design of gene-specific primers, and full length cDNAs for five different UGTs were obtained (Additional file 2A-C). CL5227 was 1.2 kb while CL809, CL8584, RP131, and RP250 were all ~1.5 kb (Additional file 2C). The unique UGT sequences were classified as belonging to four families and five sub-families as per the nomenclature of the International Union of Biochemistry and Molecular Biology and the IUPAC-IUBMB joint committee responsible for UDP-glycosyltransferases  and designated UGT74S1 (CL809), UGT94H1 (CL5227), UGT89B3 (CL8584), UGT74T1 (RP131) and UGT712B1 (RP250). Their sequences were submitted to GenBank under accession numbers JX011632 to JX011636.
UGT structural gene organization
PSPG motif characterization
Tissue-specific in silicoEST analysis of UGTs
A BLASTn search against the flax EST database that includes libraries from 13 different tissues revealed a higher level of expression in embryo and seed coat (Additional file 4). UGT712B1 expression was exclusively detected in the globular and heart stage embryos (GE and HE) whereas UGT94H1 was expressed in the torpedo (TE) and cotyledon stage embryos (CE), as well as in the torpedo stage seed coat (TC) (Additional file 4). UGT74S1, UGT74T1 and UGT89B3 were found exclusively in globular (GC) and torpedo stage seed coat (TC). UGT74S1 and UGT74T1 were the most abundant with 25 EST hits each in the TC EST library.
Quantitative expression of UGTs and PLR in developing flax seed, leaf and stem tissues
SDG lignan profiling
Heterologous expression of flax UGTs and enzyme activities
Enzyme assays and reactions conditions
UGT74S1 biochemical parameters
Because UGT94H1, UGT89B3, UGT74T1 and UGT712B1 did not glycosylate SECO into SMG, further tests were conducted to determine if they were involved in the glucosylation of SMG to form SDG. Since SMG is not commercially available, SDG was hydrolyzed to SMG . Using this SMG as a substrate, the five UGTs were assayed. But again, only UGT74S1 showed a peak corresponding to SDG retention time (data not shown). Therefore, UGT89B3, UGT74T1, UGT712B1, and UGT94H1 appeared not to be involved in SDG lignan glycosylation and their biochemical function remains to be elucidated. Thus, UGT74S1 was the only flax UGT cloned and identified in this study that used SECO as a substrate, first producing SMG and then SDG in a sequential manner. Its truncated version TrUGT74S1 was also assayed using the optimal conditions set for UGT74S1 and was also unable to glucosylate SECO (Additional file 6).
UGT74S1 kinetic parameters
By reacting UGT74S1 with SECO at pH 8.0 and 30°C, the catalytic efficiency (kcat) for SDG production was determined to be 0.89 sec−1. The estimated apparent Km values toward SECO and UDP-glucose for SDG production were determined to be 79 and 1188 μM, respectively.
UGTs are a large and complex family of enzymes that catalyze glycosidic bond formation. To get a better understanding of UGTs that may play a role in the glycosylation process of flax SDG lignan, we undertook the cloning and characterization of flax UGTs. We identified and characterized five flax full length UGTs, namely UGT74S1, UGT94H1, UGT89B3, UGT74T1, and UGT712B1. We found that UGT74S1 and UGT94H1 were highly expressed in developing seed and their expression was coordinated with that of PLR, the first-step lignan biosynthetic gene , and well correlated with the SDG lignan biosynthesis patterns in seed. By expressing each of the five UGTs and reacting the purified proteins with SECO and UDP-glucose, only UGT74S1 produced both SMG and SDG metabolites. To our knowledge, this is the first demonstration linking any flax UGT gene to SDG lignan biosynthesis.
The International Union of Biochemistry and Molecular Biology and IUPAC-IUBMB joint committee responsible for UDP-glycosyltransferase  classified the five UGTs into four families and five sub-families, representing five distinct genes. In the course of this study, Barvkar et al.  probed the recently released flax genome (; Deyholos, http://www.linum.ca) and reported 137 flax UGTs including homologs to our reported UGT74S1 (CL809), UGT94H1 (CL5227), UGT89B3 (CL8584) and UGT712B1 (RP250). These were not, however, characterized with regards to their functionality towards aglycones. Moreover, TrUGT74S1 (JN088324.1; ) was 50 amino acids shorter than UGT74S1 described herein (Additional file 5). We provided convincing evidence that TrUGT74S1 is unable to glucosylate SECO into SDG, and is thereby not functional (Additional file 6). The 50 amino acids missing in TrUGT74S1 seem to be essential for glucosyltransferase activity.
The UGTs described in this study differed in their structural organization, primary sequence, and in their PSPG motifs. Coding sequence variation among plant UGT family 1 members is generally high, varying from less than 35% to more than 95% overall identity , with the C-terminal regions that contain the PSPG box being more conserved . Although well conserved, diversity within the PSPG motif of the five flax UGT genes was revealed. At the structural level, one of the UGTs had no introns while the remaining four had one intron each, which varied in size from 71 to 739 bp. In Arabidopsis, more than half of the UGTs have no introns  and those with introns were much smaller (~100 bp), a difference somewhat proportional to the genome size differences of ~370 Mb for flax and 135 Mb for Arabidopsis. Differences were also observed in the spliced coding sequence (CDS) sizes (379 to 476 amino acids), further emphasizing the diversity within the UGT family and in agreement with its recent origin hypothesis [22, 23].
Although UGT family 1 is a very diverse gene superfamily, its members are usually classified based on their sequence identity  and the presence of the conserved PSPG motif  that includes key conserved residues for substrate recognition and catalysis . The UGTs described herein all possessed the conserved 44 amino acid PSPG motif and the two peptide sequences, WAPQV and HCGWNS, present in 95% of all β-group GTs analysed to date . Amino acid variations were nonetheless observed (italized positions) in these two short peptide motifs as well as in the C-terminal of the PSPG-box . Sugar donor specificity has been attributed to the PSPG box . For example, substitution of tryptophan (W) at position 355 (position 22 of PSPG) for arginine (R) sufficed to modify the sugar donor specificity from UPD-glucose to UDP galacturonic acid in Lamiale. The domain involved in the recognition and binding of the diverse aglycone substrates is purported to be located towards the N-terminal end, whereas the C-terminal region encodes a domain involved in binding the nucleotide sugar substrate .
The optimal enzyme conditions (pH, temperature, cofactors) for UGT74S1 were established and fall within the range of the majority of UDP-glycosyltransferases . UGT74S1 was found to be sensitive to increased ionic strength of metal ions as reported for other UDP-glycosyltransferases [39, 40]. The UGT74S1 apparent Km for UDP–glucose was higher than that for SECO, and fall in the Km ranges previously reported [41–43]. The catalytic efficiency (kcat) of UGT74S1 for SECO was close to that of UGT71A9 and UGT94D1 reported by Noguchi et al. . None of the flavonoid or phenolic acid aglycone substrates tested in this study served as good substrates for UGT74S1.
Taken together, we have cloned five UGTs from flax seeds and demonstrated through a comprehensive multi-approach analysis that UGT74S1 was a functional enzyme capable of converting SECO into SDG. Our results suggest that UGT74S1 is involved in secoisolariciresinol glucosylation in planta to form flax SDG lignan. The findings shed more clarity in flax lignan biochemistry and provide the necessary background to conduct site directed mutagenesis studies.
Flax plants (Linum usitatissimum L. cv AC McDuff) were grown at AAFC Harrington farm (Harrington, PEI, Canada) in the 2008 to 2011 growing seasons. Plants were grown in four replications each year. At anthesis, referred to as 0 days after anthesis (0 DAA), individual flowers were tagged. Developing bolls were harvested at 0, 8, 16, 24, 32 DAA and at maturity and immediately frozen in liquid nitrogen as previously described in Arabidopsis , soybean  and flax [46, 47]. The 0 DAA samples consisted of ovaries free of other flower tissues, whereas the other boll samples (8–32 DAA and maturity) contained seeds at different developmental stages (Additional file 7). At the flowering stage, young leaf and stem tissues were similarly collected. Developing bolls, leaves and stems were stored at −80°C until use.
Before RNA isolation, ovules (0 DAA) and developing seeds were first extracted from the bolls. Total RNA was isolated using Trizol (Invitrogen, Carlsbad, ON, Canada) as previously described . RNA samples were further purified using the Invitrogen PureLink™ RNA Mini kit (Invitrogen, Mississauga, ON, Canada) as per manufacturer’s instructions, quantified by spectrophotometry, and the quality was verified by agarose gel eletrophoresis and the Experion RNA analyzer (BioRad, Missisauga, ON, Canada).
Library mining and UGT cloning
The flax NAPGEN EST database (Plant Biotechnology Institute, NRC, Saskatoon) was mined using the keywords UGT, glucosyltranferase and glycosyltranferase. A total of 893 UGT hits were found amongst 178,656 ESTs. For primer design, we retained members of UGT subclasses 71 (7 hits), 88 (3 hits) and miscellaneous (7 hits). A set of 19 flax-specific and one degenerated primer pairs were designed (Additional file 8).
Total RNA (2 μg) from all developmental stages was used as template to create the cDNA using the first strand cDNA synthesis kit (Invitrogen, Mississauga, ON, Canada) following manufacturer’s instructions. After treatment with 2 U RNAse H (Invitrogen), the cDNA samples were diluted 10-fold and 1 μL was used as template. Each of the 20 primer pairs (Additional file 8) was used in PCR reactions consisting of an initial denaturation 94°C for 2 min followed by 35 cycles of 94°C for 30 s, 60-63°C for 30 s, and 72°C for 60 s prior to a final extension at 72°C for 10 min. Aliquots of 10 μL of the PCR products were resolved on 1% agarose gels stained with ethidium bromide. The amplified fragments were purified using the QIAquick gel extraction kit (Qiagen) for direct sequencing and for TOPO cloning (Invitrogen) in E. coli prior to sequencing.
The identities of the obtained partial sequences were confirmed by BLASTx against the NCBI non-redundant protein sequence (nr) database using a cut off value of 1e−30. The relationship between the partial sequences was inferred by a phylogenetic consensus tree constructed using UPGMA method with 1000 bootstrap replicates as implemented in MEGA4 .
To clone the full length UGTs, 5′ and 3′ gene specific primers (GSP) and nested gene specific RACE PCR primers were designed from representative sequences of each group observed in the consensus tree (Additional file 1) and were used in 5′ and 3′ cDNA end amplification reactions. Briefly, using the Gene Racer kit (Invitrogen, Mississauga, ON, Canada), the purified total RNA was dephosphorylated using a calf intestinal phosphatase, and decapped with a tobacco acid pyrophosphatase. The RNA oligos were ligated to the decapped mRNA by T4 RNA ligase (Invitrogen, CA, USA) before reverse transcription of mRNA using oligo-dT primers. The 5′ and 3′ RACE PCR reactions were carried out using eight pairs of GSP and nested primers (Additional file 9) following the kit’s specifications. The expected 5′ and 3′ RACE PCR products of the putative UGTs CL809, CL5227, CL8584, RP131, RP250 were gel-purified, cloned in TOPO 4.0 vector (Invitrogen, Mississauga, ON, Canada) and sequenced using M13 forward and reverse primers. New primer sets containing restriction sites compatible with the multiple cloning site of pYES2/NT C plasmid vector (Invitrogen, Mississauga, ON, Canada) were designed from the 5′ and 3′ ends (Additional file 10) for the amplification of the full length cDNAs (Additional file 2). The amplified full length cDNAs were gel-purified, restriction digested and similarly cloned into pYES2/NT C. The cDNA corresponding to one of the UGT clones reported by Bavkar et al.  (accession JN088324.1) was also cloned as described above. The plasmids carrying the full length cDNA clones were sequenced using T7 promoter primer x(5′-TAATACGACTCACTATAGGG-3′) and CYC1 reverse primer (5′-GCGTGAATGTAAGCGTGAC-3′).
UGT structural gene organization
To characterize the structural organization of the flax genomic DNA corresponding to each of the five UGTs, a BLASTn search within the flax sequence assembly (http://www.linum.ca) was performed to identify the 5′ and 3′ untranslated regions (UTRs), and the intron and exon structure of the coding regions. The PROSITE scan tool of the ExPASy web interface (http://www.expasy.org) was used to determine the position of the conserved motifs characteristic of plant UGTs such as the PSPG box.
In silicoanalysis of UGTs
To characterize the relative abundance of the cloned UGTs, an in silico EST analysis was performed. The five full length UGT sequences were compared to 13 flax tissue-specific EST libraries (globular embryo, heart embryo, torpedo embryo, cotyledon embryo, mature embryo, pooled endosperm, globular stage seed coat, torpedo stage seed coat, etiolated seedling, leaves, stem, stem peel and mature flower) previously described  and the number of EST hits corresponding to each query UGT in each library was recorded and plotted.
UGT real time gene expression analysis
To assess the gene transcript expression levels of the putative cloned UGTs in developing flax seed, leaf and stem tissue, real-time PCR primers were designed from the five flax UGTs, one PLR and one ribosomal (EU307117) RNA sequence (Additional file 11). The rRNA primers were used for data normalization. Total RNA was extracted from three separate biological replicates for each seed developmental stage (0, 8, 16, 24, 32 DAA, and mature seed). First strand cDNA was obtained as described earlier. The cDNA samples were quantified by spectrophotometry or Qubit (Invitrogen) and diluted to 100 ng/μL. Real-time PCR reactions were performed using the SYBR Green PCR Master Mix (BioRad Laboratories, Canada) on a CFX96 Real Time system (BioRad). For each sample, three biological and three technical replicates, for a total of 9 data points, were obtained. The 25 μL Real Time amplification reactions consisted of 1x SYBR Green Master Mix, 300 nM of each primer, 100 ng of first strand synthesis cDNA obtained from ovaries (0 DAA), developing seeds (8, 16, 24, 32 DAA), mature seeds, leaves, stems and water controls. Real-time PCR reactions were performed as follows: denaturation at 95°C for 10 min followed by 40 cycles of 95°C for 30 s, 60°C for 30 s. Following the final amplification cycle, a melting dissociation curve was generated to ensure specificity of the primers and to confirm the uniqueness of the amplification product. The output data was determined following the 2-∆∆CT method described by Livak and Schmittgen  and it is reported as fold changes of relative expression.
SDG lignan profiling in developing flax seeds
To assess the SDG lignan biosynthesis in developing flax seeds, 250 mg of ovary or seed at six developmental stages was used as starting material following modifications to a protocol described by Popova et al. . Developing flax seed tissue was ground to a fine powder in liquid nitrogen using mortar and pestle. The powder (200 mg) was transferred into a glass centrifuge vial and defatted with 2 mL hexane (1:10 w/v) on a Wrist Action Shaker (Burrell Scientific, PA, USA) for 2 h at room temperature. After centrifugation at 1500 rpm for 15 min, the supernatant was discarded. The pellet was rinsed with 2 mL hexane, centrifuged and air-dried for 15 min. The defatted material was extracted with 2 mL of 70% (v/v) methanol/water at 55°C for 2 h using rotation in an oven, with intermittent manual shaking 2–3 times. A final vigorous shaking was performed for 15 min on a Wrist action shaker (Burrell Scientific, Pittsburgh, PA, USA) at room temperature. The samples were centrifuged at 1500 rpm for 15 min and the supernatant (S1) was collected in new capped vial. The residue was rinsed again with 0.5 mL 70% methanol, centrifuged and the supernatant (S2) was collected and pooled with S1. The total supernatant volume was recorded before hydrolysis. The combined samples (S1 + S2) were hydrolysed for 1 h at 60°C with 0.5 N NaOH at a ratio of 3:5 (v/v). After hydrolysis, the samples were immediately neutralized using 0.5 N HCl at a ratio of 0.4 mL for every 0.5 mL extract. The hydrolysate was cooled and purified via solid phase extraction using 10 mL Waters HLB columns (Waters, Mississauga, ON, Canada). The eluted lignan fractions were collected in glass vials and dried using a rotary evaporator (Heidolph instrument Gamborg, Germany). The dried material was dissolved in methanol:water (50:50), filtered and injected for UPLC-MS analysis using a commercially available SDG standard (Chromadex, Irvine, CA, USA) as reference.
An Acquity H-Class, quaternary pump UPLC system (Waters, Mississauga, ON, Canada) equipped with in-line degassing, diode array detector (DAD), robotic autosampler, sample and column temperature controls and tandem quad mass spectrometer (TQD) was used for lignan profiling analysis. A ternary solvent system for UPLC-MS analysis consisting of water, acetonitrile and 10% formic acid in water was used for UPLC-MS analysis. UV–vis spectra were recorded from 210–600 nm, and the MS was run in ESI mode, 3000 V capillary voltage, in scanning mode from 100–2000 a.m.u., with a fragmentation setting of 150 V, 13.0 L/min carrier gas (N2) flow at 350°C and 60 psi to ensure identity of the profiled metabolites. The post-hydrolysis SDG lignan peak was identified and quantitated through comparison (UV–VIS absorption, retention time) to a commercial standard. Other phenolic compounds, including hydroxycinnamic acids liberated by the base hydrolysis were present but were not quantified. A standard curve for SDG was created, relating integrated peak area (mAU*s) (Y) versus concentration of SDG (mg/mL) (X). In brief, 1 mg of authentic standard was dissolved in 50% methanol and a serial dilution was created in triplicate, halving the concentration each time. The resulting standard curve was linear from 0.5 mg/mL to 0.00781 mg/mL (R2 = 0.9901) and was used to determine SDG content in relation to developmental stage (DAA). For each of the six developmental stages, three extractions and HPLC analyses were performed from three biological replicates and the values were presented as the mean of the three data points.
Heterologous expression of flax UGTs in yeast
The pYES2/NT C plasmid constructs harbouring the cDNA of the five UGTs described in this study were used to transform yeast INVSc1 strains using S.c. EasyComp transformation™ kit (Invitrogen, CA, USA). The flax UGT cDNA of Genebank accession JN088324.1  was similarly transformed for functional comparative analyses. Single transformant INVSc1 yeast colonies were inoculated into 15 mL of Saccharomyces cerevisiae minimal media without uracil (SC-U, prepared as recommended by Invitrogen) supplemented with 2% raffinose and grown for 3 days under shaking at 30°C until the OD600 reached 2.0. The culture was diluted in 50 mL of induction medium (SC-U supplemented with 1% raffinose and 2% galactose) to achieve an initial OD600 of 0.4. The culture was further incubated under shaking at 30°C for 24 hours, with 5 mL sub-sample collection at 0, 4, 8, 12 and 24 h to monitor the protein expression. The OD600 for each time point was recorded. The induced yeast cells were harvested by centrifugation at 1,500 g for 5 min at 4°C. The cells were washed using 500 μL cold sterile distilled water and centrifuged. The pellets were washed again at 4°C in 500 μL of lysis buffer (50 mM sodium phosphate, pH 7.4 supplemented with 5% glycerol and 1 mM PMSF). After centrifugation, the cells were mechanically disrupted by vortexing for 30 seconds in the presence of an equal volume of 425–600 μm acid-washed glass beads (Sigma Aldrich, Canada). After vortexing, the sample was incubated on ice for 30 seconds. The vortexing and incubation cycle was repeated 4 times to ensure complete cell lysis. The lysates were centrifuged at 18,620 g for 10 min at 4°C and the supernatant was collected. The optimum induction time for all the UGTs was monitored by western blot using equal amount of proteins and antibodies raised against the anti-ExpressTM epitope present between the 6x Histidine tag and the multiple cloning site of the construct. The polyhistidine containing recombinant proteins was purified using the ProBond™ (Invitrogen, CA, USA) purification system following manufacturer’s instruction. The purified enzymes were concentrated using 0.5 mL UltracelR-10 k Amicon membrane column (Millipore, Ireland). Protein concentrations were determined using the Bradford protein assay kit (BioRad Laboratories, Canada).
The crude and purified recombinant protein extracts obtained from the yeast cultures harboring the five different UGT cDNAs reported in this study and the one derived from JN088324.1  were reacted with different aglycone substrates including SECO (Chromadex, Irvine, CA, USA), sillibinin, quercetin, kaempferol and the phenolic acids coumaric acid, caffeic acid, sinnapic acid, cinnamic acid, and ferulic acid (Sigma Aldrich, Canada). The 100 μL reaction mixture consisted of a reaction buffer (50 mM sodium phosphate, 1 mM PMSF, 5% glycerol, pH 7.4), 280 μM aglycone substrate (acceptor for glycosylation), and 1.64 mM UDP-glucose (sugar donor) (Sigma Aldrich, Canada). The reaction mixtures were pre-incubated at 30°C for 10 min and the reactions were initiated with the addition of 50 μg of enzyme. After incubation at 30°C for 30 min, the reactions were stopped with 100 μL of 0.5% trifluoroacetic acid in acetonitrile. The reaction mixtures were purified using 0.2 μm filters (Pall Life Sciences, Mississauga, ON, Canada) to remove any particulates that might form during the reaction. The separation and identification of the reactants and products derived from the enzyme assays were carried out using a Waters H-Class Acquity UPLC system (Waters, Missisauga, ON) equipped with a TQD tandem mass spectrometer. The formation of glycosylated products was monitored by examining the masses and the principle fragments of eluted peaks via ESI–mass spectrometry. Two parallel MS2 scans were performed ranging from 120–800 a.m.u., using 15 and 45 V cone voltages. Selected ion recording (SIR) spectra were also collected to enhance the sensitivity of detection of SECO, SMG and SDG. The capillary voltage was 3 kV, the extractor set to 3 V, and RF lens at 0.1 V. Chromatographic conditions consisted of a binary gradient system composed of 3% formic acid in water (A) and acetonitrile (B), varied according to the following gradient: t0, A = 68%; t1 = 4.4 min, A = 0%; t2 = 6 min, A = 0% isocratic; t3 = 7 min, A = 68%; t4 = 8 min, A = 68% isocratic. Peaks detected at 280 nm, indicative of phenolic compounds, were validated using authentic standards (SECO and SDG) purchased from Chromadex (Chromadex, Irvine, CA, USA). A standard curve for SDG was created as detailed above. Standard purified SMG was prepared as described by .
Kinetic and biochemical characterization of UGT74S1
Ranges of pH from 6.0 to 9.0, temperature from 25°C to 50°C, enzyme concentration from 10 to 120 μg and two concentrations (1 and 10 mM) of seven metal cofactors (NaCl, KCl, MgCl2, MnCl2, CaCl2, FeSO4 and CuSO4) were tested in 100 μL reaction mixture for determining the optimal pH, temperature, enzyme concentration and metal cofactor effect on the enzyme activity. To determine the initial velocity of the recombinant UGT74S1 enzyme, a time course (5, 10, 15, 30, 45, 60 min) study using the optimum enzyme concentration and fixed excess substrate (280 μM SECO; 1.67 mM UDP-glucose) concentration was conducted at 30°C, pH 8. The linearity was maintained in assays up to 30 min at 30°C. The initial velocity of the reaction was measured at 10 min, where no more than 10% of SECO was converted to SDG at this time point. Then, the assays were carried out using various substrate concentrations (70–1400 μM SECO with UDP-glucose fixed at 1.67 mM; 0.82–6.56 mM UDP-glucose with SECO fixed at 280 μM), under optimum conditions, for 30 min, for the determination of kinetic parameters. The apparent Vmax and Km value for the glucosyl donor and acceptor substrate in the presence of 80 μg of the enzyme were determined from Lineweaver-Burk plots. The kcat was determined by dividing Vmax by the enzyme concentration.
This research is part of the Total Utilization Flax Genomics (TUFGEN) project funded by Genome Canada/Genome Prairie with financial contribution of Flax Council of Canada to BF. The authors thank David Main (Crops and Livestock Research Centre, Charlottetown, PEI), Leonardo Galindo (University of Alberta, Edmonton, AB) and Elsa Reimer and Andrzej Walichnowski (Cereal Research Centre, Winnipeg, MB) for their technical assistance and helpful proof reading.
- Dixon RA: Phytoestrogens. Annu Rev Plant Biol. 2004, 55: 225-261. 10.1146/annurev.arplant.55.031903.141729.View ArticlePubMedGoogle Scholar
- Thompson LU, Boucher BA, Liu Z, Cotterchio M, Kreiger N: Phytoestrogen content of foods consumed in Canada, including isoflavones, lignans, and coumestan. Nutr Cancer. 2006, 54: 184-201. 10.1207/s15327914nc5402_5.View ArticlePubMedGoogle Scholar
- Arroo RRJ, Androutsopoulos V, Beresford K, Ruparelia K, Surichan S, Wilsher N, Potter GA: Phytoestrogens as natural prodrugs in cancer prevention: dietary flavonoids. Phytochem Rev. 2009, 8: 375-386. 10.1007/s11101-009-9128-6.View ArticleGoogle Scholar
- Prakash D, Gupta C: Role of phytoestrogens as nutraceuticals in human health. Pharmacol Online. 2011, 1: 510-523.Google Scholar
- Pan JY, Chen SL, Yang MH, Wu J, Sinkkonen J, Zoud K: An update on lignans: natural products and synthesis. Nat Prod Rep. 2009, 26: 1251-1292. 10.1039/b910940d.View ArticlePubMedGoogle Scholar
- Noguchi A, Fukui Y, Iuchi-Okada A, Kakutani S, Satake H, Iwashita T, Nakao M, Umezawa T, Ono E: Sequential glucosylation of furofuran lignan (+)-sesaminol by Sesamum indicum UGT71A9 and UGT94D1 glucosyltransferase. Plant J. 2008, 54: 415-427. 10.1111/j.1365-313X.2008.03428.x.View ArticlePubMedGoogle Scholar
- Yousefzadi M, Sharif M, Behmanesh M, Moyano E, Bonfill M, Cusido RM, Palazon J: Podophyllotoxin: current approaches to its biotechnological production and future challenges. Eng Life Sci. 2010, 10: 281-292. 10.1002/elsc.201000027.View ArticleGoogle Scholar
- Berim A, Ebel R, Schneider B, Petersen M: UDP-glucose:(6-methoxy) podophyllotoxin 7-O-glucosyltransferase from suspension cultures of Linum nodiflorum. Phytochemistry. 2008, 69: 374-381. 10.1016/j.phytochem.2007.07.030.View ArticlePubMedGoogle Scholar
- Ono E, Kim HJ, Murata J, Morimoto K, Okazawa A, Kobayashi A, Umezawa T, Satake H: Molecular and functional characterization of novel furofuran-class lignan glucosyltransferase from Forsythia. Plant Biotechnol. 2010, 27: 317-324. 10.5511/plantbiotechnology.27.317.View ArticleGoogle Scholar
- Jhala AM, Hall LM: Flax (Linum usitatissimum L.): current uses and future applications. Aust J Basic Appl Sci. 2010, 4: 4304-4312.Google Scholar
- Rabetafika HN, Remoortel VV, Danthine S, Paquot M, Blecker C: Flaxseed proteins: food uses and health benefits. Intern J Food Sci Technol. 2011, 46: 221-228. 10.1111/j.1365-2621.2010.02477.x.View ArticleGoogle Scholar
- Buck K, Zaineddin AK, Vrieling A, Heinz J, Linseisen J, Flesch-Janys D, Chang-Claude J: Estimated enterolignans, lignan-rich foods, and fibre in relation to survival after postmenopausal breast cancer. Br J Cancer. 2011, 105: 1151-1157. 10.1038/bjc.2011.374.PubMed CentralView ArticlePubMedGoogle Scholar
- Webb AL, McCullough ML: Dietary lignans: potential role in cancer prevention. Nutr Cancer. 2005, 51: 117-131. 10.1207/s15327914nc5102_1.View ArticlePubMedGoogle Scholar
- Adolphe JL, Whiting SJ, Juurlink BHJ, Thorpe LU, Alcorn J: Health effects with consumption of the flax lignan secoisolariciresinol diglucoside. Br J Nutr. 2010, 103: 929-938. 10.1017/S0007114509992753.View ArticlePubMedGoogle Scholar
- Touré A, Xueming X: Flaxseed lignans: source, biosynthesis, metabolism, antioxidant activity, Bio-active components, and health benefits. Compr Rev Food Sci Food Saf. 2010, 9: 261-269. 10.1111/j.1541-4337.2009.00105.x.View ArticleGoogle Scholar
- Gachon CM, Langlois-Meurinne M, Saindrenan P: Plant secondary metabolism glycosyltransferases: the emerging functional analysis. Trends Plant Sci. 2005, 10: 542-549. 10.1016/j.tplants.2005.09.007.View ArticlePubMedGoogle Scholar
- Osmani SA, Bak S, Møller BL: Substrate specificity of plant UDP-dependent glycosyltransferases predicted from crystal structures and homology modeling. Phytochem. 2009, 70: 325-347. 10.1016/j.phytochem.2008.12.009.View ArticleGoogle Scholar
- Bowles D, Isayenkova J, Lim E, Poppenberger B: Glycosyltransferases: managers of small molecules. Curr Opin Plant Biol. 2005, 8: 254-263. 10.1016/j.pbi.2005.03.007.View ArticlePubMedGoogle Scholar
- Bowles D, Lim E-K: Glycosyltransferases of small molecules: their roles in plant biology. Encyclopedia of Life Sciences. 2010, Chichester: (ELS) John Wiley & Sons, Ltd, doi:10.1002/9780470015902.a0021537Google Scholar
- Wang J, Hou B: Glycosyltransferases: key players involved in the modification of plant secondary metabolites. Front Biol. 2009, 4: 39-46. 10.1007/s11515-008-0111-1.View ArticleGoogle Scholar
- Geisler-Lee J, Geisler M, Coutinho PM, Segerman B, Nishikubo N, Takahashi J, Aspeborg H, Djerbi S, Master E, Andersson-Gunneras S, Sundberg B: Poplar carbohydrate-active enzymes. Gene identification and expression analyses. Plant Physiol. 2006, 140: 946-962. 10.1104/pp.105.072652.PubMed CentralView ArticlePubMedGoogle Scholar
- Yonekura-Sakakibara K, Hanada K: An evolutionary view of functional diversity in family 1 glycosyltransferases. Plant J. 2011, 66: 182-193. 10.1111/j.1365-313X.2011.04493.x.View ArticlePubMedGoogle Scholar
- Caputi L, Malnoy M, Goremykin V, Nikiforova S, Martens S: A genome-wide phylogenetic reconstruction of family 1 UDP-glycosyltransferases revealed the expansion of the family during the adaptation of plants to life on land. Plant J. 2012, 69: 1030-1042. 10.1111/j.1365-313X.2011.04853.x.View ArticlePubMedGoogle Scholar
- Ross J, Li Y, Lim EK, Bowles DJ: Higher plant glycosyltransferases. Genome Biol. 2001, 2: 3004.1-3004.6.24.View ArticleGoogle Scholar
- Witte S, Moco S, Vervoort J, Matern U, Martens S: Recombinant expression and functional characterisation of regiospecific flavonoid glycosyltransferases from Hieracium pilosella L. Planta. 2009, 229: 1135-1146. 10.1007/s00425-009-0902-x.View ArticlePubMedGoogle Scholar
- Wang Z, Hobson N, Galindo L, Zhu S, Shi D, McDill J, Yang L, Hawkins S, Neutelings G, Datla R, Lambert G, Galbraith DW, Grassa CJ, Geraldes A, Cronk QC, Cullis C, Dash PK, Kumar PA, Cloutier S, Sharpe AG, Wong GK, Wang J, Deyholos MK: The genome of flax (Linum usitatissimum) assembled de novo from short shotgun sequence reads. Plant J. 2012, 72: 461-473. 10.1111/j.1365-313X.2012.05093.x.View ArticlePubMedGoogle Scholar
- Barvkar VT, Pardeshi VC, Kale SM, Kadoo NY, Gupta VS: Phylogenomic analysis of UDP glycosyltransferase 1 multigene family in Linum usitatissimum identified genes with varied expression patterns. BMC Genomics. 2012, 13: 175-10.1186/1471-2164-13-175. doi:10.1186/1471-2164-13-175PubMed CentralView ArticlePubMedGoogle Scholar
- Hemmati S, Cosima BI, Heimendahl V, Michael K, Wilhelm A, Schmidt TJ, Fuss E: Pinoresinol-lariciresinol reductases with opposite enantiospecificity determine the enantiomeric composition of lignans in the different organs of Linum usitatissimum L. Planta Med. 2010, 76: 928-934. 10.1055/s-0030-1250036.View ArticlePubMedGoogle Scholar
- Hano C, Martin I, Fliniaux O, Legrand B, Gutierrez L, Arroo RR, Mesnard F, Lamblin F, Lainé E: Pinoresinol–lariciresinol reductase gene expression and secoisolariciresinol diglucoside accumulation in developing flax (Linum usitatissimum) seeds. Planta. 2006, 224: 1291-1301. 10.1007/s00425-006-0308-y.View ArticlePubMedGoogle Scholar
- von Heimendahl CBI, Schafer KM, Eklund P, Sjohlm R, Schmidt TJ, Fuss E: Pinoresinol–lariciresinol reductases with different stereospecificity from Linum album and Linum usitatissimum. Phytochem. 2005, 66: 254-1263.View ArticleGoogle Scholar
- Kim HJ, Ono E, Morimoto K, Yamagaki T, Pkazawa A, Kobayashi A, Satake H: Metabolic engineering of lignan biosynthesis in Forsythia cell culture. Plant Cell Physiol. 2009, 50: 2200-2209. 10.1093/pcp/pcp156.View ArticlePubMedGoogle Scholar
- Mackenzie PI, Owens IS, Burchell B, Bock K, Bairoch A, Belanger A, Fournel-Gigleux S, Green M, Hum DW, Iyanagi T: The UDP glycosyltransferase gene superfamily: recommended nomenclature update based on evolutionary divergence. Pharmacogenetics. 1997, 7: 255-269. 10.1097/00008571-199708000-00001.View ArticlePubMedGoogle Scholar
- Li X, Yuan J-P, Xu S-P, Wang J-H, Liu X: Separation and determination of secoisolariciresinol diglucoside oligomers and their hydrolysates in the flaxseed extract by high-performance liquid chromatography. J Chromatogr A. 2008, 1185: 223-232. 10.1016/j.chroma.2008.01.066.View ArticlePubMedGoogle Scholar
- Jones P, Vogt T: Glycosyltransferases in secondary plant metabolism: tranquilizers and stimulants controllers. Planta. 2001, 213: 164-174. 10.1007/s004250000492.View ArticlePubMedGoogle Scholar
- Hansen EH, Osmani SA, Kristensen C, Møller BL, Hansen J: Substrate specificities of family 1 UGTs gained by domain swapping. Phytochem. 2009, 70: 473-482. 10.1016/j.phytochem.2009.01.013.View ArticleGoogle Scholar
- Noguchi O, Horikawa M, Fukui Y, Fukuchi-Mizutani M, Iuchi-Okada A, Ishiguro M, Kiso Y, Nakayama T, Ono E: Local differentiation of sugar donor specificity of flavonoid glycosyltransferase in Lamiales. Plant Cell. 2009, 21: 1556-1572. 10.1105/tpc.108.063826.PubMed CentralView ArticlePubMedGoogle Scholar
- Modolo LV, Escamilla-Treviño LL, Dixon RA, Wang X: Single amino acid mutations of Medicago glycosyltransferase UGT85H2 enhance activity and impart reversibility. FEBS Lett. 2009, 583: 2131-2135. 10.1016/j.febslet.2009.05.046.View ArticlePubMedGoogle Scholar
- Venglat P, Xiang D, Qiu S, Stone SL, Tibiche C, Cram D, Alting-Mees M, Nowak J, Cloutier S, Deyholos M, Bekkaoui F, Sharpe A, Wang E, Rowland G, Selvaraj G, Datla R: Gene expression analysis of flax seed development. BMC Plant Biol. 2011, 11: 74-10.1186/1471-2229-11-74.PubMed CentralView ArticlePubMedGoogle Scholar
- Owens DK, McIntosh CA: Identification, recombinant expression, and biochemical characterization of a flavonol 3-O-glucosyltransferase clone from Citrus paradisi. Phytochemistry. 2009, 70: 1382-1391. 10.1016/j.phytochem.2009.07.027.View ArticlePubMedGoogle Scholar
- Zimowski J: Specificity and properties of UDP-galactose:tomatidine galactosyltransferase from tomato leaves. Plant Sci. 1998, 136: 139-148. 10.1016/S0168-9452(98)00109-5.View ArticleGoogle Scholar
- Yamazaki M, Gong Z, Fukuchi-Mizutani M, Fukui Y, Tanaka Y, Kusumi T, Saito K: Molecular cloning and biochemical characterization of a novel anthocyanin 5-O-glycosyltransferase by mRNA differential display for plant forms regarding anthocyanin. J Biol Chem. 1999, 274: 7405-7411. 10.1074/jbc.274.11.7405.View ArticlePubMedGoogle Scholar
- Taguchi G, Ubukata T, Hayashida N, Yamamoto H, Okazaki M: Cloning and characterization of a glycosyltransferase that reacts on 7-hydroxyl group of flavonol and 3-hydroxyl group of coumarin from tobacco cells. Arch Biochem Biophys. 2003, 420: 95-102. 10.1016/j.abb.2003.09.027.View ArticlePubMedGoogle Scholar
- Noguchi A, Saito A, Homma Y, Nakao M, Sasaki N, Nishino T, Takahashi S, Nakayama T: A UDP–glucose: isoflavone 7-O-glycosyltransferase from the roots of soybean (Glycine max) seedlings. Purification, gene cloning, phylogenetics, and an implication for an alternative strategy of enzyme catalysis. J Biol Chem. 2007, 282: 23581-23590. 10.1074/jbc.M702651200.View ArticlePubMedGoogle Scholar
- Ruuska SA, Girke T, Benning C, Ohlrogge JB: Contrapuntal networks of gene expression during Arabidopsis seed filling. Plant Cell. 2002, 13: 1191-1206.View ArticleGoogle Scholar
- Dhaubhadel S, Gijzen M, Moy P, Farhangkhoee M: Transcriptome analysis reveals a critical role of CHS7 and CHS8 genes for isoflavonoid synthesis in soybean seeds. Plant Physiol. 2007, 143: 326-338.PubMed CentralView ArticlePubMedGoogle Scholar
- Vrinten P, Hu Z, Munchinsky M-A, Rowland G, Qiu X: Two FAD3 desaturase genes control the level of linolenic acid in flax seed. Plant Physiol. 2005, 139: 79-87. 10.1104/pp.105.064451.PubMed CentralView ArticlePubMedGoogle Scholar
- Fofana B, Cloutier S, Duguid S, Ching J, Rampitsch C: Gene expression of stearoyl-ACP desaturase and Δ12 fatty acid desaturase 2 is modulated during seed development of flax (Linum usitatissimum). Lipids. 2006, 41: 705-712. 10.1007/s11745-006-5021-x.View ArticlePubMedGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.View ArticlePubMedGoogle Scholar
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.View ArticlePubMedGoogle Scholar
- Popova IE, Hall C, Kubátová A: Determination of lignans in flaxseed using liquid chromatography with time-of-flight mass spectrometry. J Chromatogr A. 2009, 1216: 217-229. 10.1016/j.chroma.2008.11.063.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.