- Research article
- Open Access
Multi locus sequence typing of Chlamydiales: clonal groupings within the obligate intracellular bacteria Chlamydia trachomatis
BMC Microbiology volume 8, Article number: 42 (2008)
The obligate intracellular growing bacterium Chlamydia trachomatis causes diseases like trachoma, urogenital infection and lymphogranuloma venereum with severe morbidity. Several serovars and genotypes have been identified, but these could not be linked to clinical disease or outcome. The related Chlamydophila pneumoniae, of which no subtypes are recognized, causes respiratory infections worldwide. We developed a multi locus sequence typing (MLST) scheme to understand the population genetic structure and diversity of these species and to evaluate the association between genotype and disease.
A collection of 26 strains of C. trachomatis of different serovars and clinical presentation and 18 strains of C. pneumoniae were included in the study. For comparison, sequences of C. abortus, C. psittaci, C. caviae, C. felis, C. pecorum (Chlamydophila), C. muridarum (Chlamydia) and of Candidatus protochlamydia and Simkania negevensis were also included. Sequences of fragments (400 – 500 base pairs) from seven housekeeping genes (enoA, fumC, gatA, gidA, hemN, hlfX, oppA) were analysed. Analysis of allelic profiles by eBurst revealed three non-overlapping clonal complexes among the C. trachomatis strains, while the C. pneumoniae strains formed a single group. An UPGMA tree produced from the allelic profiles resulted in three groups of sequence types. The LGV strains grouped in a single cluster, while the urogenital strains were distributed over two separated groups, one consisted solely of strains with frequent occurring serovars (E, D and F). The distribution of the different serovars over the three groups was not consistent, suggesting exchange of serovar encoding ompA sequences. In one instance, exchange of fumC sequences between strains of different groups was observed. Cluster analyses of concatenated sequences of the Chlamydophila and Chlamydia species together with those of Candidatus Protochlamydia amoebophila and Simkania negevensis resulted in a tree identical to that obtained with 23S RNA gene sequences.
These data show that C. trachomatis and C. pneumoniae are highly uniform. The difference in genetic diversity between C. trachomatis and C. pneumoniae is in concordance with a later assimilation to the human host of the latter. Our data supports the taxonomy of the order of Chlamydiales.
Chlamydia trachomatis is the world's leading cause of preventable blindness. Also, C. trachomatis is considered the world's most common sexually transmitted bacterial pathogen. Many urogenital infections remain unnoticed, constituting a large reservoir of untreated individuals, a continuous threat for transmission of this pathogen. When not treated in time, infection with C. trachomatis can lead to infertility in women. C. trachomatis strains are discriminated by serotyping based on the antigenic difference between the major outer membrane proteins (MOMP). Nineteen serovars have been described: A, B, Ba, C (mainly seen among isolates from trachoma infections) D, Da, E, F, G, Ga, H, I, Ia, J, Ja, K, (urogenital infections) and L1, L2, L2a and L3 causing lymphogranuloma venereum (LGV). Among urogenital infections, serovars D – F are most frequently found . However, serotyping is laborious, needing culture and a large panel of antibodies [2, 3]. To overcome these drawbacks a PCR based RFLP of ompA was developed for the identification of genotypes corresponding to serovars [4–6]. Using this method genotypes were categorised into three geno-groups: the B group (B, E, D, Da, L1, L2, L2a), the C group (C, A, H, I, Ia, J, K, L3) and the intermediate group (F, G, Ga). Except for an immunological relationship between members of a group, the biological relevance of the geno-groups remains obscure.
Chlamydophila pneumoniae is a common cause of community-acquired pneumonia, bronchitis, pharyngitis and sinusitis . Although C. pneumoniae often causes mild or subclinical infections, its persistence in the host can lead to the establishment of chronic pathologies and an increasing number of reports indicate an association between persistent C. pneumoniae infections and arteriosclerosis  or coronary heart diseases [9, 10]. A robust typing scheme for C. pneumoniae is lacking.
Together with C. trachomatis, C. pneumoniae belongs to the family of Chlamydiaceae in the order of Chlamydiales. Based on phylogenetic analyses of 16S and 23S rRNA gene sequences, C. trachomatis, Chlamydia suis and Chlamydia muridarum all belong to the genus Chlamydia, while C. pneumoniae, Chlamydophila psittaci, Chlamydophila pecorum, Chlamydophila felis, Chlamydophila abortus, and Chlamydophila caviae all belong to the family of Chlamydophila [11–13]. Other family members of the order of Chlamydiales are Parachlamydiaceae and Simkaniaceae.
Currently, the typing scheme for C. trachomatis is based on epitopes in the major outer membrane protein (MOMP). Variants of this protein are subjected to selection and isolates of the same serovar may not be closely related [14, 15]. Here we present an MLST typing scheme using gene segment sequences of seven housekeeping genes. These genes were selected using the criteria that they are widely separated on the chromosome and not adjacent to putative outer membrane, secreted, or hypothetical proteins that might be under diversifying selection. In addition, each locus has a similar extent of nucleotide substitutions to ensure consistency . The results identified three sub-groupings within C. trachomatis, but no subdivision within C. pneumoniae. A phylogenetic tree based on the concatenated sequences of six of the housekeeping gene fragments is consistent with a tree based on 16S and 23S rRNA gene sequences.
MLST of C. trachomatis and C. pneumoniae
Analogous to the MLST schemes of e.g. Neisseria meningitidis  and Streptococcus pneumoniae  fragments of seven housekeeping genes scattered around the chromosome of C. trachomatis and C. pneumoniae were obtained (Tables 1 and 2). The gene order on the chromosome of both species is identical. None of the sequences of the seven different loci among the C. trachomatis strains contained gaps after alignment.
Variation among the sequences of the seven loci was very limited. In C. trachomatis the highest number (three) of synonymous substitutions was seen in oppA while the highest number (four) of non-synonymous substitutions was seen in fumC. Analogous to other MLST schemes, we assigned allele numbers to each unique allele sequence for each house-keeping gene . The number of alleles per locus varied between two and six. Most of the alleles were seen more than once. However, among the oppA alleles four unique alleles were found while among the gidA and enoA sequences one unique allele was observed.
For each isolate, the alleles at each of the seven loci define the allelic profile or sequence type (ST). Among the 26 C. trachomatis strains 15 ST's could be assigned. Analysis by eBurst revealed three non-overlapping groups or clonal complexes, consisting of related strains sharing identical alleles at six of the seven loci with at least one other member of the group (Figure 1). An UPGMA cluster analysis showed the same groups (Figure 2A). SplitsTree decomposition demonstrated that alternative routes of descent in the tree resulted in the same groupings (Figure 3). An UPGMA cluster analyses of the concatenated sequences of the seven gene fragments yielded the same groupings as when allelic profiles were used (Figure 2B), while SplitsTree decomposition analysis yielded a more simpler network but with the same groupings as with the distances matrix of allelic profiles.
Group I, with ST13 as the putative founder, defined as the ST with the most single locus variants, and ST6 as the founder of a subgroup consisted of C. trachomatis strains isolated from patients with urogenital infection (serovars D to K) as well as trachoma infections (serovars A, Ba C). The latter formed a separate branch in the UPGMA cluster analyses. Group II comprises the LGV strains (serovar L) and strain B/TW-5 (serovar B). Group III is formed by all, except one, serovar E strains and one serovar F and two serovar D strains.
Sequence variation among C. pneumoniae was far less than among C. trachomatis (Table 1). Substitutions were only seen among the sequences of gidA (1 synonymous) and enoA (1 synonymous and 1 non-synonymous). This means that all 16 strains shared identical alleles at least five of the seven loci, i.e. C. pneumoniae appeared to be highly uniform. In addition, none of the alleles in C. trachomatis and C. pneumoniae are the same.
Recombination in C. trachomatis
In group II, the allelic profile of ST2 (B/TW-5) differs from that of ST11 (serovar L strains) at one locus (Table 3). It shares the fumC allele with the majority of the strains in group I and II, indicating exchanges of the fumC sequences between a strain with genotype ST11 and a strain with genotype other than ST11 or ST5 (with fumC allele different from that of ST1 and all other ST's; Table 3), resulting the B/TW-5 strain with genotype ST2. The difference in fumC sequences are three substitutions in an 89 bp region, albeit that all three appear to be non-synonymous (Figure 4).
In addition, while ST11 strains are serovar L, ST 2 is serovar B, indication exchange of ompA (encoding MOMP, defining the serovar type) sequences between a serovar B strain and a serovar L strain. Other indications of recombination between different C. trachomatis genotypes and exchange of ompA sequences might be inferred from the position of the only serovar E (serovar E11A, ST6) in group I, while all other serovar E strains cluster in group III.
MLSA based phylogeny of Chlamydiales
The oppA sequences of C. pneumoniae contained several indels when compared to the sequences of C. trachomatis and other species of Chlamydiales. All genomes of the Chlamydiales contain multiple copies of oppA genes. In each genome, these copies are highly homologous, but vary between the different species, making selection of the right oppA copy from these genome sequences indecisive.
An insert was observed in the enoA sequences of Candidatus protochlamydia and Simkania negevensis. Also, the hemN sequences of these strains contained small indels as compared to the hemN sequences of the other members of the Chlamydiales. Small indels were also observed among the hlfX sequences of C. abortus, C. caviae, C. felis, C. psittaci, Candidatus protochlamydia and Simkania negevensis. Recently, multilocus sequence analysis (MLSA) was introduced to study relatedness of closely related species [19, 20]. In this analysis the sequences of multi locus housekeeping fragments are concatenated and used in cluster analysis. Phylogenetic analysis of the Chlamydiales by Neighbour-Joining method of the aligned concatenated sequences of the housekeeping gene fragments, except that of oppA, resulted in a tree (Fig 5A) comparable to that obtained with 16S rRNA gene and 23S rRNA gene (Fig 5B) sequences .
To assess the population structure of C. trachomatis and C. pneumoniae sequences of fragments of seven housekeeping genes, obtained from 26 C. trachomatis strains and 18 C. pneumoniae strains, were analysed. C. pneumoniae appeared to be highly uniform. Among the C. trachomatis strains three very coherent clonal complexes were observed, consisting of strains sharing identical alleles of at least 6 of the 7 loci with one other member of the group. C. pneumoniae appeared to be highly uniform.
Recently, an MLST scheme has been published, in which five target regions were selected based on their relatively high variability as compared to the rest of the genome and analysed. In addition, these targets were not widely separated on the genome . That typing scheme was intended to be highly discriminatory and to be applied in contact tracing.
In the present study, 7 housekeeping genes were chosen, which were widely separated on the chromosome and not adjacent to putative outer membrane, secreted, or hypothetical proteins that might be under diversifying selection. In addition, each locus has a similar level of variation in terms of nucleotide substitutions to ensure consistency . Fifteen sequence types were found among 26 C. trachomatis isolates (0.6 STs per isolate). Many organisms show more diversity, i.e. more sequence types per isolate, but a significant number of organisms shows comparable or less sequence types per isolate . In these organisms MLST has been applied for strain characterization and epidemiological surveillance (e.g. Listeria monocytogenes, 29 STs among 62 isolates; Acinetobacter baumannii, 20 STs among 49 isolates), population structure analyses (e.g. Porphyromonas gingivalis, 50 STs in 59 isolates) and evolutionary analyses studies (e.g. Batrachochytrium dendrobatidis, 10 STs in 35 isolates) . The here presented MLST scheme for C. trachomatis may be similarly used. Our results show three clonal complexes among C. trachomatis of which one is associated with LGV.
A phylogeny tree based on the concatenated sequences of 6 loci resulted in a tree consistent with that of obtained when 16S rRNA and 23S rRNA genes were used in the phylogeny analyses . This approach, using the concatenated sequences to study the relationships among strains of similar species was recently termed multilocus sequence analysis (MLSA)  and has successfully been applied to other species .
C. trachomatis show limited variation; the average number of synonymous substitutions in C. trachomatis is in the same order as that in Yersinia pseudotuberculosis [22–24]. In contrast, the average number of synonymous substitutions in C. pneumoniae is even smaller, comparable to that in Vibrio sonnei and Yersinia pestis, but larger than in Mycobacterium tuberculosis [23–25]. This may suggest that both species C. trachomatis and C. pneumoniae are evolutionarily young or recently past severe bottle necks .
Three clonal complexes were seen among the 26 strains of C. trachomatis; each group includes isolates that differ at only one locus from at least one other isolates within the group. Singeltons, differing at two or more loci from all other isolates were not observed. Our data provided some evidence of recombination, e.g. exchange of the MOMP (serovar determining) encoding ompA sequences and of fumC sequences. Discongruence between ompA and the main part of the genome has also been observed by Gomes and colleagues and Brunelle and Sensabaugh [15, 27]. In addition, earlier reports of mosaic ompA gene structures indicated that ompA or parts of ompA do exchange between C. trachomatis strains [28–30].
Brunelle and Sensabaugh observed recombination in ompA genes, pmpE genes and pmpH genes, but not in the remainder of the genome . Recently published data by Gomes and colleagues suggested frequent recombination in C. trachomatis, albeit that this recombination occurred at hotspot near or in ompA and pmp genes . Here we demonstrated in at least one instance recombination in or near the housekeeping fumC. The allelic profile of ST2 was identical to that of ST11 with the exception of the fumC allele. The fumC allele of ST2 was identical to that of the majority of the other C. trachomatis strains, while that of the ST11 strains differed at three positions within 87 nucleotides, suggesting uptake and recombination of (a part of) fumC sequence by an ST11 genotype C. trachomatis resulting in ST2 genotype. It is unlikely that pmp or ompA sequences are involved in the exchange of fumC sequences, since the nearest pmp genes are 54 Kbp upstream (pmpD, cta0884) and 23 Kbp upstream (pmpE, cta0949) of fumC.
The three clonal complexes or groups are partly associated with tissue tropism. All LGV causing strains group together in group II. The urogenital strains and ocular strains are distributed over two groups, albeit that the ocular strains group together with the less frequent occurring urogenital strains (serovar H to K). In addition, the trachoma strains form a separated branch within group I. The more frequently occurring urogenital strains formed the separate group III. High frequency occurring genotypes may be linked with symptomatic infection, but in a study among woman with urogenital C. trachomatis infections serovar E and F strains were equally isolated from patients with symptoms and from patients without symptoms . Hence, host factors may determine disease outcome. Alternatively, the high frequency occurring genotypes may be associated with higher transmission rates.
The C. pneumoniae population is highly uniform, while that of C. trachomatis shows three clonal complexes based on an MLST scheme of 7 housekeeping genes. More clonal groups may be identified when more strains will be analysed with this scheme. The difference in genetic diversity between C. trachomatis and C. pneumoniae is in concordance with a later assimilation to the human host of the latter.
Twenty-four C. trachomatis strains were used, consisting of reference strains: (A/Sa-1, H/UW-4, I', D-, Ba/Apache-2, C/UW-1, Da/mt-566, I/UW-12, K/UW-31, G/IOL-238, J/UW-36, L1/440-L, L2/434-B, L3/404-L, B/TW5, D/IC-CAL-8, D', F/MRC-301 and E/DK-20 [32–34] and additional new isolates from patients, E4a, E7a, E10a, E11a, and E12a [see additional file 1]. The serovar of these isolates was confirmed by RFLP [32, 35]. We also tested 14 reference strains of C. pneumoniae (CM-1, IOL-207, PS-32, AR-338, BAL-16, CWL-011, CWL-50, GRO-21, H-12, K-7, NWL-1, UZG-1, 2023, 2043)  [see additional file 1]. In addition, we also analyzed C. trachomatis strains D/UW-3/CX (Accession no. AE001273) , A/HAR-13 (Accession no. NC_007429)  and C. pneumoniae strain CWL-029 (Accession no. NC_000922) , AR-39 (Accession no. AE002161) , TW-183 (Accession no. AE017160) and J138 (Accession no. BA000008)  whose complete genome sequences are available in the database. The collection of C. trachomatis and C. pneumoniae strains represents all known serovars and were from patients with various disease outcomes (see additional file 1). In addition, the corresponding sequences from C. muridarum Nigg. (Accession no. AE002160) (Chlamydia)  and of C. abortus S26/3 (Accession no. NC_004552) , C. caviae GPIC (Accession no. AE015925) , C. felis Fe/C-56 (Accession no. AP006861)  and of Candidatus protochlamydia amoebophila UWE25 (Accession no. BX908798) (Parachlamydiaceae)  were obtained from publicly accessible databases. The corresponding sequences from C. pecorum E58 (McNutt), C. psittaci 6BC and Simkania negevensis (Simkaniaceae) were obtained from TIGR's unfinished genomes database.
DNA, genes, PCR products and sequences
DNA was extracted from elementary bodies from cultures of C. trachomatis or C. pneumoniae according to Boom et al . Fragments of 7 genes, i.e. gatA, oppA3, hflX, gidA, enoA, hemN and fumC encoding aspartyl/glutamyl-tRNA amidotransferase subunit A, oligopeptide-binding protein, GTP-binding protein, tRNA (uracil-5-)-methyltransferase, enolase, coproporphyrinogen III oxidase and fumarate hydratase, respectively were amplified using the oligonucleotide primers shown in Table 4. Amplification primers were designed based on the genome sequence of A/HAR-13 to yield amplicons that were short enough to obtain complete double stranded sequences in two single sequence runs. Each sequence run was performed from a different PCR amplicon and sequence traces were obtained with ABI Big-dyes and an ABI 3730 sequencer.
Phylogenetic and other analyses
The number of synonymous and non-synonymous substitutions per site was determined using DnaSP 4.0 . For C. trachomatis, unique sequences were assigned allele numbers using the Non-redundant databases (NRDB) program . Allele profile data were analysed in eBurst to define clonal complexes or groups [47, 48]. Groups were defined as sets of related strains containing pairs of strains that share at least six identical alleles at the seven loci.
A distance matrix in Nexus format was generated from the set of allelic profiles using SplitsTree . This file was then used for phylogenetic analyses in SplitsTree 4.0 , both by generating an UPGMA tree and by SplitsTree decomposition analyses. Decomposition analysis depicts all the shortest pathways linking sequences, including those that produce an interconnected network.
Phylogenetic evolutionary analyses of the sequences of the different members of Chlamydiales were conducted using MEGA version 3.1 .
Spaargaren J, Verhaest I, Mooij S, Smit C, Fennema HS, Coutinho RA, Salvador Pena A, Morre SA: Analysis of Chlamydia trachomatis. serovar distribution changes in the Netherlands (1986–2002). Sex Transm Infect. 2004, 80: 151-152. 10.1136/sti.2003.006395.
Wang SP, Kuo CC, Barnes RC, Stephens RS, Grayston JT: Immunotyping of Chlamydia trachomatis with monoclonal antibodies. J Infect Dis. 1985, 152: 791-800.
Ossewaarde JM, Rieffe M, de Vries A, Derksen-Nawrocki RP, Hooft HJ, van Doornum GJ, van Loon AM: Comparison of two panels of monoclonal antibodies for determination of Chlamydia trachomatis serovars. J Clin Microbiol. 1994, 32: 2968-2974.
Frost EH, Deslandes S, Veilleux S, Bourgaux-Ramoisy D: Typing Chlamydia trachomatis by detection of restriction fragment length polymorphism in the gene encoding the major outer membrane protein. J Infect Dis. 1991, 163: 1103-1107.
Sayada C, Denamur E, Orfila J, Catalan F, Elion J: Rapid genotyping of the Chlamydia trachomatis major outer membrane protein by the polymerase chain reaction. FEMS Microbiol Lett. 1991, 67: 73-78. 10.1111/j.1574-6968.1991.tb04392.x.
Rodriguez P, Vekris A, de Barbeyrac B, Dutilh B, Bonnet J, Bebear C: Typing of Chlamydia trachomatis by restriction endonuclease analysis of the amplified major outer membrane protein gene. J Clin Microbiol. 1991, 29: 1132-1136.
Kuo CC, Jackson LA, Campbell LA, Grayston JT: Chlamydia pneumoniae (TWAR). Clin Microbiol Rev. 1995, 8: 451-461.
Campbell LA, Kuo CC: Chlamydia pneumoniae and atherosclerosis. Semin Respir Infect. 2003, 18: 48-54. 10.1053/srin.2003.50006.
Mussa FF, Chai H, Wang X, Yao Q, Lumsden AB, Chen C: Chlamydia pneumoniae and vascular disease : an update. J Vasc Surg. 2006, 43: 1301-7. 10.1016/j.jvs.2006.02.050.
Grayston JT: Background and current knowledge of Chlamydia pneumoniae and atherosclerosis. J Infect Dis. 2000, 181: S402-S410. 10.1086/315596.
Everett KDE, Bush RM, Andersen AA: Emended description of the order Chlamydiales, proposal of ParaChlamydiaceae fam. nov. and Simkaniaceae fam. nov., each containing one monotypic genus, revised taxonomy of the family Chlamydiaceae, including a new genus and five new species, and standards for the identification of organisms. Int J Syst Bacteriol. 1999, 49: 415-440.
Bush RM, Everett KD: Molecular evolution of the Chlamydiaceae. Int J Syst Evol Microbiol. 2001, 51: 203-220.
Herrmann B, Pettersson B, Everett KD, Mikkelsen NE, Kirsebom LA: Characterization of the rnpB gene and RNase P RNA in the order Chlamydiales. Int J Syst Evol Microbiol. 2000, 50: 149-58.
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.
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.
Maiden MCJ: Multilocus sequence typing of bacteria. Annu Rev Microbiol. 2006, 60: 561-588. 10.1146/annurev.micro.59.030804.121325.
Maiden MC, Bygraves JA, Feil E, Morelli G, Russell JE, Urwin R, Zhang Q, Zhou J, Zurth K, Caugant DA, Feavers IM, Achtman M, Spratt BG: Multilocus sequence typing: a portable approach to the identification of clones within populations of pathogenic microorganisms. Proc Natl Acad Sci USA. 1998, 95: 3140-3145. 10.1073/pnas.95.6.3140.
Enright MC, Spratt BG: A multilocus sequence typing scheme for Streptococcus pneumoniae : identification of clones associated with serious invasive disease. Microbiology. 1998, 144: 3049-60.
Gevers D, Cohan FM, Lawrence JG, Spratt BG, Coenye T, Feil EJ, Stackebrandt E, Van de Peer Y, Vandamme P, Thompson FL, Swings J: Opinion: Re-evaluating prokaryotic species. Nat Rev Microbiol. 2005, 3: 733-739. 10.1038/nrmicro1236.
Hanage WP, Spratt BG, Turner KM, Fraser C: Modelling bacterial speciation. Philos Trans R Soc Lond B Biol Sci. 2006, 361: 2039-2044. 10.1098/rstb.2006.1926.
Klint M, Fuxelius HH, Goldkuhl RR, Skarin H, Rutemark C, Andersson SG, Persson K, Herrmann B: High-Resolution Genotyping of Chlamydia trachomatis Strains by Multilocus Sequence Analysis. J Clin Microbiol. 2007, 45: 1410-1414. 10.1128/JCM.02301-06.
Achtman M, Zurth K, Morelli G, Torrea G, Guiyoule A, Carniel E: Yersinia pestis, the cause of plague, is a recently emerged clone of Yersinia pseudotuberculosis. Proc Natl Acad Sci USA. 1999, 96: 14043-14048. 10.1073/pnas.96.24.14043.
Achtman M: Population structure of pathogenic bacteria revisited. Int J Med Microbiol. 2004, 294: 67-73. 10.1016/j.ijmm.2004.06.028.
Achtman M, Morelli G, Zhu P, Wirth T, Diehl I, Kusecek B, Vogler AJ, Wagner DM, Allender CJ, Easterday WR, Chenal-Francisque V, Worsham P, Thomson NR, Parkhill J, Lindler LE, Carniel E, Keim P: Microevolution and history of the plague bacillus, Yersinia pestis. Proc Natl Acad Sci USA. 2004, 101: 17837-17842. 10.1073/pnas.0408026101.
Sreevatsan S, Pan X, Stockbauer KE, Connell ND, Kreiswirth BN, Whittam TS, Musser JM: Restricted structural gene polymorphism in the Mycobacterium tuberculosis complex indicates evolutionarily recent global dissemination. Proc Natl Acad Sci USA. 1997, 94: 9869-9874. 10.1073/pnas.94.18.9869.
Rattei T, Ott S, Gutacker M, Rupp J, Maass M, Schreiber S, Solbach W, Wirth T, Gieffers J: Genetic diversity of the obligate intracellular bacterium Chlamydophila pneumoniae by genome-wide analysis of single nucleotide polymorphisms : evidence for highly clonal population structure. BMC Genomics. 2007, 8: 355-346. 10.1186/1471-2164-8-355.
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.
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.
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.
Hayes LJ, Yearsley P, Treharne JD, Ballard RA, Fehler GH, Ward ME: Evidence for naturally occurring recombination in the gene encoding the major outer membrane protein of lymphogranuloma venereum isolates of Chlamydia trachomatis. Infect Immun. 1994, 62: 5659-5663.
Pannekoek Y, van der Ende A, Eijk PP, van Marle J, de Witte MA, Ossewaarde JM, van den Brule AJ, Morre SA, Dankert J: Normal IncA expression and fusogenicity of inclusions in Chlamydia trachomatis isolates with the incA I47T mutation. Infect Immun. 2001, 69: 4654-4656. 10.1128/IAI.69.7.4654-4656.2001.
Lan J, Walboomers JM, Roosendaal R, van Doornum GJ, MacLaren DM, Meijer CJ, van den Brule AJ: Direct detection and genotyping of Chlamydia trachomatis in cervical scrapes by using polymerase chain reaction and restriction fragment length polymorphism analysis. J Clin Microbiol. 1993, 31: 1060-1065.
Meijer A, Morre SA, van den Brule AJ, Savelkoul PH, Ossewaarde JM: Genomic relatedness of Chlamydia isolates determined by amplified fragment length polymorphism analysis. J Bacteriol. 1999, 181: 4469-4475.
Jurstrand M, Falk L, Fredlund H, Lindberg M, Olcén P, Andersson S, Persson K, Albert J, Bäckman A: Characterization of Chlamydia trachomatis omp1 genotypes among sexually transmitted disease patients in Sweden. J Clin Microbiol. 2001, 39: 3915-3919. 10.1128/JCM.39.11.3915-3919.2001.
Morré SA, Moes R, Van Valkengoed I, Boeke JP, van Eijk JT, Meijer CJ, Van den Brule AJ: Genotyping of Chlamydia trachomatis in urine specimens will facilitate large epidemiological studies. J Clin Microbiol. 1998, 36: 3077-3078.
Stephens RS, Kalman S, Lammel CJ, Fan J, Marathe R, Aravind L, Mitchell WP, 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.
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.
Kalman S, Mitchell W, Marathe R, Lammel C, Fan J, Hyman RW, Olinger L, Grimwood J, Davis RW, Stephens RS: Comparative genomes of Chlamydia pneumoniae and C. trachomatis. Nat Genet. 1999, 21: 385-389. 10.1038/7716.
Read T, Brunham R, Shen C, Gill S, Heidelberg J, White O, Hickey E, 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 S, Eisen J, Fraser C: Genome sequences of Chlamydia trachomatis MoPn and Chlamydia pneumoniae AR39. Nucleic Acids Res. 2000, 28: 1397-1406. 10.1093/nar/28.6.1397.
Shirai M, Hirakawa H, Kimoto M, Tabuchi M, Kishi F, Ouchi K, Shiba T, Ishii K, Hattori M, Kuhara S, Nakazawa T: Comparison of whole genome sequences of Chlamydia pneumoniae J138 from Japan and CWL029 from USA. Nucleic Acids Res. 2000, 28: 2311-23144. 10.1093/nar/28.12.2311.
Thomson NR, Yeats C, Bell K, Holden MT, Bentley SD, Livingstone M, Cerdeno-Tarraga 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.
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.
Azuma Y, Hirakawa H, Yamashita A, Cai Y, Rahman MA, Suzuki H, Mitaku S, Toh H, Goto S, Murakami T, Sugi K, Hayashi H, Fukushi H, Hattori M, Kuhara S, Shirai M: Genome sequence of the cat pathogen, Chlamydophila felis. DNA Res. 2006, 13: 15-23. 10.1093/dnares/dsi027.
Boom R, Sol CJ, Salimans MM, Jansen CL, Wertheim-van Dillen PM, van der Noordaa J: Rapid and simple method for purification of nucleic acids. J Clin Microbiol. 1990, 28: 495-503.
Rozas J, Sanches-DelBarrio JC, Messeguer X, Rozas R: DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003, 19: 2496-2497. 10.1093/bioinformatics/btg359.
Feil EJ, Li BC, Aanensen DM, Hanage WP, Spratt BG: eBURST: inferring patterns of evolutionary descent among clusters of related bacterial genotypes from multilocus sequence typing data. J Bacteriol. 2004, 186: 1518-1530. 10.1128/JB.186.5.1518-1530.2004.
Spratt BG, Hanage WP, Li B, Aanensen DM, Feil EJ: Displaying the relatedness among isolates of bacterial species – the eBURST approach. FEMS Microbiol Lett. 2004, 241: 129-134. 10.1016/j.femsle.2004.11.015.
Huson DH, Bryant D: Application of phylogenetic networks in evolutionary studies. Mol Biol Evol. 2006, 23: 254-267. 10.1093/molbev/msj030.
Kumar S, Tamura K, Nei M: MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment. Briefings in Bioinformatics. 2004, 5: 150-163. 10.1093/bib/5.2.150.
Mr. Paul Eijk is acknowledged for his technical assistance. The authors like to thank dr. Mark Achtman for suggestions and critical comments on the manuscript. The preliminary sequence data of C. pecorum E58, C. psittaci 6BC and Simkania negevensis were obtained from The Institute for Genomic Research. Sequencing of Simkania negevensis was accomplished with support from TIGR. This research was partly funded by the Sixth Framework Programme of the European Commission, Proposal/Contract no.: 512061 (Network of Excellence 'European Virtual Institute for Functional Genomics of Bacterial Pathogens' ).
YP participated in the design and coordination of the study and helped to draft the manuscript. GM participated in the design of the study and carried out the sequencing. BK participated in the design of the study and carried out the sequencing. SAM participated in the design of the study helped to draft the manuscript. JMO participated in the design and coordination of the study. AAL carried out sequencing. AvdE participated in the design and coordination of the study, did the analyses and interpretation of the sequence data and helped to draft the manuscript. All authors read and approved the final manuscript.