OPEN ACCESS ARTICLE

This Article Free via Open Access: OA OA Abstract Full Text (PDF) Supplemental Data OA All Versions of this Article: 20/8/2160    most recent tpc.108.058040v2 tpc.108.058040v1 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 PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via CrossRef Google Scholar Articles by Yonekura-Sakakibara, K. Articles by Saito, K. PubMed PubMed Citation Articles by Yonekura-Sakakibara, K. Articles by Saito, K. Agricola Articles by Yonekura-Sakakibara, K. Articles by Saito, K.

Comprehensive Flavonol Profiling and Transcriptome Coexpression Analysis Leading to Decoding Gene–Metabolite Correlations in Arabidopsis^[W],[OA]

^a RIKEN Plant Science Center, Tsurumi-ku, Yokohama 230-0045, Japan
^b Graduate School of Pharmaceutical Sciences, Chiba University, Inage-ku, Chiba 263-8522, Japan

³ Address correspondence to ksaito{at}faculty.chiba-u.jp.

	ABSTRACT

TOP ABSTRACT INTRODUCTION RESULTS DISCUSSION METHODS REFERENCES

 
To complete the metabolic map for an entire class of compounds, it is essential to identify gene–metabolite correlations of a metabolic pathway. We used liquid chromatography–mass spectrometry (LC-MS) to identify the flavonoids produced by Arabidopsis thaliana wild-type and flavonoid biosynthetic mutant lines. The structures of 15 newly identified and eight known flavonols were deduced by LC-MS profiling of these mutants. Candidate genes presumably involved in the flavonoid pathway were delimited by transcriptome coexpression network analysis using public databases, leading to the detailed analysis of two flavonoid pathway genes, UGT78D3 (At5g17030) and RHM1 (At1g78570). The levels of flavonol 3-O-arabinosides were reduced in ugt78d3 knockdown mutants, suggesting that UGT78D3 is a flavonol arabinosyltransferase. Recombinant UGT78D3 protein could convert quercetin to quercetin 3-O-arabinoside. The strict substrate specificity of UGT78D3 for flavonol aglycones and UDP-arabinose indicate that UGT78D3 is a flavonol arabinosyltransferase. A comparison of flavonol profile in RHM knockout mutants indicated that RHM1 plays a major role in supplying UDP-rhamnose for flavonol modification. The rate of flavonol 3-O-glycosylation is more affected than those of 7-O-glycosylation by the supply of UDP-rhamnose. The precise identification of flavonoids in conjunction with transcriptomics thus led to the identification of a gene function and a more complete understanding of a plant metabolic network.

	INTRODUCTION

TOP ABSTRACT INTRODUCTION RESULTS DISCUSSION METHODS REFERENCES

 
Plants employ many diverse metabolic pathways to produce >200,000 compounds (Dixon and Strack, 2003). Completion of the Arabidopsis thaliana and rice (Oryza sativa) genome sequencing projects has allowed (1) the annotation of genes involved in metabolic pathways based on homology and (2) genome-wide identification of whole metabolic pathways. Assigning functions to each of the genes that correlate with metabolites in a pathway is a prerequisite for developing a comprehensive understanding of complicated metabolic pathways. Because of its small genome and the development of extensive genetic, molecular, and physiological tools, Arabidopsis may be the only plant in which the complete network of secondary metabolism can currently be characterized within a single plant species.

The flavonoids comprise one of the major secondary metabolite groups, with >7000 known compounds (Harborne et al., 1999; Andersen and Markham, 2006). Formation of the basic flavonoid skeleton has been well studied in terms of molecular biology and natural products chemistry. Flavonol and anthocyanidin regulatory and biosynthetic pathways have been characterized, and the corresponding genes have been isolated from various plants (Figure 1 ; Andersen and Markham, 2006; Tanaka and Filippa, 2006). However, the pathways for sequential modification, such as glycosylation, acylation, and methylation, are still relatively unexplored even though modification produces a huge chemical diversity and is essential for the stable accumulation of flavonoids. A wide variety of flavonoids have been isolated and identified from specific organs, such as flowers or fruits (Andersen and Markham, 2006), but there has been no comprehensive description of flavonoid metabolism in the organs of a single plant species. Although the flavonoids are important as flower pigments, UV-B protectants, phytoalexins, signaling molecules, and regulators of auxin transport (Dooner et al., 1991; Dixon and Paiva, 1995), the precise relationship between structure and function is also largely unclear.

View larger version (14K): [in this window] [in a new window] [as a PowerPoint slide]   Figure 1. The Arabidopsis Flavonoid Biosynthetic Pathway. 4CL, 4-coumarate:CoA ligase; CHS, chalcone synthase; CHI, chalcone isomerase; F3H, flavanone 3-hydroxylase; F3'H, flavonoid 3'-hydroxylase; FLS, flavonol synthase; OMT1, O-methyltransferase; UGT, UDP-dependent glycosyltransferase; DFR, dihydroflavonol 4-reductase; LDOX/ANS, leucoanthocyanidin dioxygenase/anthocyanidin synthase; AAT, anthocyanin acyltransferase.

Because of the rich research resources available for Arabidopsis, it is the best candidate for cataloging the entire flavonoid pathway. The integration of metabolome analysis with transcriptome analysis has been successfully used in Arabidopsis to decipher gene functions in biological processes (Hirai et al., 2005, 2007; Tohge et al., 2005). Transcriptome coexpression analysis is also a powerful tool for the efficient targeting of genes within a multigene family (Yonekura-Sakakibara et al., 2007; Saito et al., 2008). Use of these novel tools in combination with reverse genetics and a biochemical approach using recombinant proteins can provide robust evidence from which a whole pathway can be constructed.

In this study, we attempted to complete the identification of flavonols in each of the organs and flavonol annotation using flavonoid mutant lines. Taken together with previous reports, a total of 32 flavonol end or intermediate products were listed. Thirty of 32 flavonols were detected, and the structures of 19 compounds were annotated. Transcriptome coexpression analysis using known flavonoid genes as queries suggested 139 genes that may be involved in flavonoid metabolism. Finally, integration of comprehensive flavonol identification/annotation with transcriptome coexpression analysis identified a gene encoding flavonol 3-O-arabinosyltransferase and revealed the physiological role of a gene coding for UDP-rhamnose synthase in flavonoid biosynthesis.

	RESULTS

TOP ABSTRACT INTRODUCTION RESULTS DISCUSSION METHODS REFERENCES

 
Comprehensive Flavonol Profiling in Wild-Type and tt4 Mutants
To catalog all of the flavonol derivatives in Arabidopsis, flavonol profiling in the flowers, leaves, stems, and roots of the wild type and tt4 mutants were performed using ultraperformance liquid chromatography (UPLC)–photodiode array (PDA)–electrospray ionization (ESI)/quadrupole time-of-flight (Q-TOF)/MS (Figure 2A ). The metabolites detected only in the wild type should be flavonoid derivatives because TT4 encodes chalcone synthase, the first committed enzyme in flavonoid biosynthesis, and no flavonoid derivatives were detected in the tt4 mutant (Shirley et al., 1995; von Roepenack-Lahaye et al., 2004).

View larger version (21K): [in this window] [in a new window] [as a PowerPoint slide]   Figure 2. UPLC-PDA-MS Analyses of the Extracts from Leaves, Flowers, Roots, and Stems of Wild-Type (Col-0) and tt4 Plants. (A) UPLC-PDA and mass chromatograms of aqueous methanol extracts of the Arabidopsis wild type and the tt4 knockout mutant. Absorbance at 320 nm was used for the detection of flavonols. Labels correspond to compounds shown in Figure 3. (B) Principal component analyses of the data obtained by LC-MS. Proportions of the first and second components are in parentheses. Wild-type results are in red, and those from the tt4 mutant are in blue.

View larger version (50K): [in this window] [in a new window] [as a PowerPoint slide]   Figure 3. Flavonol Glycosides in Arabidopsis. Asterisks indicate that the compounds were identified based on a comparison of retention times and UV/mass spectra of the standards used in this study. R1=H, kaempferol; R1=OH, quercetin; R1=OMe, isorhamnetin. The presence of a compound in the indicated tissue is denoted by a gray box.

The MS-detectable peaks in the wild type and tt4 were compared to extract all flavonoid-related peaks (see Supplemental Figure 1 and Supplemental Data Set 1 online). Principal component analysis was conducted with 1804 MS peaks in leaves, 1864 in flowers, 1576 in roots, and 2138 in stems, indicating the presence of clear clusters according to the wild-type and tt4 genotypes (Figure 2B). The wild-type–specific peaks (53 peaks in leaves, 322 in flowers, 152 in roots, and 75 in stems) were picked up using a high PC1 loading value, and most of the peaks are due to the ions of flavonols (see Supplemental Figure 1 online). For example, 42 of 53 peaks are three major kaempferol glycosides in leaves (f1 to f3). A total of 30 flavonol derivatives were detected in Arabidopsis (four in leaves, 25 in flowers, 14 in roots, and four in stems). Flowers, in which most flavonol derivatives were detected, were used for further analysis.

Peak Annotation by Comparison with Flavonoid Biosynthetic Gene Knockout Mutants
Twelve flavonols (f1 to f8, f17, f18, f20, and f24) in Arabidopsis have been identified by UV spectroscopy, MS, ¹H-NMR, and ¹³C-NMR spectroscopy or annotated using UPLC-PDA-ESI/Q-TOF/MS (Veit and Pauli, 1999; Jones et al., 2003; Tohge et al., 2005, 2007; Kerhoas et al., 2006). We identified two additional flavonols (f16 and f28) by comparing retention times and mass-to-charge ratio (m/z) values with the standard compounds (Table 1 ; see Supplemental Data Set 2 online). To annotate the remaining 18 uncharacterized flavonols, mutants lacking CHS, F3'H, flavonol 3'-O-methyltransferase (OMT1), flavonoid 3-O-glucosyltransferase (UGT78D2), flavonol 3-O-rhamnosyltransferase (UGT78D1), flavonol 7-O-rhamnosyltransferase (UGT89C1), or anthocyanin 5-O-glucosyltransferase (UGT75C1) were used for flavonol analysis. Of the 14 identified (f1 to f8, f16 to f18, f20, f24, and f28), flavonol derivatives missing from each mutant were well correlated with the loss of gene function (Figure 4 , Table 1). No quercetin and isorhamnetin conjugates (f5, f6, f8, f16, f24, and f28) accumulated in tt7, which lacks F3'H. No isorhamnetin conjugates (f28) were detected in omt1 mutants. The mutants ugt78d1 and ugt89c1 had no conjugates of flavonol 3-O-rhamnoside (f1, f5, f17, f20, and f24) and flavonol 7-O-rhamnoside (f1 to f3, f5, f6, and f8), respectively. Interestingly, in the ugt78d2 mutant, accumulation of the kaempferol/quercetin 3-O-glucosides (f2, f3, f6, and f8) was drastically reduced but detectable. Other kaempferol/quercetin 3-O-glucosides (f16 and f18) also decreased. This result is in agreement with the finding that cyanidin derivatives were detected at low levels in the ugt78d2 mutant (Tohge et al., 2005). On the other hand, accumulation of the isorhamnetin 3-O-glucosides (f28) was similar or slightly increased. These data indicate that another enzyme is involved in 3-O-glucosylation of isorhamnetin. Flavonol profiles in the flowers of the wild type and ugt75c1 were essentially identical, demonstrating that UGT75C1 is not involved in flavonol glycosylation at any detectable level.

View larger version (36K): [in this window] [in a new window] [as a PowerPoint slide]   Figure 4. UPLC-PDA-MS Analyses of the Extracts from Flowers of Wild-Type (Col-0) and Flavonoid-Defective Mutants. Absorbance at 320 nm was used for detection. A5GlcT, anthocyanin 5-O-glucosyltransferase; F3RhaT, flavonol 3-O-rhamnosyltransferase; F7RhaT, flavonol 7-O-rhamnosyltransferase; Fd3GlcT, flavonoid 3-O-glucosyltransferase. Labels correspond to compounds shown in Figure 3. The labels in parentheses indicate compounds at <5% in the peak.

New Flavonoid Biosynthetic Genes Deduced by Transcriptome Coexpression Analysis
Transcriptome coexpression analysis (Saito et al., 2008) was used to identify all of the flavonoid-related genes and was executed with the 24 genes encoding flavonoid biosynthetic enzymes or transcription factors listed in Table 2 as query. The lowest correlation coefficient (r) between query genes, except for the flavonol 7-O-glucosyltransferase gene (UGT73C6) and OMT1 is 0.527. UGT73C6 and OMT1 were excluded because UGT73C6 may cross-hybridize with UGT73C5, a brassinosteroid-related UGT, and the OMT1 coexpression profile suggests that it is under the regulation of lignin biosynthesis but not of flavonoid biosynthesis. The coexpression relationships of all Arabidopsis genes exhibiting correlation coefficients > 0.525 (r > 0.525, all data sets version 3) with the query genes are shown in Supplemental Figure 2 online and Table 3 . The cutoff value was set according to Aoki et al. (2007). The clusters of coexpressed flavonoid genes were classified into two major groups. One cluster contains the general flavonoid/flavonol subgroup; another cluster represents the anthocyanin subgroup. UGT73C6 and OMT1 did not belong to either of these clusters. The 139 genes were highly correlated with the general flavonoid-related or anthocyanin-related genes (see Supplemental Figure 2 online). Twenty-four genes had more than two correlations with flavonoid-related genes. One-third of the genes listed in Table 3 were coregulated with flavonoid-related genes by flavonoid-specific MYBs (PAP1, MYB11, MYB12, and MYB111; Tohge et al., 2005, Stracke et al., 2007). The genes encoding embryo-specific protein (At5g62210), 3-keto-acyl-CoA thiolase (At5g48880), cryptochrome 3 (At5g24850), UGT84A1 (At4g15480), cinnamoyl-CoA reductase-like protein (At2g23910), early light-induced protein (At3g22840), TT8 (At4g09820), and glutathione S-transferase (GST; At1g10370) were regulated in a manner similar to the flavonoid biosynthetic genes by high irradiance and blue light in a CRYPTOCHROME1 (CRY1) photoreceptor and/or LONG HYPOCOTYL5 (HY5) transcription factor-dependent manner (Kleine et al., 2007). Other genes were also regulated by HY5 under UV-B radiation, much like the flavonoid biosynthetic genes (Table 3; Brown et al., 2005). For further analysis, we focused on UGT78D3 as a putative flavonol glycosyltransferase and RHM1, which was selected from the genes highly correlated with the flavonol biosynthetic genes.

To test our hypothesis, T-DNA insertion mutant SALK_114099, designated as ugt78d3, was used for reverse genetics analysis. T-DNA was inserted between –122 and –100 bp upstream from the start codon of UGT78D3 (Figure 5A ). UGT78D3 transcripts were abundant in floral buds but present at low levels in siliques and stems (Figure 5B). Real-time PCR showed that the accumulation level of UGT78D3 transcripts in flowers of the homozygous ugt78d3 mutant decreased to 30% of those in the wild type (Figure 5C). Flavonol target analysis revealed that three flavonol 3-O-pentoside conjugates (f19, f25, and f29) were under the detection limit in flowers of the ugt78d3 mutant (Figure 5D). This flavonol 3-O-pentoside conjugate-deficient phenotype was complemented by the constitutive expression of UGT78D3 cDNA under the control of the cauliflower mosaic virus 35S promoter (Figure 5D), indicating that UGT78D3 encodes a flavonol 3-O-pentosyltransferase.

View larger version (20K): [in this window] [in a new window] [as a PowerPoint slide]   Figure 5. A T-DNA Insertion Mutant of UGT78D3. (A) Schematic representation of UGT78D3 with a T-DNA insertion mutant used in this work. The thick black line indicates coding sequence. Numbers indicate the position of the T-DNA insertion. The gray box adjacent to T-DNA indicates the region containing the UGT78D3 promoter from –77 bp to +4 bp. LB, left border; RB, right border. (B) Real-time PCR analysis of UGT78D3 transcripts in organs of the Arabidopsis wild type (Col-0). Error bars represent SD of three technical replicates per sample. (C) Real-time PCR analysis of UGT78D3 transcripts in wild-type (Col-0) and ugt78d3 mutant flowers. Error bars represent SD of three technical replicates per sample. (D) Extracted fragment mass chromatograms of aqueous methanol extracts from flowers of the wild type (Col-0), the ugt78d3 mutant, and ugt78d3 complemented with p35S:UGT78D3. The extracted fragment mass ions m/z 565, m/z 581, and m/z 595 correspond to compounds f19, f25, and f29 in Figure 3, respectively.

View larger version (15K): [in this window] [in a new window] [as a PowerPoint slide]   Figure 6. UPLC Analyses of the Reaction Products of UGT78D3 Recombinant Protein. Elution profiles of reaction products of GST protein (control) or GST-fused UGT78D3 protein (UGT78D3) and the standards (quercetin 3-O-arabinoside) are shown.

View larger version (27K): [in this window] [in a new window] [as a PowerPoint slide]   Figure 7. UPLC-MS Analyses of the Reaction Products of UGT78D3 or UGT89C1 Recombinant Protein. Extracted fragment mass chromatograms of aqueous methanol extracts from flowers of wild type (Col-0) ([A] and [F]), the reaction products of UGT89C1 with kaempferol 3-O-arabinoside (B), (A) coeluted with (B) (C), the reaction products of UGT78D3 with kaempferol 7-O-rhamnoside (D), (A) coeluted with (D) (E), the reaction products of UGT89C1 with quercetin 3-O-arabinoside (G), and (F) coeluted with (G) (H). The extracted fragment mass ions m/z 565 ([A] to [E]) and m/z 581 ([F] to [H]) correspond to compounds f19 and f25 in Figure 3, respectively. Kae 3-Ara, kaempferol 3-O-arabinoside; Kae 7-Rha, kaempferol 7-O-rhamnoside; Que 3-Ara, quercetin 3-O-arabinoside.

View larger version (36K): [in this window] [in a new window] [as a PowerPoint slide]   Figure 8. T-DNA Insertion Mutants of RHM1 and RHM2, and UPLC-MS Analyses of the rhm Mutant Lines. (A) Real-time PCR analysis of RHMs transcripts in Arabidopsis organs. Error bars represent SD of three technical replicates per sample. (B) Schematic representation of RHM1 and RHM2 with the mutations used in this work. rol1-1, a nonsense mutation at position 318; rol1-2, a missense mutation changing an Arg at position 283 to a Lys; rhm2-1 and rhm2-3, T-DNA insertion mutants. LB/RB, left/right borders. Numbers indicate the position of T-DNA insertion. (C) UPLC-PDA-MS analyses of the extracts from leaves and flowers in the wild type (Col-0) and the rhm mutants. Absorbance at 320 nm was used for detection. Labels correspond to compounds shown in Figure 3. The labels in parentheses indicate compounds at <5% in the peak. (D) The ratio of flavonol derivatives to total flavonol in leaves of the wild type (Col-0) and the mutants. Error bars represent SD for five independent experiments. (E) The ratio of flavonol derivatives to total flavonol in flowers of the wild type (Col-0) and the mutants. Error bars represent SD for five independent experiments. N.D., not determined.

	DISCUSSION

TOP ABSTRACT INTRODUCTION RESULTS DISCUSSION METHODS REFERENCES

 
The Power of Detailed Targeted Flavonol Analysis and Transcriptome Coexpression Analysis for Functional Genomics
Here, we used a novel strategy for functional genomics in Arabidopsis using flavonol as a case study (Figure 9 ). Coupling the detailed analysis of secondary metabolites with coexpression can be used as a more general method for identifying the biosynthetic pathways of plant compounds associated with specific organs and developmental stages as well as their genetic components. The use of T-DNA mutations in genes encoding flavonoid biosynthetic enzymes to provide comprehensive annotation of the flavonols in Arabidopsis is effective and can provide an efficient method for elucidating flavonol structures. The integration of metabolic profiling and transcriptome coexpression could efficiently narrow the pool of candidate genes. Using these methods, a flavonol 3-O-arabinosyltransferase gene was found, and the relationship between flavonols and UDP-rhamnose synthesis has emerged. After listing the candidates potentially involved in the reactions or regulation of a metabolic pathway, the specific identification of a gene function could be accomplished using a traditional reverse genetics approach with knockout mutants and/or by biochemical characterization with the recombinant proteins, as has been exemplified in this study. The above strategy will be much improved as large metabolome data sets follow transcriptome data sets into to the public domain.

View larger version (30K): [in this window] [in a new window] [as a PowerPoint slide]   Figure 9. Proposed Comprehensive Flavonol Pathway in Arabidopsis. Each number in blue corresponds to the compounds shown in Figure 3. The compounds detected in Arabidopsis are shown in a blue background. Dotted lines indicate proposed but unidentified pathways. F3AraT, flavonol 3-O-arabinosyltransferase; Fd3GlcT, flavonoid 3-O-glucosyltransferase; F3RhaT, flavonol 3-O-rhamnosyltransferase; F7GlcT, flavonol 7-O-glucosyltransferase; F7RhaT, flavonol 7-O-rhamnosyltransferase; Ara, arabinose; DeoxyHex, deoxyhexose; Glc, glucose; Hex, hexose; Rha, rhamnose.

Transcriptome Coexpression Analysis Delimits Genes Related to Flavonoid Metabolism
Transcriptome coexpression analyses suggest other candidate genes that might be related to flavonoid metabolism. Eight of 24 candidate genes were regulated in PAP1-overexpressing or in MYB11/MYB12/MYB111 triple knockout plants. CHI_put (At5g05270) may play a supplementary role to the primary CHI (At3g55120). The coordinated expression of CHI_put with CHI tends to support this hypothesis. UGT79B1 belongs to a cluster containing the anthocyanin biosynthetic group. Phylogenetic trees of flavonoid UGTs suggest that UGT79B1 encodes anthocyanin 3-O-glucoside 2''-xylosyltransferase (Tohge et al., 2005). These data suggest that CHI_put and UGT79B1 are involved in the flavonoid pathway. UGT84A1 (At4g15480) has 1-O-glucosyltransferase activity for phenylpropanoids such as p-coumaric acid and sinapic acid (Lim et al., 2001). Sinapoylmalate converted from 1-O-sinapoylglucose by malate sinapoyltransferase (SNG1), is a UV-B protective compound much as the flavonoids are (Ruegger and Chapple, 2001). Under UV stress conditions, accumulation of sinapoylmalate was higher in tt4 mutants than in the wild type (Kusano et al., 2007), suggesting that sinapoylmalate biosynthesis acts as a backup system for flavonoid metabolism under stress conditions. The positive correlation between MYB12 and UGT84A1 (r = 0.781) and the negative correlation of MYB12/UGT84A1 and SNG1 (r = –0.357/–0.323) tends to support this hypothesis. 1-O-Sinapoylglucose is also a substrate for anthocyanin sinapoyltransferase (Fraser et al., 2007), and anthocyanins accumulate with UV exposure and high sucrose stress. Arabidopsis may thus coregulate 1-O-sinapoylglucose with the flavonoids for prompt responses to stress. Cinnamoyl-CoA reductase–related protein (At2g23910) also may be involved in phenylpropanoid biosynthesis. Although coexpression does not necessarily indicate a functional relationship (Stuart et al., 2003), we believe that the significance of transcriptome coexpression analysis for functional gene identification has been established, at least for secondary metabolism.

A Flavonol Glycosyltransferase:Flavonol 3-O-Arabinosyltransferase
We demonstrated that UGT78D3 encodes a flavonol 3-O-arabinosyltransferase. The UGT78D3 function in planta was established by the identification of f19 and f25 as kaempferol/quercetin 3-O-arabinoside-7-O-rhamnoside, respectively. The conformation of substrates (kaempferol/quercetin 3-O--L-arabinosides and UDP-β-L-rhamnose) and previous data are consistent with the hypothesis that f19 and f25 are kaempferol/quercetin 3-O--L-arabinoside-7-O--L-rhamnosides, respectively (Mulinacci et al., 1995; Sakar et al., 2005). Compared with UGT78D1 and UGT78D2, UGT78D3 is weakly correlated with genes in the flavonoid metabolic pathway, which is consistent with the anatomical distribution of flavonol glycosides. The flavonol 3-O-arabinosides make up only a minor fraction of flavonols or are under the detection limit in organs except for flowers, whereas the flavonol 3-O-glucosides and 3-O-rhamnosides are distributed in all tested organs as major flavonols. Among the 6850 flavonoids registered in the Flavonoid Viewer database (http://www.metabolome.jp/software/FlavonoidViewer/viewer), 2.9% are arabinosylated flavonoids (M. Arita, personal communication). To find genes involved in such minor metabolite biosynthesis, an in-depth analysis of gene-to-metabolite correlation is essential.

Functional identification of three Arabidopsis flavonol 3-O-glycosyltransferases (3GlyTs) in the UGT78D subfamily with distinct UDP-sugar specificity allowed us to deduce the amino acid residues involved in substrate recognition. The three-dimensional structure of grape (Vitis vinifera) flavonoid 3-O-glucosyltransferase (3GlcT) was determined, and a structural analysis suggests the presence of several key residues that interact with UDP-sugar and the flavonoid backbone (Offen et al., 2006). The amino acids Gln-375, Asp-374, and Thr-141 have been proposed to interact with hydroxyl groups at the C-2 and C-3, C-3, and C-4, and C-6 positions of the glucose moiety of UDP-glucose, respectively, but Gln-375 and Thr-141 are not conserved in the Arabidopsis 3GlyTs (see Supplemental Figure 3 online). His residues corresponding to Gln-375 are conserved among the galactosyltransferases (Kubo et al., 2004). A His residue was also found as the corresponding residue in UGT78D3. The configurations of hydroxyl groups at C-2, C-3, and C-4 positions are the same in arabinose and galactose. However, UGT78D3 could not use UDP-galactose. Asp-374 in grape flavonoid 3GlcT is conserved in the three Arabidopsis 3GlyTs, although the configurations of hydroxyl groups at the C-3 and C-4 positions are different (see Supplemental Figure 3 online). Given this disparity, the residues involved in sugar donor specificity cannot be ascribed to a single amino acid residue.

RHM1 Is Involved in Flavonol Biosynthesis
Rhamnosylation is a major glycosylation of flavonols in Arabidopsis. Twenty-one of 32 detected flavonol derivatives in wild-type Arabidopsis are rhamnosides (Figure 3). A comparison of flavonol profile indicates that RHM1 plays a major role in supplying UDP-rhamnose for flavonol biosynthesis. RHM1 transcripts are accumulated as the major RHM species in most of the tested Arabidopsis organs. The expression pattern of RHM1 was well coordinated with the known flavonoid biosynthetic genes. The RHM1 correlation coefficients for flavonol rhamnosyltransferase (UGT89C1, r = 0.67; UGT78D1, r = 0.51) are higher than those for flavonol glucosyltransferase (UGT78D2, r = 0.45; UGT78D3, r = 0.32). In leaves, glycosylation at the C-3 position was strongly affected by the restriction of UDP-rhamnose supply because subsequent glycosylations at the C-7 or C-2'' positions for major flavonols are rhamnosylation only. In flowers of rol1s, additional flavonol 3-O-glucosides were detected, suggesting that more UDP-rhamnose may be necessary in flowers and glycosylation at the C-7 position is also affected in the severe limitation of UDP-rhamnose. In fact, major six flavonol contents in flowers are 10- to 70-fold higher than in leaves (Yonekura-Sakakibara et al., 2007). The role of RHM1 for UDP-rhamnose supply to flavonols was also supported by a recent report with young shoots of the rol1 mutants (Ringli et al., 2008) and our data with RHM1 knockdown mutants (SALK_027926 and SALK_143589; see Supplemental Figure 4 online).

The extent of the RMH3 contributions is unclear, and it may be committed to other pathways. RHM2/MUM4 is involved in the synthesis of pectinaceous rhamnogalacturonan I in seed mucilage (Usadel et al., 2004; Western et al., 2004). Interestingly, RHM3 has a low correlation coefficient with RHM1 (r = 0.46, all data sets version 3). The molecular evolution of RHMs may explain this correlation. The Arabidopsis genome contains a number of large duplicated chromosomal segments. Blanc et al. (2003) proposed that the Arabidopsis lineage underwent at least two distinct episodes (old and recent) of duplication. According to Blanc et al. on their website (http://wolfe.gen.tcd.ie/athal/dup) for exploring paralogous blocks identified within the Arabidopsis genome, RHM1 is the ancestral gene form. RHM3 was duplicated and diverged from RHM1 during the old duplication, and RHM2/MUM4 emerged from RHM3 by some recent duplication (see Supplemental Figure 5 online). Comprehensive analysis of rhamnose-containing metabolites using RHM knockout mutants will be helpful to clarify the additional functions of RHM1, RHM2/MUM4, and RHM3.

	METHODS

TOP ABSTRACT INTRODUCTION RESULTS DISCUSSION METHODS REFERENCES

 
Plant Materials
Arabidopsis thaliana (accession Columbia-0; Lehle Seeds) was used as the wild type in this study. The mutants ugt75c1, ugt78d1, ugt78d2, ugt89c1, rol1-1, rol1-2, and rhm2-1 were described previously (Jones et al., 2003; Usadel et al., 2004; Tohge et al., 2005; Diet et al., 2006; Yonekura-Sakakibara et al., 2007). The T-DNA–inserted knockout mutants of omt1 (CS25167) and tt7 (CS6509) were obtained from the ABRC. The tt4 (C1) mutant was obtained from tt mutant lines induced by ion beam irradiation of Arabidopsis (Shikazono et al., 2003). The T-DNA–inserted mutant of Arabidopsis, lines SALK_114099 for UGT78D3 and SALK_085051 (rhm2-3, a homozygous knockout line) for RHM2, were obtained from the Salk Institute. T-DNA insertion lines were screened by PCR using specific primers for UGT78D3, RHM2 ,and T-DNA: UGT78D3f, UGT78D3r, RHM2f, RHM2r, LBa1, and RBa1 (see Supplemental Table 1 online). PCR products were sequenced to determine the exact insertion points.

Plants were cultured on Murashige and Skoog–agar medium containing 1% sucrose (Valvekens et al., 1988) in a growth chamber at 22°C with 16 h/ 8 h light and dark cycles for 18 d or in a greenhouse at 22°C with 16 h/ 8 h light/dark for 4 weeks. Plants were harvested, immediately frozen with liquid nitrogen, and stored at –30°C until use. Five biological replicates were used for flavonoid profiling.

Chemicals
Flavonoid standards, including quercetin 3-O-arabinoside, were purchased from Extrasynthese and Analyticon. Kaempferol 3-O--L-rhamnoside-7-O--L-rhamnoside, kaempferol 3-O--L-rhamnosyl(12)-β-D-glucoside-7-O--L-rhamnoside, quercetin 3-O--L-rhamnoside-7-O--L-rhamnoside, and quercetin 3-O--L-rhamnosyl(12)-β-D-glucoside-7-O--L-rhamnoside were from our laboratory collection. UDP-β-L-arabinose and UDP--D-xylose were purchased from CarboSource Services (supported in part by National Science Foundation–Plant Cell Wall Biosynthesis Research Network Grant 0090281).

Flavonol Profiling by UPLC-PDA-ESI/Q-TOF/MS
Flavonol analyses were performed in quintuplicate. Frozen leaves were homogenized in 5 µL of extraction solvent (methanol: CH₃COOH:H₂O = 9:1:10, 0.02 mM naringenin-7-O-glucoside) per milligram of fresh weight of tissue in a mixer mill (MM300; Retsch) for 5 min at 30 Hz. After centrifugation at 12,000g, the supernatants were immediately used for flavonoid analysis.

For flavonol profiling, a Waters Acquity UPLC system (Waters) fitted with a Q-TOF Premier mass spectrometer (Micromass MS Technologies) was used. UPLC was performed on a UPLC phenyl C18 column (2.1 mm x 100 mm; Waters) at a flow rate of 0.5 mL/min at 35°C. Compounds were separated with a linear elution gradient with solvent A (0.1% formic acid in water) and solvent B (0.1% formic acid in acetonitrile) from 0 min, 100% A to 10 min, 40% B. PDA was used for the detection of UV-visible absorption in the range of 210 to 500 nm. The TOF mass analyzer was used for the detection of flavonol glycosides [M + H]⁺ and the peak of fragment ions in a positive ion scanning mode with the following setting: desolvation temperature, 400°C with a nitrogen gas flow of 600 L/h; capillary spray, 3.0 kV; source temperature, 150°C; and cone voltage of 10 V for flavonoid glycosides [M + H]⁺ or 30 V for fragment ions. For comprehensive flavonol profiling, nonprocessed MS data were converted to NetCDF format by MassLynx software (Micromass MS Technologies). Data analyses including principal component analysis were performed by the Phenomenome Profiler (Phenomenome Discoveries). Flavonol glycoside standards were used for the identification of the peaks in the plant extracts based on retention times, UV-visible absorption spectra, and mass fragmentation by tandem MS analysis. Other flavonol peaks were annotated by comparing their UV-visible absorption spectra, elution times, m/z values, and MS² fragmentation patterns with 85 reference flavonoid compounds and the reported data (Tohge et al., 2005, 2007; Routaboul et al., 2006; Yonekura-Sakakibara et al., 2007). The mass spectrum data of standard compounds (see Supplemental Data Set 2 online) were recorded in the MASSBANK Database (http://www.massbank.jp/index-e.html).

Transcriptome Coexpression Analysis
Coexpression analyses were performed using a Coexpression Gene Search algorithm on the RIKEN PRIMe website (http://prime.psc.riken.jp/; Akiyama et al., 2008) as described previously (Yonekura-Sakakibara et al., 2007; Saito et al., 2008). Twenty-four genes (listed in Table 2) were used as core query genes. The genes exhibiting positive correlations (r > 0.525) against the flavonoid biosynthetic genes were extracted using ATTED-II (Obayashi et al., 2007) based on all data sets version 3 in AtGenExpress. The coexpression graph was depicted by the Pajek program.

Evaluation of T-DNA Insertion Mutants
For complementation tests, full-length UGT78D3 cDNA was amplified by PCR using Arabidopsis flower cDNA as a template with primers UGT78D3-GWf and UGT78D3-GWr (see Supplemental Table 1 online). The Gateway system was used for construction of the binary vectors for Arabidopsis transformation. cDNA was cloned into the pCR8/GW/TOPO vector (Invitrogen) as an entry vector and sequenced to confirm the absence of PCR errors. pGWB2 was used as a destination vector. Transformation into Agrobacterium tumefaciens and Arabidopsis and the selection of transformants were performed as described previously (Yonekura-Sakakibara et al., 2007).

Arabidopsis plants were grown at 22°C under 16 h/8 h light and dark cycles in a plant growth room for 5 weeks. The leaves and flowers of plants were harvested, immediately frozen with liquid nitrogen, and stored at –30°C until use for metabolite profiles.

Quantitative Real-Time PCR
RNA Extraction and cDNA synthesis were performed as described previously (Yonekura-Sakakibara et al., 2004). The developmental stage of each organ used for analysis was as described previously (Yonekura-Sakakibara et al., 2007). Accumulation levels of UGT78D3, RHM1, RHM2, and RHM3 transcripts were analyzed by real-time PCR with an ABI PRISM 7500 real-time PCR system (Applied Biosystems), monitoring amplification with SYBR-Green I dye (Applied Biosystems). The primers UGT78D3-RTf, UGT78D3-RTr, RHM1-RTf, RHM1-RTr, RHM2-RTf, RHM2-RTr, RHM3-RTf, and RHM3-RTr (see Supplemental Table 1 online) were designed using Primer Express software (Applied Biosystems) and checked for specific product formation by a dissociation program. In each case, plasmid DNA containing the corresponding gene was used as a template to generate a calibration curve. Real-time PCR was performed in triplicate on a single biological sample for each genotype.

Production of Recombinant UGT78D3 Protein and Assay of Glycosyltransferase
Full-length UGT78D3 cDNA was amplified and cloned into pCR8/GW/TOPO (Invitrogen) as described above. To construct the protein expression vector, the resultant plasmid was amplified by PCR with the primers UGT78D3-pET41bf and UGT78D3-pET41br (see Supplemental Table 1 online), the PCR product was digested with StuI and XhoI, and ligated with StuI-XhoI-digested pET-41b(+) vector (Novagen). The sequence of the resultant plasmid, pKYS342, was confirmed. Escherichia coli strain BL21star(DE3) was used as a host for expression. Transformed cells were cultivated at 37°C until A₆₀₀ reached 0.5. After the addition of isopropyl-β-D-thiogalactopyranoside to a final concentration of 1 mM, cells were cultured at 25°C for 4 h. The cells were collected, and the protein was purified as a GST fusion according to the manufacturer's instructions (Qiagen).

The standard enzyme assay reaction mixture (final volume of 50 µL) consisted of 50 mM HEPES-KOH, pH 7.5, 150 µM flavonoid substrate, and 500 µM UDP-sugar. The mixture was preincubated at 30°C for 2 min, and the reaction was started by the addition of enzyme. Reactions were stopped after 0, 4, 8, 12, or 30 min of incubation at 30°C by the addition of 50 µL ice-cold 0.5% (v/v) trifluoroacetic acid/methanol for flavonols or 50 µL ice-cold 1.0% (v/v) HCl/methanol for anthocyanidins and anthocyanins, and the supernatant was recovered by centrifugation at 12,000g for 3 min. Flavonoids in the resultant solution were analyzed as described above.

Preparation and NMR Analysis of Kaempferol 3-O-Arabinoside-7-O-Rhamnoside and Quercetin 3-O-Arabinoside-7-O-Rhamnoside
A reaction mixture of UGT89C1 with kaempferol 3-O--L-arabinoside and UDP-β-L-rhamnose as substrates (total 300 mL) was concentrated and chromatographed on COSMOSIL 75C₁₈-OPN (Nacalai Tesque) and stepwise eluted with water, 10% aqueous CH₃CN, 30% aqueous CH₃CN, 50% aqueous CH₃CN, 70% aqueous CH₃CN, and 100% methanol. Kaempferol glycosides in the fractions of 30% aqueous CH₃CN, 50% aqueous CH₃CN, 70% aqueous CH₃CN, and 100% methanol were further separated by preparative HPLC using a C18 reversed-phase column (10 mm x150 mm I.D., particle size 5 µm, Inertsil ODS-EP; GL Science) using 30% aqueous CH₃CN with flow rates of 2.0 mL/min at 30°C. Kaempferol 3-O-L-arabinoside-7-O-L-rhamnoside (1.8 mg) was obtained and the structure was supported by ¹H-NMR (600 MHz, methanol-d₄): 5.16 (1H, d, J = 6.6 Hz, arabinose H-1), 5.52 (1H, br. s, rhamnose H-1), 6.43 (1H, d, J = 1.8 Hz, H-6), 6.73 (1H, d, J = 1.2 Hz, H-8), 6.86 (2H, d, J = 9.0 Hz, H-3'and H-5'), 8.06 (2H, d, J = 8.4 Hz, H-2' and H-6'), ROE correlation between the anomeric proton of rhamnose (rhamnose H-1, 5.52), and the protons at C-6 (H-6, 6.43) and C-8 position (H-8, 6.73), ESI-TOF-MS (positive ion mode): m/z 565.1554 [M+H]⁺, calculated for C₂₆H₂₉O₁₄, 565.1557, ESI-TOF-MS/MS (positive ion mode, collision energy 30 eV): m/z (relative intensity) 433.1145 [M+H-arabinose]⁺ (39), 287.0605 [M+H-arabinose-rhamnose]⁺ (100), and UV _max nm: 345 (on flow).

A large-scale enzymatic product of UGT89C1 with quercetin 3-O--L-arabinoside and UDP-β-L-rhamnose was isolated by a similar procedure except for the mobile phase (20% aqueous CH₃CN) of preparative HPLC step. Quercetin 3-O-L-arabinoside-7-O-L-rhamnoside (2.1 mg) was obtained and the structure was supported by ¹H-NMR (600 MHz, methanol-d₄): 5.22 (1H, d, J = 6.0 Hz, arabinose H-1), 5.56 (1H, br. s, rhamnose H-1), 6.46 (1H, br. s, H-6), 6.75 (1H, br. s, H-8), 6.88 (1H, d, J = 7.8 Hz, H-5'), 7.61 (1H, br. d, J = 8.4 Hz, H-6'), 7.76 (1H, br. s, H-2'). ROE (irradiation 5.56 [rhamnose H-1]): H-6 ( 6.46) and H-8 ( 6.75), ROE correlation between the anomeric proton of rhamnose (rhamnose H-1, 5.52) and the protons at C-6 (H-6, 6.43) and C-8 position (H-8, 6.73), ESI-TOF-MS (positive ion mode): m/z 581.1496 [M+H]⁺, calculated for C₂₆H₂₉O₁₅, 581.1506. ESI-TOF-MS/MS (positive ion mode, collision energy 30 eV): m/z (relative intensity) 449.1086 [M+H-arabinose]⁺ (41), 303.0554 [M+H-arabinose-rhamnose]⁺ (100), and UV _max nm:340 (on flow). Monitoring was accomplished by PDA. ¹H-NMR spectra were recorded on a JNM ECP-600 spectrometer (JEOL).

Accession Numbers
Sequence data from this article can be found in the Arabidopsis Genome Initiative database under accession numbers TT4 (At5g13930), TT7 (At5g07990), UGT75C1 (At4g14090), UGT78D1 (At1g30530), UGT78D2 (At5g17050), UGT78D3 (At5g17030), UGT89C1 (At1g06000), OMT1 (At5g54160), RHM1 (At1g78570), RHM2 (At1g53500), and RHM3 (At3g14790).

Supplemental Data
The following materials are available in the online version of this article.

Supplemental Figure 1. Venn Diagrams of MS Peaks Detected in Various Organs of Col-0 and tt4.

Supplemental Figure 2. Coexpression Relationships of Genes Involved in the Flavonoid Pathway in Arabidopsis.

Supplemental Figure 3. Multiple Alignment of the Deduced Amino Acid Sequences of Flavonoid 3-O-Glucosyltransferases from Arabidopsis and Grape (Vitis vinifera).

Supplemental Figure 4. T-DNA Insertion Mutants of RHM1 and UPLC/MS Analyses of the rhm1 Mutant Lines.

Supplemental Figure 5. Schematic Representation of Duplicated Regions Containing RHMs in the Arabidopsis Genome.

Supplemental Table 1. Primers Used in This Study.

Supplemental Data Set 1. MS-Detectable Peaks in the Wild Type and tt4.

Supplemental Data Set 2. Characteristics of Standard Compounds.

	Acknowledgments

 
We thank C. Ringli (University of Zurich) for the rol1-1 and rol1-2 mutants, M. Pauly (Michigan State University) for the rhm2-1 mutant, S. Kitamura (Japan Atomic Energy Research Institute) for the tt4 mutant, T. Nakagawa (Shimane University) for the pGWB2 vector, Y. Sekiyama (RIKEN) for helpful comments on NMR analysis, T. Obayashi (University of Tokyo), H. Suzuki, and N. Sakurai (Kazusa DNA Research Institute) for advice on coexpression analysis, and K. Akiyama (RIKEN) for excellent technical support for using the RIKEN-PRIMe website. This study was supported in part by grants-in-aid from the Ministry of Education, Culture, Sports, Science, and Technology of Japan and the Uehara Memorial Foundation, Japan.

	Footnotes

 
¹ These authors contributed equally to this work.

² Current address: Max Planck Institute of Molecular Plant Physiology, Am Mühlenberg 1, 14476 Potsdam Golm, Germany.

The author responsible for distribution of materials integral to the findings presented in this article in accordance with the policy described in the Instructions for Authors (www.plantcell.org) is: Kazuki Saito (ksaito{at}faculty.chiba-u.jp).

^[W] Online version contains Web-only data.

^[OA] Open Access articles can be viewed online without a subscription.

www.plantcell.org/cgi/doi/10.1105/tpc.108.058040

Received January 15, 2008; Revision received July 8, 2008. accepted August 10, 2008.

	REFERENCES

TOP ABSTRACT INTRODUCTION RESULTS DISCUSSION METHODS REFERENCES

 
Akiyama, K., Chikayama, E., Yuasa, H., Shimada, Y., Tohge, T., Shinozaki, K., Hirai, M.Y., Sakurai, T., Kikuchi, J., and Saito, K. (2008). PRIMe: A Web site that assembles tools for metabolomics and transcriptomics. In Silico Biol. 8: 0027.

Andersen, Ø.M., and Markham, K.R. (2006). Flavonoids: Chemistry, Biochemistry, and Applications. (Boca Raton, FL: CRC Taylor & Francis).

Aoki, K., Ogata, Y., and Shibata, D. (2007). Approaches for extracting practical information from gene co-expression networks in plant biology. Plant Cell Physiol. 48: 381–390.[Abstract/Free Full Text]

Blanc, G., Hokamp, K., and Wolfe, K.H. (2003). A recent polyploidy superimposed on older large-scale duplications in the Arabidopsis genome. Genome Res. 13: 137–144.[Abstract/Free Full Text]

Bloor, S.J., and Abrahams, S. (2002). The structure of the major anthocyanin in Arabidopsis thaliana. Phytochemistry 59: 343–346.[CrossRef][ISI][Medline]

Brown, B.A., Cloix, C., Jiang, G.H., Kaiserli, E., Herzyk, P., Kliebenstein, D.J., and Jenkins, G.I. (2005). A UV-B-specific signaling component orchestrates plant UV protection. Proc. Natl. Acad. Sci. USA 102: 18225–18230.[Abstract/Free Full Text]

D'Auria, J.C., and Gershenzon, J. (2005). The secondary metabolism of Arabidopsis thaliana: Growing like a weed. Curr. Opin. Plant Biol. 8: 308–316.[CrossRef][ISI][Medline]

Diet, A., Link, B., Seifert, G.J., Schellenberg, B., Wagner, U., Pauly, M., Reiter, W.D., and Ringli, C. (2006). The Arabidopsis root hair cell wall formation mutant lrx1 is suppressed by mutations in the RHM1 gene encoding a UDP-L-rhamnose synthase. Plant Cell 18: 1630–1641.[Abstract/Free Full Text]

Dixon, R.A., and Paiva, N.L. (1995). Stress-induced phenylpropanoid metabolism. Plant Cell 7: 1085–1097.[CrossRef][ISI][Medline]

Dixon, R.A., and Strack, D. (2003). Phytochemistry meets genome analysis, and beyond. Phytochemistry 62: 815–816.[CrossRef][ISI][Medline]

Dooner, H.K., Robbins, T.P., and Jorgensen, R.A. (1991). Genetic and developmental control of anthocyanin biosynthesis. Annu. Rev. Genet. 25: 173–199.[CrossRef][ISI][Medline]

Fraser, C.M., Thompson, M.G., Shirley, A.M., Ralph, J., Schoenherr, J.A., Sinlapadech, T., Hall, M.C., and Chapple, C. (2007). Related Arabidopsis serine carboxypeptidase-like sinapoylglucose acyltransferases display distinct but overlapping substrate specificities. Plant Physiol. 144: 1986–1999.[Abstract/Free Full Text]

Graham, T.L. (1998). Flavonoid and flavonol glycoside metabolism in Arabidopsis. Plant Physiol. Biochem. 36: 135–144.[CrossRef][ISI]

Harborne, J.B., Baxter, H., and Moss, G.P. (1999). Phytochemical Dictionary: A Handbook of Bioactive Compounds from Plants, 2nd ed. (Londong: Taylor & Francis).

Hirai, M.Y., et al. (2005). Elucidation of gene-to-gene and metabolite-to-gene networks in Arabidopsis by integration of metabolomics and transcriptomics. J. Biol. Chem. 280: 25590–25595.[Abstract/Free Full Text]

Hirai, M.Y., et al. (2007). Omics-based identification of Arabidopsis Myb transcription factors regulating aliphatic glucosinolate biosynthesis. Proc. Natl. Acad. Sci. USA 104: 6478–6483.[Abstract/Free Full Text]

Jones, P., Messner, B., Nakajima, J., Schaffner, A.R., and Saito, K. (2003). UGT73C6 and UGT78D1, glycosyltransferases involved in flavonol glycoside biosynthesis in Arabidopsis thaliana. J. Biol. Chem. 278: 43910–43918.[Abstract/Free Full Text]

Kerhoas, L., Aouak, D., Cingoz, A., Routaboul, J.M., Lepiniec, L., Einhorn, J., and Birlirakis, N. (2006). Structural characterization of the major flavonoid glycosides from Arabidopsis thaliana seeds. J. Agric. Food Chem. 54: 6603–6612.[CrossRef][ISI][Medline]

Kleine, T., Kindgren, P., Benedict, C., Hendrickson, L., and Strand, A. (2007). Genome-wide gene expression analysis reveals a critical role for CRYPTOCHROME1 in the response of Arabidopsis to high irradiance. Plant Physiol. 144: 1391–1406.[Abstract/Free Full Text]

Kubo, A., Arai, Y., Nagashima, S., and Yoshikawa, T. (2004). Alteration of sugar donor specificities of plant glycosyltransferases by a single point mutation. Arch. Biochem. Biophys. 429: 198–203.[CrossRef][ISI][Medline]

Kusano, M., Fukushima, A., Arita, M., Jonsson, P., Moritz, T., Kobayashi, M., Hayashi, N., Tohge, T., and Saito, K. (2007). Unbiased characterization of genotype-depedent metabolic regulations by metabolomic