Previous Article | Next Article ![]()
Journal of Bacteriology, January 2006, p. 275-286, Vol. 188, No. 1
0021-9193/06/$08.00+0 doi:10.1128/JB.188.1.275-286.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Center for Immunobiology and Vaccine Research, Children's Hospital Oakland Research Institute, 5700 Martin Luther King Jr. Way, Oakland, California 94609,1 Centro de Bacteriologia, Instituto Nacional de Saúde, Av. Padre Cruz, 1649-016 Lisboa, Portugal,2 T-10 Theoretical Biology and Biophysics, MS-K710 Los Alamos National Laboratory, Los Alamos, New Mexico 87545,3 University of California at San Francisco School of Medicine, San Francisco, California 941434
Received 25 July 2005/ Accepted 6 October 2005
|
|
|---|
|
|
|---|
Seven chlamydial genomes have been sequenced over the last 7 years, revealing a polymorphic membrane protein (Pmp) family unique to Chlamydiales that likely plays an important role in chlamydial biology and disease pathogenesis (46). There are 17 pmp genes identified for Chlamydophila caviae (formerly Chlamydia psittaci; strain GPIC) (36), 21 for Chlamydia pneumoniae (strain CWL029) (21), and 9 for Chlamydia trachomatis (serovar D/UW-3 and mouse pneumonitis strain Nigg) (35, 42). The pmp genes are predicted to be localized in the outer membrane, due to the presence of a C-terminal phenylalanine, a cleavable peptide leader sequence (15, 18, 39, 44), and the fact that all Pmps have been detected by proteomics as membrane constituents, at least for serovar L2 (40, 41, 45). The C. trachomatis Pmps account for 3.15% of the organism's coding capacity (18) and are immunogenic for humans (16, 17, 20, 23). By protein structure analysis, they are predicted to be autotransporters (19, 47) that mediate the translocation of the N terminus to the bacterial surface (19, 47). They also contain multiple GGAI motifs, which have been associated in other organisms with adhesion to the host cell (18).
We recently performed statistical modeling and phylogenetic analyses of 18 reference serovars and genovariant Ja for pmpC (15) and included the partial sequences of pmpE, pmpH, and pmpI for 15 serovars that were available from another study (43). pmpC and pmpH were found to group serovars according to cell tropism properties, including a distinct clade for the most prevalent serovars, E and F. pmpC also revealed remnants of two putative insertion sequence (IS)-like elements differentially present in some serovars, suggesting a role in genome plasticity (15). Because of these findings, the goal of the present study was to sequence the remaining five pmp genes (pmpA, -B, -D, -F, and -G) for all 18 C. trachomatis serovars and genovariant Ja to completion and provide a detailed genomics and molecular analysis, including phylogenetics of the entire pmp family. We also identified specific protein regions that differentiate disease groups, as these may be responsible for significant tissue tropism properties of the organism.
|
|
|---|
Serovar confirmation of each reference strain was indirectly performed by ompA sequencing and further nucleotide comparison with available GenBank data as previously described (15). Briefly, DNA was extracted and amplified using 1x OptiBuffer (Bioline, London, United Kingdom), 2.8 mM MgCl2, 200 µM deoxynucleoside triphosphates (Promega, Madison, WI), 25 pmol of each primer, 1.5 U Bio-X-Act Proof-reading DNA Polymerase (Bioline), and 5 µl of strain DNA, in a final volume of 25 µl. PCR was run in a GeneAmp PCR System 9700 (Applied Biosystems, Foster City, CA) with the following thermocycling profile: 1 min at 95°C; 1 min at 50°C; 30 cycles, each consisting of 1 min at 70°C, 30 s at 95°C, and 30 s at 50°C; and a final elongation step of 10 min at 70°C. The amplicons were gel purified and sequenced as previously described (15) using the BigDye Terminator Cycle Sequencing kit, version 1.1, according to the manufacturer's instructions (Applied Biosystems) in an ABI 3700 automated sequencer (Applied Biosystems). The sequences were compared with available ompA chlamydial sequences from GenBank using MegAlign software (DNASTAR, Madison, WI) to identity and confirm each ompA genotype.
Amplification and sequencing of pmpA, pmpB, pmpD, pmpF, and pmpG for each serovar. PCR primers for pmpA, pmpB, pmpD, pmpF, and pmpG (Table 1) were designed with Primer Select software (DNASTAR) based on already available GenBank sequences of the respective pmp genes and intergenomic regions for serovar D (UW-3; GenBank no. AE001273). Due to the large size of pmpB (5,256 bp) and pmpD (4,596 bp), two PCR primer pairs (Table 1) were designed for each to generate two overlapping amplicons for the entire gene. PCR reagents and thermocycling profiles were the same as described for ompA, except for the primers (Table 1), annealing temperature, and extension time: 57°C and 4 min for PmpA-1/PmpA-2, 57°C and 3 min 20 s for PmpB-1/PmpB-2 and PmpB-3/PmpB-4, 56°C and 2 min for PmpD-5/PmpD-6, and 56°C and 4 min for PmpD-9/PmpD-10, PmpF-1/PmpF-2, and PmpG-5/PmpG-6. Amplicons were visualized in ethidium bromide-stained 0.8% agarose gels, purified, and sequenced as described above for ompA.
|
View this table: [in a new window] |
TABLE 1. Oligonucleotide primers used for PCR and sequencing
|
Amplification and sequencing of pmpE, pmpH, and pmpI for serovars Da, Ia, and L2a and genovariant Ja.
Partial sequences of pmpE, pmpH, and pmpI were previously available for 15 serovars (excluding Da, Ia, and L2a and genovariant Ja) (43). These sequences did not include
80 bases at the N terminus and
20 bases at the C terminus for pmpE and pmpI. We had previously sequenced the N terminus of pmpH for the 15 serovars (15). All serovars, including Da, Ia, and L2a, and genovariant Ja were sequenced to completion for these three pmp genes as described above. In addition, prior to the GenBank release of the partial pmpI sequence, our laboratory had sequenced pmpI for all 18 serovars and Ja. Therefore, our sequences for pmpI were used in this study.
PCR reagents and thermocycling profiles were the same as described above except for primers (Table 1), annealing temperature, and extension time: 57°C and 3 min 30 s for PmpE-1/PmpE-2, 57°C and 4 min 30 s for PmpH-3/PmpH-4, and 56°C and 3 min 30 s for PmpI-13/PmpI-14. The sequencing strategy was identical to that described above.
Genomic analysis of the pmp family for all known serovars and genovariant Ja. The genomic analysis of all pmp genes was performed using Molecular Evolutionary Genetics Analysis software, version 2.1 (MEGA; Institute of Molecular Evolutionary Genetics, Pennsylvania State University; http://www.megasoftware.net) as previously described (15). Pairwise comparisons were performed between individual serovar sequences, between group means where sequences were clustered according to disease properties (ocular, genital, and LGV), and between serovar sequences within each of these groups. Briefly, mean genetic distances and pairwise comparisons were based on the number of nucleotide differences including both transitions and transversions where only valid sites were compared (gaps were excluded). The genetic distance evaluation was also done using the p-distance method that calibrates the obtained differences relative to the total number of sites under comparison that varied from 2,637 for pmpI to 5,283 for pmpC. A similar strategy was applied to the respective proteins where distances were evaluated based on all amino acid changes and respective protein lengths. For the protein analysis, the total number of valid sites under comparison varied from 878 for PmpI to 1,758 for PmpC.
To study the molecular evolution of the pmp family, a bootstrap test of phylogeny was conducted where neighbor-joining trees (38) were created using the Kimura two-parameter model (22), which takes into account transitional and transversional substitution rates while assuming that the four nucleotide frequencies are the same and that rates of substitution do not vary among sites (30). Phylogenetic reconstructions were also done at the protein level using the gamma distance model, which takes into account the inequality of the substitution rates among sites (30).
Another molecular evolution study was performed where each member of the pmp family was studied by the Nei-Gojobori method (29), where the ratio of nonsynonymous (dN) to synonymous (dS) substitutions was calculated. The p-distance option was used to normalize the nonsynonymous and synonymous differences against the number of potential nonsynonymous and synonymous sites (the latter is much smaller than the former, as substitutions occurring in the first two bases of each codon usually yield a residue substitution). For short genetic distances (e.g., on the order of 10% or less) as for our data, p-distances are nearly indistinguishable from other genetic distances. Also, considering the different protein regions that fit the autotransporter model proposed for some Pmps (see the introduction), this method was also applied separately to the 3' end of each pmp (corresponding to the 40-kDa C-terminal domain) and to the 5' half (defined by us to involve all conserved GGAI motifs) to evaluate possible evolutionary differences between these two major protein domains. Significant differences in mean dN and dS substitutions were determined by comparing 95% confidence intervals.
For all bioinformatics tests, a standard error (SE) computation was performed by bootstrap analysis with 1,000 replications.
A search for specific pmp regions that phylogenetically define tropism differences among serovars was also conducted: pmp genes with phylogenetic reconstructions showing serovars belonging to a specific disease group clustered in an evolutionary branch that was distant from the other serovars were subjected to a SimPlot software analysis (http://sray.med.som.jhmi.edu/SCRoftware/simplot/) (24). In SimPlot, serovars were grouped according to disease properties, followed by nucleotide pairwise comparison between the mean of the three groups using the Jukes-Cantor (one-parameter) correction for superimposed substitutions (gaps were included), yielding a graphic representation of nucleotide changes relative to their gene positions.
Analysis of protein features of the pmp family. Based on either presumptive (15, 18, 42) or experimentally demonstrated (28, 40, 41, 45) membrane location of Pmps for C. trachomatis (in most studies, only one serovar was used), several features that are characteristic of membrane proteins were studied for the entire Pmp family for all 18 serovars and genovariant Ja, namely, the presence of a cleavable peptidase signal, the presence of a C-terminal phenylalanine, and the distribution of tryptophan residues. Briefly, the first was evaluated using SignalIP 3.0 server software (1, 31) for prediction of cleavage sites and a signal peptide based on a combination of several artificial neural networks and hidden Markov models (32). The second and the third were determined by using EditSeq software (DNASTAR).
A search for cysteine residues was also performed, since they appeared to be nonrandomly distributed for some Pmps previously studied (15); in addition, disulfide bonds can hinder autotransporter function (19). The distribution of the tetrapeptide GGAI motif, which is characteristic of the Pmp family, throughout the protein was also examined. These two analyses were performed with EditSeq software (DNASTAR).
Nucleotide sequence accession numbers. The sequences of C. trachomatis pmp genes determined in this study were submitted to GenBank under accession numbers AY884090 through AY884108 (for pmpA), AY884109 through AY884127 (for pmpB), AY299408 through AY299426 (for pmpD), AY967735 through AY967738 (for pmpE), AY887644 through AY887662 (for pmpF), AY967739 through AY967757 (for pmpG), AY967759 through AY967761 (for pmpH), and AY299427 through AY299445 (for pmpI).
|
|
|---|
The mean genetic distance between all C. trachomatis serovars (Table 2) varied from the least polymorphic at 4.0 (SE, 1.2) (0.1%) nucleotides (nt) and 2.4 (SE, 0.9) (0.2%) amino acid differences for pmpA, to the most polymorphic at 215.5 (SE, 8.1) (7.0%) nucleotide and 71.4 (SE, 4.9) (6.9%) amino acid differences for pmpF (Fig. 1). The analysis of the genetic distance within each disease group showed that the LGV serovars were the most closely related for all pmp genes (except for pmpB and pmpE) with, in general,
0.1% nucleotide differences; four Pmps (PmpA, PmpD, PmpF, and PmpI) had 100% similarity among LGV serovars (Table 2). The urogenital group was the most heterogeneous, in particular for pmpF, which showed 3.3% nucleotide and amino acid differences among urogenital serovars.
|
View this table: [in a new window] |
TABLE 2. Average number of replacements observed and percent genetic distances of the nine pmp genes of C. trachomatis for all 18 serovars and the genovariant Jaa
|
![]() View larger version (19K): [in a new window] |
FIG. 1. Mean genetic distance within each of the nine C. trachomatis pmp genes based on the average p-distance between all possible pairs of sequences (different serovars) for the same pmp gene. Minimum and maximum values represent lower and upper limits of the 95% confidence interval of the estimate (vertical bar), while values plotted at the horizontal level represent the mean estimate. Absolute and percent values for the nucleotide and amino acid mean genetic distance for each pmp gene are also presented.
|
All genes except pmpA, pmpC, and pmpI revealed small deletion events for some serovars, varying from 1 for pmpD and pmpG to 21 codon deletions for pmpH (Table 3). Furthermore, pmpB had an insertion of 60 bp corresponding to a 20-amino-acid sequence at nt position 2427 for all serovars except for serovar Ia. This insertion included an imperfect direct target repeat (DTR) of 8 bp in the C terminus (2480CT[T/C]CAGCA2487) also present in the region preceding the insertion. As extensively discussed previously for pmpC, where we discovered putative DTRs for two IS-like elements (15), we think that this insertion might reflect a remnant fragment of an IS-like element, where the DTRs were created by duplication of the target sequence upon insertion of the IS element. A BLAST search revealed no homology with any gene or gene segments.
|
View this table: [in a new window] |
TABLE 3. Deletion events for pmp genes among C. trachomatis reference serovars
|
![]() View larger version (7K): [in a new window] |
FIG. 2. Phylogenetic reconstruction of the nucleotide (A) and amino acid (B) sequences showing the evolutionary history of pmpF by neighbor-joining tree topologies, based on distance estimates using a Kimura two-parameter model for substitution events (nucleotides) and the gamma distance model (amino acids). These reconstructions are based on pmpF sequences of the 18 serovars and the genovariant Ja of C. trachomatis. Branch lengths are proportional to distances between serovars. The values at the nodes are the bootstrap confidence levels representing the percentage of 1,000 bootstrap replicates for which the serovar to the right was separated from the others.
|
![]() View larger version (7K): [in a new window] |
FIG. 3. Phylogenetic reconstruction of the nucleotide (A) and amino acid (B) sequences showing the evolutionary history of pmpD by neighbor-joining tree topologies, based on distance estimates using a Kimura two-parameter model for substitution events (nucleotides) and the gamma distance model (amino acids). These reconstructions are based on pmpD sequences of the 18 serovars and the genovariant Ja of C. trachomatis. Branch lengths are proportional to distances between serovars. The values at the nodes are the bootstrap confidence levels representing the percentage of 1,000 bootstrap replicates for which the serovar to the right was separated from the others.
|
![]() View larger version (14K): [in a new window] |
FIG. 4. Mean genetic distance (p-distance) within the C. trachomatis genital group (11 serovars) with ( ) and without () serovars E and F, for pmpB to -D and pmpH to -I.
|
99% for nucleotide trees grouping serovar Da with ocular strains in pmpE through pmpI but with urogenital strains for pmpB and pmpC implied a P value of 0.01 x 0.01 = 0.0001 and that both groupings were not spurious. For third-codon position trees (not shown), which largely rule out any effects due to protein evolutionary pressure, the bootstrap support for these groupings was always at least 95%, implying a P value of 0.0025.
![]() View larger version (15K): [in a new window] |
FIG. 5. Similarity of pmp sequences to those from other ocular and urogenital strains. For the more variable pmp genes, the ocular strains (red) are well separated from the urogenital ones (blue asterisks). The exception is the urogenital serovar Da (green asterisks), which spans the two clusters, providing evidence for recombination. The lower-right Da points represent pmp genes pmpE, -F, and -H, each near the corresponding pmp genes from ocular strains. The upper-left Da points represent pmp genes pmpB and pmpC. Urogenital sequences from the more variable pmp genes (i.e., pmpE, -F, and -H) are not shown because they fall far off the scale to the left. Note that the percent identity is expressed as the fraction identity.
|
![]() View larger version (17K): [in a new window] |
FIG. 6. Top and middle graphs show mean synonymous and nonsynonymous mutation substitutions, respectively, for each of the nine pmp genes, based on the method of Nei-Gojobori (p-distance model). Minimum and maximum values represent lower and upper limits of the 95% confidence interval of the estimate (vertical bar), while values plotted at horizontal level (diamonds) represent the mean estimates. The graph on the bottom shows the nonsynonymous-to-synonymous mutation ratio for each pmp gene, based on the mean estimates shown in the prior two graphs.
|
1050 to 1900) (Fig. 7A). For pmpF (Fig. 7B), several main gene regions could be observed: one where the LGV group separated from the ocular and other urogenital serovars (from nt
300 to 1000), two regions where there was a clear separation of all groups (from nt
1250 to 1500 and
nt 2050 to 2250), and another two regions where the segregation was seen for the ocular group (from nt
1000 to 1200 and
nt 1600 to 2000). For pmpH (Fig. 7C), the ocular group was distinct from the other groups at nt 350 to 650 and nt 1600 to 1900; the LGV group had lower similarity to the other groups from nt 1100 to 1550 and from nt 2500 until the end of the gene. A clear separation of the three disease groups was found at nt 1900 to 2350. For PmpE, PmpF, and PmpH, these clusters of point mutations corresponded to clusters of nonconservative amino acid changes in the protein alignment (data not shown).
![]() View larger version (21K): [in a new window] |
FIG. 7. SimPlot graphs for pmpE (A), pmpF (B), and pmpH (C) of nucleotide similarity (using the Kimura two-parameter method) between the three disease groups, where each group represents the mean genetic distance of the respective serovars. Genital and LGV groups are compared to the ocular group. Ocular group, serovars A/SA-1, B/TW-5, Ba/Apache-2, and C/TW-3; urogenital group, urogenital group of serovars as described in Materials and Methods.
|
The number and distribution of the conserved cysteine residues varied from nine to 26 cysteines for PmpA and PmpD, respectively (see Fig. S7 in the supplemental data). The distribution of the GGAI motifs that were conserved for all serovars varied from 2 for PmpH to 11 for PmpD, respectively, and were almost exclusively located in the N-terminal half of the protein (see Fig. S7 in the supplemental data). With few exceptions, they seemed to be grouped in clusters. Furthermore, for PmpD, PmpE, PmpF, PmpG, PmpH, and PmpI, the distribution of these motifs seemed nonrandom, as
60% of them are adjacent or within 10 nt of the cysteine residues.
|
|
|---|
1.1 Mbp, about 25% of Escherichia coli's genome) is thought to have undergone progressive genome reduction from a larger precursor organism (33), while maintaining the ability to infect different tissue types and cause human diseases of major public health significance. It is thus presumed that this intracellular pathogen retains only the most important genes for survival and propagation. Consequently, the presence of a nine-member polymorphic membrane protein family that is unique to Chlamydiales and that represents 3.2% of the coding capacity of the genome (18) suggests that these proteins likely play a role in diversifying exposed antigens for escape from host immune surveillance (33). Our results support this theory. We found extensive polymorphisms among serovars for many of the pmp genes, the most notable of which was pmpF, where the mean genetic distance among all C. trachomatis serovars was 215.5 nt (7.0% change) (Table 2). Considering the variability in pmp and Pmp sequences, we expected to find an evolutionary divergence for serovars with different cellular tropism. In support of this and in contrast to other chlamydial membrane proteins that have been studied, six of the nine pmp genes presented phylogenetic reconstructions where serovars were clustered according to disease properties. Further, genetic distances were significant in distinguishing the ocular, urogenital, and LGV disease groups for pmpF and pmpH (Table 2).
Unexpectedly, the urogenital serovar Da showed an evolutionary trend similar to the ocular serovars for the pmp cluster from pmpE to pmpI. However, the high similarity with the ocular group within this large genomic region (
15,000 bp) and the remarkable genetic distance between the ocular and the urogenital groups for pmpE to pmpI (Table 2; Fig. 5) point to a recombination event instead of the accumulation of point mutations throughout evolution. Also, the phylogenetic analysis based on third-codon position trees that rule out any effects due to evolutionary pressure support these groupings (P = 0.0025). Thus, assuming that a partial arsenal of ocular Pmps can be effective, it appears that serovar Da has acquired a genetic composition that could hypothetically expand its host cell range. We know that urogenital serovars can cause a self-limited conjunctivitis (9) and we previously described ompA mosaics of B/D and Ba/D, the former of which was isolated from a patient with trachoma in Tunisia (8), while the latter was isolated from a patient with a sexually transmitted disease in the United States (26). While these occurrences do not directly support host cell expansion, they do support recombination among strains with different cell appetences. These collective data lend further support to the plasticity of the genome and the evidence for adaptive responses to environmental or host pressures.
We previously described two putative IS-like elements for pmpC (15). As for pmpC, the IS-like element identified for pmpB might reflect the prior presence of a complete insertion with a nearly precise excision phenomenon, leaving the imperfect downstream DTR and remnants of the IS at the site of the insertion (3, 6, 25). The differential presence among serovars of remnants of this putative IS-like element for pmpB (as for pmpC) may indicate either a differential insertion event among the serovars or a selective excision upon insertion. A complete deletion of the IS-like element may also have occurred for some serovars, a phenomenon that involves the two DTRs and DNA replication slippage (2, 37). It has been suggested that DTRs may lead to sequence deletions, as observed in this paper, and genome reductive evolution of intracellular bacteria (37), as has been documented for C. trachomatis (35). Considering the evidence for a high frequency of deletions occurring between two direct repeats in close proximity (2, 34), we would predict the complete deletion of these IS remnants during pmp evolution.
Phylogenetic analyses revealed that the pmp genes have evolved quite heterogeneously. pmpB showed the highest dN/dS ratio (1.59 [SE, 0.45]) (Table 2), indicating an evolutionary trend toward protein modification, which may occur via host immune pressure or horizontal gene transfer. In contrast, the most polymorphic pmp genes, pmpF and pmpH, encode the most evolutionary conserved proteins (with dN/dS values of 0.23 [SE, 0.02] and 0.14 [SE, 0.02], respectively) (Table 2). PmpF and PmpH may therefore play an important functional role where the occurrence of amino acid substitutions may be disadvantageous, while not being especially exposed to immune pressure. On the other hand, the phylogenetic clustering of serovars by disease properties and the remarkable genetic distance between serovars from different disease groups (Table 2) for pmpF and pmpH suggest an evolutionary adaptation of serovars to a specific tissue. Thus, we speculate that the biologic role of PmpF and PmpH may be associated with adhesion and entry into the host cell, where amino acid substitutions could be disadvantageous, by altering essential structural and functional characteristics. This is supported by the recent finding that antibodies specific to C. pneumoniae Pmp21 blocked chlamydial infectivity in the epithelial cells, which suggests a role for pmp genes in the adhesion process (50).
We also found that for five of the nine pmp genes (pmpB, -C, -D, -H, and -I), serovars E and F diverged separately from the other urogenital serovars and were closely related in a distant clade. This is interesting because these two serovars are the most prevalent among sexually transmitted diseases worldwide. The findings are distinct from ompA, where serovars E and F are evolutionarily closer to other serovars (27). The influence of these serovars on the heterogeneity of the urogenital group is illustrated in Fig. 4, where the exclusion of serovars E and F from this group for the five pmp genes dramatically decreased the mean genetic distance within the urogenital disease group. It is possible that these pmp genes confer a biological advantage for serovars E and F over the other serovars that may facilitate infectivity and transmission. Globally, as the majority of genes among the pmp family present phylogenetic reconstructions that largely segregate serovars according to disease properties, it indicates that Pmps may be involved in modulation of tissue tropism by either structural or functional constraints. While there is still a lack of data on the conformation and structure of Pmps, we have previously shown diverse immunoreactivity of patient sera to recombinant PmpC and heterogeneous transcriptome profiles for pmpC for serovars representing the three disease groups (16).
We also examined the distribution of polymorphisms in pmp genes that might segregate at least one disease group. For pmpE, pmpF, and pmpH, it was possible to observe a nonrandom distribution of point mutations with clustering in specific regions that clearly distinguished disease groups. pmpE contained one large gene region that distinguished the ocular serovars from all others; five other discrete regions were observed for pmpF and pmpH that differentiated at least one disease group (Fig. 7). For all these regions, clusters of point mutations also corresponded to smaller clusters of residue substitutions, which suggests that, if Pmps are hypothetically involved in tropism for chlamydial strains, these specific peptide segments on Pmps would certainly be responsible for those biological differences, either by the modification of the protein conformation structure or by alteration of its function. Unfortunately, the lack of similarity between Pmps and other proteins for which the three-dimensional structure has already been determined makes any prediction of hypothetical changes in protein folding due to these specific regions highly speculative.
Only recently has the translation and membrane localization of all Pmps been demonstrated by proteomics for C. trachomatis (serovar L2) (41), although previous studies had detected some Pmps (also for L2) (40, 45). The membrane location of Pmps was supported by our results, as demonstrated by a cleavable peptidase signal at the N-terminal domain (except for PmpA), a C-terminal phenylalanine, and tryptophan residues that were almost exclusively clustered at the C terminus for all serovars and Ja (see Fig. S6 in the supplemental material). Also, some Pmps, including PmpA, had high numbers of cysteine residues (PmpD; 26 residues), which are usually involved in the intra- and intermolecular interactions of membrane proteins through disulfide bonds. Although PmpA did not contain a signal peptide, the lack of one can be overcome by the presence of an anchor sequence (which is less well defined than the signal peptide) in a protein region that spans the C-terminal domain of the membrane and initiates translocation. Further, PmpA was the only Pmp reported by Skipp et al. (41) to localize to the reticulate body and not the elementary body, suggesting that it may be an important reticulate body membrane constituent but not under the influence of immune modulation.
The GGAI motifs were generally conserved among all serovars, were almost exclusively located in the N-terminal half of the protein, and in general were grouped in clusters (see Fig. S7 in the supplemental material). For PmpD, -E, -F, -G, -H, and -I, they were close or adjacent to cysteine residues, indicating a nonrandom distribution. In the autotransporter model proposed for chlamydial Pmps (19, 47), the ß-barrel structure of the C-terminal domain forms a pore in the membrane, mediating the translocation of the N terminus to the outer membrane and becoming exposed in bacterial surface. This theoretical model was experimentally demonstrated for C. pneumoniae Pmp21 (ortholog of C. trachomatis PmpD) (50). In view of this, the N-terminal location of the GGAI motifs, in close proximity to cysteine residues which are responsible for the highly disulfide cross-linked proteins of the outer membrane complex of C. trachomatis, suggests both a relevant biological role for the GGAI motifs and a hypothetical cooperative action with the cysteine residues for structural or functional constraints of the passenger domain. Considering this model, we have also separately studied the molecular evolution of both the N-terminal (passenger) domain and the
40 kDa C-terminal domain for each Pmp; curiously, PmpE (but none of the other Pmps) showed a remarkable difference of dN/dS between those two domains (1.21 [SE, 0.61] and 0.22 [SE, 0.05], respectively). Accordingly, it seems that the N-passenger domain of PmpE has been undergoing amino acid changes throughout evolution, which supports the above-described theoretical model in which this domain is exposed in the bacterial surface and subjected to immune pressure. In contrast, the C-terminal region shows clear evolutionary conservation, which supports conserved function without immune pressure for this domain.
The comparative genetic, genomic, and molecular analyses performed in this study suggest a possible structural or functional role that may vary among Pmps in promoting antigenic polymorphisms on the bacterial surface and/or promoting a set of different receptors that may be involved in host cell adhesion and differential tissue tropism. Additional genetic and genome sequencing, especially of clinical isolates, will be required to tease out the effects of each genetic event described above.
Supplemental material for this article may be found at http://jb.asm.org/. ![]()
|
|
|---|
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»