Plant Cell Illumina
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


First published online June 13, 2008; 10.1105/tpc.108.060541

The Plant Cell 20:1603-1622 (2008)
© 2008 American Society of Plant Biologists

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplemental Data
Right arrow All Versions of this Article:
20/6/1603    most recent
tpc.108.060541v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via CrossRef
Right arrow Citing Articles via Web of Science (5)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Lee, D. W.
Right arrow Articles by Hwang, I.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Lee, D. W.
Right arrow Articles by Hwang, I.
Agricola
Right arrow Articles by Lee, D. W.
Right arrow Articles by Hwang, I.

Arabidopsis Nuclear-Encoded Plastid Transit Peptides Contain Multiple Sequence Subgroups with Distinctive Chloroplast-Targeting Sequence Motifs[W]

Dong Wook Leea,1, Jong Kyoung Kimb,1, Sumin Leea, Seungjin Choib, Sanguk Kimc,2 and Inhwan Hwanga,2

a Laboratory of Cellular Systems Biology, Division of Molecular and Life Sciences, POSTECH, Pohang 790-784, Korea
b Department of Computer Science, POSTECH, Pohang 790-784, Korea
c Laboratory of Structural Bioinformatics, Division of Molecular and Life Sciences, POSTECH, Pohang 790-784, Korea

2 Address correspondence to sukim{at}postech.ac.kr or ihhwang{at}postech.ac.kr.


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 RESULTS
 DISCUSSION
 METHODS
 REFERENCES
 
The N-terminal transit peptides of nuclear-encoded plastid proteins are necessary and sufficient for their import into plastids, but the information encoded by these transit peptides remains elusive, as they have a high sequence diversity and lack consensus sequences or common sequence motifs. Here, we investigated the sequence information contained in transit peptides. Hierarchical clustering on transit peptides of 208 plastid proteins showed that the transit peptide sequences are grouped to multiple sequence subgroups. We selected representative proteins from seven of these multiple subgroups and confirmed that their transit peptide sequences are highly dissimilar. Protein import experiments revealed that each protein contained transit peptide–specific sequence motifs critical for protein import into chloroplasts. Bioinformatics analysis identified sequence motifs that were conserved among members of the identified subgroups. The sequence motifs identified by the two independent approaches were nearly identical or significantly overlapped. Furthermore, the accuracy of predicting a chloroplast protein was greatly increased by grouping the transit peptides into multiple sequence subgroups. Based on these data, we propose that the transit peptides are composed of multiple sequence subgroups that contain distinctive sequence motifs for chloroplast targeting.


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 RESULTS
 DISCUSSION
 METHODS
 REFERENCES
 
Nuclear-encoded plastid proteins contain N-terminal transit peptide sequences with all the information necessary for their import into plastids (von Heijne et al., 1989Go; Bruce, 2000Go; Emanuelsson and von Heijne, 2001Go; Schein et al., 2001Go; Zhang and Glaser, 2002Go; Kessler and Schnell, 2004Go; Schleiff and Soll, 2005Go). The acquisition of transit peptides by plastid proteins is one of the most critical events in plastid evolution (McFadden, 1999Go), making the transit peptide's evolutionary origin of significant interest. Various studies have probed how transit peptides may have evolved during plastid endosymbiosis and have proposed that plastid proteins acquire N-terminal transit peptides via exon shuffling, alternative splicing, or gene duplication (Gantt et al., 1991Go; Arimura et al., 1999Go; McFadden, 1999Go). The origin of the transit peptide has also been linked to the origin of the import channel (Reumann et al., 1999Go). Despite extensive research on the topic, the evolutionary origin of the transit peptide and how it became linked to nuclear-encoded plastid proteins remain refractory to investigation.

Another important question is the nature of the information encoded by the transit peptide. This question has been addressed by many different approaches (von Heijne et al., 1989Go; Rensink et al., 1998Go; Wienk et al., 1999Go). One approach is to dissect the sequence information encoded by transit peptides and elucidate how this information affects plastid protein delivery into chloroplasts. Such sequence information can be determined by identifying the sequence motifs or domains involved in interactions with import factors, such as import receptors, heat shock proteins, and guidance proteins, that play critical roles in different steps of the import process. Such analysis has shown that transit peptides are composed of multiple domains (Rensink et al., 1998Go; Gutensohn et al., 2000Go; Rial et al., 2000Go; Hinnah et al., 2002Go; Jarvis and Soll, 2002Go; Becker et al., 2004Go; Kessler and Schnell, 2004Go; Smith et al., 2004Go). Recently, Lee et al. (2006)Go demonstrated that the small subunit of the ribulose-1,5-bis-phosphate carboxylase/oxygenase complex (RbcS) contains multiple sequence motifs that play important roles in various steps of protein translocation through the plastid envelope membrane.

Another approach used to address the information carried by transit peptides is to biochemically study the transit peptide structure. Transit peptides are largely unstructured in aqueous solution (Wienk et al., 1999Go). However, in aqueous solutions containing detergent micelles or in membrane mimetic solutions such as trifluoroethane, transit peptides form {alpha}-helical structures (Wienk et al., 1999Go, 2000Go). These results suggest that transit peptides form helical structures during their association with membranes, such as chloroplast envelope membranes. In addition, bioinformatics has been used to analyze the amino acid sequences of all known transit peptides (the transit peptidome) in several studies, which revealed that there is a high preference for hydroxylated amino acids (i.e., Ser, Thr, and Pro) and a lack of acidic amino acids in the transit peptidome, and a high degree of hydrophobicity at the transit peptide N terminus (von Heijne et al., 1989Go; Zhang and Glaser, 2002Go; Bhushan et al., 2006Go). However, results obtained from a few prototypical transit peptides do not generally extend to other transit peptides in the transit peptidome, and the identification of transit peptide functional sequences remains elusive.

Here, we analyzed transit peptide sequences in detail by two different approaches, in vivo protein import experiments using Ala substitution mutants and a bioinformatics-based prediction approach. These combined approaches provide evidence that transit peptides can be grouped into multiple subgroups with distinct sequence motifs.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 RESULTS
 DISCUSSION
 METHODS
 REFERENCES
 
The Transit Peptidome Is Composed of Highly Diverse Sequences
Early attempts to identify transit peptide homology blocks from the transit peptidome failed, due to the high degree of sequence variability among transit peptides (Emanuelsson et al., 1999Go; Bruce, 2000Go; Schein et al., 2001Go). One possibility is that this failure arises from a lack of identified sequence motifs that could be used to define consensus sequences or common sequence motifs. We recently identified sequence motifs from the RbcS transit peptide (Lee et al., 2006Go). We attempted to define any common sequence motifs from the transit peptidome. Similar to previous studies, we were unable to identify any consensus sequences within the transit peptidome. Instead, transit peptides from only a small number of proteins contained identical or similar sequence motifs to those of RbcS. These results, together with the high sequence variability of plastid transit peptides, raise the possibility that the transit peptidome does not contain any consensus sequence motifs, but rather, different transit peptides contain distinct sequence motifs.

To test this hypothesis, we performed a hierarchical clustering on the transit peptidome of Arabidopsis thaliana, which consisted of 208 plastid proteins that have been experimentally confirmed as chloroplast proteins (see Methods). The hierarchical clustering yielded a dendrogram showing that the transit peptide sequences are highly dissimilar (see Supplemental Figure 1 online). We chose seven proteins from this dendrogram whose transit peptide sequences were highly different from each other: RbcS, chlorophyll a/b binding protein (Cab), biotin carboxyl carrier protein (BCCP), protochlorophyllide oxidoreductase A (PORA), DnaJ-J8, tocopherol cyclase (TOCC), and ferredoxin-dependent glutamate synthase 2 (GLU2). When these transit peptides were compared, they displayed highly diverse amino acid sequences (Figure 1A ). Their low sequence identity and insignificant P values of the global alignment scores when aligned with each other (Figure 1B) suggested that these transit peptides are markedly distant from each other at the amino acid sequence level.


Figure 1
View larger version (41K):
[in this window]
[in a new window]
[as a PowerPoint slide]
 
Figure 1. In Vivo Chloroplast Targeting of GFP Fusion Proteins in Protoplasts.

(A) The amino acid sequences and the cleavage site predictions. The first 80 residues of RbcS, BCCP, Cab, DnaJ-J8, PORA, TOCC, and GLU2 are shown. The transit peptides of seven representative proteins were analyzed by the ChloroP algorithm to predict the cleavage sites, which are indicated by arrowheads. Underlining in the Cab transit peptide indicates a region where actual processing may occur.

(B) Sequence identity and corresponding P values for seven representative transit peptides. Sequence identity was calculated with the Needleman-Wunsch algorithm after aligning two sequences. The gap open and extension penalty was 5, and the substitution matrix was BLOSUM 50. Statistical significance of the global alignment scores was assessed by calculating P values, where scores of the alignments to 10,000 randomly permuted sequences were approximated by Gaussian distribution.

(C) GFP reporter protein localization in protoplasts. Protoplasts were transformed with the GFP construct indicated, and GFP localization was examined by fluorescence microscopy. Bars = 20 µm.

(D) Protein gel blot analysis of the reporter protein. Protoplasts were transformed with the GFP construct indicated together with red fluorescent protein (RFP), which was used as a control for soluble protein contamination in the purified chloroplast fraction. Intact chloroplasts were purified from transformed protoplast lysates by Percoll gradient and were analyzed by protein gel blots using anti-GFP and anti-RFP antibodies. Total protoplast extracts were included in the analysis. Immunoblots were stained with Coomassie blue, and the large subunit of the ribulose-1,5-bis-phosphate carboxylase/oxygenase complex (RbcL) was used as a loading control. Pre, precursor forms; Pro, processed mature form; T, total extracts; CH, purified chloroplasts.

(E) Confirmation of predicted cleavage sites. Based on prediction and the apparent sizes of processed reporter proteins, we prepared transit peptide deletion mutants lacking the N-terminal region (black in [A]). For Cab, the processed Cab:GFP apparent size was much smaller than that predicted. We therefore estimated the cleavage site based on apparent size. For in vitro translation, an ATG initiation codon was added to the PCR products' N terminus, and the deletion mutants were translated using wheat germ extracts. In vitro–translated products and processed forms were separated together. P, protoplast extracts; T, in vitro translation.

 
Transit Peptides of the Representative Proteins Contain Highly Specific Sequence Motifs
To examine the significance of the identified sequence diversity with respect to protein import into chloroplasts, we generated the reporter proteins Cab:GFP, BCCP:GFP, DnaJ-J8:GFP, PORA:GFP, GLU2:GFP, and TOCC:GFP, by fusing the N-terminal regions (67 or 80 amino acids) of Cab, BCCP, DnaJ-J8, PORA, GLU2, and TOCC to green fluorescent protein (GFP), respectively. The N-terminal regions used for the construction of GFP reporter proteins included the processing site predicted by ChloroP plus its downstream 8 to 34–amino acid region (Figure 1A). Fluorescence microscopy revealed that these reporter proteins were targeted to chloroplasts (Figure 1C), as was previously observed for RbcS-nt:GFP (Hand et al., 1989Go; Reinbothe et al., 1995Go; Kohler et al., 1997Go; Lee et al., 2006Go). To confirm the import of these reporter proteins into chloroplasts, protein extracts from the transformed protoplasts were analyzed by protein gel blotting using an anti-GFP antibody. In all six transit peptides, >95% of the reporter proteins were proteolytically processed, resulting in the production of smaller protein species. The processed forms copurified with chloroplasts (Figure 1D), confirming that the N-terminal regions of the six proteins were sufficient to function as transit peptides for targeting proteins to chloroplasts.

We next determined the approximate processing site in these transit peptides. We initially employed the prediction algorithm ChloroP to predict the cleavage sites in the transit peptides. The predicted cleavage sites were located between amino acid positions 54 to 55, 35 to 36, 46 to 47, 61 to 62, 53 to 54, 72 to 73, and 47 to 48 in RbcS, Cab, DnaJ-J8, BCCP, PORA, GLU2, and TOCC transit peptides, respectively (arrowheads, Figure 1A, arrowheads). It was previously shown that processing of the RbcS transit peptide occurs at the cleavage site predicted here (Clark and Lamppa, 1991Go). For DnaJ-J8, BCCP, PORA, GLU2, and TOCC, the apparent sizes of the processed forms agreed with those predicted by ChloroP.

To confirm that processing occurred at the predicted cleavage sites, we generated transit peptide deletion mutants lacking the N-terminal region upstream of the cleavage site. In contrast with the other transit peptides, the processed form of Cab reporter proteins migrated faster than the predicted one. Thus, based on the apparent size of the Cab processed form, we generated a deletion mutant that lacked 55 amino acids from the N-terminal region (Figure 1A, arrow). These six deletion constructs were then translated in vitro and separated along with their corresponding processed forms on SDS-polyacrylamide gels. The deletion mutants of DnaJ-J8, BCCP, PORA, GLU2, and TOCC migrated to the same positions as their corresponding processed forms (Figure 1E), indicating that the cleavage occurred at or near the predicted amino acid positions. By contrast, the processed form of Cab reporter proteins in protoplasts migrated slightly faster than the in vitro–translated deletion mutant, indicating that processing occurred further downstream of the site used to generate the Cab deletion mutant. Thus, this processing likely occurred at a position in the underlined region (Figure 1A).

We next identified sequence motifs in individual transit peptides, beginning with the Cab transit peptide. The Cab transit peptide was divided into 10–amino acid segments (T1 to T6), and the segments were replaced one at a time with Ala residues to generate single-T mutants (Figure 2A ) as described previously (Lee et al., 2006Go). GFP constructs with N-terminal mutant transit peptides were introduced into protoplasts, and GFP import into the chloroplasts was examined by fluorescence microscopy. For all of the proteins tested, the T1A mutants showed undetectably low expression levels and are discussed separately below. The mutants, Cab[T2A]:GFP, and Cab[T6A]:GFP displayed a GFP pattern similar to that of Cab:GFP (Figure 2B, panels b and f), indicating that the Ala substitution in T2 and T6 did not affect the ability of the Cab transit peptide to target a protein to chloroplasts. By contrast, Cab[T3A]:GFP and Cab[T5A]:GFP showed strong GFP signals in the cytosol (Figure 2B, panels c and e), indicating that these proteins were not imported into chloroplasts.


Figure 2
View larger version (58K):
[in this window]
[in a new window]
[as a PowerPoint slide]
 
Figure 2. Identification of Sequence Motifs in Cab Transit Peptides.

(A) Sequences of wild-type transit peptides and their corresponding Ala substitution mutants.

(B) and (C) Targeting of GFP reporter protein. Protoplasts were transformed with the construct indicated, and reporter protein localization was directly examined by fluorescence microscopy (green = GFP fluorescence; red = chlorophyll fluorescence; yellow = overlap of signals). In panels (a) and (d) of (B), GFP reporter protein localization was detected by immunostaining with anti-GFP antibodies, as no GFP signals were observed directly. Although Cab:GFP had strong GFP signals in chloroplasts (Figure 1C, panel a), protoplasts transformed with Cab:GFP (a) were immunostained with anti-GFP antibody as a control for Cab[T4A]:GFP (d). Proteins from transformed protoplasts were analyzed by protein gel blots using an anti-GFP antibody. Pre, precursor form; Pro, processed mature form; I, intermediate form.

 
Since Cab[T4A]:GFP did not display any detectable GFP signal, protoplasts transformed with Cab[T4A]:GFP were immunostained with an anti-GFP antibody(Figure 2B, panel d). Strong signals were detected as a ring pattern, reminiscent of localization to the envelope membrane as previously seen for some mutants of the RbcS transit peptide (Lee et al., 2006Go), indicating that Cab[T4A]:GFP was not imported into chloroplasts, but rather was trapped at the envelope membrane in a unfolded form. These results demonstrate that the T3, T4, and T5 segments contain important sequence motifs. To confirm this, protein extracts were prepared from transformed protoplasts and analyzed by protein gel blotting using an anti-GFP antibody (Jin et al., 2001Go; Kim et al., 2001Go). Consistent with image analyses, T3A, T4A, and T5A displayed severe defects in protein import into chloroplasts. The majority of T3A:GFP and T5A:GFP accumulated in precursor forms, whereas T4A:GFP generated intermediates with no processed form (Figure 2B, panels i to k). These results indicate that the Cab transit peptide segments T3, T4, and T5 contain sequence motifs critical for the import process.

To further delineate sequence motifs in these regions, the first or second half of each T segment was restored to wild-type sequence in the background of T3A, T4A, or T5A. Restoration of the amino acids KSKF in T3A completely rescued chloroplast import (Figure 2C, panels a and d), and restoration of the sequences PLPN or GNVGR in T4A almost completely rescued the targeting efficiency to wild-type levels (Figure 2C, panels b, e, and h). Restoration of the sequence IRMAA in T5A also restored targeting efficiency to >90% of that of the wild type (Figure 2C, panels c and f). Together, these results strongly suggest that the sequence motifs KSKF (T3), PNPL (T4), GNVGR (T4), and IRMAA (T5) are critical for protein import into chloroplasts and that the two halves of the T4 segment are functionally redundant. These sequence motifs are different from those in the RbcS transit peptide (Lee et al., 2006Go).

A similar approach was employed to identify DnaJ-J8 transit peptide sequence motifs. Protoplasts transformed with single-T mutants were directly observed with fluorescence microscopy. Except for DnaJ-J8[T6A]:GFP, which produced a staining pattern similar to that of the wild type, all the single T-mutants displayed strong GFP signals in the cytoplasm. These results suggest that critical sequence motifs are dispersed throughout the T2 to T5 segments of the DnaJ-J8 transit peptide (Figure 3B ). To confirm this, protein extracts from transformed protoplasts were analyzed by protein gel blotting using an anti-GFP antibody. The mutants produced varying amounts of precursors. Consistent with the image analysis, all of the single-T mutants, except T6A, were severely handicapped in delivering proteins to chloroplasts. T2A, T3A, T4A, and T5A import efficiencies were <5 to 10% that of the wild type (Figure 3B, panel f). The T2A:GFP and T4A:GFP proteolytic products (Figure 3B, asterisks) migrated slightly faster than the mature processed wild-type form.


Figure 3
View larger version (58K):
[in this window]
[in a new window]
[as a PowerPoint slide]
 
Figure 3. In Vivo Targeting of GFP Fusion Proteins by Various Mutant Forms of the DnaJ-J8 Transit Peptide.

(A) Sequences of wild-type and Ala substitution mutant forms of the DnaJ-J8 transit peptide.

(B) and (C) Reporter protein targeting. Protoplasts were transformed with the constructs indicated and analyzed as described in Figure 2. We confirmed that the protein species indicated by asterisks are imported into chloroplasts (see Supplemental Figure 2A online). Pre, precursor forms; Pro, processed form; I, intermediate form. Bars = 20 µm.

 
To test whether the processed forms of T2A:GFP and T4A:GFP were also imported into chloroplasts, intact chloroplasts were purified from gently lysed protoplast extracts on a Percoll gradient and analyzed by protein gel blotting using an anti-GFP antibody. Processed T2A:GFP and T4A:GFP forms copurified with chloroplasts, suggesting that they had been imported (see Supplemental Figure 2A online). In addition, the upper bands of T2A:GFP and T4A:GFP migrated slightly faster than their wild-type precursors. To verify that the upper bands of these mutants represented precursors and not proteolytic products, they were translated in vitro using wheat germ extract, and their motilities were compared with those produced in protoplasts. The migration patterns of the in vitro–translated mutant proteins were identical to their corresponding protoplast forms, indicating that the size differences in the single-T mutant precursors were likely due to Ala substitution (see Supplemental Figure 3A online). These results indicate that T2, T3, T4, and T5 of the DnaJ-J8 transit peptide contain sequence motifs critical for protein import.

DnaJ-J8 sequence motifs were further delineated into half-T regions (Figure 3A). In T2A, restoration of either the first or second half of T2 improved targeting efficiency to the wild-type level (Figure 3C, panels a and e). Restoration of SSSSS in T3A rescued targeting efficiency to >90% of that of the wild type (Figure 3C, panels b and f). Restoration of either the first or second half of T4 or T5 was sufficient to rescue T4A and T5A import efficiencies to >90% of that of the wild type (Figure 3C, panels c, d, g, and h). These results indicate that SSSSS (T3), NSRRK (T4), NTKML (T4), NRSKV (T5), and VCSSS (T5) are critical sequence motifs for protein import and that GFSGL and PGSSF of T2 also play an additional role in import. The results also strongly suggest that GFSGL and PGSSF in T2, NSRRK and NTKML in T4, and NRSKV and VCSSS in T5 are functionally redundant to each other.

It has been shown that there is a high preference for hydroxylated amino acids, such as Ser, Thr, and Pro, in the transit peptidome (von Heijne et al., 1989Go; Zhang and Glaser, 2002Go; Bhushan et al., 2006Go). Interestingly, the T3 sequence motif of DnaJ-J8 consists only of Ser residues, and T5 contains three Ser residues out of its 5–amino acid residues. This result may contribute to the idea that Ser is abundant in the transit peptide (Zhang and Glaser, 2002Go).

Sets of single- and double-T mutants were also generated for the BCCP transit peptide. These mutants were transformed into protoplasts, and GFP staining patterns were observed directly with fluorescence microscopy. Of the single-T mutants, only T3A:GFP and T4A:GFP displayed weak GFP signals in the cytoplasm together with strong GFP signals in chloroplasts (Figure 4B , panels c and d), indicating that the import efficiency of each of the single-T mutants was comparable to the wild type. To confirm these findings, protein extracts from transformed protoplasts were analyzed by protein gel blotting using an anti-GFP antibody. Consistent with the image analysis of single-T mutants, only T3A and T4A displayed minor defects in protein import (Figure 4B, panels j and k). T3A and T4A import efficiencies were quantified at three different time points. During these periods, T3A targeting efficiency increased from 65 to 78%, and T4A targeting efficiency increased from 83 to 90% (Figure 4E, panels a and b). These results strongly suggest a high degree of redundancy in the sequence motifs of the BCCP transit peptide.


Figure 4
View larger version (46K):
[in this window]
[in a new window]
[as a PowerPoint slide]
 
Figure 4. In Vivo Targeting of GFP Fusion Proteins by Various Mutant Forms of the BCCP Transit Peptide.

(A) and (C) Sequences of wild-type transit peptide and their corresponding Ala substitution mutants.

(B) and (D) Targeting of GFP reporter proteins. Reporter proteins were analyzed as described in Figure 2. Pre, precursor form; Pro, processed form; I, intermediate form. Bars = 20 µm.

(E) Quantification of import efficiency. Protein extracts prepared from protoplasts transformed with the constructs indicated at the time points indicated were analyzed by protein gel blotting using an anti-GFP antibody. The intensity of the precursors and processed forms was measured using image processing software connected to LAS3000 (Fuji Photo Film). The intensity was obtained at various exposure times to obtain a linear range of signal intensity. For import efficiency, the intensity of processed forms was divided by the sum of the intensities of the precursors and processed forms.

 
A small proportion of an intermediate form of T6A was detected, suggesting that T6A may have a minor defect in crossing the envelope membrane (Figure 4B, panel m). Of the double-T mutants (Figure 4C), T35A, T36A, and T46A displayed severe import defects (Figure 4D, panels a to c, h, i, and k). In particular, their targeting efficiencies were reduced to 38, 43, and 27%, respectively, at 8 h after transformation (Figure 4E, panels c to e). Their targeting efficiencies gradually increased to 65 to 75% with time. These results indicate that the T3, T4, T5, and T6 segments of BCCP contained critical sequence motifs and that sequence motifs in between T3 and T5 or T6, and between T4 and T6 were functionally redundant.

T segments in the double-T mutants were divided into 5–amino acid regions, and each region was restored in the background of double-T mutants (Figure 5A ). Of the four T35A restoration mutants, T35A+FPIQN and T35A+AKPKL displayed targeting efficiencies similar to T5A and T3A, respectively (Figure 5B, panels a, d, g, and h), indicating that FPIQN (T3) and AKPKL (T5) are critical for protein import. Of the four T36A restoration mutants, T36A+FPIQN and T36A+PSRSS import efficiencies were similar to those of T6A and T3A, respectively, while the T36A+YPVVK import efficiency was similar to that of T3A (Figure 5B, panels b, e, i, and j). These results indicated that FPIQN (T3), PSRSS (T6), and YPVVK (T6) are critical for protein import. Of the four T46A restoration mutants, T46A+RSRRV and T46A+SFRLS import efficiencies were ~90% of that of the wild-type peptide, while the T46A+YPVVK targeting efficiency was similar to that of T4A, indicating that RSRRV (T4), SFRLS (T4), and YPVVK (T6) are critical for protein import and that at least one of them is required for efficient targeting to chloroplasts (Figure 5B, panels c, f, k, and l). These results also suggest that the sequence motifs of the BCCP transit peptide have a high degree of functional redundancy.


Figure 5
View larger version (44K):
[in this window]
[in a new window]
[as a PowerPoint slide]
 
Figure 5. Identification of Sequence Motifs in the BCCP Transit Peptide.

(A) Sequences of wild-type transit peptide and their corresponding Ala substitution mutants.

(B) Targeting of GFP reporter proteins. Reporter proteins were analyzed as described in Figure 2. Pre, precursor form; Pro, processed form; I, intermediate form. Bars = 20 µm.

 
Sequence motifs in the PORA transit peptide were defined using a similar approach. Of the single-T mutants, T5A:GFP produced strong GFP signals in the cytoplasm with almost no signal in chloroplasts (Figure 6B , panel d). GFP signals in T2A:GFP, T3A:GFP, and T6A:GFP were primarily detected in the chloroplasts, with weak signals in the cytoplasm (Figure 6B, panels a, b, and e). Targeting efficiency was determined by protein gel blot analysis using an anti-GFP antibody. The targeting efficiency of T5A was <5% of that of the wild-type transit peptide, whereas T2A, T3A, and T6A retained 40 to 60% of the wild-type targeting efficiency (Figure 6B, panel f). The proteolytic products (Figure 6B, asterisks) of T2A, T4A, and T6A migrated slightly faster than the processed form of the wild type.


Figure 6
View larger version (55K):
[in this window]
[in a new window]
[as a PowerPoint slide]
 
Figure 6. In Vivo Targeting of GFP Fusion Proteins by Various Mutant Forms of the PORA Transit Peptide.

(A) Sequences of the PORA transit peptide mutants.

(B) and (C) Targeting of GFP reporter proteins. Reporter proteins were analyzed as described in Figure 2. We confirmed that the protein species indicated by asterisks are imported into chloroplasts (see Supplemental Figure 2B online). Pre, precursor form; Pro, processed mature form; I, intermediate form. Bars = 20 µm.

 
We examined whether these protein species were also imported into chloroplasts. Chloroplasts were purified from gently lysed protoplast extracts and analyzed by protein gel blotting using an anti-GFP antibody. The proteolytic proteins copurified with chloroplasts, suggesting that they were imported into chloroplasts (see Supplemental Figure 2B online). In addition, as observed with DnaJ-J8, the upper bands of the T2A, T4A, and T6A mutants migrated faster than the wild type (Figure 6B, panel f). To verify that this was also due to Ala substitution, these mutants were translated in vitro using wheat germ extract and their motilities were compared with those expressed in protoplasts. The in vitro–translated proteins produced identical migration patterns to those expressed in protoplasts, confirming that the upper migrating protein species were mutant precursors (see Supplemental Figure 3B online).

In the single-T mutants, the first or second half regions of the T segments were restored (Figure 6A). T2A+SAFSV, T3A+SSSFK, and T6A+SLRCK displayed targeting efficiencies similar to the wild type (Figure 6C, panels a, b, d to f, and h). Restoration of either the first or second half of T5A was equally effective in restoring targeting efficiency to about the wild-type level (Figure 6C, panels c and g). Together, these results suggest that SAFSV (T2), SSSFK (T3), EQSKA (T5), DFVSS (T5), and SLRCK (T6) of the PORA transit peptide are critical sequence motifs for protein import. Furthermore, EQSKA and DFVSS of T5 may be functionally redundant.

To define sequence motifs in the GLU2 transit peptide, we generated single T Ala substitution mutants (Figure 7A ). Of the single-T mutants, T3A:GFP and T5A:GFP produced strong GFP signals in the cytoplasm in addition to the chloroplasts (Figure 7B, panels b and d). GFP signals in T2A:GFP, T4A:GFP, and T6A:GFP were primarily detected in the chloroplasts (Figure 7B, panels a, c, and e). Targeting efficiency was determined by protein gel blot analysis using an anti-GFP antibody. The targeting efficiency of T3A and T5A was ~50% that of the wild-type transit peptide. In addition, T4A also produced ~20% of precursor, whereas T2A and T6A were targeted as efficiently as the wild type (Figure 7B, panel f). These results indicated that T3 and T5 contain critical sequence motifs and the motif in T4 may play a minor role in import into the chloroplast.


Figure 7
View larger version (54K):
[in this window]
[in a new window]
[as a PowerPoint slide]
 
Figure 7. In Vivo Targeting of GFP Fusion Proteins by Various Mutant Forms of the GLU2 Transit Peptide.

(A) Sequences of the GLU2 transit peptide mutants.

(B) and (C) Targeting of GFP reporter proteins. Reporter proteins were analyzed as described in Figure 2. Pre, precursor form; Pro, processed mature form. Bars = 20 µm.

 
In T3A, T4A, and T5A mutants, the first or second half regions of the T segments were restored (Figure 7A). In these mutants, restoration of either the first or second half region was equally effective in restoring the targeting efficiency to approximately wild-type levels (Figure 7C, panels a to g). Together, these results suggest that SAKLS (T3), STKTI (T3), ISKGT (T5), and KRRNE (T5) of GLU2 are critical sequence motifs for protein import. In addition, FSVDF (T4) and VRSYC (T4) also play a role in protein import. Finally, the lack of a sequence motif in T2 and T6 was confirmed by the T26A mutant (Figure 7C, panels d and h).

Finally, to define sequence motifs in the TOCC transit peptide initially, single T Ala substitution mutants were used in protein import experiments (Figure 8A ). Of the single-T mutants, T3A:GFP, T4A:GFP, and T5A:GFP produced strong GFP signals in the cytoplasm in addition to the chloroplasts (Figure 8B, panels b to d), whereas GFP signals of T2A:GFP and T6A:GFP were primarily detected in the chloroplasts (Figure 8B, panels a and e). Targeting efficiency was determined by protein gel blot analysis using an anti-GFP antibody. The targeting efficiency of T3A, T4A, and T5A was ~60, 40, and 85% that of the wild-type transit peptide, respectively, indicating that the sequence motif in T4 is most critical in protein import into the chloroplasts. By contrast, T2A and T6A were targeted as efficiently as the wild type (Figure 8B, panel f). These results indicated that T3 and T4 contain critical sequence motifs and that the motif in T5 may play a minor role in import into the chloroplast.


Figure 8
View larger version (50K):
[in this window]
[in a new window]
[as a PowerPoint slide]
 
Figure 8. In Vivo Targeting of GFP Fusion Proteins by Various Mutant Forms of the TOCC Transit Peptide.

(A) Sequences of the TOCC transit peptide mutants.

(B) and (C) Targeting of GFP reporter proteins. Reporter proteins were analyzed as described in Figure 2. Pre, precursor form; Pro, processed mature form. Bars = 20 µm.

 
In T3A, T4A, and T5A mutants, the first or second half regions of the T segments were restored (Figure 8A). In the T3A mutant, restoration of either the first or second half region restored the targeting efficiency >90% that of the wild-type level (Figure 8C, panel d). In the cases of T4A and T5A mutants, restoration of the first half region restored the targeting efficiency to >95% that of the wild type (Figure 8C, panels e and f). Together, these results suggest that RPVSP (T3), LTRSL (T3), VPFRS (T4), and RSISR (T5) of TOCC are critical sequence motifs for protein import.

The Hydrophobic Domain of Some T1 Regions Plays a Role in the Efficient Targeting of Proteins to Chloroplasts
For all of the fusions characterized above, T1A:GFP (containing Ala substitutions from the second to 10th amino acid positions) was always expressed at undetectable levels, discussed below. Hydrophobic amino acid residues in the T1 segment are critical for protein import into chloroplasts (Emanuelsson et al., 1999Go; Lee et al., 2006Go). Therefore, we generated point mutants by replacing two or three hydrophobic amino acid residues in T1 with Ala residues (Figure 9A ). DnaJ-J8[T1-hy]:GFP, PORA[T1-hy]:GFP, and TOCC[T1-hy]:GFP targeting efficiencies were reduced to ~50% those of wild-type transit peptides (Figure 9B, panels d, e, and g), similar to the RbcS transit peptide (Figure 9B, panel a). By contrast, Cab[T1-hy]:GFP, BCCP[T1-hy]:GFP, and GLU2[T1-hy] displayed a targeting efficiency similar to their wild types (Figure 9B, panels b, c, and f), suggesting that other motifs in the transit peptide of these proteins may compensate for the loss of these hydrophobic residues.


Figure 9
View larger version (25K):
[in this window]
[in a new window]
[as a PowerPoint slide]
 
Figure 9. The Hydrophobic Domain of Some T1 Regions Plays a Role in the Efficient Targeting of Proteins to Chloroplasts.

(A) Sequences of the transit peptide mutants.

(B) Protein gel blot analysis of the GFP reporter fusion proteins.

 
Sequence Motif Prediction of Plastid Transit Peptides
We also asked whether the sequence motifs of these seven highly dissimilar transit peptides could be predicted by sequence information alone. We developed a motif discovery algorithm to identify a set of motif segments discriminating between plastid and nonplastid proteins (see Supplemental Methods online). Because we hypothesize that plastid transit sequences comprise multiple subgroups with distinct sequence motifs, this algorithm attempts to find sequence motifs that occur more frequently in a subgroup of plastid transit sequences compared with the N termini of nonplastid proteins. Using this motif discovery algorithm, we tested whether we could find sequence motifs from transit peptides. We used 208 Arabidopsis plastid proteins that have been experimentally shown to localize to the plastid. The list of these 208 plastid proteins with their predicted sequence motifs is available in Supplemental Data Set 1 online. From these data, we used the previously identified seven highly dissimilar transit peptides and compared their sequence motifs obtained by the motif discovery algorithm with those obtained by our import experiments above. The sequence motifs identified by the two independent approaches were nearly identical or significantly overlapping (Figure 10A ). We calculated P values to evaluate the overlap between the experimentally determined motifs and the predicted ones and found that all of the predictions were statistically significant (P < 0.1) except RbcS and GLU2, suggesting that our algorithm adequately captured the sequence features of plastid transit peptides (Figure 10B). In addition, these results raised the possibility that the sequence motifs identified by transit peptide functional analysis are conserved among members of the subgroups characterized by these sequence motifs.


Figure 10
View larger version (33K):
[in this window]
[in a new window]
[as a PowerPoint slide]
 
Figure 10. Transit Peptide Sequence Motifs Can Be Predicted Based Solely on Sequence Information.

(A) Motif alignment identified by two different approaches. Transit peptide sequences of RbcS, BCCP, Cab, DnaJ-J8, PORA, TOCC, and GLU2 are shown. Sequence motifs characterized experimentally in this study are in red, and the predicted sequence motifs are in blue. The predicted motifs were obtained by minimizing the upper bound of the general misclassification error of discriminating 208 plastid transit peptides from 778 nonplastid proteins of Arabidopsis.

(B) Significance of the overlap between motifs characterized in vivo and those predicted by modeling. The P value was calculated under the null hypothesis, which can be modeled with the hypergeometric distribution. It is unlikely that the predicted motifs are unrelated to the motifs characterized in vivo (P value < 0.1).

(C) The percentage of transit peptides classified into one of the seven sequence groups. The percentage of transit peptides classified into one of the seven groups was based on the identification of sequence motifs at two different P values ({alpha} = 0.01 and 0.05). The transit peptides not classified into one of the seven groups are represented as "None."

 
In the analysis, we initially used the first 80 amino acid residues of each transit protein in our motif discovery algorithm. It was therefore likely that some sequences include residues from mature proteins. To exclude this possibility, we performed motif prediction again, using only amino acid residues from the N terminus to the ChloroP-predicted cleavage site. This caused the performance of our prediction algorithm to degrade significantly, with poor alignment (insignificant P values) between the experimentally determined and predicted motifs (see Supplemental Figure 4 online).

One possible explanation for our failure to predict sequence motifs in this manner is that the analyzed sequences may not have been sufficiently long to reliably produce any conserved sequence motifs. In the cases of ChloroP and TargetP, 100–amino acid residues of the N-terminal regions were used for reliable transit peptide predictions (Emanuelsson et al., 1999Go, 2000Go). This suggests that the transit peptide region downstream of the potential cleavage site may also include information necessary for chloroplast targeting. Consistent with this, critical sequence motifs in the RbcS transit peptide are found in the region downstream of the cleavage site and the mature portion increases targeting efficiency (Comai et al., 1988Go; Lee et al., 2006Go).

To further confirm this possibility, we generated new reporter constructs by fusing GFP to the N-terminal upstream regions of the predicted cleavage sites from Cab, BCCP, PORA, and DnaJ-J8 (see Supplemental Figure 5A online). These GFP reporter constructs, Cab-cs:GFP, BCCP-cs:GFP, DnaJ-J8-cs:GFP, and PORA-cs:GFP, were introduced into protoplasts and examined for their targeting to chloroplasts by imaging and protein gel blot analyses. BCCP-cs:GFP produced primarily processed forms with a minor portion of precursors. Its ratio of processed forms to precursors was nearly identical to that of BCCP-nt:GFP (see Supplemental Figure 5C online), indicating that the N-terminal portion upstream of the cleavage site is sufficient for targeting. By contrast, for other constructs, the precursor levels were greatly increased compared with their respective full-length transit peptides with levels equal to or greater than those of processed forms. In particular, PORA-cs:GFP primarily produced precursors and only a small amount of processed forms. These results demonstrate that in the case of the four transit peptides analyzed here, the upstream regions of the cleavage site are not sufficient for efficient protein import into chloroplasts. This is consistent with the observation that the downstream region of the cleavage site also contains important sequence motifs (Comai et al., 1988Go; Lee et al., 2006Go).

To confirm that the processed forms were imported into chloroplasts, we purified chloroplasts from gently lysed protoplast extracts by Percoll gradient and analyzed the chloroplast fractions by protein gel blotting analysis using an anti-GFP antibody. The processed forms but not the precursors were detected in the chloroplast fractions (see Supplemental Figure 5D online), strongly suggesting that in most transit peptides the region downstream of the cleavage site is necessary for efficient protein import into chloroplasts.

We next examined the 208 transit peptides to determine whether the experimentally identified sequence motifs were also present in other transit peptides. A transit peptide was assumed to have a particular sequence motif and was assigned to the corresponding sequence group if the P value of our proposed alignment score with dual gap penalties between the transit peptide and sequence motifs was lower than a specified cutoff value ({alpha} = 0.01 or 0.05) (see Supplemental Methods online). Of the 208 transit peptides, 48 (23.1%, {alpha} = 0.01) were assigned to at least one of the RbcS, Cab, BCCP, DnaJ-J8, PORA, TOCC, and GLU2 sequence groups (Figure 10C). The RbcS sequence subgroup had the most members (22 transit peptides), while the TOCC sequence subgroup had the least (two transit peptides). When the cutoff P value was set at 0.05, the number of transit peptides increased in all seven groups, with 107 (51.4%) being assigned to at least one and with the number of members in the DnaJ-J8 sequence group increasing the most (from eight transit peptides to 30). These results supported our hypothesis that the transit peptidome is composed of multiple subgroups with distinct sequence motifs. The predicted transit peptide sequences of the seven groups and their P values ({alpha} = 0.05) and motif alignments are available in Supplemental Data Set 2 online. We found that approximately half of the 208 transit peptides could be assigned into one of the seven subgroups; however, a large number of transit peptides were still unassigned to any of the subgroups. These results again confirmed our initial hypothesis that the transit peptide sequences are highly dissimilar, indicating that there are more transit peptide groups yet to be discovered.

Sequence Motif–Based Approach Increases the Accuracy of Chloroplast Transit Peptide Prediction
We were interested in whether the sequence motifs predicted by our motif discovery algorithm could be used to improve prediction of plastid transit peptides. Previously, we reported several methods for extracting the desired features of protein sequences and predicted the subcellular localization of eukaryotic proteins with an empirical kernel map using the scores of pairwise alignments (Kim et al., 2006aGo, 2006bGo). We employed a similar approach to predict plastid transit peptides but explicitly incorporated sequence motifs into our prediction model (see Methods).

The prediction power of our method (MultiP) was compared with four other methods, including ChloroP (Emanuelsson et al., 1999Go), PCLR (Schein et al., 2001Go), TargetP (Emanuelsson et al., 2000Go), and Predotar (Small et al., 2004Go). To obtain an accurate assessment, our method was trained and tested using the same training and test set that were used in each study (see Supplemental Methods online). The resulting three independent comparisons allowed us to remove the effect of the size of training data, which is known to be the main determinant to improved performance (Millar et al., 2006Go). Our method showed significantly improved performance compared with ChloroP or PCLR (Table 1 ). The observed accuracy and Matthew's correlation coefficient (see Supplemental Methods online) were significantly higher than those for ChloroP and PCLR and were largely due to a reduced number of false positives (increased specificity). We then compared our method with TargetP and Predotar trained on large training sets and found that our method also showed a better performance against these methods (Table 1). Specifically, our method provided both high sensitivity and specificity compared with TargetP and Predotar. Our program is available from the website http://sbi.postech.ac.kr/MultiP/.


View this table:
[in this window]
[in a new window]

 
Table 1. Performance Comparison of Our Algorithm (MultiP) with ChloroP, PCLR, TargetP, and Predotar in Predicting Chloroplast Transit Peptides

 

    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 RESULTS
 DISCUSSION
 METHODS
 REFERENCES
 
Here, we demonstrated that the transit peptides of nuclear-encoded plastid proteins can be clustered into multiple subgroups that contain subgroup-specific distinctive sequence motifs. These subgroup-specific sequence motifs were identified by two independent approaches. The first approach was based on the results of protoplast protein import experiments, in which we used Ala substitution mutants (Lee et al., 2006Go) to identify critical sequence motifs for protein import into chloroplasts. The second approach used bioinformatics-based predictions, in which highly conserved sequence motifs were identified among members of each subgroup. These two approaches uncovered nearly identical or significantly overlapping sequence motifs that turned out to be highly subgroup specific. However, we also found that a large number of transit peptides were still unassigned to any of these subgroups. Since our grouping method relied upon sequence alignment with seven arbitrarily selected transit peptides displaying amino acid sequences very different from each other, the selected representative transit peptide for each subgroup may not be the best candidate for clustering. It is also important to note that the highly diverse transit peptidome can be composed of additional subgroups.

These results demonstrate not only that transit peptides contain highly conserved sequence motifs but also that the conserved sequence motifs play critical roles in protein import into chloroplasts. This is in contrast with the current notion that transit peptides do not contain any conserved sequence motifs (Emanuelsson et al., 1999Go; Bruce, 2000Go; Schein et al., 2001Go). Previous attempts to identify conserved sequence motifs may have been hampered by the fact that these sequences are not characteristic of all transit peptides but are confined to specific subgroups. Therefore, the limited success of previous works in predicting chloroplast transit peptides, characterized by low specificity, can be explained by the lack of knowledge of the existence of multiple transit peptide subgroups (Emanuelsson et al., 1999Go; Schein et al., 2001Go). Specifically, the specificity of our method in predicting chloroplast transit peptides was significantly better than that of ChloroP and PCLR. Furthermore, compared with TargetP and Predotar, our method provided well-balanced sensitivity and specificity by exploiting the fact that the transit peptidome is composed of multiple subgroups with distinctive sequence motifs.

Detailed analysis of these representative transit peptides revealed that they display distinct characteristics. RbcS transit peptides contain a large number of sequence motifs throughout the entire transit peptide region. In addition, these sequence motifs display a high degree of functional redundancy, despite a lack of amino acid sequence similarity (Lee et al., 2006Go). Similarly, the BCCP group displayed a high degree of functional redundancy for the sequence motifs; all of the single-T Ala substitution mutants displayed >90% of the wild-type targeting efficiency, and at least two T regions in the transit peptide had to be simultaneously substituted with Ala residues to observe >50% inhibition of chloroplast import. By contrast, Cab, DnaJ-J8, PORA, GLU2, and TOCC transit peptides contained unique sequence motifs, as single-T Ala substitution mutants of these peptides showed severe defects in protein import into chloroplasts. These results raise the possibility that the organization, with respect to redundancy and relative order, of the transit peptide sequence motifs may be different in different peptides. However, it is still not clear if the functional redundancy of sequence motifs in a transit peptide plays a role in determining its targeting efficiency.

In addition, analysis of the amino acid composition of the critical sequence motifs yielded a different picture from that obtained from total transit peptide sequences. It is generally accepted that hydroxylated amino acid residues, such as Ser and Thr, are well represented in transit peptides (Bruce, 2000Go). However, except for the sequence motifs in DnaJ-J8, a large number of sequence motifs identified in the representative transit peptides did not display any clear overrepresentation of hydroxylated amino acids. In the case of DnaJ-J8, the four critical sequence motifs contained substantial numbers of Ser residues. In fact, one of the sequence motifs was composed entirely of Ser residues.

The exact role of these sequence motifs was not directly addressed in this study. Depending on the individual sequence motif, Ala substitution produced characteristic features in the extent of precursor processing and/or production of intermediate forms, similar to the RbcS transit peptide (Lee et al., 2006Go). Thus, it is possible that these sequence motifs are recognized at different steps during the import process by specific protein factors (Bruce, 2000Go; Jarvis and Soll, 2002Go; Kessler and Schnell, 2004Go; Smith et al., 2004Go). Another possibility is that certain sequence motifs are necessary to adopt specific structures, such as {alpha}-helices, which may be required for association with chloroplast membranes (Wienk et al., 1999Go, 2000Go).

The presence of multiple subgroups in the transit peptidome raises the question of how transit peptides belonging to different subgroups are recognized by import receptors. One possibility is that there exist multiple types of import receptors. This hypothesis is supported by results showing that the import receptors Toc159 and Toc33 are responsible for the import of photosynthetic proteins, whereas other receptors, such as Toc132 and Toc34, are needed for the import of nonphotosynthetic ones (Bédard and Jarvis, 2005Go). However, this explanation is not entirely satisfactory, as two transit peptides of Cab and RbcS, which are both involved in photosynthesis, belong to different subgroups. Alternatively, import receptors such as Toc159 might be able to recognize more than one sequence group. For instance, Sec24, a component of COPII vesicles that is involved in cargo recognition during anterograde trafficking from the endoplasmic reticulum to the Golgi complex, can accommodate more than one type of cargo through its multiple cargo binding sites for different endoplasmic reticulum exit signals (Miller et al., 2003Go). Additional experiments will be necessary to confirm this hypothesis.


    METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 RESULTS
 DISCUSSION
 METHODS
 REFERENCES
 
Growth of Plants
Arabidopsis thaliana (Columbia ecotype) was grown in soil at 20 to 25°C in a greenhouse with a 16 h/8 h light/dark cycle or on Murashige and Skoog plates in a growth chamber at 20°C. Leaf tissues were harvested from 2-week-old plants and used immediately for protoplast isolation.

PCR-Based Mutagenesis and Plasmid Construction
The Ala substitution mutant constructs, including the single-T and double-T mutants, were generated using standard molecular cloning techniques as described previously (Lee et al., 2006Go). The primer sequences used to prepare the constructs are shown in Supplemental Table 1 online.

Transient Expression in Protoplasts
The plasmids used in protoplast transformations were purified using Qiagen columns according to the manufacturer's protocol. DNA was introduced into Arabidopsis protoplasts prepared from leaf tissues by polyethylene glycol–mediated transformation (Jin et al., 2001Go).

Images were obtained using a cooled CCD camera and a Zeiss Axioplan fluorescence microscope. The filter sets used were XF116 (exciter, 474AF20; dichroic, 500DRLP; emitter, 510AF23) and XF33/E (exciter, 535DF35; dichroic, 570DRLP; emitter, 605DF50) (Omega) for GFP/FITC and RFP/TRITC, respectively. The data were processed using Adobe Photoshop software, and the images were rendered in pseudocolor.

Protein Gel Blot Analysis and Signal Intensity Quantification
Protein gel blot analysis was performed using protein extracts from transformed protoplasts as described previously (Jin et al., 2001Go; Kim et al., 2001Go). Briefly, protoplasts were lysed by brief sonication in lysis buffer (10 mM Tris-HCl, pH 7.5, 150 mM NaCl, 2.0 mM EDTA, and 0.2 mM PMSF) and subjected to centrifugation at 10,000g at 4°C for 10 min to remove debris. Proteins in the supernatants were separated by SDS-PAGE. Immunoreactive proteins were visualized using enhanced chemiluminescence (ECL kit; Amersham Pharmacia Biotech), and images were obtained using a LAS3000 image capture system (Fuji Photo Film). The immunoblots were quantified by measuring the intensity of the protein bands with LAS3000 software. Immunoblot images were obtained at different exposure times. Only images whose band intensities were within the linear range between intensity increase and exposure time were selected. For import efficiency, the intensity of processed forms was divided by the sum of the intensities of the precursors and processed forms.

In Vitro Translation in Wheat Germ Extracts
The templates used for in vitro translation were generated by PCR. The forward primer consisted of the T7 promoter, transcription start site, and 18 bases complementary to the 5' region of the template. The reverse primer consisted of 18 bases complementary to the 3' region of the nos terminator. One microgram of PCR product was transcribed and translated using wheat germ extracts (TNT T7 coupled wheat germ extract system; Promega) according to the manufacturer's manual. Five microliters of in vitro–translated proteins were mixed with 25 µL of 1x sample buffer and boiled for 5 min, followed by immunoblotting with an anti-GFP antibody.

Chloroplast Purification
To isolate chloroplasts, protoplasts were gently lysed in ice-cold HMS buffer (330 mM sorbitol, 50 mM HEPEs/KOH, pH 7.6, and 3 mM MgCl2). The lysed mixtures were overlaid onto two silica sol gradients (Percoll) (2 mL of 80% [v/v] Percoll in 330 mM sorbitol, 50 mM HEPES/KOH, pH 7.6; 4 mL of 40% [v/v] Percoll in 330 mM sorbitol, 50 mM HEPES/KOH, pH 7.6) and centrifuged for 5 min at 3000g. After centrifugation, intact chloroplasts were carefully collected, transferred to 10 mL of HMS buffer, and pelleted by centrifugation for 5 min at 3000g. Chloroplast extracts were prepared from the isolated chloroplasts and used for immunoblotting with an anti-GFP antibody.

Global Alignment with Dual Gap Penalties
We developed a novel global alignment algorithm between a protein sequence and sequence motifs that was based on dynamic programming, guaranteeing that the optimal alignment maintained motif block structures. The key concept was to allow dual gap penalties such that a relatively higher gap penalty was given within a motif block than between blocks. We denoted two protein sequences of lengths m and n by a = a1a2...am and b = b1b2...bn, where ai and bj represent amino acids of protein sequence a and b, respectively. The second sequence b that represents a motif set comprising p motif blocks was constructed by concatenating p motif blocks in order. For example, the sequence motifs of RbcS are represented by the sequence MLLKSSFPRKCMQVWKKFET. Therefore, sequence b could also be denoted b = Formula, where Formula represents the first motif block of b (i.e., Formula). Note that jp = n. We then constructed a matrix Formula whose (i, j) element, F(i, j), is the score of the best alignment between a1:i and b1:j. Initializing F(1,1) = 0, the matrix F was built up recursively, filled from top left to bottom right, through the following recursion:

Formula
or

Formula
where s(ai,bj) is the score for the match between ai and bj, and g1, g2 are dual gap penalties such that g1 > g2. For boundary conditions, we defined F(I + 1,1) = –ig2 and F(1, j + 1) = –jg2 for the first column and first row of F. We identified a global alignment by tracing the choices that led to the final value of F(m + 1, n + 1). In this study, the substitution matrix was BLOSUM 50 and the dual gap penalties g1, g2 were set to 5 and 0, respectively.

Prediction of Sequence Motifs in Plastid Transit Peptides
Discriminative motif discovery takes a query sequence as input and predicts putative sequence motifs. The sequences from plastid transit peptide peptides (positive set) and the N termini of nonplastid proteins (negative set) were used as a training data set (see Supplemental Figure 6 online). Our model was designed to identify an ordered set of sequence motifs that maximally discriminated between the positive and negative sets. The model had two steps: sequentially adding or deleting a motif segment into a motif set. The selection criterion was the upper bound of misclassification error derived from the Bayesian decision theory (Duda et al., 2001Go).

In the forward motif selection step, the model selected each motif segment with a length of five amino acids by scanning the query sequence. It then converted each motif sequence into the three-dimensional feature vector. The first two components of the feature vector were the scores of the global and local alignments between the query sequence and the sequence of the feature vector (Needleman and Wunsch, 1970Go; Smith and Waterman, 1981Go). The third component was the score of a global alignment with dual gap penalties between the chosen motif segment and the sequence of feature vector. All the training data points were then projected into the one-dimensional space using Fisher's linear discriminant. We modeled the class-conditional densities of the projected data using normal distributions and derived the upper bound of the misclassification error. After obtaining error bounds of all possible motif segments, we chose the one with minimum error bound. We continued this forward motif selection until the number of motif segments in the motif set reached a predefined number (6).

Next, the backward motif elimination step deleted unlikely motif segments chosen from forward motif selection. Specifically, motif segments that produced the smallest error bound value were eliminated at each iteration. Iteration was stopped when the smallest value of the error bound was larger than the previous one.

Prediction of Plastid Transit Peptides
Plastid transit peptide prediction had three steps (see Supplemental Figure 7 online). First, a clustering module partitioned the plastid transit peptides in the training data into groups with common sequence motifs. Second, a feature extraction module converted each sequence of the training data into a corresponding feature vector based on the representative sequences chosen from the clustering module. Finally, a classification module learned a support vector machine (SVM) classifier where the input vector was that obtained from the feature extraction module. The learned SVM classifier takes an input sequence and predicts plastid transit peptides.

In the clustering step, we started with a randomly selected plastid transit peptide and predicted its sequence motifs through a model for discriminative motif discovery. We used motif alignment scores to group sequences sharing similar motifs. We assumed that sequences belonged to the same group when the P value of its motif alignment score was below a predefined cutoff (P < 0.001). We then selected another plastid transit peptide with a maximum P value from the positive set where the grouped sequences had been removed. We repeated these processes until the positive set was empty. The output of the clustering module was a set of representative sequences with predicted sequence motifs.

The feature extraction step converted each sequence into a corresponding feature vector with the empirical kernel map technique using alignment scores (Kim et al., 2006bGo). The feature vector components were the alignment scores between the input and representative sequences. To further enhance prediction performance, we also randomly selected a set of representative sequences from the negative set of the training data.

In the last step, an SVM classifier learned to predict the binary class label of an input sequence (Cristianini and Shawe-Taylor, 2000Go). We used the linear kernel function where the input vector was from the feature extraction module. For testing, the input sequence was first converted into the corresponding feature vector through the empirical kernel map (see Supplemental Methods online).

Preparation of 986 Arabidopsis Protein Sequences
Sequences of all nuclear-encoded proteins of Arabidopsis were downloaded from the SWISS-PROT database, release 51.4 (UniProt Consortium, 2007Go), with the exception of Cab, which was retrieved from the GenBank database (Benson et al., 2007Go). Chloroplast transit peptides were obtained from literature. Sequences confirmed by the SUBA database (Heazlewood et al., 2005Go) were also included, where only plastid proteins supported from GFP analysis were considered. Nonchloroplast proteins were extracted by collecting all protein sequences that were experimentally verified in the CC and FT line. All sequences containing ambiguous amino acids, such as B, Z, or X, or with a length of ≤80 amino acids were excluded. The resulting data set consisted of 208 chloroplast transit peptides and 778 nonchloroplast protein sequences, for a total of 986 protein sequences.

Accession Numbers
Sequence data from this article can be found in the Arabidopsis Genome Initiative or GenBank/EMBL data libraries under the following accession numbers: AAG40356 (RbcS), 2208465A (BCCP), NP_001078288 (Cab), NP_178207 (DnaJ-J8), NP_200230 (PORA), Q94FY7 (TOCC), and NP_181655 (GLU2).

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

Supplemental Figure 1. A Dendogram Constructed by Hierarchical Clustering of 208 Chloroplast Transit Peptides in Arabidopsis thaliana.
Supplemental Figure 2. Fractionation Analysis of Selected Alanine Substitution Mutants on a Percoll Gradient.
Supplemental Figure 3. In Vitro Translation of Alanine Substitution Mutants.
Supplemental Figure 4. The Performance of the Motif Discovery Algorithm Was Degraded by Removal of Sequences Downstream of Predicted Cleavage Sites.
Supplemental Figure 5. GFP Reporter Constructs with the N-Terminal Transit Peptide Region Upstream of the Cleavage Sites Were Not Efficiently Imported into Chloroplasts.
Supplemental Figure 6. Plastid Transit Peptide Motif Prediction Scheme.
Supplemental Figure 7. Overview of Plastid Transit Peptide Prediction.
Supplemental Figure 8. Gaussian Distribution Assessment of Alignment Score Significance.
Supplemental Table 1. The Nucleotide Sequences of Primers Used to Construct Various Alanine Substitution Mutants.
Supplemental Methods. Assessment of Alignment Significance, Description of Prediction of Plastid Transit Peptide Sequence Motifs, Description of Plastid Transit Peptide Prediction, Performance Evaluation of Predicting Chloroplast Transit Peptides, and Performance Comparisons of Predicting Chloroplast Transit Peptides.
Supplemental Data Set 1. List of 208 Plastid Proteins with Their Predicted Sequence Motifs.
Supplemental Data Set 2. The Alignments between Plastid Transit Peptides and the Sequence Motifs of the Seven Representative Plastid Proteins.


    Acknowledgments
 
This work was supported by grants from the Systems Bio-Dynamics Center (R15-2004-033-05002-0), the Ministry of Education, Science, and Technology (I.H.), the Agricultural Research and Planning Center of the Ministry of Agriculture, Fishery, and Food (I.H.), BioGreen 21 Program (20070401-034-026-008-03-00), Rural Development Administration (I.H.), the POSTECH Core Research Program (I.H.), the Korea Science and Engineering Foundation (R15-2004-033-07002-0) (S.K.), and the Korea Research Foundation (KRF-2006-311-C00503) (S.K.) and by a Microsoft Research Asia fellowship (J.K.K.).


    Footnotes
 
1 These authors contributed equally to this work. Back

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: Inhwan Hwang (ihhwang{at}postech.ac.kr).

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

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

Received May 5, 2008; Revision received May 5, 2008. accepted May 23, 2008.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 RESULTS
 DISCUSSION
 METHODS
 REFERENCES
 
Arimura, S., Takusagawa, S., Hatano, S., Nakazono, M., Hirai, A., and Tsutsumi, N. (1999). A novel plant nuclear gene encoding chloroplast ribosomal protein S9 has a transit peptide related to that of rice chloroplast ribosomal protein L12. FEBS Lett. 450: 231–234.[CrossRef][Web of Science][Medline]

Becker, T., Jelic, M., Vojta, A., Radunz, A., Soll, J., and Schleiff, E. (2004). Preprotein recognition by the Toc complex. EMBO J. 23: 520–530.[CrossRef][Web of Science][Medline]

Bédard, J., and Jarvis, P. (2005). Recognition and envelope translocation of chloroplast preproteins. J. Exp. Bot. 56: 2287–2320.[Abstract/Free Full Text]

Benson, D.A., Karsch-Mizrachi, I., Lipman, D.J., Ostell, J., and Wheeler, D.L. (2007). GenBank. Nucleic Acids Res. 35: D21–D25.[Abstract/Free Full Text]

Bhushan, S., Kuhn, C., Berglund, A.K., Roth, C., and Glaser, E. (2006). The role of the N-terminal domain of chloroplast targeting peptides in organellar protein import and miss-sorting. FEBS Lett. 580: 3966–3972.[CrossRef][Web of Science][Medline]

Bruce, B.D. (2000). Chloroplast transit peptides: Structure, function and evolution. Trends Cell Biol. 10: 440–447.[CrossRef][Web of Science][Medline]

Clark, S.E., and Lamppa, G.K. (1991). Determinants for cleavage of the chlorophyll a/b binding protein precursor: A requirement for a basic residue that is not universal for chloroplast imported proteins. J. Cell Biol. 114: 681–688.[Abstract/Free Full Text]

Comai, L., Larson-Kelly, N., Kiser, J., Mau, C.J., Pokalsky, A.R., Shewmaker, C.K., McBride, K., Jones, A., and Stalker, D.M. (1988). Chloroplast transport of a ribulose bisphosphate carboxylase small subunit-5-enolpyruvyl 3-phosphoshikimate synthase chimeric protein requires part of the mature small subunit in addition to the transit peptide. J Biol Chem. 263: 15104–15109.[Abstract/Free Full Text]

Cristianini, N., and Shawe-Taylor, J. (2000). An Introduction to Support Vector Machines and Other Kernel-Based Learning Methods. (Cambridge, UK: Cambridge University Press).

Duda, R.O., Hart, P.E., and Stork, D.G. (2001). Pattern Classification. (New York: John Wiley & Sons).

Emanuelsson, O., Nielsen, H., Brunka, S., and von Heijne, G. (2000). Predicting subcellular localization of proteins based on their N-terminal amino acid sequence. J. Mol. Biol. 300: 1005–1016.[CrossRef][Web of Science][Medline]

Emanuelsson, O., Nielsen, H., and von Heijne, G. (1999). ChloroP, a neural network-based method for predicting chloroplast transit peptides and their cleavage sites. Protein Sci. 8: 978–984.[Web of Science][Medline]

Emanuelsson, O., and von Heijne, G. (2001). Prediction of organellar targeting signals. Biochim. Biophys. Acta 1541: 114–119.[Medline]

Gantt, J.S., Baldauf, S.L., Calie, P.J., Weeden, N.F., and Palmer, J.D. (1991). Transfer of rpl22 to the nucleus greatly preceded its loss from the chloroplast and involved the gain of an intron. EMBO J. 10: 3073–3078.[Web of Science][Medline]

Gutensohn, M., Schulz, B., Nicolay, P., and Flugge, U.I. (2000). Functional analysis of the two Arabidopsis homologues of Toc34, a component of the chloroplast protein import apparatus. Plant J. 23: 771–783.[CrossRef][Web of Science][Medline]

Hand, J.M., Szabo, L.J., Vasconcelos, A.C., and Cashmore, A.R. (1989). The transit peptide of a chloroplast thylakoid membrane protein is functionally equivalent to a stromal-targeting sequence. EMBO J. 8: 3195–3206.[Web of Science][Medline]

Heazlewood, Z.L., Tonti-Filippini, J., Verboom, R.E., and Millar, A.H. (2005). Combining experimental and predicted datasets for determinination of the subcellular location of proteins in Arabidopsis. Plant Physiol. 139: 598–609.[Abstract/Free Full Text]

Hinnah, S.C., Wagner, R., Sveshnikova, N., Harrer, R., and Soll, J. (2002). The chloroplast protein import channel Toc75: Pore properties and interaction with transit peptides. Biophys. J. 83: 899–911.[Web of Science][Medline]

Jarvis, P., and Soll, J. (2002). Toc, tic, and chloroplast protein import. Biochim. Biophys. Acta 1590: 177–189.[Medline]

Jin, J.B., Kim, Y.A., Kim, S.J., Lee, S.H., Kim, D.H., Cheong, G.W., and Hwang, I. (2001). A new dynamin-like protein, ADL6, is involved in trafficking from the trans-Golgi network to the central vacuole in Arabidopsis. Plant Cell 13: 1511–1526.[Abstract/Free Full Text]

Kessler, F., and Schnell, D.J. (2004). Chloroplast protein import: Solve the GTPase riddle for entry. Trends Cell Biol. 14: 334–338.[CrossRef][Web of Science][Medline]

Kim, D.H., Eu, Y.J., Yoo, C.M., Kim, Y.W., Pih, K.T., Jin, J.B., Kim, S.J., Stenmark, H., and Hwang, I. (2001). Trafficking of phosphatidylinositol 3-phosphate from the trans-Golgi network to the lumen of the central vacuole in plant cells. Plant Cell 13: 287–301.[Abstract/Free Full Text]

Kim, J.K., Bang, S.Y., and Choi, S. (2006b). Sequence driven features for prediction of subcellular localization of proteins. Pattern Recognit. 39: 2301–2311.[CrossRef][Web of Science]

Kim, J.K., Raghava, G.P.S., Bang, S.Y., and Choi, S. (2006a). Prediction of subcellular localization of proteins using pairwise sequence alignment and support vector machine. Pattern Recognit. Lett. 27: 996–1001.[CrossRef]

Kohler, R.H., Cao, J., Zipfel, W.R., Webb, W.W., and Hanson, M.R. (1997). Exchange of protein molecules through connections between higher plant plastids. Science 276: 2039–2042.[Abstract/Free Full Text]

Lee, D.W., Lee, S., Lee, G.J., Lee, K.H., Kim, S., Cheong, G.W., and Hwang, I. (2006). Functional characterization of sequence motifs in the transit peptide of Arabidopsis small subunit of rubisco. Plant Physiol. 140: 466–483.[Abstract/Free Full Text]

McFadden, G.I. (1999). Endosymbiosis and evolution of the plant cell. Curr. Opin. Plant Biol. 2: 513–519.[CrossRef][Web of Science][Medline]

Millar, A.H., Whelan, J., and Small, I. (2006). Recent surprises in protein targeting to mitochondria and plastids. Curr. Opin. Plant Biol. 9: 610–615.[CrossRef][Web of Science][Medline]

Miller, E.A., Beilharz, T.H., Malkus, P.N., Lee, M.C., Hamamoto, S., Orci, L., and Schekman, R. (2003). Multiple cargo binding sites on the COPII subunit Sec24p ensure capture of diverse membrane proteins into transport vesicles. Cell 114: 497–509.[CrossRef][Web of Science][Medline]

Needleman, S.B., and Wunsch, C.D. (1970). A general method applicable to the search for similarities in the amino acid sequence of two proteins. J. Mol. Biol. 48: 443–453.[CrossRef][Web of Science][Medline]

Reinbothe, S., Runge, S., Reinbothe, C., van Cleve, B., and Apel, K. (1995). Substrate-dependent transport of the NADPH:protochlorophyllide oxidoreductase into isolated plastids. Plant Cell 7: 161–172.[Abstract]

Rensink, W.A., Pilon, M., and Weisbeek, P. (1998). Domains of a transit sequence required for in vivo import in Arabidopsis chloroplasts. Plant Physiol. 118: 691–699.[Abstract/Free Full Text]

Reumann, S., Davila-Aponte, J., and Keegstra, K. (1999). The evolutionary origin of the protein-translocating channel of chloroplastic envelope membranes: Identification of a cyanobacterial homolog. Proc. Natl. Acad. Sci. USA 96: 784–789.[Abstract/Free Full Text]

Rial, D.V., Arakaki, A.K., and Ceccarelli, E.A. (2000). Interaction of the targeting sequence of chloroplast precursors with Hsp70 molecular chaperones. Eur. J. Biochem. 267: 6239–6248.[Web of Science][Medline]

Schein, A.I., Kissinger, J.C., and Ungar, L.H. (2001). Chloroplast transit peptide prediction: a peek inside the black box. Nucleic Acids Res. 29: E82.[CrossRef][Medline]

Schleiff, E., and Soll, J. (2005). Membrane protein insertion: Mixing eukaryotic and prokaryotic concepts. EMBO Rep. 6: 1023–1027.[CrossRef][Web of Science][Medline]

Small, I., Peeters, N., Legeai, F., and Lurin, C. (2004). Predotar: A tool for rapidly screening proteomes for N-terminal targeting sequences. Proteomics 4: 1581–1590.[CrossRef][Web of Science][Medline]

Smith, M.D., Rounds, C.M., Wang, F., Chen, K., Afitlhile, M., and Schnell, D.J. (2004). atToc159 is a selective transit peptide receptor for the import of nucleus-encoded chloroplast proteins. J. Cell Biol. 165: 323–334.[Abstract/Free Full Text]

Smith, T.F., and Waterman, M.S. (1981). Identification of common molecular subsequences. J. Mol. Biol. 147: 195–197.[CrossRef][Web of Science][Medline]

UniProt Consortium (2007). The universal protein resource (UniProt). Nucleic Acids Res. 35: D193–D197.[Abstract/Free Full Text]

von Heijne, G., Steppuhn, J., and Herrmann, R.G. (1989). Domain structure of mitochondrial and chloroplast targeting peptides. Eur. J. Biochem. 180: 535–545.[Web of Science][Medline]

Wienk, H.L., Czisch, M., and de Kruijff, B. (1999). The structural flexibility of the preferredoxin transit peptide. FEBS Lett. 453: 318–326.[CrossRef][Web of Science][Medline]

Wienk, H.L., Wechselberger, R.W., Czisch, M., and de Kruijff, B. (2000). Structure, dynamics, and insertion of a chloroplast targeting peptide in mixed micelles. Biochemistry 39: 8219–8282.[CrossRef][Web of Science][Medline]

Zhang, X.P., and Glaser, E. (2002). Interaction of plant mitochondrial and chloroplast signal peptides with the Hsp70 molecular chaperone. Trends Plant Sci. 7: 14–21.[CrossRef][Web of Science][Medline]




This article has been cited by other articles:


Home page
Plant CellHome page
S. Lee, D. W. Lee, Y. Lee, U. Mayer, Y.-D. Stierhof, S. Lee, G. Jurgens, and I. Hwang
Heat Shock Protein Cognate 70-4 and an E3 Ubiquitin Ligase, CHIP, Mediate Plastid-Destined Precursor Degradation through the Ubiquitin-26S Proteasome System in Arabidopsis
PLANT CELL, December 1, 2009; 21(12): 3984 - 4001.
[Abstract] [Full Text] [PDF]


Home page
Plant CellHome page
A. Kitajima, S. Asatsuma, H. Okada, Y. Hamada, K. Kaneko, Y. Nanjo, Y. Kawagoe, K. Toyooka, K. Matsuoka, M. Takeuchi, et al.
The Rice {alpha}-Amylase Glycoprotein Is Targeted from the Golgi Apparatus through the Secretory Pathway to the Plastids
PLANT CELL, September 1, 2009; 21(9): 2844 - 2858.
[Abstract] [Full Text] [PDF]


Home page
Plant Physiol.Home page
D. W. Lee, S. Lee, Y. J. Oh, and I. Hwang
Multiple Sequence Motifs in the Rubisco Small Subunit Transit Peptide Independently Contribute to Toc159-Dependent Import of Proteins into Chloroplasts
Plant Physiology, September 1, 2009; 151(1): 129 - 141.
[Abstract] [Full Text] [PDF]


Home page
RNAHome page
D. Karcher and R. Bock
Identification of the chloroplast adenosine-to-inosine tRNA editing enzyme
RNA, July 1, 2009; 15(7): 1251 - 1257.
[Abstract] [Full Text] [PDF]


Home page
Plant Physiol.Home page
S. Huang, N. L. Taylor, J. Whelan, and A. H. Millar
Refining the Definition of Plant Mitochondrial Presequences through Analysis of Sorting Signals, N-Terminal Modifications, and Cleavage Motifs
Plant Physiology, July 1, 2009; 150(3): 1272 - 1285.
[Abstract] [Full Text] [PDF]


Home page
Plant CellHome page
S. Kikuchi, M. Oishi, Y. Hirabayashi, D. W. Lee, I. Hwang, and M. Nakai
A 1-Megadalton Translocation Complex Containing Tic20 and Tic21 Mediates Chloroplast Protein Import at the Inner Envelope Membrane
PLANT CELL, June 1, 2009; 21(6): 1781 - 1797.
[Abstract] [Full Text] [PDF]


Home page
Plant CellHome page
A. Jamai, P. A. Salome, S. H. Schilling, A. P.M. Weber, and C. R. McClung
Arabidopsis Photorespiratory Serine Hydroxymethyltransferase Activity Requires the Mitochondrial Accumulation of Ferredoxin-Dependent Glutamate Synthase
PLANT CELL, February 1, 2009; 21(2): 595 - 606.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplemental Data
Right arrow All Versions of this Article:
20/6/1603    most recent
tpc.108.060541v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via CrossRef
Right arrow Citing Articles via Web of Science (5)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Lee, D. W.
Right arrow Articles by Hwang, I.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Lee, D. W.
Right arrow Articles by Hwang, I.
Agricola
Right arrow Articles by Lee, D. W.
Right arrow Articles by Hwang, I.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
ASPB Publications THE PLANT CELL PLANT PHYSIOLOGY
Copyright © 2008 by the American Society of Plant Biologists