Previous Article | Next Article ![]()
Journal of Bacteriology, February 2006, p. 1405-1410, Vol. 188, No. 4
0021-9193/06/$08.00+0 doi:10.1128/JB.188.4.1405-1410.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Institute of Genetics, Martin-Luther-University, 06099 Halle, Germany
Received 27 September 2005/ Accepted 21 November 2005
|
|
|---|
-helical nature of the TTS signal. This is the first report of a surface structure from a plant-pathogenic bacterium that evolved under the constraint of positive selection and hints to the evolutionary adaptation of this extracellular appendage to avoid recognition by the plant defense surveillance system. |
|
|---|
Our laboratory studies the plant pathogen Xanthomonas campestris pathovar vesicatoria, which is the causal agent of bacterial spot disease in pepper and tomato plants. The TTS system of X. campestris pv. vesicatoria is encoded by the 23-kb hrp (hypersensitive response and pathogenicity) gene cluster, which is essential to cause disease in susceptible plants as well as to induce the hypersensitive response in resistant plants (9). A key component of the X. campestris pv. vesicatoria pathogenicity machinery is the Hrp pilus, which is thought to serve as a conduit for protein translocation between the bacterium and the plant host cell (30). The Hrp pilus is a filamentous structure that extends up to several micrometers from the bacterial cell surface and is mainly composed of a large number of identical pilin subunits, the 10-kDa HrpE protein, which is unique for the genus Xanthomonas (30). HrpE itself is a TTS substrate and harbors a secretion signal in the first 17 N-terminal amino acids (29). Sequence alignment and secondary structure predictions of five hrpE sequences showed that the very-C terminus is entirely conserved and predominantly
-helical, whereas the N terminus is hypervariable and mostly devoid of regular secondary structure. Reporter protein fusions and pentapeptide insertions demonstrated that HrpE contains two functional domains, a surface-exposed domain, including the TTS signal, and the polymerization domain (29). The sequence hypervariability of the surface-exposed domain prompted us to analyze the evolutionary rate of 26 hrpE sequences, representing Xanthomonas strains belonging to 14 different pathovars. Our results reveal different modes of selection acting on different portions of the protein. The N terminus of HrpE undergoes positive selection, whereas the C terminus is subjected to purifying selection. Since the positively selected parts are exposed at the surface of the assembled Hrp pilus, the plant pathogen Xanthomonas thus generates surface diversity, most probably to evade recognition by the plant perception and defense systems. Furthermore, we propose an amphipathic
-helix as the TTS signal of HrpE which is on one side highly conserved and polar and on the opposite side positively selected and hydrophobic.
|
|
|---|
Sequence data. The coding regions of 21 hrpE alleles were amplified using PCR primers derived from conserved regions of hrpD6 (hrpEup, 5'-CAACATCGTCGCGCACATCGC) and hpaB (hrpEdown, 5'-TCGAATCGCGCGCTGCTCAT) located upstream and downstream of the hrpE coding sequence, respectively. For colony PCR, a single colony of each strain was used as the source of template DNA. PCR conditions were an initial 95°C step for 300 s, followed by 35 cycles of 95°C for 30 s, 58°C for 30 s, and 72°C for 40 s, and a final step at 72°C for 10 min. PCR products were cloned into pCR 2.1-TOPO and transformed into DH10b (Invitrogen, Carlsbad, CA). Five clones of each hrpE allele were pooled and sequenced using vector-specific primers (TOPO TA cloning kit; Invitrogen).
Additionally, the hrpE nucleotide sequences of X. campestris pv. vesicatoria strain 75-3 (GenBank accession number [AC] AF056246), X. oryzae pv. oryzae strain KACC 10331 (AC AE013598), X. axonopodis pv. glycines strain 8ra (AC AF499777), X. axonopodis pv. citri strain 306 (AC AE011665), and X. campestris pv. campestris strain ATCC 33913 (AC AE012221) were used for bioinformatic analyses.
Sequence analysis. The nucleotide sequences were translated and aligned using CLUSTAL X (26). The multiple alignment was then manually edited using the BIOEDIT program (11). The amino acid alignments were transposed back to the nucleotide sequence level to gain a codon-based alignment. Mean diversity estimates within sequences were performed with the MEGA 2 software, version 2.1 (15), using the modified Nei-Gojobori method with the Jukes-Cantor corrected distance and the complete-deletion-of-gaps option. Unweighted pair group method with averages (UPGMA) trees were constructed using PAUP 3.1.1 (25). Distances were determined by using the Kimura two-parameter model. Support for nodes was assessed with a bootstrap confidence level using 1,000 replicates.
Calculation of synonymous and nonsynonymous mean pairwise diversities.
The method of Comeron, as implemented in the K-Estimator program, version 6.1 (5), was applied to test whether positive or purifying selection has taken place on hrpE. K-Estimator takes into account the divergence values, number of nucleotides or codons, the transition/transversion substitution ratio (
/ß), the amino acid composition, and the G+C content at the third codon position. The confidence intervals for Ka and Ks were obtained by Monte Carlo simulations (5). Synonymous mean pairwise diversity per synonymous site (Ks) and nonsynonymous mean pairwise diversity per nonsynonymous site (Ka) were calculated with the complete-deletion-of-gaps option and a sliding window of 15 nucleotides with a step size of 9 nucleotides.
Nucleotide sequence accession numbers. The nucleotide sequences for the 21 hrpE alleles were submitted to GenBank (accession numbers DQ286679 to DQ286699).
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Sequence variations of five genes among five Xanthomonas strainsa
|
![]() View larger version (27K): [in a new window] |
FIG. 1. UPGMA tree of hrpE sequences from 26 Xanthomonas strains. Numbers on branches indicate the number of times (percent) that the node was supported by 1,000 replicates of the bootstrap analysis (only values greater than 80% are shown).
|
-helix (30). A helical wheel plot revealed that the conserved polar residues are located on one side of the helix, whereas the opposite side consists of hydrophobic amino acids (Fig. 3). These findings strongly suggest that the HrpE secretion signal is formed by an amphipathic
-helix.
![]() View larger version (52K): [in a new window] |
FIG. 2. Nonredundant amino acid sequence alignment of HrpE. HrpE sequences of group I and II xanthomonads were aligned using CLUSTAL X (26). Residues which are identical in at least 90% of the sequences are shown in white on a black background. Proposed functional domains are indicated on top of the alignment (29).
|
![]() View larger version (38K): [in a new window] |
FIG. 3. Type III secretion signal of the X. campestris pv. vesicatoria 85-10 HrpE in a helical wheel representation. Conserved amino acids and hydrophobic residues are shaded in gray and black, respectively. The N-terminal proline residue is oriented to the bottom.
|
Modes of selection operating on hrpE. Using the enlarged set of hrpE sequences, we investigated whether certain parts of HrpE show different modes of selection. A sliding window analysis of Ka and Ks values was performed. Since the first amino acid of HrpE from group I and the first four amino acids of HrpE from group II give no substantial information, they were not considered. Windows in which Ka exceeded Ks (one tailed Z-test, P < 0.05) (Fig. 4A) were found in the N-terminal and central regions of the protein. This finding shows that diversifying selection had occurred at these regions. In contrast, Ks exceeded Ka in windows 98 and 215 to 266 (one tailed Z-test, P < 0.05 [Fig. 4A]; the number indicates the central nucleotide of the window), indicating that these regions evolved under purifying selection.
![]() View larger version (25K): [in a new window] |
FIG. 4. Mean pairwise diversity plots of Ks and Ka values. The plots were constructed for windows of 15 nucleotides sliding along the hrpE sequence alignment. The first codon of group I and the first four codons of group II sequences are not included. (A) Comparison of all hrpE sequences. (B) Group I sequence comparison. (C) Group II sequence comparison. Filled boxes indicate the windows (step size of 9 nucleotides) in which nonsynonymous differences significantly exceed synonymous differences (one-tailed Z test, P < 0.05). White boxes indicate positions in the sequence where the reverse condition applies (Ks > Ka; P < 0.05) or sequences which are identical. Solid line, mean Ka; dashed line, mean Ks.
|
The positively selected positions are surface exposed. Certain modes of selection act often highly localized within genes. For many plant resistance proteins, plant chitinases, and the outer membrane porin PorB of Neisseria meningitidis, positive selection has been detected primarily for solvent-exposed amino acids (3, 7, 28). Since solvent-exposed parts of a protein are in general rather hydrophilic and core regions of proteins are more hydrophobic, we performed hydrophobicity plot analyses of groups I and II and compared the localizations of positively selected and hydrophilic regions. As shown in Fig. 5, the HrpE regions which are positively selected in group I or II correspond to hydrophilic regions. Of special interest is the region between amino acids 25 and 35, which is hydrophilic in group I and hydrophobic in group II. This different physico-chemical character is perfectly reflected by different modes of selection, indicating that this region might be surface exposed and therefore positively selected in group I. Notably, for the C-terminal hydrophilic region an indication of positive selection is missing. This finding may be related to the proposed function of this domain in the polymerization process (29).
![]() View larger version (20K): [in a new window] |
FIG. 5. Positively selected parts of HrpE are surface exposed. Mean hydrophobicity plots of group I (dashed line) and group II (solid line) HrpE proteins were calculated over a sliding window of five amino acid residues, using the Kyte and Doolittle hydrophobicity scale (16). Values below 0.4 are considered hydrophilic. Positively selected parts in groups I and II are indicated by open and filled boxes, respectively.
|
|
|
|---|
-helical secondary structure, and sequence hypervariability (14). Interestingly, sequence hypervariability is not only observed between different species but also between pathovars of the same species (12). In this study, we analyzed the sequence variability of the Xanthomonas Hrp pilus subunit HrpE, showed that the observed sequence hypervariability is caused by positive selection, and speculated about the biological significance of hypervariability. A sliding window analysis of 18 hrpE alleles revealed several regions in the N-terminal two-thirds in which Ka exceeds Ks, thus indicating that these regions are positively selected. In contrast, the far C terminus is almost identical or subjected to purifying selection in all Xanthomonas sequences. Positive selection has been reported as a beacon for proteins involved in genetic conflict and has been documented in genes that encode surface proteins of pathogen (2, 28) and host defense systems, such as the human major histocompatibility complex (13), plant chitinases (3), and plant resistance genes (7). The correlation between positive selection and host-pathogen interactions is particularly strong. For example, a GenBank survey uncovered remarkably few sequences (0.45%) evolving under positive selection, but more than half of these sequences were involved in host-pathogen interactions (8). Only six examples of bacterial surface proteins which have undergone positive selection have been described, but all of them originated from animal pathogens (27). This is the first report of a surface protein from a plant pathogen that is subjected to positive selection. This points to an unexpected analogy between animal and plant pathogenic bacteria in their strategy to avoid recognition by the host.
Despite the lack of an adaptive immune system, plants possess recognition capacities for several microbial compounds often indispensable for the microbial lifestyle. This set of general elicitors are referred to as pathogen-associated molecular patterns (PAMPs) and include the eubacterial flagellin, lipopolysaccharide of gram-negative bacteria, fungal cell wall-derived glucans, chitins, mannans, and surface proteins of microbial pathogens (20). An Hrp pilus of 500 nm in length would consist of approximately 1,000 HrpE subunits. As an abundant key component of the TTS system which is crucial for pathogenicity, HrpE would be a prime candidate to act as a PAMP. However, there is no report of a plant species or cultivar which has evolved a mechanism to detect HrpE or Hrp pilins from other plant pathogens as a PAMP.
Diversification of the Hrp pilus surface by positive selection acting on surface-exposed parts of the protein would make perfect sense in order to conceal the attack by Xanthomonas spp. Such a surface-limited diversification is a well-known strategy of human pathogenic bacteria (2, 28). Indeed, in both phylogenetic groups of Xanthomonas a remarkable correlation between hydrophilic regions and positive selection was found. Notably, only the first four hydrophilic stretches of the HrpE sequence were positively selected, whereas the C-terminal hydrophilic stretch was not. These findings are in perfect agreement with the proposed domain structure of HrpE (29). In this previous work it was found that the N-terminal two-thirds of the protein are tolerant to pentapeptide insertions, whereas the C-terminal third is not. Moreover, we isolated four dominant-negative insertions in the C-terminal third, thus supporting our hypothesis that this domain is conserved due to its function in polymerization. In line with this model, the C-terminal hydrophilic stretch would be hidden within the assembled Hrp pilus, and varying this sequence would give no selective advantage. In conclusion, the N-terminal surface-exposed domain functions as a camouflage suit for the highly conserved assembly domain at the C terminus.
An alternative mechanism to achieve sequence variation is based on a high frequency of recombination, as has been reported for components of the outer membrane or cell surface appendages of mammalian pathogens, such as the outer membrane protein PorB and the pilus protein PilE of N. meningitidis and the flagellin FliC of pathogenic E. coli (2, 22, 28). However, we did not observe discrepancies in the phylogenies between the full-length sequence and the separately analyzed N-terminal and C-terminal halves, indicating that recombination is not responsible for the diversification of hrpE. This finding is supported by the absence of a second copy of hrpE-related sequences in the completely sequenced genomes of five Xanthomonas strains. Interestingly however, the plant pathogens P. syringae and R. solanacearum possess a second copy of their Hrp pilin genes. It has been speculated that these second copies might act as capping or scaffolding proteins (12). Alternatively, one might speculate about a possible role of these extra genes in surface variation.
The very-N terminus of HrpE was shown to serve as a TTS signal (29). Despite 15 years of intense research, and even with the three-dimensional structures of the N termini of TTS substrates at hand, the nature of the TTS signal is still enigmatic (24). The fact that the N termini of these substrates do not share any obvious consensus sequence at the primary sequence level makes it difficult to understand how TTS substrates are recognized and secreted even by heterologous TTS machineries (24). Comparisons of Yersinia TTS substrates indicated that a more general property serves as secretion signal. Yersinial secretion signals consist of an amphipathic series of randomly alternating polar and hydrophobic residues, and even synthetic amphipathic signals consisting of serine and isoleucine were able to mediate secretion of a reporter protein (18). Similar biophysical features were identified in TTS substrates from the plant pathogen P. syringae (21). However, amphipathic N-terminal sequences derived from cytosolic components of the yersinial Yop regulon did not promote secretion of the reporter protein YopE (17). These findings imply that, in addition to the amphipathic character, there must be a sequence dependence on N-terminal secretion signals that is currently not understood.
While most studies on TTS signals have compared the N termini of distantly related or unrelated TTS substrates or have introduced point mutations into the N terminus of TTS substrates (21, 24), we analyzed a set of naturally occurring HrpE variants. From previous work it is known that the HrpE TTS signal is located in the first 17 amino acids (29). Since all HrpE N termini were obtained from pathogens, it is reasonable to assume that they are all functional TTS signals. The trade-off between sequence diversification (to avoid host recognition) and functional conservation (to allow secretion) operating at the HrpE N terminus offered a unique advantage to identify residues and principles that are of crucial importance for function as a TTS signal. The
-helix that would be formed by the amphipathic N-terminal sequence displays a polar and a hydrophobic side. In addition to the conserved amphipathic character of the
-helix, the five amino acid residues at the polar side seem to be critical for function, as indicated by their strong conservation. Probably, these two features are crucial for function of the HrpE TTS signal. To what extent this rule may apply to other TTS substrates remains to be clarified in future studies.
E.W. was supported by a GradFG fellowship from Sachsen-Anhalt.
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»