- Open Access
Chlamydia trachomatisdiversity viewed as a tissue-specific coevolutionary arms race
Genome Biology volume 9, Article number: R153 (2008)
The genomes of pathogens are thought to have evolved under selective pressure provided by the host in a coevolutionary arms race (the 'Red Queen's Hypothesis'). Traditionally, adaptation by pathogens is thought to rely not on whole chromosome dynamics but on gain/loss of specific genes, yielding differential abilities to infect distinct tissues. Thus, it is not known whether distinct host organs differently shape the genome of the same pathogen. We tested this hypothesis using Chlamydia trachomatis as model species, looking at 15 serovars that infect different organs: eyes, genitalia and lymph nodes.
We analyzed over 51,000 base pairs from all serovars using various phylogenetic approaches and a non-phylogenetic indel-based algorithm to study the evolution of individual and concatenated loci. This survey comprised about 33% of all single nucleotide polymorphisms in C. trachomatis chromosomes. We present a model in which genome evolution indeed correlates with the cell type (epithelial versus lymph cells) and organ (eyes versus genitalia) that a serovar infects, illustrating an adaptation to physiologically distinct niches, and discarding genetic drift as the dominant evolutionary driving force. We show that radiation of serovars occurred primarily by accumulation of single nucleotide polymorphisms in intergenomic regions, housekeeping genes, and genes encoding hypothetical and cell envelope proteins. Furthermore, serovar evolution also correlates with ecological success, as the two most successful serovars showed a parallel evolution.
We identified a single nucleotide polymorphism-based tissue-specific arms race for strains in the same species, reflecting global chromosomal dynamics. Studying such tissue-specific arms race scenarios is crucial for understanding pathogen-host interactions during the course of infectious diseases, in order to dissect pathogen biology and develop preventive and therapeutic strategies.
When two species interact with each other, such as a pathogen and human, a never-ending reciprocal and dynamic adaptation process takes place. Whereas the 'goal' of the human being is to try to avoid, solve or minimize the infection, the 'goal' of the pathogen is to deal with this constant host environmental and immune pressure, through genomic evolutionary changes, in order to win this arms race [1–4]. Typically, genome evolution within same-species strains of a pathogen has been studied mainly in the light of horizontal gene transfer (HGT) at specific chromosome loci [5, 6], as for Escherichia coli [7, 8], Staphylococcus aureus , Streptococcus pyogenes , Salmonella enterica , Shigella flexneri , and Pseudomonas syringae . An extreme example is provided by the well-studied E. coli, where strains K-12 and O157 differ by more than 1 million base pairs , and same-serovar strains were found to present profound differences in gene content [13, 14]. Globally, these targeted HGT events reflect different pathoadaptation processes for microrganisms with reversible genome size-plasticity; depending on the transitory 'cassette-genes' carried at any specific time, the pathogenecity or ability of these microrganisms to infect different tissues may vary . Thus, generally, these processes rely on gain/loss of virulence/colonization factors rather than reflect whole chromosomal dynamics, the evaluation of which remains complex. Indeed, assessment of tissue-specific adaptive evolution at the whole genome level demands that same-species strains of a pathogen specifically and non-transitorily infect different tissues. Therefore, on behalf of the arms race theory assumed by the evolutionary Red Queen's Hypothesis [15, 16], one question arises: do distinct host organs differently shape the genome of the same pathogen? No microrganism is more suitable than Chlamydia trachomatis, the most prevalent sexually transmitted bacterial pathogen worldwide, to test this hypothesis, as the species comprises several serovars with a wide range of specific human tissue tropism. This pathogen is mainly classified into 15 serovars based on the differential immunoreactivity of the major outer membrane protein (MOMP), constituting three disease groups : serovars A-C and Ba are commonly associated with ocular trachoma; serovars D-K infect the epithelial cells of genitalia and are normally found in non-invasive sexually transmitted infections (where serovar E represents about one-third of all infections, and together with serovar F constitute up to 50% of them); serovars L1-L3 are also sexually transmitted but are invasive and disseminate into the local lymph nodes causing lymphogranuloma venereum (LGV). However, in the context of this classification system, the evaluation of adaptive evolution becomes enigmatic because there is no correlation between it and C. trachomatis tropism nor with the ecological success of the different serovars, as strains with different organ specificities are placed within the same classification group.
As occurred for Mycobacterium leprae , Rickettsia prowazekii , and the aphid endosymbiont Buchnera aphidicola , the first stages of Chlamydia evolution consisted of a massive genome reduction upon becoming an obligate intracellular parasite [21, 22]. However, comparative genomics over the few currently fully sequenced C. trachomatis genomes [20, 23–25] revealed that gene decay is not involved in the more recent evolutionary stages. Indeed, contrary to most pathogens, the core- and the pan-genome  of this microrganism are near identical, indicating that the factors involved in the differential organ specificity among serovars are not acquired by gene transfer .
To evaluate if distinct arms races occur between different infected human organs and this pathogen's serovars, we performed high-scale concatenation-based phylogenomics, using about one-third of all chromosome single nucleotide polymorphisms (SNPs). So far, in contrast to the ocular group, only one strain from the epithelial-genital and LGV groups has been fully sequenced [20, 23–25], making our multiple-loci scrutiny of all 15 serovars the ideal tool to track the evolutionary diversity of a microrganism characterized for its distinct infection niches. Here, we show a matchless model of SNP-based adaptive evolution of same-species strains to each infected cell-type and organ that relies on whole chromosome evolutionary dynamics, unlike previous reports for other pathogens focused on specific gene gain/loss.
Evaluation of the degree of polymorphism for the selected loci
Considering that the strain radiation yielding the present-day chlamydial serovars likely occurred over millions of years , the use of reference strains is an accurate strategy as they were isolated only a few decades ago. Thus, in this evolutionary survey, we used the traditional reference strains that represent all 15 C. trachomatis serovars. We selected 51 polymorphic loci (approximately 51,000 bp) dispersed throughout the chromosome (Figure 1; Additional data file 1) that represent the following loci categories: 16 intergenomic regions (IGRs); 16 genes encoding cell envelope proteins (CEPs); 13 housekeeping genes (HKs); and 6 genes encoding hypothetical or unclassified proteins (HPs) (Additional data file 2). In order to evaluate the degree of polymorphism of these loci in comparison with the whole chromosome, we used the data generated from two of the five fully sequenced genomes, A/Har13 (ocular)  and D/UW3 (epithelial-genital) . We observed in the studied 51 loci a global mutation rate 14.3-fold higher than in the remaining chromosome regions (Fisher's exact test, P < 0.001). Moreover, we found 1,099 SNPs in these 51 loci between A/Har13 and D/UW3, which is greater than 200-fold more than what has been studied to date through concatenation , and comprises about 33% of the whole chromosome SNPs, indicating that our results could be scaled up to the full-chromosome level.
Additionally, a global overview of GC content revealed a mean value for all loci categories (data not shown) that is similar to the total mean GC content of approximately 41% observed for the fully sequenced genomes [21, 23–25] with a standard deviation of 2.9%, which is not indicative of any putative HGT event.
Correlation of individual loci with tissue-specific strain radiation
We used phylogenomics to correlate each individual locus with tissue-specific strain radiation. Only four (25.0%) CEPs (incD, incE, pmpF and pmpH) and one (6.3%) IGR (incD/incE) comprehensively grouped the strains according to their cell-type/organ appetence (that is, revealed a larger evolutionary distance between strains with different niche appetencies than between strains infecting the same niche; Figure 2a). This clustering seems to be associated with loci revealing a higher p-distance-based polymorphism (Mann-Whitney P = 0.025). A full segregation by cell-type/organ appetence was not seen for most of the remaining CEPs due to the heterogeneity among the genital strains, where serovars E and F frequently form a separate cluster for 62.5% of CEPs (Figure 2a). Globally, 77.6% of loci belonging to different functional categories grouped strains that invade the lymph nodes as an individual cluster (LGV cluster), and the clustering of strains infecting the ocular tissue (ocular cluster) was also frequent. As above, we identified a significant association between a higher absolute number of SNPs and both the occurrence of a LGV cluster and an ocular cluster for each locus (Mann-Whitney P = 0.037 and P = 0.045, respectively). Interestingly, from the loci that better illustrate adaptation to lymph nodes, 80% of HPs and 53% of CEPs, compared with only 29% of HKs, show >50% non-synonymous SNPs (Figure 2b). Considering the DNA replication process, all SNPs on one strand that may imply strain segregation will also have the same impact on the other DNA strand. However, from the 51 loci that we used, only 4 pairs of loci overlap and the overlapping region never exceeds 10 bp (data not shown), which makes this effect negligible. Overall, these results suggest that the distinct genetic variability of strains infecting a specific cell-type/organ likely reflects an evolutionary adaptation process.
By performing intra-locus analysis, we observed that three HPs (CT049, CT144 and CT622) and two IGRs (rs2/ompA and ompA/pbpB) revealed distinct domains in which SNPs are concentrated, instead of being randomly distributed, and are associated with strains that infect a specific cell-type/organ (Figure 3). For these HPs, the SNP domains correspond to clusters of amino acid changes in the protein sequence (data not shown), mirroring the previous findings for some polymorphic membrane protein genes . Unfortunately, there is no assigned role for these open reading frames, which rules out any speculation about the functional implications of these specific clustered amino acid alterations. Nevertheless, this tissue-specific amino acid clustering points to a targeted fixation of mutations that may reflect the host-pathogen specific interaction within each organ.
Genomic analysis of the concatenated loci
We evaluated the nucleotide sequence variation in each concatenated loci category (Table 1). We highlight the multi-loci concatenation approach as a powerful tool to generate robust phylogenomic inferences, even when individual loci have evolved with different substitution patterns [29–31]. Overall, the HPs exhibit the highest number of variable sites (10.3%), whereas the HKs are the least variable (3.3%), which is supported by the mean p-distance values. Curiously, the IGRs show polymorphism similar to the CEPs. Globally, concatenation of all 51 loci yielded a 'super' sequence of up to 51,074 bp for each of the 15 reference strains, showing a mean of 1,032.1 (standard error (SE) 17.2) nucleotide differences.
Evolutionary history of C. trachomatis
Due to the speed and efficiency of the neighbor joining (NJ) method in inferring large phylogenies [32, 33], we used this approach on concatenated data. The NJ phylogenies inferred from the four concatenated loci categories (Additional data file 3) are consistent with most of the respective individual loci trees. Although only the CEP category clearly segregates strains by the disease they cause, the other categories show a notable segregation of at least one disease group, suggesting that heterogeneous loci categories are involved in the arms race process. The global phylogenetic tree presented in Figure 4 (where each taxon is represented by about 50,000 bp) reveals the putative final picture of C. trachomatis's evolution, showing strain grouping according to the cell-type (epithelial and lymph cells) and organ (eyes and genitalia) that they infect. These distinct segregations are supported by maximum bootstrap values (99-100%) in the nodes that separate disease groups, reinforcing that the targeted and distinct fixation of nucleotide changes on strains infecting a specific cell-type/organ are likely adaptive and barely the consequence of genetic drift. In fact, the genetic distance matrix (Table 2) shows that all strains that preferentially infect the eyes revealed only 0.27% (SE 0.02%) differences among them, but shows a mean genetic distance 7.4- and 11.2-fold higher (corresponding to 983 (SE 20) and 1,484 (SE 42) nucleotides) to strains infecting the epithelial-genital and lymph node tissues, respectively. Also, the LGV strains differ by only 69 (SE 8) nucleotides, whereas their distance to the epithelial-genital strains is 1,226 (SE 34) nucleotides. A separate main branch involving all epithelial-genital strains was not comprehensively seen for any individual loci (except for the CEPs pmpF and pmpH; data not shown) due to the separation of E and F strains. Indeed, the latter has a mean genetic distance of 673 (SE 16) nucleotides to the other epithelial-genital strains (Table 2). Similar NJ tree topologies were obtained for the three models used to estimate evolutionary distances (Kimura 2-parameter (K2P), Jukes-Cantor or Tamura-Nei) as well as for the maximum parsimony method (data not shown), with only slight variations in the bootstrap values, which supports the robustness of these distinct arms race scenarios.
We also highlight the loci that most contribute to the final tree topology (Figure 4), as they may be relevant for the evolutionary adaptation to each specific niche. Among these loci, we have found either highly conserved or polymorphic loci for strains infecting the same cell-type/organ. The former may represent a step forward in the evolutionary process by revealing the final stages  of this tissue-specific adaptive evolution, while the latter may also be involved in pathogenic differences between strains infecting the same tissue . The most extreme case is given by the CEP pmpF, where all the strains that infect the lymph nodes are 100% similar but show a mean distance of 312 and 421 SNPs to strains infecting the epithelial-genital and ocular tissues, respectively. In contrast, the epithelial-genital strains reveal up to 129 SNPs among them (data not shown). Although less markedly, CT049 is polymorphic among the LGV strains but near 100% identical among the ocular strains.
Additionally, we identified loci that do not seem to have influenced adaptation to each niche, since they generate an incongruent strain-radiation (Table 3), and whose polymorphism may thus be a consequence of genetic drift. However, previous results have demonstrated the involvement of some of these loci (CT622, tsf, rs2 and pbpB) in the pathogenesis of trachoma . As expected because of the serovar multiplicity, the epithelial-genital group revealed a higher number of polymorphic loci, and, overall, these loci belong to different categories. In contrast, strains infecting the lymph nodes constitute the most homogeneous group.
Impact of small insertions/deletions (indels) on tissue-specific strain radiation
In order to have a more complete picture of the evolution of the serovars, we studied the chromosomal occurrence of small insertion/deletion (indel) events, which are non-phylogenetic parameters. We observed 84 small indel events (from 1-43 bp) inside the global concatenated loci for all strains, which mainly occurred within the IGR and CEP categories (Additional data file 4). None of these events was found to disrupt the coding sequence of the respective loci, indicating the absence of gene decay in the studied regions.
For the global concatenated data, we estimated the evolutionary distances using the indel-based parameter γ , which computes the number of gap nucleotides per nucleotide site between those sequences, while SNPs are not considered. The γ-distances (Figure 5a) are highly concordant with phylogenomic analyses, showing high heterogeneity within the epithelial-genital strains, and remarkable homogeneity among the LGV strains. Also, they revealed a segregation of strains by their cell-type/organ appetencies, which supports the tissue-specific arms race scenario.
Evolutionary inferences on the ecological success
Analysis of the global phylogenetic tree (Figure 4) also shows that the two most prevalent genital serovars worldwide, E and F, are closely related and separated from the other epithelial-genital strains. This segregation is observed for the majority of loci, with the exception of the HPs (Figure 2a). From all these loci, 70% of CEPs show an amino acid replacement for >50% of SNPs, compared to only 20% of HKs (Figure 2b). Curiously, the most remarkable segregation of E and F was seen for two IGRs (rs2/ompA and yfh0_1/parB) and three HKs (karG, tsf and rs2) (Figure 4). Furthermore, for the still unclassified protein gene CT622 and for the IGR rs2/ompA, we observed a non-random distribution of SNPs that are present in serovars E and F but not in the other epithelial-genital strains (Figure 3c,d). Finally, the mean γ-distance from any epithelial-genital strain to serovar E or F was from 3.4-fold (between G and E/F) to 4.7-fold (between I and E/F) higher than the distance between E and F (Figure 5b), which supports this close relationship between the two most ecologically successful serovars.
We have hypothesized that distinct arms races may occur inside the same host when the same pathogen is able to infect different organs. In contrast to free living bacteria, where HGT is strongly associated with a pathogen's adaptive evolution [3, 5–11], Chlamydia has been characterized by genetic isolation and, while cumulative studies suggest that HGT has almost certainly occurred in Chlamydiaceae [35–37], there is no report to date of transferable mobile elements in C. trachomatis. Here, we demonstrate that C. trachomatis strains that preferentially infect the eyes, the epithelial-genital cells or the lymph nodes present a distinct evolutionary pattern likely illustrating a SNP-based tissue-specific arms race.
In order to develop a more compelling argument for a causal link between genome profile and cell/organ appetence, the use of genetic modification and especially the use of animal models are appealing approaches. However, C. trachomatis is genetically non-tractable and, except for the cynomolgus monkey (accurate for studying the trachoma pathology) , no suitable animal model exists for the three types of C. trachomatis disease. Also, there is no in vitro model, such as cell culture, that mirrors the chlamydial infection in vivo, and it has been previously demonstrated that intensive serial passaging of chlamydial strains yielded no mutations on the most variable chlamydial gene (ompA) . Furthermore, it would be inconceivable that these approaches could represent millions of years of chlamydial evolution.
It is believed that the LGV biovar was the first to diverge from a common C. trachomatis ancestor when new primate hosts evolved after the dinosaur extinction, whereas separation of genital and ocular serovars might have occurred with the appearance of early humanoid primate hosts . The skill to colonize different organs and cell-types likely developed through indel events and SNP accumulation on virulence/colonization factors. So far, chlamydial putative virulence factors, such as the type III effector tarp , the cytotoxin gene , and especially the tryptophan operon [40, 41], are the best candidates for providing that skill. In particular, while the first of these factors differentiates the LGV strains from the other groups, the other two differentiate the strains colonizing the genitalia from the strains colonizing other niches. For example, it was clearly demonstrated that only strains possessing a functional trpBA operon are able to colonize the genital tract . With respect to type III effectors, although their role in C. trachomatis tropism is not clear, it was shown that evolutionary genetic diversification of the type III effector HopZ family, via horizontal transfer, had clear implications for Pseudomonas syringae host specificity . However, none of the chlamydial putative virulence factors fully explain the existence of the three major tropism groups made up from the different serovars. Also, the putative emergence of tissue-specific adhesins cannot be discarded.
With regard to our results (Figure 4), strain radiation within each disease group likely occurred because of accumulation of mutations throughout the chromosome caused by environmental and immune pressure in each niche, giving rise to the contemporary serovars. Within the genitalia, the higher serovar multiplicity and radiation of epithelial-genital strains compared to the LGV strains would be unexpected in the light of the earlier evolutionary divergence of the latter . However, besides the different host immune responses in those niches, the epithelial-genitalia environment presents pH and hormonal fluctuations that are variable among individuals, and also an abundant nutrient-competing flora, which could have strongly influenced the evolutionary pathway of the infecting strains. In support of this, nutrient-competing flora were shown to be a major factor in the successful pathoadaptation of Salmonella enterica serovar Typhimurium to the intestinal tract, as the inflammatory process induced by this pathogen was shown to make a negative impact on mainly the other colonizing microrganisms and, thus, a positive impact on its arms race with the host .
Globally, we have observed that the loci that most contribute to strain segregation by cell-type/organ are spread throughout the chromosome (Figure 1) and belong to different functional categories, suggesting that this dynamic evolutionary adaptation is a general trait of the entire genome. Whereas the contribution of CEPs is likely associated with putative structural, antigenic or host-adhesion roles, no assumption can be made for the HPs. However, we found that HPs were the most variable among the serovars, with an overall polymorphism 2.2-fold higher than the CEPs (Table 1), which suggests a higher involvement in chromosomal dynamics. With respect to IGRs, we speculate that their contribution to strain segregation may be associated with recombination events that may promote genetic variability, as we recently described . Nevertheless, the high variability of IGRs was surprising, as they commonly involve regulatory regions that are expected to be conserved; thus, the existence of random genetic drift may also be considered for IGRs. Finally, although the HKs are involved in strain segregation, the vast majority of them showed <50% non-synonymous mutations (Figure 2b), which is consistent with their role in essential biological functions.
It is known that in populations without HGT and with bottlenecks, as is the case for C. trachomatis, random genetic drift can play a major role in evolution, being responsible for the fixation of unfavorable mutations . However, our results suggest that chlamydial strain segregation according to tropism properties occurred mainly through an adaptive evolutionary process and not through dominant genetic drift. Several arguments point in this direction: the statistical association found between most polymorphic loci (number of SNPs/loci and p-distance/loci) and the strain clustering according to their tissue specificity; Chlamydiae presents a relatively high ratio of non-synonymous to synonymous changes when compared, for example, to E. coli and Buchnera , further supported by our findings where the majority of HPs and CEPs involved in the segregation of the LGV strains showed >50% non-synonymous SNPs (Figure 2b); for at least eight loci (CT049, CT144, CT622, pmpE, pmpF, pmpH, rs2/ompA IGR and ompA/pbpB IGR), we observed a non-random fixation of SNPs exclusive of same niche-infecting strains (Figure 3), corresponding to specific clusters of amino acid changes in coding sequences; the extremely robust global phylogenetic tree with maximum bootstrap support (99-100%) in the branch nodes where strains are separated by their cell-type/organ specificity (Figure 4); 20 out of the 22 loci that contribute to the segregation of strains that preferentially infect the eyes are also involved in the segregation of strains that colonize the lymph nodes (Figure 4) by presenting a dissimilar and specific SNP pattern; and finally, the well-known differences in environmental and immune pressure as well as competing flora and physiological specificities between ocular, epithelial-genital and lymph node tissues.
Within all the loci that are more likely to be involved in the adaptive evolution to each specific niche, we have found either highly conserved or polymorphic loci among strains infecting the same cell-type/organ (Figure 4), where the most remarkable examples are pmpF and CT049 (see Results). We hypothesize that pmpF and CT049 may be good representatives of a final stage of the adaptive evolution to the lymph nodes and the eyes, respectively, considering their extreme conservation among the corresponding strains. On the other hand, these genes may be responsible for pathogenic differences among epithelial-genital and LGV strains, respectively, based on their strong polymorphism among the corresponding strains. While PmpF has been implicated as a potential target for the host immune response, as it contains several putative major histocompatibility epitopes , biological information for CT049 is lacking.
Additionally, we found several loci that are polymorphic among strains infecting the same cell-type/organ that seem not to have been involved in the adaptation to each niche, but which may have been involved in the pathogenesis of trachoma, genital infections or LGV disease (Table 3). Indeed, 4 of these loci (CT622, tsf, rs2 and pbpB) belong to a pool of 22 genes that are responsible for profound differences in virulence among two C. trachomatis ocular strains in nonhuman primates .
Interestingly, we also observed a clear evolutionary co-segregation of the two most ecologically successful serovars (E and F). This is intriguing as there is a 15% difference between them in the gene coding for the major antigen (the major outer membrane protein (MOMP)), which constitutes about 60% of the membrane dry-weight  and is a putative cytoadhesin . Although it is not known why serovars E and F are the most prevalent worldwide, their ecological success seems not to be associated with intracellular multiplication rate , indicating that it is likely defined at the host cell adhesion and entry steps. However, the existence of E/F specific virulence factors or adhesins cannot be addressed in this study. Even so, tarp is the unique virulence factor that distinguishes serovar E from the other epithelial-genital serovars (including F), as it presents fewer repeat motifs in the 5' region , but its phenotypic consequences are not known. Moreover, a more successful host immune evasion could also be speculated for serovars E and F considering the well-known different antigenic profile among epithelial-genital serovars .
Regarding the loci that most markedly contribute to the segregation of serovars E and F, we highlight the IGRs tsf, rs2 and rs2/ompA (Figure 4). The first two of these may be involved in hypothetical differences in strain growth , while the last involves the regulatory region of rs2. This IGR includes specific domains where most SNPs are exclusive of strains E and F (Figure 3d), suggesting a potential impact on the rs2 regulation and, thus, on strain growth. Also, the IGR rs2/ompA is a recombination hotspot for the generation of mosaic structures within chlamydial strains , suggesting that recombination may contribute to the ecological success of the two serovars. However, as most SNPs of the CEPs involved in the E/F segregation confer amino acid replacements (Figure 2b), we suggest that the positive selection for the membrane proteins may also be a driving force for the E/F evolutionary divergence, likely through antigenic variability.
It is not surprising that bacterial populations that evolved in different ecological niches have different profiles of genetic variability. However, contrary to all previous reports for other pathogens focused on HGT events and gene decay, we present evidence of SNP-based, tissue-specific evolutionary adaptation relying on whole chromosome dynamics, as a consequence of the occurrence of dissimilar arms races between the pathogen and diverse host organs. Answering the proverbial question of 'which came first' (tropism or SNPs), the scenario presented here suggests that while some SNPs, on very few and specific loci, are likely responsible for tropism differences, the vast majority of SNPs throughout the chromosome are a consequence of different tissue tropisms and are expected to be involved in maintaining organ appetence, as per the Red Queen's Hypothesis. Mirroring bacterial virulence , we present evidence that a 'one size fits all' approach cannot be applied to adaptive evolution. This phenomenon is illustrated by a pathogen believed to infect 140 million people, where the incidence rate can be as high as 30% among adolescent females . We believe that grasping a pathogen's genetic trends with regard to its interaction with the host will be an essential tool in deciphering the molecular genetic aspects of infectious diseases.
Materials and methods
Culture of C. trachomatisreference strains
We used the most common reference strains representing the 15 C. trachomatis serovars: A/Har13, B/TW5, Ba/Apache2, C/TW3, D/UW3, E/Bour, F/IC-Cal3, G/UW57, H/UW4, I/UW12, J/UW36, K/UW31, L1/440, L2/434 and L3/404. McCoy cell culture of all strains plated in T-25 cm2 flasks was performed as previously described . At 48-72 h post-infection, elementary bodies were harvested, and DNA was extracted using QIAamp® DNA Mini Kit (Qiagen, Valencia, CA, USA) according to the manufacturer's instructions. Serovar confirmation of each reference strain was performed using ompA genotyping with BLAST comparison of the available GenBank sequences.
Selection of loci
A GenBank search was performed to look for genomic regions that had been sequenced for at least one C. trachomatis reference strain from each of the three disease groups. Up to 93 loci were found, comprising about 84,000 bp of the chromosome, and involving IGRs, HKs, HPs and CEPs. Only non-constant loci were selected (51 of the 93; Figure 1; Additional data file 2) for sequencing the other reference strains if their sequences were not available yet. Automated sequencing was performed as previously described . The DNA sequence data have been deposited in a public database ([GenBank: EU239694-EU239702], [GenBank:EU239705-EU239712], and [GenBank:EU247618-EU247753]). Primer sequences are given in Additional data file 5. For all strains, five types of concatenated sequences were created in a head-to-tail fashion: one for each loci category (IGRs, HKs, HPs and CEPs) and a global concatenated sequence involving all loci (approximately 50,000 bp for each taxon).
We used data from the fully sequenced genomes A/Har13 and D/UW3 for this evaluation. Thus, considering the 3,354 SNPs identified between these two genomes , we evaluated whether 1,099 SNPs restricted to the 51,074 bp analyzed in this study are overrepresented relative to the 2,255 SNPs found in the rest of the chromosome. We framed this as a contingency table (Table 4) with a restricted sequence of 1,519,042 bp for each strain (corresponding to the length of the D/UW3 chromosome), and we estimated P-values using the Fisher's exact test as well as the odds ratios with a 95% confidence interval.
For all individual loci and concatenated sequences, alignments of all strains were generated using LaserGene (DNASTAR, Madison, WI, USA) and MEGA 3.1 . MEGA 3.1 was also used to create matrices of pairwise comparisons and to estimate the number of variable sites, the number of parsimony informative sites and overall mean genetic distances. The pairwise-deletion option was chosen to remove all sites containing missing data or alignment gaps from all distance estimations, only when the need arose and not prior to the analysis.
In order to search for distinct regions that may be associated with strains belonging to a specific disease group, SimPlot 3.5.1  was used on all 51 loci. For each similarity plot, serovars were grouped according to the cell-type/organ that they infect, and nucleotide pairwise distances were calculated using the K2P method (gaps excluded; ts/tv of 2.0) in a sliding window size of 160 bp moved across the alignment in a step size of 10 bp. Additionally, for all loci where serovars E and F clustered apart from the other epithelial-genital serovars was observed, a SimPlot analysis was also performed to evaluate if the E/F nucleotide differences compared to the other genital serovars were clustered in specific domains of each locus.
Prior to the phylogenetic reconstructions, and in order to select the appropriate evolutionary models, we evaluated the homogeneity of substitution patterns between sequences by calculating the Monte Carlo test-based Disparity Index per site . This gives the probability of rejecting the null hypothesis that sequences have evolved with the same pattern of substitution. The NJ method  was used with K2P , Jukes-Cantor  and Tamura-Nei  models to generate phylogenies. For the concatenated sequences, in order to examine the accuracy of the major conclusions reached from the NJ analysis, trees were also constructed under the maximum parsimony criterion , using the max-min branch-and-bound algorithm.
Considering that recombination disturbs a phylogenetic signal since the two parts of the recombined region may have different evolutionary histories [61, 62], one locus (ompA) was excluded from the phylogenetic concatenated analysis, as its highly recombinant nature has already been demonstrated [48, 63]. The use of outgroup sequences was discarded in the present study because no rooted trees were needed to achieve the objectives defined above. Also, the most suitable strain for use as an outgroup, C. muridarum (MoPn strain), has several loci that vary greatly in size and diverge from those in the C. trachomatis strains, which would entail the removal of a huge portion of the sequences being analyzed.
We used a non-phylogenetic method for estimating the evolutionary distance between each pair of homologous DNA sequences, which is given by the parameter γ :
γ = -2logeP
P = n xy /√n x n y
where n xy is the number of nucleotides shared by the two sequences, and n x and n y are the number of nucleotides of each sequence. For comparative purposes, we used the same set of loci as for the phylogenetic concatenated analyses, that is, we excluded the recombinant ompA gene. The γ variability was estimated by Monte Carlo using the alignments of each individual locus through the statistical platform R 2.5.1 . Each time, 20 loci were randomly selected with replacement and γ-distances were calculated by repeating this procedure 50 times.
The statistical association between genetic and phylogenetic variables was performed using the ANOVA test by comparing groups' population means. We considered as genetic variables the overall mean values of percent GC content, p-distance and absolute SNPs obtained for each of the selected loci. The phylogenetic variables were: clustering of strains according to tropism properties; co-segregation of E/F strains; segregation of a LGV cluster or an ocular cluster; and the 'weight' of each locus in the final concatenated tree. The homogeneity of variances was tested using the Levene's test. Whenever the hypothesis of homogeneity of variances was rejected, the non-parametric Mann-Whitney test was used to compare distributions among groups. A P-value of 0.05 or less was considered significant.
Additional data files
The following additional data are available with the online version of this paper. Additional data file 1 is a figure showing the overall mean genetic distances among all 15 C. trachomatis serovars for the 51 loci. Additional data file 2 is a table listing the cellular roles of the 51 loci. Additional data file 3 is a figure showing C. trachomatis's evolutionary history by loci category. Additional data file 4 is a table listing the indel events found among the 15 C. trachomatis serovars for all loci. Additional data file 5 is a table listing the primers used for PCR and sequencing of selected loci.
cell envelope protein gene
horizontal gene transfer
hypothetical or unclassified protein gene
single nucleotide polymorphism.
Bush RM: Predicting adaptive evolution. Nat Rev Genet. 2001, 2: 387-392. 10.1038/35072023.
Woolhouse ME, Webster JP, Domingo E, Charlesworth B, Levin BR: Biological and biomedical implications of the co-evolution of pathogens and their hosts. Nat Genet. 2002, 32: 569-577. 10.1038/ng1202-569.
Ma W, Dong FFT, Stavrinides J, Guttman DS: Type III effector diversification via both pathoadaptation and horizontal transfer in response to a coevolutionary arms race. PLoS Genet. 2006, 2: e209-10.1371/journal.pgen.0020209.
Grant AJ, Restif O, McKinley TJ, Sheppard M, Maskell DJ, Mastroeni P: Modelling within-host spatiotemporal dynamics of invasive bacterial disease. PLoS Biol. 2008, 6: e74-10.1371/journal.pbio.0060074.
Ochman H, Moran NA: Genes lost and genes found: evolution of bacterial pathogenesis and symbiosis. Science. 2001, 292: 1096-1099. 10.1126/science.1058543.
Pallen MJ, Wren BW: Bacterial pathogenomics. Nature. 2007, 449: 835-842. 10.1038/nature06248.
Groisman EA, Ochman H: Pathogenicity islands: bacterial evolution in quantum leaps. Cell. 1996, 87: 791-794. 10.1016/S0092-8674(00)81985-6.
Ohnishi M, Kurokawa K, Hayashi T: Diversification of Escherichia coli genomes: are bacteriophages the major contributors?. Trends Microbiol. 2001, 9: 481-485. 10.1016/S0966-842X(01)02173-4.
Brussow H, Canchaya C, Hardt WD: Phages and the evolution of bacterial pathogens: from genomic rearrangements to lysogenic conversion. Microbiol Mol Biol Rev. 2004, 68: 560-602. 10.1128/MMBR.68.3.560-602.2004.
Groisman EA, Ochman H: How Salmonella became a pathogen. Trends Microbiol. 1997, 5: 343-349. 10.1016/S0966-842X(97)01099-8.
West NP, Sansonetti P, Mounier J, Exley RM, Parsot C, Guadagnini S, Prévost MC, Prochnicka-Chalufour A, Delepierre M, Tanguy M, Tang CM: Optimization of virulence functions through glucosylation of Shigella LPS. Science. 2005, 307: 1313-1317. 10.1126/science.1108472.
Hayashi T, Makino K, Ohnishi M, Kurokawa K, Ishii K, Yokoyama K, Han CG, Ohtsubo E, Nakayama K, Murata T, Tanaka M, Tobe T, Iida T, Takami H, Honda T, Sasakawa C, Ogasawara N, Yasunaga T, Kuhara S, Shiba T, Hattori M, Shinagawa H: Complete genome sequence of enterohemorrhagic Escherichia coli O157:H7 and genomic comparison with a laboratory strain K-12. DNA Res. 2001, 8: 11-22. 10.1093/dnares/8.1.11.
Ohnishi M, Terajima J, Kurokawa K, Nakayama K, Murata T, Tamura K, Ogura Y, Watanabe H, Hayashi T: Genomic diversity of enterohemorrhagic Escherichia coli O157 revealed by whole genome PCR scanning. Proc Natl Acad Sci USA. 2002, 99: 17043-17048. 10.1073/pnas.262441699.
Ogura Y, Ooka T, Asadulghani , Terajima J, Nougayrède JP, Kurokawa K, Tashiro K, Tobe T, Nakayama K, Kuhara S, Oswald E, Watanabe H, Hayashi T: Extensive genomic diversity and selective conservation of virulence-determinants in enterohemorrhagic Escherichia coli strains of O157 and non-O157 serotypes. Genome Biol. 2007, 8: R138-10.1186/gb-2007-8-7-r138.
Van Valen L: A new evolutionary law. Evol Theory. 1973, 1: 1-30.
Bell G: The Master Piece of Nature: The Evolution and Genetics of Sexuallity. 1982, Berkeley: University of California Press
Fields PI, Barnes RC: The Genus Chlamydia. The Prokaryotes. Edited by: Balows A, Truper HG, Dworkin M, Harder W, Schleifer KH. 1992, New York: Springer-Verlag, 3691-3709.
Cole ST, Eiglmeier K, Parkhill J, James KD, Thomson NR, Wheeler PR, Honoré N, Garnier T, Churcher C, Harris D, Mungall K, Basham D, Brown D, Chillingworth T, Connor R, Davies RM, Devlin K, Duthoy S, Feltwell T, Fraser A, Hamlin N, Holroyd S, Hornsby T, Jagels K, Lacroix C, Maclean J, Moule S, Murphy L, Oliver K, Quail MA, et al: Massive gene decay in the leprosy bacillus. Nature. 2001, 409: 1007-1011. 10.1038/35059006.
Andersson SG, Zomorodipour A, Andersson JO, Sicheritz-Pontén T, Alsmark UC, Podowski RM, Näslund AK, Eriksson AS, Winkler HH, Kurland CG: The genome sequence of Rickettsia prowazekii and the origin of mitochondria. Nature. 1998, 396: 133-140. 10.1038/24094.
Pérez-Brocal V, Gil R, Ramos S, Lamelas A, Postigo M, Michelena JM, Silva FJ, Moya A, Latorre A: A small microbial genome: the end of a long symbiotic relationship?. Science. 2006, 314: 312-313. 10.1126/science.1130441.
Stephens RS, Kalman S, Lammel C, Fan J, Marathe R, Aravind L, Mitchell W, Olinger L, Tatusov RL, Zhao Q, Koonin EV, Davis RW: Genome sequence of an obligate intracellular pathogen of humans: Chlamydia trachomatis. Science. 1998, 282: 754-759. 10.1126/science.282.5389.754.
Zomorodipour A, Andersson SG: Obligate intracellular parasites: Rickettsia prowazekii and Chlamydia trachomatis. FEBS Lett. 1999, 452: 11-15. 10.1016/S0014-5793(99)00563-3.
Carlson JH, Porcella SF, McClarty G, Caldwell HD: Comparative genomic analysis of Chlamydia trachomatis oculotropic and genitotropic strains. Infect Immun. 2005, 73: 6407-6418. 10.1128/IAI.73.10.6407-6418.2005.
Thomson NR, Holden MT, Carder C, Lennard N, Lockey SJ, Marsh P, Skipp P, O'Connor CD, Goodhead I, Norbertzcak H, Harris B, Ormond D, Rance R, Quail MA, Parkhill J, Stephens RS, Clarke IN: Chlamydia trachomatis: genome sequence analysis of lymphogranuloma venereum isolates. Genome Res. 2007, 18: 161-171. 10.1101/gr.7020108.
Kari L, Whitmire WM, Carlson JH, Crane DD, Reveneau N, Nelson DE, Mabey DC, Bailey RL, Holland MJ, McClarty G, Caldwell HD: Pathogenic diversity among Chlamydia trachomatis ocular strains in nonhuman primates is affected by subtle genomic variations. J Infect Dis. 2008, 197: 449-456. 10.1086/525285.
Stephens RS: Chlamydiae and evolution: A billion years and counting. Chlamydial Infections: Proceedings of the 10th International Symposium on Human Chlamydial Infections: 16-21 June, 2002; Antalya, Turkey. Edited by: Schachter J, Christiansen G, Clarke IN, Kattenboeck B, Kuo CC, Rank RG, Ridgway GL, Saikku P, Stamm WE, Stephens RS, Summersgill JT, Timms P, Wyrick PB. 2002, San Francisco: International Chlamydia Symposium, 3-12.
Brunelle BW, Sensabaugh GF: The ompA gene in Chlamydia trachomatis differs in phylogeny and rate of evolution from other regions of the genome. Infect Immun. 2006, 74: 578-585. 10.1128/IAI.74.1.578-585.2006.
Gomes JP, Nunes A, Bruno WJ, Borrego MJ, Florindo C, Dean D: Polymorphisms in the nine polymorphic membrane proteins of Chlamydia trachomatis across all serovars: evidence for serovar Da recombination and correlation with tissue tropism. J Bacteriol. 2006, 188: 275-286. 10.1128/JB.188.1.275-286.2006.
Rokas A, Williams BL, King N, Carroll SB: Genome-scale approaches to resolving incongruence in molecular phylogenies. Nature. 2003, 425: 798-804. 10.1038/nature02053.
Wolf YI, Rogozin IB, Koonin EV: Coelomata and not Ecdysozoa: evidence from genome-wide phylogenetic analysis. Genome Res. 2004, 14: 29-36. 10.1101/gr.1347404.
Gadagkar SR, Rosenberg MS, Kumar S: Inferring species phylogenies from multiple genes: concatenated sequence tree versus consensus gene tree. J Exp Zoolog B Mol Dev Evol. 2005, 304: 64-74. 10.1002/jez.b.21026.
Kumar S, Gadagkar SR: Efficiency of the neighbour-joining method in reconstructing deep and shallow evolutionary relationships in large phylogenies. J Mol Evol. 2000, 51: 544-553.
Tamura K, Nei M, Kumar S: Prospects for inferring very large phylogenies by using the neighbour-joining method. Proc Natl Acad Sci USA. 2004, 101: 11030-11035. 10.1073/pnas.0404206101.
Tajima F, Nei M: Estimation of evolutionary distance between nucleotide sequences. Mol Biol Evol. 1984, 1: 269-285.
Read TD, Myers GS, Brunham RC, Nelson WC, Paulsen IT, Heidelberg J, Holtzapple E, Khouri H, Federova NB, Carty HA, Umayam LA, Haft DH, Peterson J, Beanan MJ, White O, Salzberg SL, Hsia RC, McClarty G, Rank RG, Bavoil PM, Fraser CM: Genome sequence of Chlamydophila caviae (Chlamydia psittaci GPIC): examining the role of niche-specific genes in the evolution of the Chlamydiaceae. Nucleic Acids Res. 2003, 31: 2134-2147. 10.1093/nar/gkg321.
Dugan J, Rockey DD, Jones L, Andersen AA: Tetracycline resistance in Chlamydia suis mediated by genomic islands inserted into the chlamydial inv-like gene. Antimicrob Agents Chemother. 2004, 48: 3989-3995. 10.1128/AAC.48.10.3989-3995.2004.
Thomson NR, Yeats C, Bell K, Holden MT, Bentley SD, Livingstone M, Cerdeño-Tárraga AM, Harris B, Doggett J, Ormond D, Mungall K, Clarke K, Feltwell T, Hance Z, Sanders M, Quail MA, Price C, Barrell BG, Parkhill J, Longbottom D: The Chlamydophila abortus genome sequence reveals an array of variable proteins that contribute to interspecies variation. Genome Res. 2005, 15: 629-640. 10.1101/gr.3684805.
Stothard DR, Pol Van Der B, Smith NJ, Jones RB: Effect of serial passage in tissue culture on sequence of omp1 from Chlamydia trachomatis clinical isolates. J Clin Microbiol. 1998, 36: 3686-3688.
Carlson JH, Hughes S, Hogan D, Cieplak G, Sturdevant DE, McClarty G, Caldwell HD, Belland RJ: Polymorphisms in the Chlamydia trachomatis cytotoxin locus associated with ocular and genital isolates. Infect Immun. 2004, 72: 7063-7072. 10.1128/IAI.72.12.7063-7072.2004.
Fehlner-Gardiner C, Roshick C, Carlson JH, Hughes S, Belland RJ, Caldwell HD, McClarty G: Molecular basis defining human Chlamydia trachomatis tissue tropism. A possible role for tryptophan synthase. J Biol Chem. 2002, 277: 26893-26903. 10.1074/jbc.M203937200.
Caldwell HD, Wood H, Crane D, Bailey R, Jones RB, Mabey D, Maclean I, Mohammed Z, Peeling R, Roshick C, Schachter J, Solomon AW, Stamm WE, Suchland RJ, Taylor L, West SK, Quinn TC, Belland RJ, McClarty G: Polymorphisms in Chlamydia trachomatis tryptophan synthase genes differentiate between genital and ocular isolates. J Clin Invest. 2003, 111: 1757-1769.
Stecher B, Robbiani R, Walker AW, Westendorf AM, Barthel M, Kremer M, Chaffron S, Macpherson AJ, Buer J, Parkhill J, Dougan G, von Mering C, Hardt WD: Salmonella enterica serovar typhimurium exploits inflammation to compete with the intestinal microbiota. PLoS Biol. 2007, 5: 2177-2189. 10.1371/journal.pbio.0050244.
Gomes JP, Bruno WJ, Nunes A, Santos N, Florindo C, Borrego MJ, Dean D: Evolution of Chlamydia trachomatis diversity occurs by widespread interstrain recombination involving hotspots. Genome Res. 2007, 17: 50-60. 10.1101/gr.5674706.
Suzuki DT, Griffiths AJF, Miller JH, Lewontin RC: An Introduction to Genetic Analysis. 1989, New York: WH Freeman and Company
Caldwell HD, Kromhout J, Schachter J: Purification and partial characterization of the major outer membrane protein of Chlamydia trachomatis. Infect Immun. 1981, 31: 1161-1176.
Su H, Raymond L, Rockey DD, Fischer E, Hackstadt T, Caldwell HD: A recombinant Chlamydia trachomatis major outer membrane protein binds to heparan sulfate receptors on epithelial cells. Proc Natl Acad Sci USA. 1996, 93: 11143-11148. 10.1073/pnas.93.20.11143.
Gomes JP, Borrego MJ, Atik B, Santo I, Azevedo J, Brito de Sá A, Nogueira P, Dean D: Correlating Chlamydia trachomatis infectious load with genital ecological success and disease pathogenesis. Microbes Infect. 2006, 8: 16-26. 10.1016/j.micinf.2005.05.014.
Brunham R, Yang C, Maclean I, Kimani J, Maitha G, Plummer F: Chlamydia trachomatis from individuals in a sexually transmitted disease core group exhibit frequent sequence variation in the major outer membrane protein (omp1) gene. J Clin Invest. 1994, 94: 458-463. 10.1172/JCI117347.
World Health Organization. [http://www.who.int]
Catry MA, Borrego MJ, Cardoso J, Azevedo J, Santo I: Comparison of the Amplicor Chlamydia trachomatis test and cell culture for the detection of urogenital chlamydial infections. Genitourin Med. 1995, 71: 247-250.
MEGA 3.1 software. [http://www.megasoftware.net]
SimPlot 3.5.1 software. [http://sray.med.som.jhmi.edu/SCRoftware/]
Read TD, Brunham RC, Shen C, Gill SR, Heidelberg JF, White O, Hickey EK, Peterson J, Utterback T, Berry K, Bass S, Linher K, Weidman J, Khouri H, Craven B, Bowman C, Dodson R, Gwinn M, Nelson W, DeBoy R, Kolonay J, McClarty G, Salzberg SL, Eisen J, Fraser CM: Genome sequences of Chlamydia trachomatis MoPn and Chlamydia pneumoniae AR39. Nucleic Acids Res. 2000, 28: 1397-1406. 10.1093/nar/28.6.1397.
SWAAP 1.0.2 software. [http://www.bacteriamuseum.org/SWAAP/SwaapPage.htm]
Kumar S, Gadagkar SR: Disparity Index: a simple statistic to measure and test the homogeneity of substitution patterns between molecular sequences. Genetics. 2001, 158: 1321-1327.
Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4: 406-425.
Kimura M: A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980, 16: 111-120. 10.1007/BF01731581.
Jukes TH, Cantor CR: Evolution of protein molecules. Mammalian Protein Metabolism. Edited by: Munro HN. 1969, New York:Academic Press, 21-132.
Tamura K, Nei M: Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol. 1993, 10: 512-526.
Nei M, Kumar S: Phylogenetic inference: Maximum Parsimony methods. Molecular Evolution and Phylogenetics. Edited by: Nei M, Kumar S. 2000, New York: Oxford University Press, 115-145.
Feil EJ, Holmes EC, Bessen DE, Chan MS, Day NP, Enright MC, Goldstein R, Hood DW, Kalia A, Moore CE, Zhou J, Spratt BG: Recombination within natural populations of pathogenic bacteria: short-term empirical estimates and long-term phylogenetic consequences. Proc Natl Acad Sci USA. 2001, 98: 182-187. 10.1073/pnas.98.1.182.
Spratt BG, Hanage WP, Feil EJ: The relative contributions of recombination and point mutation to the diversification of bacterial clones. Curr Opin Microbiol. 2001, 4: 602-606. 10.1016/S1369-5274(00)00257-5.
Millman KL, Tavare S, Dean D: Recombination in the ompA gene but not the omcB gene of Chlamydia contributes to serovar-specific differences in tissue tropism, immune surveillance, and persistence of the organism. J Bacteriol. 2001, 183: 5997-6008. 10.1128/JB.183.20.5997-6008.2001.
Statistical platform R 2.5.1. [http://www.r-project.org/]
The Institute for Genomic Research (TIGR). [http://cmr.tigr.org]
We would like to thank Dr Brendan Wren and Dr Bush for useful comments, critical discussions and reading of the manuscript prior to publication. This work was supported by grants from Fundação para a Ciência e a Tecnologia (FCT) (PTDC/BIA-BCM/71117/2006) and Comissão de Fomento da Investigação em Cuidados de Saúde (n°112/2007) to JPG. AN is the recipient of a PhD Grant (SFRH/BD/25651/2005) from FCT.
AN performed the experimental work, analyzed and interpreted the data, performed the statistical analysis, and wrote the manuscript. PJN performed the statistical analysis. MJB contributed to the experimental work. JPG designed the project, obtained the funding, interpreted the data, and wrote the manuscript.
Electronic supplementary material
About this article
Cite this article
Nunes, A., Nogueira, P.J., Borrego, M.J. et al. Chlamydia trachomatisdiversity viewed as a tissue-specific coevolutionary arms race. Genome Biol 9, R153 (2008) doi:10.1186/gb-2008-9-10-r153
- Horizontal Gene Transfer
- Additional Data File
- Chlamydia Trachomatis
- Neighbor Join