|
|
||||||||
|
First published online June 3, 2005; 10.1105/tpc.105.030767 © 2005 American Society of Plant Biologists
CINNAMYL ALCOHOL DEHYDROGENASE-C and -D Are the Primary Genes Involved in Lignin Biosynthesis in the Floral Stem of Arabidopsis
|
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Cinnamyl alcohol dehydrogenase (CAD) is a specialized enzyme involved in the reduction of cinnamaldehydes into cinnamyl alcohols, the last step of monolignol biosynthesis before oxidative polymerization in the cell wall. CAD genes and their proteins in different taxa display distinct features depending on their phylogenetic origin. Duplication of CAD genes has only been observed in angiosperms (Knight et al., 1992
; Brill et al., 1999
). In gymnosperms, CAD genes are reported to be monogenic despite the fact that different alleles could be detected in some species (Galliano et al., 1993b
; MacKay et al., 1995
). In addition, conifer CAD proteins are believed to be highly specific for the reduction of coniferaldehyde (Kutsuki et al., 1982
; Galliano et al., 1993a
), whereas angiosperm CAD proteins have been shown to have a significant affinity for both coniferaldehyde and sinapaldehyde (Mansell et al., 1974
; Grima-Pettenati et al., 1994
; Hawkins and Boudet, 1994
; Brill et al., 1999
). These characteristics, combined with the fact that angiosperm lignin displays mainly both G and S units, whereas lignin from most gymnosperms contains predominantly G units, has led some authors to hypothesize that specialized CAD proteins could be specifically involved in either sinapyl or coniferyl alcohol biosynthesis (Hawkins and Boudet, 1994
). This hypothesis has been supported by the identification of a novel enzyme in aspen (Populus tremuloides), sinapyl alcohol dehydrogenase (SAD), which was shown to be specifically involved in the reduction of sinapaldehyde in vitro (Li et al., 2001
). These authors suggested that all the previously characterized CAD proteins were involved only in the reduction of coniferaldehyde. Together with aldehyde-5-hydroxylase and aldehyde 5-O-methyl transferase, SAD is proposed to be the third enzyme involved in sinapyl alcohol biosynthesis. The biological role of SAD could be of practical importance because it may regulate the channeling of S units during lignification. Indeed, modulation of the S:G ratio in transgenic lines could potentially lead to advantageous alterations in wood properties, such as reduced resistance to kraft delignification (reviewed in Baucher et al., 1998
, 2003
; Mellerowicz et al., 2001
).
Approaches using mutants and transgenic plants have previously been shown to be of considerable use for studying the impact of metabolic disruption upon secondary cell wall structure. For instance, the cell wall properties of the naturally CAD-deficient mutants brown midrib1 (bm1) of maize (Zea mays) and cad-n1 of loblolly pine (Pinus taeda) have been intensively studied (Lapierre et al., 2000
; Marita et al., 2003
). bm1 showed a 20% decrease in lignin content with no alteration in the S:G ratio, although lignin-associated aldehydes were detected (Halpin et al., 1998
). The loblolly pine CAD mutant was shown to have a very slight reduction in lignin content (<5%) but to also contain an unusually high content of lignin-associated coniferaldehyde and dihydroconiferyl alcohols. The production of downregulated CAD plants via antisense strategies has been performed in various species, including alfalfa (Medicago sativa), poplar (Populus spp), tobacco (Nicotiana tabacum), and, most recently, fescue (Festuca arundinacea) (Halpin et al., 1994
; Baucher et al., 1996
, 1999
; Chen et al., 2003
). Lignin analysis of these transgenic plants led to controversial interpretations of the potential role of CAD (Li et al., 2001
; Anterola and Lewis, 2002
). Indeed, lignin content was found to be impacted only modestly, if at all. Moreover, even though the S:G ratio was found reduced in some studies, there was a general lack of consistency (Halpin et al., 1994
; Higuchi et al., 1994
; Baucher et al., 1996
; Stewart et al., 1997
; Chen et al., 2003
), a result that could be related to incomplete disruption of CAD expression. In addition, the potential for functional redundancy, as for example produced by duplicate CAD genes, could lead to further ineffectiveness of antisense strategies. Transcriptomic tools such as large-scale microarrays in lignin-altered plants could help confirm or reject such hypotheses by revealing the extent of transcriptional modification in these transgenic plants.
A better understanding of the complexity of gene families involved in lignification has been provided by the complete genome sequencing and annotation of rice (Oryza sativa) and Arabidopsis thaliana (Sasaki and Sederoff, 2003
). For example, Tobias and Chowk (2005)
have recently shown that there are 12 CAD genes in rice, consistent with duplication of many of the genes coding for enzymes involved in the phenylpropanoid pathway of Arabidopsis, although the precise role of these duplicates in lignin metabolism remains to be elucidated (Costa et al., 2003
; Goujon et al., 2003
; Raes et al., 2003
). Previously, we surveyed the complete CAD gene family in Arabidopsis (Sibout et al., 2003
), which led to the identification of nine CAD-like proteins distributed into four different classes based on their amino acid similarity. CAD-C (At3g19450) and CAD-D (At4g34230) belong to one class that is highly similar to other well-characterized angiosperm CAD proteins (eucalyptus [Eucalyptus gunnii] CAD2, alfalfa CAD2, aspen CAD, and tobacco CAD14) (reviewed in Baucher et al., 2003
). The second class contains CAD-A (At4g34970), CAD-B1 (At4g37980), and CAD-B2 (At4g37990), which are most closely related to the poplar SAD identified by Li et al. (2001)
. The third class contains CAD1 (At4g34930), CAD-E (At2g21730), and CAD-F (At2g21890), which are most closely related to the alfalfa CAD2 identified by Brill et al. (1999)
. The fourth class contains CAD-G (At1g72680), which does not have well-identified homologs. Importantly, two mutants (cad-c and cad-d) were also characterized with decreased CAD activities, although only cad-d produced a modest reduction in lignin content along with a substantial decrease in conventional S lignin (Sibout et al., 2003
).
In this study, we show that the double cad-c cad-d mutant in Arabidopsis produces the strongest phenotype achieved to date in any CAD-disrupted dicotyledonous angiosperm. This mutant possesses a highly condensed lignin enriched in aldehydes. Thioacidolysis results also show that conventional H, G, and S subunits are drastically altered. Using Affymetrix Arabidopsis 24K GeneChip and real-time PCR, we showed that this double mutation alters the expression of several genes, including CAD genes. Moreover, Fourier transform infrared (FTIR) spectroscopy analysis clearly showed that the metabolic impact of CAD mutations affects both fiber and xylem cell wall composition. Finally, the complementation of the double mutant by various CAD genes from different species suggests different abilities of these genes/proteins to produce syringyl-lignin moieties but does not indicate the requirement for any specific SAD gene.
| RESULTS |
|---|
|
|
|---|
|
|
All these observations suggested that the lignin of the double mutant is reduced within various organs and tissues and that the S unit content is reduced below that detectable by Mäule staining. It is also evident that the polysaccharides of the cell wall fibers in the double mutant were more accessible to Safranin/Alcian Blue staining analysis as compared with the wild type, which may be due to lower lignin content. These features prompted us to perform additional analyses to better define the lignin characteristics of the double mutant.
Cell Wall Modifications of the cad-c cad-d Mutant
The previous histological analysis showed that interfascicular fiber elements and xylem bundles can be distinguished easily within transverse sections, despite the small size of Arabidopsis stems. These tissues are of interest because lignin composition is substantially different: xylem cell walls are mainly lignified with G units, whereas fiber lignin is composed of both G and S units. We performed localized FTIR analysis on these specialized tissues and analyzed spectra modifications within the range of 1790 to 839 cm1 in both wild-type and double mutant plants.
First, we compared the absorbance of xylem and fiber tissues from wild-type plants. Although the spectra from the xylem and fibers in wild-type plants were similar, the magnitude of absorbance was often different at specific wavelengths (Figure 3A). For example, the difference in peak height at 1743 cm1 may be due to the relative concentration of hemicellulose polymer in the fibers as compared with the xylem (Figure 3A). Owen and Thomas (1989)
suggested that strong absorption at that wavelength is likely produced by the stretching of the free carbonyl group of the hemicellulose polymer. Similarly, variation in absorbance in the range of 1190 to 850 cm1 is characteristic of changes in polysaccharide structure or content such as cellulose (Mouille et al., 2003
).
|
In a previous study, cad-c was not found to have a reduced Klason lignin content, whereas cad-d produced a very modest reduction in lignin content compared with wild-type plants (Sibout et al., 2003
). By contrast, the double mutant showed a drastic reduction (40%) in Klason lignin content (Table 1). We also quantified the susceptibility of the double mutant to cellulase using a method adapted from Rexen (1977)
. We observed a 57% dry mass loss in the double mutant after cellulase treatment but only a 29% dry mass loss in the wild type. The H, G, and S lignin-derived units released by thioacidolysis were also analyzed by gas chromatographymass spectrometry (GC-MS). When calculated on the basis of the Klason lignin content, the yield of lignin-derived monomers was drastically reduced in the mutant (Table 1). This result clearly shows a severe reduction in conventional units specifically involved in ß-O-4 linkage and suggests that the lignin is considerably enriched in condensed bonds that resist the thioacidolysis procedure. Besides the conventional H, G, and S thioacidolysis monomers, we found relatively large amounts of indene derivatives (Table 1). These lignin-derived monomers specifically originate from ß-O-4linked sinapaldehyde and coniferaldehyde units, as shown in recent articles (Kim et al., 2002
; Lapierre et al., 2004
). In addition, the dithioketal derivatives of coniferaldehyde, vanillin, and syringaldehyde were also found to be released in relatively larger amounts from the thioacidolysis of the extract-free double mutant sample (data not shown), suggesting that these aldehydes were incorporated in the mutant lignin by ether bonds at C4OH. The reduced thioacidolysis S:G ratio of the conventional S and G units (0.38 in the wild type versus 0.08 in the mutant, Table 1) suggests that CAD mutations more specifically affect the formation of conventional S units than that of G units. In addition to the determination of the lignin-derived monomers, the amount of ferulic and sinapic acids released by thioacidolysis was found to be higher in the mutant (Table 1).
|
|
24,000 genes). Our first goal was to use the technology as a tool to explore gene expression profiles in the double mutant transcriptome. In particular, we investigated in detail the expression profile of CAD homologs as well as other genes upstream in the phenylpropanoid pathway. Four biological repetitions of RNA extracted from the bases of stems (which share differentiated interfascicular fibers and xylem) of wild-type and mutant plants were used for the microarray analysis. The complete array data for each sample is presented in Supplemental Tables 1 to 8 online. A statistical treatment of eight ATH1 GeneChips (unpaired Student's t test, P = 0.05) using Affymetrix Microarray Suite (MAS 5.0) software detected 5789 genes in the double mutant that had significantly different expression compared with the wild type (see Supplemental Table 9 online). Two screen parameters were chosen to avoid little or unreliable evidence of differential expression: probe set variations between mutant and wild-type hybridization were tagged "significant" when the change in raw signal was higher than 200 units with a minimum twofold change. A preliminary analysis was done to validate the stringency of our screen parameters and to reveal the reproducibility between ATH1 GeneChip hybridization experiments. Therefore, scatterplots of the average signal intensity of each probe set from two control ATH1 GeneChips were compared with those of two other control ATH1 GeneChips. Different combinations were done with all control or all cad-c cad-d ATH1 DNA chips (see Supplemental Figure 1 online). After application of our screen parameters to the 5789 genes with altered expression by computing all ATH1 GeneChip data (four wild-type versus four mutant ATH1 GeneChips), we found 546 genes with altered expression in the double mutant. Means of signals with standard errors, fold changes, and the Student's t test P value of each probe set are presented in Supplemental Table 10 online.
A few small differences in signal were observed for the genes involved in monolignol biosynthesis, from Phe to the cinnamaldehydes (Table 3). For example, signals corresponding to Phenylalanine ammonia-lyase 3 (PAL3; At5g04230) and 4-Coumarate-CoA ligase 3 (4CL3; At1g65060) suggest that these genes are overexpressed 2-fold (Student's t test, P = 0.0087) and 2.6-fold (Student's t test, P = 0.0092), respectively, relative to the wild type. Furthermore, cinnamyl CoA reductase 1 (CCR1; At1g15950) and ferulate 5-hydroxylase (F5H; At4g36220) seem to be underexpressed 1.7-fold (Student's t test, P = 0.0004) and 1.8-fold (Student's t test, P = 0.014), respectively, relative to the wild type. However, low hybridization signal (<200) and/or weak change in expression (less than twofold) were not relevant to our screen parameters and were consequently not considered in the analysis. Contrastingly, transcripts of CAD-G were found to be 5.5 times more abundant in the double mutant. In addition, five laccase-like and two peroxidase-like genes were found to be downregulated in cad-c cad-d. Several other genes encoding proteins known to be involved in cell wall biosynthesis also had altered expression, including arabinogalactans, pectin esterases, proline-rich protein (PRP), and glycine-rich protein (GRP) (Table 3). It is noteworthy that among several glucosyl transferases with altered expression in the mutant (see Supplemental Table 10 online), the IRX3 (At5g17420) and IRX5 (At5g44030) genes that encode cellulose synthase subunits showed at least fourfold lower expression than the wild type (Table 3). Nevertheless, we did not observe an alteration in the transcript levels of the UDP glycosyl transferase genes (see Supplemental Table 2 online), which are specifically involved in the glycosylation of sinapyl and coniferyl alcohol or related acids (Lim et al., 2001
). Interestingly, distinct increases in transcript abundance were observed for 10 genes encoding glutathione S-transferase (GST), some being increased >10-fold (At2g29480), 29-fold (At2g29490), and 40-fold (At1g17170) (see Supplemental Table 10 online). The 1-aminocyclopropane-1-carboxylic acid oxidase, which is involved in ethylene biosynthesis, exhibited a twofold increase in transcript abundance in the double mutant relative to the wild type. This enzyme is known to play a role in wood formation and also in stress response (for review, see Andersson-Gunnerås et al., 2003
). The double mutant also showed alterations in the expression level of several transcription factors (WRKY, MYB, HD-ZIP, BZIP, and MADS box). We did not detect, however, any significant alteration in the expression of PAP1, SHATTERPROOF1 or 2, BREVIPEDICELLUS, or MYB61, all of which have been suggested to be involved in lignin biosynthesis (Borevitz et al., 2000
; Liljegren et al., 2000
; Mele et al., 2003
; Newman et al., 2004
).
|
|
Complementation of cad-c cad-d with CAD Genes from Forest Trees
An obvious way to demonstrate that the physiological features observed in the double mutant were really due to the tagged mutations in both CAD genes was to rescue the double mutant with CAD-D and CAD-C expression. In addition, to broaden our study in the context of recently published data, we attempted to complement the double mutant with three CADs known to be phylogenetically distinct and supposedly displaying different substrate affinities. The first gene was CAD (Z19568) from Populus deltoides, which has previously been used in antisense strategies in poplar (Baucher et al., 1996
). P. deltoides CAD is 97.5% identical at the amino acid level to P. tremuloides CAD, the isoform suggested to be involved in the production of guaiacyl units in aspen by Li et al. (2001)
. The second gene was CAD (AJ868574) from spruce (Picea abies), and the third, SAD1 (AY830131), was cloned in our laboratory from Populus tremula x tremuloides. P. tremula x tremuloides SAD1 is 99.2% identical at the amino acid level to P. tremuloides SAD, a protein found to have high affinity with sinapaldehyde in vitro (Li et al., 2001
). Both a constitutive expression system (35S Cauliflower mosaic virus promoter) and a tissue-specific expression system using the Arabidopsis CAD-D promoter were used to perform complementation experiments. Because the results obtained using both approaches were similar, we present here only the lignin features of the complemented cad-c cad-d lines constructed using the Arabidopsis CAD-D promoter.
First, we selected 16 transgenic lines per type of molecular construct expressing Arabidopsis CAD-C, Arabidopsis CAD-D, P. deltoides CAD, or P. abies CAD and named these chimeric lines ChimAtCAD-C, ChimAtCAD-D, ChimPdCAD, and ChimPaCAD, respectively. Most of the 16 lines for each construct reacted to Mäule staining, allowing positive correlation with the disappearance of the natural typical reddish color observed in the double mutant (Figure 5A). Moreover, typical phenotypes observed in cad-c cad-d, such as bending of stems, were not present in those lines (data not shown). Among 25 ChimPttSAD1 lines transformed with the P. tremula x tremuloides SAD1 gene, only six reacted very weakly to Mäule staining. Moreover, all ChimPttSAD1 lines still showed a reddish-brown color in lignified tissues typical of the lignin-deficient mutant (Figure 5A).
|
Except for the ChimPttSAD1 lines, most of the constructs were able to restore the G and S yield to levels similar to the wild type (Figure 5B). Nevertheless, the S:G ratios were not equally restored between constructs (Figure 5C). ChimAtCAD-D and ChimPdCAD lines showed a complete recovery of the proportion of S and G thioacidolysis main monomers. The S:G ratios were found similar to those of the wild-type plants (0.34 and 0.36 for each wild-type line). In addition, the unusually high level of indene monomers released by thioacidolysis of the double mutant lignin was reduced drastically in those lines, although still detectable. ChimAtCAD-C lines showed S:G ratios of 0.27 and 0.29, consistent with the S:G ratio observed in the Arabidopsis cad-d mutant (Sibout et al., 2003
) deficient in CAD-D but not CAD-C transcripts. In agreement with this finding, indenes corresponding to coniferaldehydes were found significantly lower in ChimAtCAD-C than indenes corresponding to sinapaldehydes that were detected at a level comparable to that of the mutant. Similar results were observed with ChimPaCAD; however, the S:G ratio was further reduced (0.20 and 0.21 for each line). Total yield of conventional thioacidolysis monomers in the ChimPttSAD1 lines was found increased relative to the double mutant but still 80% lower than the wild type. Surprisingly, the S:G ratio was restored to levels close to that of the wild type (0.31 and 0.34), suggesting that these lines could produce S units (Figure 5C). Nevertheless, yield of syringyl units was still 82% lower than that of the wild type, and indene monomers were still present in those lines at a level comparable to that of the mutant. P. tremula x tremuloides SAD1 mRNAs were detected in those lines (data not shown), and a protein blot analysis using a specific antibody allowed the detection of P. tremula x tremuloides SAD1 (see Supplemental Figure 2 online).
Thioacidolysis analysis of four other ChimpttSAD1 lines found positive in the Mäule test revealed similar results (data not shown). Finally, results obtained with P. tremula x tremuloides SAD1 were later confirmed with transgenic lines expressing the P. tremuloides SAD gene kindly provided by V. Chiang and L. Li (data not shown).
| DISCUSSION |
|---|
|
|
|---|
The cad-c cad-d Double Mutation Affects Plant Growth and Xylem/Fiber Cell Walls
The cad-c cad-d mutant has a delay in both vegetative growth and bolting, combined with reduced size. Similar growth characteristics have previously been observed in irx (irregular xylem) mutants (Turner and Somerville, 1997
), such as irx3, which is mutated in a cellulose synthase gene, or irx4, which has incorrect splicing of CCR transcripts (Jones et al., 2001
) and a strong decrease in lignin content. In addition to altered growth characteristics, irx mutants are distinguished by stems that tend to recline, a trait also observed with the cad-c cad-d double mutant. Moreover, xylem elements in the double mutant were found to be distorted and collapsed, whereas interfascicular fibers were not, features that were also observed in irx mutants. Consequently, it is likely that the collapse of the functioning xylem elements in the cad-c cad-d double mutant is also due to the negative pressure generated by water in this tissue, as was deduced for irx by Turner and Somerville (1997)
. This alteration in the conducting vessels could impact water and nutrient movements to the upper part of the stem and reduce the stem size of the double mutant. This trait might make the double mutant susceptible to environmental conditions. Another possible explanation for the delay in growth and shorter stems in this mutant could be the high content of accumulated soluble aldehydes that have been associated with the coloration of lignin. Such coloration has often been observed in transgenic or mutant plants with altered expression of caffeic acid O-methyl transferase, caffeoyl CoA O-methyl transferase, CCR, or CAD (reviewed in Baucher et al., 2003
). Cinnamaldehydes have been shown to be growth inhibitors of bacterial and fungal strains (Rutten and Gocke, 1988
; Utama et al., 2002
), with recent studies reporting cytostatic properties of cinnamaldehyde derivatives by inducing mitotic arrest in mammalian cells (Jeong et al., 2003
). Brown et al. (2001)
have previously suggested some phenolics to be endogenous modulators of auxin transport. Similar impacts of aldehydes on growth of the double mutant should not be excluded.
Fibers also play an important role in the upright growth habit of Arabidopsis stems, as shown by studies of the ifl1 mutant, in which fibers are lacking (Zhong et al., 1997
). The strong tendency of cad-c cad-d stems to collapse further suggested that physical properties of the fibers could be altered. Whereas most CAD genes and/or proteins have been detected in both the xylem and sclerenchyma cells of different species (Feuillet et al., 1995
; Hawkins et al., 1997
), Li et al. (2001)
showed, using a stringent in situ hybridization technique, that P. tremuloides SAD expression is restricted to the fibers and colocalized with G-S units but not with the xylem elements that display only G units. Because we suspected cad-c cad-d to be deficient not only in G, but also in G-S units, we took advantage of FTIR spectroscopy to determine whether both the xylem and fibers were affected within the double mutant without compromising the integrity of these tissues. Xylem bundles can be easily distinguished from the interfascicular fibers within an Arabidopsis stem, facilitating the spectral properties of the tissues of the mutant to be compared with that of the wild type. Statistical analysis of FTIR spectra showed unambiguously that cell walls of both lignified tissues were impacted in the cad-c cad-d stem. It is important to note that the results of FTIR analysis were strongly consistent with those of the thioacidolysis analysis, allowing FTIR to be further developed to detect Arabidopsis mutants differentially altered in fibers and/or vessels. The impact of the CAD deficiency on both fiber and xylem tissues is in accordance with the expression patterns of Arabidopsis CAD-C and Arabidopsis CAD-D (Sibout et al., 2003
). It is also interesting to note that the impact of the mutations is different between the two tissue types. This distinction is likely due to differences in cell wall composition between the fibers and xylem elements as suggested in Figure 3A. Limited knowledge exists on the chemical composition of the cell walls of Arabidopsis fiber and xylem tissues, and new techniques of dissections, such as laser capture coupled with chemical cell wall analysis, could help considerably to increase our understanding.
The inability of not only the fibers but also the xylem cell walls of the double mutant to be stained by the Mäule procedure as they are in the wild-type plants is astonishing. This result could be attributed to lignin units under detectable levels and particularly to the structure of the lignin polymer in the double mutant. Further experiments are currently underway to explain this phenomenon.
Alteration in the Transcript Levels of Cell WallAssociated Genes
Although perturbations of lignin-related genes such as CCR and C3H could induce severe phenotypes (Jones et al., 2001
; Franke et al., 2002
), most CAD upregulated or downregulated single mutants or transgenics did not result in a strong phenotype and/or modifications in lignin content. Baucher et al. (2003)
explained this phenomenon by the ability of CAD-modified plants to incorporate aldehyde precursors into lignin to compensate for the reduced availability of corresponding monolignols. Because the cad-c cad-d double mutant produced the strongest impact on lignin described to date in CAD downregulated plants, we were interested in determining whether other genes of the phenylpropanoid pathway were altered in this mutant and contributed to its phenotype. We took advantage of the Affymetrix Arabidopsis ATH1-24K GeneChip to analyze alterations of whole genome expression in the lignified tissues of the mutant compared with the wild type.
The consistent expression level of the majority of the genes encoding enzymes positioned upstream of the CAD step in both lines suggested that the double mutation does not induce substantial feedback on their transcription in the tissues harvested. This finding is in accordance with the unchanged expression of most transcription factors previously shown to promote lignification (Borevitz et al., 2000
; Liljegren et al., 2000
; Mele et al., 2003
; Newman et al., 2004
). By contrast, alterations in transcript abundance of other genes encoding enzymes downstream from the CAD step, such as laccases and peroxidases, suggest different regulation for those genes compared with genes involved in monolignol biosynthesis. For example, altered transcript levels of laccase and peroxidase genes could be related to a decrease in the amount of conventional monomer availability for oxidative polymerization. However, microarray analyses were only performed at one time point during the growth and development of the Arabidopsis stem. Therefore, it is possible that expression of these genes might be impacted in the double mutant at other developmental stages.
Alteration of transcription levels of other genes involved in cell wall biosynthesis, such as cellulose synthases, arabinogalactans, pectinacetylesterases, glycosyl hydrolases, and UDP-glucose transferases, suggests that the cell wallrelated phenotype of the double mutant goes beyond modified lignin deposition, potentially affecting synthesis of cell wall polysaccharides. However, additional chemical analysis is needed to substantiate this possibility. The correlation of GRP transcript abundance with the decrease in PRP transcripts, as well as the decrease in IRX3 and IRX5 transcripts, which are specific to secondary cell wall biosynthesis, suggests a delay in secondary cell wall formation in this double mutant. Indeed, under normal conditions, GRP abundance in the primary wall is negatively correlated with lignification (Condit, 1993
). It has been suggested that GRP may be part of a structural mechanism stabilizing conductive elements to prevent water flow between adjacent protoxylem (Ringli et al., 2001
), and the increased level of GRP transcripts in the double mutant may thus reflect activation of a compensatory mechanism in response to a mechanical deficiency.
Interestingly, transcripts related to detoxification, such as GST proteins, are also overrepresented in the double mutant. GST catalyzes the conjugation of a variety of hydrophobic, electrophilic, and cytotoxic substrates that are often xenobiotics. Aldehyde derivatives that accumulate in the double mutant are also members of a group suggested to be GST substrates, based upon a common chemical signal named the Michael acceptor, consisting of a carboncarbon double bond adjacent to an electron-withdrawing group (Marrs, 1996
). GST-encoding genes can be induced by their own substrates (Marrs, 1996
). One could therefore speculate that a detoxifying process is activated in stem tissues of the double mutant in response to aldehyde accumulation.
Neither microarray nor real-time PCR analyses revealed a significant decrease in transcript levels for any of the Arabidopsis CAD genes in the double mutant other than the expected CAD-C and CAD-D. These findings suggest that the other CAD members are not associated with the dramatic phenotype of cad-c cad-d. Notably, CAD-A and CAD1 transcripts were slightly more abundant in cad-c, cad-d, and cad-c cad-d mutants, and this increase in transcript level seems to be correlated with the severity of lignin phenotypes in these mutants (Sibout et al., 2003
; this work). One could therefore speculate that a CAD compensatory mechanism may take place in the double mutant. This mechanism could be responsible to some extent for the biosynthesis of the conventional monolignols still weakly detectable in the double mutant during thioacidolysis despite the lack of significantly detectable in vitro CAD activity of CAD1, CAD-A, and CAD-G (Kim et al., 2004
). However, ß-glucuronidase promoter expression analysis for CAD-G confirmed a weak but significant expression pattern in lignified tissues of wild-type stems and also a strong expression pattern in fiber tissues of the double mutant stems (data not shown). Finally, a robust analysis of CAD gene transcription in wild-type and cad mutant plants was obtained by absolute quantification with real-time PCR. This method allowed detection of the magnitude of CAD gene expression in a lignified tissuethe base of the Arabidopsis stem. Thus, CAD-B1 and -B2 were shown to be very poorly expressed, which would support an insignificant role in stem lignification.
Precursors of Syringyl and Guaiacyl Units Are Synthesized by CAD-C and -D
Our results suggest that both CAD-C and -D are involved in the synthesis of syringyl and guaiacyl precursors in Arabidopsis for the following reasons. First, both of the single mutants were affected in sinapyl and coniferyl alcohol activities, with cad-d being more deficient in sinapyl alcohol activity than cad-c (Sibout et al., 2003
). Second, the slight decrease in lignin-derived thioacidolysis monomers observed in the single mutant cad-d, combined with a dramatic decrease in the double mutant, highlights the importance of the redundant expression pattern observed previously for the two genes (Sibout et al., 2003
), validating the synergistic roles for supplying G and S precursors. Third, the three types of lignin units (G, S, and H) are drastically affected in the double mutant, resulting in the production of a highly condensed lignin in which sinapaldehyde and coniferaldehyde are incorporated through a ß-O-4 bonding pattern. Finally, additional confirmation of a complementary role for CAD-D and CAD-C in G- and S-unit biosynthesis was provided by the ability of each gene to complement the double mutant when expressed under the control of the CAD-D promoter.
That two genes should be involved in both coniferyl and sinapyl alcohol biosynthesis is intriguing. The possibility that they are the result of a gene duplication event (Tavares et al., 2000
) that led to the evolution of additional functions is supported by expression analysis that showed that CAD-C has a broader pattern of expression when compared with CAD-D (Sibout et al., 2003
). A range of aldehyde substrates are potentially available to CAD-C in the Arabidopsis leaf (for review, see Kirch et al., 2004
). CAD-C may also be able to modulate the pools of ferulate and sinapate in Arabidopsis leaves because an aldehyde dehydrogenase REF1 has been discovered recently (Nair et al., 2004
). The authors showed that cell walllinked ferulate esters were decreased in the ref1 plants. The slight amount of ferulic and sinapic acids released by thioacidolysis of extract-free stems of the double mutant is consistent and suggests that aldehydes accumulated in the mutant because of the CAD deficiency could be oxidized by REF1 to produce their corresponding acids. The high raw signal detected for REF1 in our microarray analysis of Arabidopsis stems (see Supplemental Table 9 online) is consistent with this hypothesis.
cad-c cad-d as a Tool to Study CAD/SAD Genes from Woody Plants
The absence of S units in gymnosperm lignin has long been attributed to a lack of the enzymatic function (such as 5-hydroxylase or 5-O-methyl transferase activity) required for the transformation of coniferaldehyde into sinapyl alcohol. CAD proteins isolated from gymnosperms showed very low sinapyl alcohol activity in vitro, whereas most angiosperm CAD proteins displayed suitable activity for both sinapyl and coniferyl alcohols, consistent with the lignin structure and the taxonomy of these plants. The ability of the spruce CAD to produce S units in the double mutant was surprising at first. However, the low but bona fide sinapaldehyde affinity often detected in gymnosperm CAD proteins and the amounts of syringyl moieties found in some conifers support our results (Fullerton and Franisch, 1983
; Lewis and Yamamoto, 1990
). Our data clearly show that the absence of S units in some conifers (at least in spruce) is not due to the CAD gene. Therefore, introduction of angiosperm CAD genes into transgenic conifers to produce enriched S lignins for kraft pulping improvement may not be necessary, although it may be useful for increasing syringyl yield.
It is interesting to note that whereas the tissue-specific expression of Arabidopsis CAD-D restored the lignin phenotype of the mutant with a recovery of the S:G thioacidolysis ratio and decrease of aldehydes, CAD-C overexpression only partially restored the S:G ratio, suggesting slight differences in substrate affinity between these enzymes. In addition, restored lignin in ChimAtCAD-C lines is comparable to lignin of ChimPaCAD lines, suggesting similar protein properties between them, consistent with results obtained by Kim et al. (2004)
, who found that Arabidopsis CAD-C protein showed higher affinity to coniferaldehyde than to sinapaldehyde.
Recently, Lapierre et al. (2004)
reported a detailed analysis of structural traits of lignins in CAD-deficient poplars. They showed that the lignin of their CAD-deficient poplar lines had fewer syringyl and ß-O-4 units and displayed a slight but significant level of sinapaldehyde. The amount of coniferaldehyde units, however, was not enhanced. In our work, we show that the same gene complements cad-c cad-d, completely restoring the level of conventional G units in the lignin, but also that of the conventional S units. These observations suggest a bifunctional activity for this protein. These results, however, are not consistent with the role proposed by Li et al. (2001)
for P. tremuloides CAD, which is homologous to P. deltoides CAD, which is supposed to be involved in coniferyl alcohol biosynthesis only. Our study also shows that the role of SAD genes in lignification requires further investigation. The ability of P. tremula x tremuloides SAD1 to restore the thioacidolysis S:G ratio in ChimPttSAD1 lines confirms the existence of an aldehyde reduction activity of the SAD protein in accordance with in vitro assays shown by Li et al. (2001)
and by Bomati and Noel (2005)
. Nevertheless, the very low yield in syringyl units and the amount of sinapaldehyde remaining in ChimPttSAD1 lines compared with other restored lines (including ChimPaCAD) do not support an essential role for this gene in the incorporation of S units into the constitutive lignin polymer. Nevertheless, the similarity of SAD to defense-related proteins suggests that this gene might be involved in synthesis of syringyl units or monolignol-like compounds in very specific conditions in poplar (Bomati and Noel, 2005
).
In conclusion, the unusual incorporation of sinapaldehydes and coniferaldehydes in the lignin of the double mutant and the significant sinapaldehyde amounts remaining in lines overexpressing spruce CAD supports an aldehyde pathway (Osakabe et al., 1999
) rather than an alcohol pathway for the production of sinapyl alcohol in Arabidopsis, despite the fact that an alcohol pathway may also be present (Humphreys et al., 1999
). Also, accumulation of ferulic and sinapic acids in the double mutant is consistent with our model (Figure 6) and is supported by the recent description of an aldehyde dehydrogenase (REF1; Nair et al., 2004
). However, our work does not exclude involvement of at least a third enzyme displaying very weak CAD activity responsible for the low content of conventional units (mainly guaiacyl units) remaining within the lignin of the double mutant. Thus, our work shows that the last step of monolignol biosynthesis in Arabidopsis is under the control of a multigenic process. CAD-D and CAD-C, which display a synergistic role in reducing coniferaldehyde and sinapaldehyde into corresponding alcohols, are the two main protagonists of the last reduction step in monolignol biosynthesis in Arabidopsis stems.
|
| METHODS |
|---|
|
|
|---|
For RNA isolation and lignin and histological analyses of mature stems, each line was grown in the same growth chamber with random arrangement for 8 weeks at 23°C under 12 h of light per day and low moisture. Plants were grown as described by Chaffey et al. (2002)
for inducing hypocotyl secondary growth.
Protein Gel Blot Analysis
Protein gel blot analysis was conducted as described by Sibout et al. (2003)
. Briefly, total protein extracts were obtained by homogenization of mature stems in 100 mM Tris-HCl, pH 7.5, containing 10% polyvinylpolypyrrolidone (Sigma-Aldrich, St. Louis, MO) and 10 µM ß-mercaptoethanol. Protein samples (15 µg) were heated at 95°C for 5 min and loaded on a 15% acrylamide SDS-PAGE with a 12% resolving gel using a Bio-Rad Protean II apparatus (Hercules, CA). Proteins were transferred onto 0.45-µm nitrocellulose membrane (Amersham Biosciences, Piscataway, NJ) by electroblotting. The polyclonal antibodies directed against Nicotiana tabacum xylem CAD2 were used at a 1:800 dilution. Blots were developed using the Western Light Plus kit (Applied Biosystems, Foster City, CA) according to the manufacturer's instructions.
Arabidopsis Transformation
Binary vectors were introduced into the Agrobacterium tumefaciens strain C58pMP90 (Koncz and Schell, 1986
) by electroporation. Plants were transformed using the flower infiltration protocol (Bechtold and Pelletier, 1998
). T1 transgenic plants were selected on medium containing hygromycin (50 mg/L) or kanamycin (100 mg/L).
Histology
Fresh hand-cut sections were subjected to histochemical analysis or immediately observed under a Zeiss Axioskop microscope (Jena, Germany). Wiesner staining was performed by incubating sections in 1% phloroglucinol in ethanol:water (7:3) with 30% HCl. Mäule staining was performed by first incubating sections in KMnO4. After 10 min, sections were washed and acidified with HCl for 1 min, washed again, and then incubated in NaHCO3. For detection of carbohydrates and phenolics in the hypocotyls, sections were incubated 15 min in a Safranin O/Alcian Blue 1:1 (v/v) solution and then washed before observation under UV light.
RNA Extraction and Real-Time PCR
Because cad-c cad-d is delayed in growth and bolting and to assure uniformity of samples, the first emerging stems (5 to 10 cm long) were removed by cutting at the base and not used for analysis. This removal of the first stem allowed us to standardize the growth of all secondary stems that emerged subsequently at the same time for all lines. The development of the first 20 cm of stem was identical between mutants and controls under our growth conditions. At this stage of growth, the base of the stems (first 5 cm) was harvested at random from individual plants, three of which were pooled to constitute one sample, thereby providing enough material for analysis. Cauline leaves and secondary buds were systematically removed from stem tissues. For each line (the wild type, cad- d, cad-c, and cad-c cad-d), five samples were harvested. Samples were immediately frozen in liquid nitrogen and stored at 80°C until extraction.
RNA from each sample was extracted using the RNeasy extraction kit (Qiagen, Valencia, CA) according to the manufacturer's recommendations, with an extra DNase I treatment. RNA samples were quantified using a Beckman spectrophotometer (Fullerton, CA). Reverse transcriptase (Superscript II; Invitrogen, Carlsbad, CA) reactions were performed using 1 µg of total RNA in a 20-µL volume according to the manufacturer's instructions. These reactions were done at the same time for all samples. The reverse-transcribed cDNA was then treated with RNase H and diluted fivefold before PCR analysis.
Real-time PCR analyses were performed using the Quantitec SYBR green system (Qiagen) and the Opticon PCR machine (MJ Research, Waltham, MA). Data were treated using the Opticon 2 software provided by MJ Research. Because Arabidopsis CAD genes belong to a small family, oligonucleotide primers were designed to take advantage of the unique sequences on the 3' untranslated cDNA regions and had a Tm close to 70°C (see Supplemental Table 12 online). PCR amplification was done in two steps: DNA denaturation at 94°C for 10 s and elongation at 62°C for 2 min. Fluorescence was evaluated at the end of the 2-min elongation. PCR reactions were maintained for 45 cycles. Absolute quantity of transcripts was calculated using standard curves as described by Rutledge and Côté (2003)
.
Microarray Analysis
Samples were harvested and prepared as described for real-time PCR analysis, except that 10 bases of stems of individual plants were pooled to provide enough material for one sample. Four samples were harvested for each line (cad-c cad-d and control). Ten micrograms of total RNA from each biological sample were used for each ATH1 GeneChip hybridization (Arabidopsis Genome ATH1 Array). Complementary RNA was prepared as described in the Affymetrix ATH1 GeneChip Expression Analysis Technical Manual (Affymetrix, Santa Clara, CA; available at http://genomequebec.mcgill.ca). Data generated from Affymetrix Microarray Suite software (MAS5.0) was downloaded in Excel format for each single ATH1 GeneChip (see Supplemental Tables 1 to 8 online) from the Affymetrix Web site to verify the level of confidence (P value) for the probe set of interest. Each signal value was assigned a detection call of P (present, indicating that the signal is of high quality), A (absent, indicating that the signal is generally too low and/or unreliable), or M (marginal, indicating that the hybridization signal is negligible [no transcript detected]). The unpaired Student's t test (P = 0.05) was performed on the data from the eight ATH1 GeneChip samples using MAS5.0 software to select a first set of genes potentially upregulated or downregulated in the mutant. False positives were eliminated by subtracting probe set signals that showed a difference lower than 200 and a fold change of less than two.
Transformation Vectors and Cloning Procedures
The binary plasmid pCAMBIA-1390-HygR (Cambia, Canberra, Australia) was modified and used for this study. A plant expression system was used to express specifically different CAD genes with the goal to complement the cad-c cad-d mutant. In each vector, the different CAD genes were under the control of the Arabidopsis CAD-D promoter and the nopaline synthase gene terminator. The CAD-D promoter, including the 5' untranslated region, was amplified from the previous construction pCAD-D:ß-glucuronidase (see Sibout et al., 2003
) using primers EcoRIpcad3const 5'-GGAATTCGAAATTCTCCACTCGTAGCTCTTCGTTCTG-3' and SpeIpcad3const 5'-GACTAGTTCTTTCTTCTTTCTTATCTTGATCTGCTGC-3'. The PCR product was cloned using the pGEM-T Easy vector cloning kit (Promega, Madison, WI). Selected clones were sequenced to verify the absence of sequence errors and digested to insert the amplicon between the EcoRI-SpeI sites of the binary vector pCAMBIA 1390. This new vector, named CHIMpCAMBIA 1390, was used for all subsequent constructions.
cDNAs corresponding to Populus tremula x tremuloides SAD1 (AY850131), Picea abies CAD (AJ868574), and Populus deltoides CAD (Z19568) were cloned from P. tremula x tremuloides, P. abies, and P. deltoides, respectively. P. abies CAD was a kind gift from M.H. Walter (Leibniz Institute of Plant Biochemistry, Halle, Germany), whereas cDNAs for Arabidopsis CAD-D (At4g34230) and Arabidopsis CAD-C (At3g19450) were cloned from stem RNA using RT-PCR (ecotype Wassilewskija). Besides the SpeI site inserted into both primers and used to amplify each cDNA, five adenosines were added between the SpeI site and the ATG start codon of each gene (only the translated region, not the 5' and 3' untranslated regions, was amplified for each gene). Resulting amplicons were cloned into the pGEM-T easy vector system, sequenced, digested with SpeI, and inserted into the binary vector CHIMpCAMBIA 1390. Clones with the appropriate sense of open reading frames were selected after sequencing.
FTIR Spectroscopy
Stems from 3-week-old plants were used f