Search tips
Search criteria 


Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
PLoS One. 2010; 5(12): e14179.
Published online 2010 December 2. doi:  10.1371/journal.pone.0014179
PMCID: PMC2996290

Multi Locus Sequence Typing of Chlamydia Reveals an Association between Chlamydia psittaci Genotypes and Host Species

Wenjun Li, Editor


Chlamydia comprises a group of obligate intracellular bacterial parasites responsible for a variety of diseases in humans and animals, including several zoonoses. Chlamydia trachomatis causes diseases such as trachoma, urogenital infection and lymphogranuloma venereum with severe morbidity. Chlamydia pneumoniae is a common cause of community-acquired respiratory tract infections. Chlamydia psittaci, causing zoonotic pneumonia in humans, is usually hosted by birds, while Chlamydia abortus, causing abortion and fetal death in mammals, including humans, is mainly hosted by goats and sheep. We used multi-locus sequence typing to asses the population structure of Chlamydia. In total, 132 Chlamydia isolates were analyzed, including 60 C. trachomatis, 18 C. pneumoniae, 16 C. abortus, 34 C. psittaci and one of each of C. pecorum, C. caviae, C. muridarum and C. felis. Cluster analyses utilizing the Neighbour-Joining algorithm with the maximum composite likelihood model of concatenated sequences of 7 housekeeping fragments showed that C. psittaci 84/2334 isolated from a parrot grouped together with the C. abortus isolates from goats and sheep. Cluster analyses of the individual alleles showed that in all instances C. psittaci 84/2334 formed one group with C. abortus. Moving 84/2334 from the C. psittaci group to the C. abortus group resulted in a significant increase in the number of fixed differences and elimination of the number of shared mutations between C. psittaci and C. abortus. C. psittaci M56 from a muskrat branched separately from the main group of C. psittaci isolates. C. psittaci genotypes appeared to be associated with host species. The phylogentic tree of C. psittaci did not follow that of its host bird species, suggesting host species jumps. In conclusion, we report for the first time an association between C. psittaci genotypes with host species.


Chlamydia comprises a group of obligate intracellular bacterial parasites responsible for a variety of diseases in humans and animals, including several zoonoses. It was proposed in 1999 that the single genus of Chlamydia should be reassigned into two genera, Chlamydia and Chlamydophila, based on clustering analyses of the 16S rRNA and 23S rRNA genes [1], which has not been widely accepted by the chlamydial research community. Recently however, reversion to the single genus Chlamydia was recommended [2], with the Chlamydia nomenclature used here.

Chlamydia trachomatis can cause diseases with severe morbidity, such as trachoma, urogenital infection and lymphogranuloma venereum [3][5]. Several serovars and genotypes have been identified, but which have not been linked to disease or clinical outcome [6], [7]. Chlamydia pneumoniae is a common cause of community-acquired pneumonia, bronchitis, pharyngitis and sinusitis [8]. Although C. pneumoniae often causes mild or subclinical infections, its persistence in the host can lead to the establishment of chronic pathologies and has been implicated with arteriosclerosis [9] and coronary heart diseases [10], [11]. Chlamydia psittaci which can cause zoonotic pneumonia in humans are usually hosted by birds [12][16]. Transmission of C. psittaci from birds to humans is frequently reported and veterinarians, poultry farmers, bird breeders and pet shopkeepers are in particular at high risk [17][27].

Chlamydia abortus has been associated with abortion and fetal death in mammals, including humans, and is hosted by goats, sheep and less frequently by cattle, horses and pigs [28]. The microorganism has also been associated with pneumonia, conjunctivitis, arthritis as well as epididymitis and has been isolated from the faeces of healthy sheep and goats [29][32]. Pregnant women are at risk when exposed to animals infected with C. abortus and may suffer severe infections, including spontaneous abortion [33][36].

Multi-locus sequence typing (MLST) based on the partial sequences of seven housekeeping genes was first used to evaluate the population structure of Neisseria meningitidis [37]. Previously, we developed an MLST scheme to understand the population genetic structure of C. trachomatis and C. pneumoniae and the diversity of these species and to evaluate the association between genotype and disease [7]. In the present study, the MLST scheme was used to evaluate the population structure of C. psittaci and C. abortus. Results of cluster analyses of concatenated sequences of the 7 housekeeping fragments expanded and validated the proposed typing system for all chlamydial species. The results indicated that the 7 housekeeping fragments used in our study were likely representative for the whole genome sequence. C. psittaci genotypes were associated with their host species; the C. psittaci phylogenetic tree however, did not follow that of its host bird species. Furthermore, C. psittaci 84/2334, formerly considered as the missing link between C. abortus and C. psittaci, was clearly typed as C. abortus.


Population structure of Chlamydia

Among 132 Chlamydia isolates 44 sequence types (STs) were identified; 19 STs among 60 C. trachomatis isolates, 4 STs among 18 C. pneumoniae isolates, 12 STs among 34 C. psittaci isolates and 4 among 16 C. abortus isolates (Table S1). Previously, it was shown that in C. trachomatis the homologous recombination rate is low [38]. We determined the index of association for the different Chlamydia species according to Haubold [39]. Significant linkage disequilibrium was detected for C. psittaci, C. abortus, C. trachomatis and C. pneumoniae (Table S2). In addition, we tested the sequences for evidence of recombination using the maximum chi square [40]. Recombination events were not detected in C. pneumoniae and C. abortus sequences. In C. trachomatis, 5 putative recombination events were detected between a pair of hflX alleles (putative recombination site at position 422) and 4 oppA alleles (putative recombination sites at positions 500 and 506), respectively. In C. psittaci, we detected one putative recombination event in hemN (at position 426) (Table S3). Together these data suggest that the role of recombination in diversity of Chlamydia species is low and that Chlamydia species are clonal. In phylogenetic analyses with clonal bacterial species, where mutation is more important than recombination, it is preferable to use concatenated sequences of the MLST loci rather than allelic profiles, because the magnitude of changes between alleles is lost in allelic profiles [41]. This method was subsequently used to study relatedness of closely related species [42], [43], [44]. Phylogenetic analysis of all 132 strains of Chlamydia using the Neighbour-Joining algorithm with the maximum composite likelihood model of the single 3120 base pairs (bp) sequence of the aligned concatenated loci resulted in a tree comparable to that obtained with 16S rRNA gene and 23S rRNA gene sequences [1] and to the previous reported tree based on concatenated sequences of 6 MLST loci [7] (Fig. 1). In the present tree, three main groups are identified among 60 C. trachomatis isolates, consistent with earlier reported results of analyses of 26 C. trachomatis isolates (not shown) [7]. In addition, C. pneumoniae LPCoLN isolated from koala branched separately from C. pneumoniae from patients, consistent with the results of Myers and colleagues [45]. C. psittaci M56 from a muskrat branched separately from the main cluster of C. psittaci isolates, while C. psittaci 84/2334 grouped together with C. abortus, suggesting that 84/2334 belongs to C. abortus.

Figure 1
Phylogenetic analyses of concatenated sequences of 7 housekeeping gene fragments of Chlamydia strains.

Diversity among Chlamydia species

Concatenated sequences of the MLST loci were used to estimate the diversity among and divergence between Chlamydia species. The diversity among C. trachomatis and C. pneumoniae is limited (Table 1). Nucleotide substitutions in protein encoding genes can be either synonymous or non-synonymous (resulting in a changed amino acid). Darwinian selection may lead to the retention of non-synonymous substitutions. The number of synonymous and non-synonymous substitutions may indicate the degree of selection operating in the population. The average number synonymous substitutions per synonymous site dS in C. trachomatis and C. pneumoniae were comparable, while the average number of nonsynonymous substitutions per nonsynonymous site dN was three times higher among C. trachomatis (Table 1). However, considering C. pneumoniae from humans only, than dS in C. pneumonia was much lower than that in C. trachomatis and dN of C. pneumonia even 7 fold lower than that of C. trachomatis, indicating that C. pneumoniae from human is much more clonal than C. trachomatis in accordance with previous observations [7].

Table 1
Diversity of Chlamydia.

The dS and dN among the 34 C. psittaci isolates was much higher than among any other Chlamydia included in this study (Table 1). However, the dS of C. psittaci without M56 and 84/2334 was considerably lower, but still higher than among C. pneumoniae, C. abortus or C. trachomatis. The dS and dN among the 16 C. abortus isolates were comparable to those among the C. pneumoniae isolates, reconfirming the homogeneity observed in this species, in accord with previous results [46]. Inclusion of C. psittaci 84/2334 to C. abortus increased the diversity among the species, which remained lower than the diversity among C. psittaci without M56 and 84/2334 (Table 1). The degree of selection can be expressed by the dN/dS ratio; a ratio dN/dS higher than 1 indicates a positive selection (altered amino acid substitutions are common), while a rate lower than 1 indicates negative selection (silent substitutions). The C. trachomatis and C. pneumoniae have the highest dN/dS, but in both species synonymous substitutions are in excess (Table 1).

Divergence between Chlamydia species

The divergences between Chlamydia species adjacent in the tree can be expressed as the average number of nucleotide substitutions per site between two species populations (DXY) and the number of fixed differences, i.e. mutations uniform within the species populations. Since most fixed differences are neutral and accumulate at a rate proportional to the mutation rate (molecular clock), the number of fixed differences are indicative for the elapsed time since two populations evolved from a common ancestor. The DXY between Chlamydia species closely positioned in the N-J tree were similar with two exceptions (Table 2). First, the divergence between C. pneumoniae and C. pecorum was twofold higher than that between other neighbours in the tree. Of note, the DXY and the number of fixed differences between C. pneumoniae from human and LPCoLN was small and consistent with the results obtained by comparisons of available C. pneumoniae genome sequences by Myers and colleagues [45]. Second, the DXY between C. psittaci and C. abortus was remarkably low and the number of fixed differences was even lower when compared to that between C. pneumoniae isolated from human and LPCoLN (Table 2). In addition, only C. abortus and C. psittaci shared polymorphisms. However, while the DXY between C. psittaci without 84/2334 and C. abortus with 84/2334 was similar to that between the whole C. psittaci and C. abortus populations, the number of fixed differences was considerably higher, but still much lower than between other neighbouring species in the tree. In addition, shared polymorphism were absent between C. psittaci without 84/2334 and C. abortus with 84/2334. The divergence between 84/2334 and C. abortus was limited and of the same magnitude as that between human C. pneumoniae and LPCoLN, supporting the notion that 84/2334 and C. abortus are one species.

Table 2
Divergence between Chlamydia species.

C. psittaci 84/2334 belongs to C. abortus

The N-J tree clearly showed that C. psittaci M56 branched separately from the main group of C. psittaci strains isolated mainly from birds, while C. psittaci 84/2334 grouped together with C. abortus isolates (Fig. 1 and and2).2). Cluster analyses of the 7 individual housekeeping gene fragments showed that with four gene fragments M56 grouped close to the remaining (except 84/2334) C. psittaci isolates. With gene fragments gatA, hflX and fumC M56 grouped neither with C. psittaci nor with C. abortus. Cluster analyses of the 7 individual housekeeping gene fragments showed that in all instances 84/2334 grouped together with C. abortus (Fig. S1). In addition, the hemN sequence of 84/2334 was identical to that of all but one C. abortus isolates. Together, with the analyses of the divergence between C. abortus and C. psittaci (Table 2) these results suggest that 84/2334 belongs to C. abortus.

Figure 2
Phylogenetic analyses of concatenated sequences of 7 housekeeping gene fragments of C. psittaci and C. abortus strains.

C. psittaci genotypes associate with host species

Using BURST ST24 and ST28 were defined as singletons, although 13 isolates from parrots/parakeets were ST24 (group I) and 8 isolates (7 from ducks and one from human) were ST28 (group II). Two clonal complexes or groups were defined: group III comprises 5 isolates from pigeons and one isolate from human and group IV was formed by the two isolates from turkey (data not shown).

Results of cluster analyses performed with concatenated sequences of MLST loci of C. psittaci excluding 84/2334 and M56 showed association between host species and C. psittaci genotype with the aforementioned four groups being recognized (Fig. 3). One isolate, C. psittaci VS225 from a parakeet, was found to lie between the group I isolates from parrots/parakeets and the group II isolates from duck. The C. psittaci WC bovine isolate did not group with any of the four main groups. Clustering was not associated with the geographic origin of the isolates or of their corresponding host species as parrots/parakeets which have their origin in three different continents, grouped together (Table S1).

Figure 3
Phylogenetic analyses of concatenated sequences of 7 housekeeping gene fragments of C. psittaci strains, excluding strain 84/2334 and M56.


In a previous study we used the present MLST scheme to analyze clonal groupings among C. trachomatis and C. pneumoniae strains [7]. A database hosted at developed and sited at the University of Oxford. The website offers a large number of ways to query the database and to further break down and export the results [47]. In the present study we used the MLST scheme to explore the population genetics of C. psittaci and C. abortus isolates. The main findings are that C. psittaci genotypes are associated with their host species and that isolate 84/2334, formerly classified as C. psittaci, is most likely typed as C. abortus.

We have shown in our previous study that an UPGMA tree produced from the allelic profiles and from concatenated allele sequences of 28 C. trachomatis isolates resulted in three groups of sequence types [7]. The urogenital strains were distributed over two separated groups; one consisted solely of strains with the frequent occurring serovars E, D and F. Strains isolated from patients with lymphogranuloma venereum (LGV strains) grouped in a single cluster, which also included C. trachomatis B/TW5, although not being identical. The close relatedness of B/TW5 was supported by results of comparisons of IncA sequences, showing that B/TW5 shares IncA polymorphisms with LGV strains, which were not found among other serovars [6]. The present study, with C. trachomatis strains supplemented with isolates from a study by Ikryannikova and colleagues [48] to a total of 60, yielded the same three groups in cluster analyses. Recently, another MLST scheme was described showing clonal groupings among C. trachomatis strains with a group consisting exclusively of LGV strains [49]. These results differ somewhat from those reported in our previous study [7].

Recently, a phylogenetic tree was reported that was based on concatenated sequences of 110 conserved protein sequences extracted from the available whole genome sequences of C. trachomatis, C. pneumoniae, C. pecorum, C. felis, C. psittaci, C. abortus and C. caviae [2], [45]. Our tree that is based on the concatenated sequences of 7 housekeeping fragments is very similar to the original tree. In addition, we showed that in our tree C. pneumoniae LPCoLN isolated from koala branched separately from the remaining human C. pneumoniae strains consistent with the results reported by Myers and colleagues who compared whole genome sequences of LPCoLN and C. pneumoniae AR39, TW-183, CWL029 and J138 isolated from human providing evidence that humans were originally infected zoonotically [45]. Together, these results indicate that the 7 housekeeping fragments used in our study are representative for the whole genome sequence.

Phylogenetic analyses of the concatenated allele sequences of C. psittaci revealed an association between C. psittaci genotype and host species. Genotype was not associated with geographic origin, as the host birds had their natural habitat in different countries and covering four different continents, yet the C. psittaci isolated from them were of a single genotype. This identified association between C. psittaci genotype and host species has not been previously observed, which may due to the difference in sequences studied. Cluster analyses of C. psittaci has previously been performed using ompA [50], [51], ompA based PCR-RFLP [52] or multi loci variable number of tandem repeats (MLVA) [53]. ompA encodes the major outer membrane protein and is subjected to host immuno-pressure and its sequence may therefore not reflect the genetic make up of Chlamydia. MLVA based on genetic variation of tandem repeats may be also less suitable to assess evolutionary history by phylogenetic analyses, because it lacks sequence information.

Recently Hackett and colleagues examined ~32 kilobases of aligned nuclear DNA sequences from 19 independent loci for 169 bird species, representing all major extant groups, and recovered a robust phylogeny from a genome-wide signal supported by multiple analytical methods [54]. In their phylogenetic tree based on Maximum Likelhood analyses ducks and turkeys group relatively close together, while pigeons and parrots are more distantly related. In contrast, our results show that C. psittaci from ducks grouped closer to those from parrots and more distantly from those isolated from turkeys, indicating that C. psittaci phylogeny did not follow that of their host species birds. Possibly, C. psittaci may have been transferred from Psittaciformes (order of birds to which parrots and parakeets belong) to Ansiformes (order of birds to which ducks belong) or vice versa. Alternatively, C. psittaci has been transferred to different bird species more than once from a source other than birds. Both scenarios imply host species jump. Incidental transmission is not likely because C. psittaci isolates from these birds and from ducks were isolated in different countries and as far as known in different years. In addition, C. psittaci M56 from a muskrat has it own branch in the phylogenetic tree obtained by using concatenated sequences of MLST loci similar as in a split network graph based on a sequence similarity matrix of ompA sequences reported by Sachse and colleagues [51]. This and the observation, that C. psittaci WC from a cow also does not group with C. psittaci from birds, might indicate that C. psittaci may have spread among mammals, allowing for the development of mammal species specific C. psittaci genotypes as well. Analyses of more C. psittaci isolates from different host species will shed more light on this interesting question. Recently Mitchell and colleagues showed that C. pneumoniae found among animals is more divers than those found among humans [55]. In addition, their results support two separate animal-to-human cross species transfer events in the evolutionary history of this pathogen. Although C. psittaci has a broad host range and numerous cases of transmissions from birds to humans have been described [16], reports of human to human transmission of C. psittaci are rare [21], [56], indicating that C. psittaci has not yet adapted to the human host.

Of note, two C. psittaci strains were isolated from human, suggesting that these were incidental zoonotic infections, without establishing C. psittaci in humans. C. psittaci CPMN isolated from human and thereafter repeatedly passaged in ferrets [57] and C. psittaci humaan E strain grouped together with strains isolated from pigeons and from ducks, respectively. Interestingly, the latter case was a turkey farmer, of whom pharyngeal and nasal swab were taken before the arrival of a new flock of turkeys arrived at the farm and repeatedly thereafter [24]. Initially, C. psittaci with ompA genotype E was cultured, but after three to six weeks ompA genotype A was found, the same ompA genotype found among the turkeys. Our data suggest that the initial infection may be acquired from waterfowl.

The 16 C. abortus isolates comprised 4 STs (5 when 84/2334 is included). Of note, strains LLG and POS have the same genotype and grouped separately from the remaining C. abortus isolates, consistent with the results of Laroucau and colleagues using MLVA [46] and immunological analysis [58]. Furthermore, MLST analyses in this study differentiated the vaccine strain 1B from its parental strain AB7, both classified as the same genotype by MLVA [46]. The 4 C. abortus genotypes were not grouped according to their host species (sheep or goat).

The much lower divergence between C. psittaci and C. abortus compared to that between the other species, suggests C. psittaci and C. abortus have diverged from a common ancestor and much more recently than the others in the tree have from their common ancestor. Alternatively, C. psittaci might have been transmitted from birds to ruminants and adapted to the new hosts as already suggested by Pudjiatmoko and colleagues [59]. In our study, host species jump by C. abortus may be indicated by C. psittaci 84/2334, which grouped close to all C. abortus strains, showing limited divergence with C. abortus with only 7 fixed nucleotide substitutions and should therefore be classified as C. abortus. Originally, C. psittaci 84/2334 was isolated from a yellow-crown amazon (parrot) and classified as C. psittaci based on its reaction with specific sera against the major outer membrane protein (MOMP). However, among 60 C. psittaci isolates from birds, strain 84/2334 had a unique ompA AluI restriction pattern, indicating that this strain differs from the majority of C. psittaci strains [60]. It was suggested that C.psittaci 84/2334 was an intermediate between C. psittaci and C. abortus based on analysis of ompA, rnpB and the rrn (part of the region between the 16S rRNA 23SrRNA genes) sequences [50]. In addition, C.psittaci 84/2334 was found to have DNA sequences that were identical to an extrachromosomal plasmid in duck C. psittaci strain N352, while extrachromosomal plasmids are not found in C. abortus strains [50]. However, trees based on full length 16S rRNA and 23S rRNA gene sequences showed that C. abortus strains grouped together with C. psittaci strains [1]. Similar results were obtained with cluster analyses based on 390 bp rnp sequences. In an N-J tree based on these sequences C. psittaci 6BC and C. abortus B577 were indistinguishable [61]. Cluster analyses of rnp sequences extracted from GenBank showed C. psittaci 84/2334 grouped closer to C. abortus isolates than to C. psittaci isolates (data not shown). Also, the difference in the 222 base pairs rrn sequence between C. psittaci strains and C. abortus strains is limited. C. psittaci 84/2334 differs at only one position with 5 C. abortus isolates of which their rrn sequences are publicly available. C. psittaci 84/2334 and the 5 C. abortus isolates differ at 3 to 4 positions with C. psittaci strains. Again, in a tree based on rrn sequences extracted from GenBank, C. psittaci 84/2334 grouped closer to C. abortus than to C. psittaci (data not shown). In addition, an extended rrn sequence of 315 bp of C. psittaci 84/2334 also differed at only one position from the C. abortus sequences [62]. In our cluster analyses we used up to 10-fold more sequences information than the single 16S RNA, 26S rRNA, rnn and rnp sequences and showed clear separation between C. psittaci isolates and C. abortus isolates with one exception: C. psittaci 84/2334 grouped together with C. abortus. Cluster analyses of the individual housekeeping fragments showed that in all instances strain 84/2334 grouped together with C. abortus. Nevertheless, for six of the seven housekeeping fragments strain 84/2334 has a unique allele. This could indicate that strain 84/2334 is unique within the group of C. abortus isolates. In addition, only C. abortus from a limited set of host (goats and sheep) have been included in the analyses. Inclusion of more C. abortus isolates from different hosts will most likely result in more C. abortus genotypes some of these maybe more related or identical to that of C. psittaci 84/2334. Ultimately, whole genome sequences of more C. psittaci strains isolated from different host including C. psittaci 84/2334 will be ideal to assess these questions.

In conclusion, C. psittaci genotypes are associated with host species and C. psittaci 84/2334 should be reclassified as C. abortus.



All 132 Chlamydia isolates (currently present in the MLST data base for Chlamydiales at were included in the study (Supplementary Table S1). This includes 26 C. trachomatis and 16 C. pneumoniae strains from a earlier study [7] and 30 C. trachomatis strains submitted to the MLST database by Dr. Ikryannikova, who kindly gave her permission to use the data [48]. Detailed analyses were performed with 34 C. psittaci strains isolated from different bird species and mammals from different geographic locations and 16 C. abortus strains. For comparison, sequences of C. caviae, C. felis, C. pecorum and C. muridarum were also included [7]. C. psittaci strains were isolated form dead animals brought to the veterinary clinic for autopsy, or from samples sent to the laboratory by veterinary practitioners for chlamydial diagnosis. C. abortus B577 (VR656) was purchased from ATCC. Greek C. abortus strains were isolated from infected placentae or aborted fetuses. C. abortus MA/231184, MB/312, ME/4004, MF/337, MD/3920, FAS, FAG, VPIG and the variant strains LLG and POS have been previously described, as well as the reference strains AB7, A22, and S26/3 [58], [46]. The vaccine strain 1B (Enzovax) was purchased from Intervet (Intervet-Hellas) and strain Krauss-15 was kindly provided by Dr Jones (Moredun Research Institute, UK). C. psittaci was cultured in Buffalo green monkey cells (ATCC CCL-26) identifying the organisms by the IMAGEN™ Chlamydia immunofluorescence staining (Dakocytomation, Denmark), as previously described [24] and C. abortus stocks were produced in McCoy cell monolayers in the continuous presence of cycloheximide (1 µg/ml) and subsequently typed with monoclonal antibodies against the major outer membrane protein (MOMP) and the polymorphic membrane proteins (Pmps) as previously described [58].

DNA, genes, PCR products and sequences

DNA extraction, PCRs and DNA sequencing were performed as previously described [7]. All alleles of the partial sequences of the 7 housekeeping genes are accessible via the MLST website for Chlamydiales ( New allele sequences have been deposited in GenBank (accession no.: HM776459-HM776511).

Phylogenetic and other analyses

Sequences of fragments from seven housekeeping genes (enoA, fumC, gatA, gidA, hemN, hlfX, oppA) were analyzed as previously described [7]. Allele numbers and genotypes were identified at Clonal complexes were identified using BURST in the START2.0 software package at[63]. Clonal complexes consisted of sequence types that shared 6 of 7 alleles with at least 1 other sequence type in the complex. The index of association for the different Chlamydia species according to Haubold [39] and sequences were tested for recombination using maximum chi square [40] with the START2.0 software package.The number of synonymous and non-synonymous substitutions per site was determined using DnaSP 4.0 [64]. Phylogenetic and molecular evolutionary analyses of the seven housekeeping fragments (individually or concatenated) were conducted using MEGA version 4 [65] or SplitsTree version 4.0 [66], generating a Neighbor-Joining (N–J) tree using the Maximum Composite Likelihood model.

Supporting Information

Table S1

All isolates used in this study according to Chlamydia species, host and ST.

(0.12 MB XLS)

Table S2

Detecting linkage disequilibrium in MLST data of Chlamydia species.

(0.03 MB DOC)

Table S3

Detection of putative recombination events in Chlamydia species.

(0.04 MB DOC)


The authors like to acknowledge stimulating discussion with Dr. A. Bart. Melanie Nguyen is acknowledged for technical assistance. Dr. W.A. Paxton is acknowledged for critical reading of the manuscript.


Competing Interests: The authors have declared that no competing interests exist.

Funding: The development of the Chlamydiales MLST web site has been funded by the Wellcome Trust. Veerle Dickx is funded by the Federal Public Service of Health, Food Chain Safety and Environment (convention RF-6177). 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’ []) and partially by the Belgian Federal Public Service of Health, Food Chain Safety and Environment (convention RF-6177). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.


1. 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. [PubMed]
2. Stephens RS, Myers G, Eppinger M, Bavoil PM. Divergence without difference: phylogenetics and taxonomy of Chlamydia resolved. FEMS Immunol Med Microbiol. 2009;55:115–119. [PubMed]
3. White JA. Manifestations and management of lymphogranuloma venereum. Curr Opin Infect Dis. 2009;22:57–66. [PubMed]
4. Burton MJ, Mabey DC. The global burden of trachoma: a review. PLoS Negl Trop Dis. 2009;3:e460. [PMC free article] [PubMed]
5. Morrison RP. New insights into a persistent problem – chlamydial infections. J Clin Invest. 2003;111:1647–9. [PMC free article] [PubMed]
6. Pannekoek Y, van der Ende A, Eijk PP, van Marle J, dd Witte MA, et al. Normal IncA expression and fusogenicity of inclusions in Chlamydia trachomatis isolates with the incA I47T mutation. Infect Immun. 2001;69:4654–6. [PMC free article] [PubMed]
7. Pannekoek Y, Morelli G, Kusecek B, Morré SA, Ossewaarde JM, et al. Multi locus sequence typing of Chlamydiales: clonal groupings within the obligate intracellular bacteria Chlamydia trachomatis. BMC Microbiol. 2008;8:42–52. [PMC free article] [PubMed]
8. Kuo CC, Jackson LA, Campbell LA, Grayston JT. Chlamydia pneumoniae (TWAR). Clin Microbiol Rev. 1995;8:451–61. [PMC free article] [PubMed]
9. Campbell LA, Kuo CC. Chlamydia pneumoniae and atherosclerosis. Semin Respir Infect. 2003;18:48–54. [PubMed]
10. Grayston JT. Background and current knowledge of Chlamydia pneumoniae and atherosclerosis. J Infect Dis. 2000;181:S402–S410. [PubMed]
11. Kern JM, Maass V, Maass M. Chlamydia pneumoniae-induced pathological signaling in the vasculature. FEMS Immunol Med Microbiol. 2009;55:131–139. [PubMed]
12. Vanrompay D, Ducatelle R, Haesebrouck F. Chlamydia psittaci infections: a review with emphasis on avian chlamydiosis. Vet Microbiol. 1995;45:93–119. [PubMed]
13. Moroney JM, Guevara R, Iverson C, Chen FM, et al. Detection of chlamydiosis in a shipment of pet birds, leading to recognition of an outbreak of clinically mild psittacosis in humans. Clinic Infect Dis. 1998;26:1425–1429. [PubMed]
14. Essig A, Zucs P, Susa M, Wasenauer G, Mamat U, et al. Diagnosis of ornithosis by cell culture and polymerase chain reaction in a patient with chronic pneumonia. Clin Infect Dis. 1995;21:1495–1497. [PubMed]
15. Harkinezhad T, Geens T, Vanrompay D. Chlamydophila psittaci infections in birds: a review with emphasis on zoonotic consequences. Vet Microbiol. 2009;135:68–77. [PubMed]
16. Beeckman D, Vanrompay D. Zoonotic Chlamydophila psittaci infections from a clinical perspective. Clin Microbiol and Infect. 2009;15:11–17. [PubMed]
17. Hinton DG, Shipley A, Galvin JW, Harkin JT, Brunton RA. Chlamydiosis in workers at a duck farm and processing plant. Aust Vet J. 1993;70:174–176. [PubMed]
18. Huminer D, Pitlik S, Kitayin D, Weissman Y, Samra Z. Prevalence of Chlamydia psittaci infection among persons who work with birds. Isr J Med Sci. 1992;28:739–741. [PubMed]
19. Huminer D, Samra Z, Weisman Y, Pitlik S. Family outbreaks of psittacosis in Israel. Lancet. 1988;2:615–618. [PubMed]
20. Kaibu H, Iida K, Ueki S, Ehara H, Shimasaki Y, et al. Psittacosis in all four members of a family in Nagasaki, Japan. Jpn J Infect Dis. 2006;59:349–350. [PubMed]
21. Saito T, Ohnishi J, Mori Y, Iinuma Y, Ichiyama S, et al. Infection by Chlamydophilia avium in an elderly couple working in a pet shop. J Clin Microbiol. 2005;43:3011–3013. [PMC free article] [PubMed]
22. Vanrompay D, Harkinezhad T, van de Walle M, Beeckman D, Van Droogenbroeck C, et al. Chlamydophila psittaci transmission from pet birds to humans. Emerg Infect Dis. 2007;13:1108–1110. [PMC free article] [PubMed]
23. Harkinezhad T, Verminnen K, Van Droogenbroeck C, Vanrompay D. Chlamydophila psittaci genotype E/B transmission from African grey parrots to humans. J Med Microbiol. 2007;56:1097–1100. [PubMed]
24. Verminnen K, Duquenne B, De Keukeleire D, Duim B, Pannekoek Y, et al. Evaluation of a Chlamydophila psittaci diagnostic platform for zoonotic risk assessment. J Clin Microbiol. 2008;46:281–285. [PMC free article] [PubMed]
25. Laroucau K, de Barbeyrac B, Vorimore F, Clerc M, Bertin C, et al. Chlamydial infections in duck farms associated with human cases of psittacosis in France. Vet Microbiol. 2009;135:82–89. [PubMed]
26. Van Droogenbroeck C, Beeckman DS, Verminnen K, Marien M, Nauwynck H, et al. Simultaneous zoonotic transmission of Chlamydophila psittaci genotypes D, F and E/B to a veterinary scientist. Vet Microbiol. 2009;135:78–81. [PubMed]
27. Heddema ER, van Hannen EJ, Duim B, de Jongh BM, Kaan JA, et al. An outbreak of psittacosis due to Chlamydophila psittaci genotype A in a veterinary teaching hospital. J Med Microbiol. 2006;55:1571–1575. [PubMed]
28. Pantchev A, Sting R, Bauerfeind R, Tyczka J, Sachse K. Detection of all Chlamydophila and Chlamydia spp. of veterinary interest using species-specific real-time PCR assays. Comp Immunol Microbiol Infect Dis. 2009 doi: 10.1016/j.cimid.2009.08.002. [PubMed]
29. Denamur E, Sayada C, Souriau A, Orfila J, Rodolakis A, et al. Restriction pattern of the major outer-membrane protein gene provides evidence for a homogeneous invasive group among ruminant isolates of Chlamydia psittaci. J Genet Microbiol. 1991;137:2525–2530. [PubMed]
30. Souriau A, Le Rouzic E, Bernard F, Rodolakis A. Differentiation of abortion-inducing and intestinal strains of Chlamydia psittaci isolated from ruminants by the microimmunofluorescence test. Vet Rec. 1993;132:217–219. [PubMed]
31. Salinas J, Souriau A, Cuello F, Rodolakis A. Antigenic diversity of ruminant Chlamydia psittaci strains demonstrated by the indirect microimmunofluorescence test with monoclonal antibodies. Vet Microbiol. 1995;43:219–226. [PubMed]
32. Salti-Montesanto V, Tsoli E, Papavassiliou P, Psarrou E, Markey BK, et al. Diagnosis of ovine enzootic abortion, using a competitive ELISA based on monoclonal antibodies againstvariable segments 1 and 2 of the major outer membrane protein of Chlamydia psittaci serotype 1. Am J Vet Res. 1997;58:228–235. [PubMed]
33. Buxton D. Potential danger to pregnant women of Chlamydia psittaci from sheep. Vet Rec. 1986;118:510–511. [PubMed]
34. Pospischil A, Thoma R, Hilbe M, Grest P, Zimmermann D, et al. Abortion in humans caused by Chlamydophila abortus (Chlamydia psittaci serovar 1). Schweiz Arch Tierheilkd. 2002;144:463–466. [PubMed]
35. Meijer A, Brandenburg A, de Vries J, Beentjes J, Roholl P, et al. Chlamydophila abortus infection in a pregnant woman associated with indirect contact with infected goats. . Eur J Clin Microbiol Infect Dis. 2004;23:487–90. [PubMed]
36. Walder G, Hotzel H, Brezinka C, Gritsch W, Tauber R, et al. An unusual cause of sepsis during pregnancy: recognizing infection with Chlamydophila abortus. Obstet Gynecol. 2005;106:1215–1217. [PubMed]
37. Maiden MC, Bygraves JA, Feil E, Morelli G, Russell JE, et al. 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. [PubMed]
38. Vos M, Didelot X. A comparison of homologous recombination rates in bacteria and archaea. ISME J. 2009;3:199–208. [PubMed]
39. Haubold B, Hudson RR. LIAN 3.0: detecting linkage disequilibrium in multilocus data. Linkage Analysis. Bioinformatics. 2000;16:847–8. [PubMed]
40. Maynard-Smith J. Analyzing the mosaic structure of genes. J Mol Evol. 1992;34:126–9. [PubMed]
41. Priest FG, Barker M, Baillie LW, Holmes EC, Maiden MC. Population structure and evolution of the Bacillus cereus group. J Bacteriol. 2004;186:7959–70. [PMC free article] [PubMed]
42. Gevers D, Cohan FM, Lawrence JG, Spratt BG, Coenye T, et al. Opinion: Re-evaluating prokaryotic species. Nat Rev Microbiol. 2005;3:733–739. [PubMed]
43. Hanage WP, Spratt BG, Turner KM, Fraser C. Modelling bacterial speciation. Philos Trans R Soc Lond B Biol Sci. 2006;361:2039–2044. [PMC free article] [PubMed]
44. Bishop CJ, Aanensen DM, Jordan GE, Kilian M, Hanage WP, et al. Assigning strains to bacterial species via the internet. BMC Biol. 2009;7:3. [PMC free article] [PubMed]
45. Myers GS, Mathews SA, Eppinger M, Mitchell C, O'Brien KK, et al. Evidence that human Chlamydia pneumoniae was zoonotically acquired. J Bacteriol. 2009;191:7225–33. [PMC free article] [PubMed]
46. Laroucau K, Vorimore F, Bertin C, Mohamad KY, Thierry S, et al. Genotyping of Chlamydophila abortus strains by multilocus VNTR analysis. Vet Microbiol. 2009;137:335–344. [PubMed]
47. Jolley KA, Chan MS, Maiden MC. mlstdbNet - distributed multi-locus sequence typing (MLST) databases. BMC Bioinformatics. 2004;5:86. [PMC free article] [PubMed]
48. Ikryannikova LN, Shkarupeta MM, Shitikov EA, Il'ina EN, Govorun VM. Comparative evaluation of new typing schemes for urogenital Chlamydia trachomatis isolates. FEMS Immunol Med Microbiol. 2010;59:188–196. [PubMed]
49. Dean D, Bruno WJ, Wan R, Gomes JP, Devignot S, et al. Predicting phenotype and emerging strains among Chlamydia trachomatis infections. Emerg Infect Dis. 2009;15:1385–94. [PMC free article] [PubMed]
50. Van Loock M, Vanrompay D, Herrmann B, Vander Stappen J, Volckaert G, et al. Missing links in the divergence of Chlamydophila abortus from Chlamydophila psittaci. Int J Syst Evol Microbiol. 2003;53:761–770. [PubMed]
51. Sachse K, Laroucau K, Hotzel H, Schubert E, Ehricht R, et al. Genotyping of Chlamydophila psittaci using a new DNA microarray assay based on sequence analysis of ompA genes. BMC Microbiol. 2008;17:8:63. [PMC free article] [PubMed]
52. Sayada C, Andersen AA, Storey C, Milon A, Eb F, et al. Usefulness of omp1 restriction mapping for avian Chlamydia psittaci isolate differentiation. Res Microbiol. 1995;146:155–165. [PubMed]
53. Laroucau K, Thierry S, Vorimore F, Blanco K, Kaleta E, et al. High resolution typing of Chlamydophila psittaci by multilocus VNTR analysis (MLVA). Infect Genet Evol. 2008;8:171–181. [PubMed]
54. Hackett SJ, Kimball RT, Reddy S, Bowie RC, Braun EL, et al. A phylogenomic study of birds reveals their evolutionary history. Science. 2008;320:1763–1768. [PubMed]
55. Mitchell CM, Hutton S, Myers GSA, Brunham R, Timms P. Chlamydia pneumoniae is genetically diverse in animals and appears to have crossed the host Barrier to Humans on (At Least) Two Occasions. PLoS Pathog. 2010;6:e1000903. doi: 10.1371/journal.ppat.1000903. [PMC free article] [PubMed]
56. Hughes C, Maharg P, Rosario P, Herrell M, Bratt D, et al. Possible nosocomial transmission of psittacosis. Infect Control Hosp Epidemiol. 1997;18:165–168. [PubMed]
57. Francis T, Magill TP. An unidentified virus producing acute meningitis and pneumonitis in experimental animals. J Exp Med. 1938;68:147–160. [PMC free article] [PubMed]
58. Vretou E, Loutrari H, Mariani L, Costelidou K, Eliades P, et al. Diversity among abortive strains of Chlamydia psittaci demonstrated by inclusion morphology, polypeptide profiles and monoclonal antibodies. Vet Microbiol. 1996;51:275–289. [PubMed]
59. Pudjiatmoko, Fukushi H, Ochiai Y, Yamaguchi T, Hirai K. Phylogenetic analysis of the genus Chlamydia based on 16S rRNA gene sequences. Int J Syst Bacteriol. 1997;47:425–31. [PubMed]
60. Vanrompay D, Butaye P, Sayada C, Ducatelle R, Haesebrouck F. Characterization of avian Chlamydia psittaci strains using omp1 restriction mapping and serovar-specific monoclonal antibodies. Res Microbiol. 1997;148:327–33. [PubMed]
61. 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;2000;50(Pt 1):149–58. [PubMed]
62. Yatsentyuk SP, Obukhov IL. Molecular Genetic Characterization of Avian Chlamydophila psittaci Isolates. Rus J Gen. 2007;43:1215–1220.
63. Jolley KA, Feil EJ, Chan MS, Maiden MC. Sequence type analysis and recombinational tests (START). Bioinformatics. 2001;17:1230–1. [PubMed]
64. Rozas J, Sanches-DelBarrio JC, Messeguer X, Rozas R. DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003;19:2496–2497. [PubMed]
65. Tamura K, Dudley J, Nei M, Kumar S. MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007;24:1596–1599. [PubMed]
66. Huson DH, Bryant D. Application of phylogenetic networks in evolutionary studies. Mol Biol Evol. 2006;23:254–267. [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science