Previous Article | Next Article ![]()
Journal of Clinical Microbiology, May 2007, p. 1410-1414, Vol. 45, No. 5
0095-1137/07/$08.00+0 doi:10.1128/JCM.02301-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Department of Clinical Microbiology, Uppsala University Hospital, Uppsala, Sweden,1 Department of Molecular Evolution, EBC, Uppsala University, Uppsala, Sweden,2 Department of Clinical Microbiology, Malmö University Hospital, Malmö, Sweden3
Received 13 November 2006/ Returned for modification 4 January 2007/ Accepted 15 February 2007
|
|
|---|
|
|
|---|
Serotyping of Chlamydia is laborious since it requires multiple passages in cell culture and the use of a large panel of monoclonal antibodies. Genotyping methods are more commonly used, and the major outer membrane protein gene, ompA, provides the best discriminatory capacity of the genes tested. Restriction fragment length polymorphism analysis of ompA is rapid, and its results show a high level of agreement with the results serotyping (9), while DNA sequencing of ompA has a higher resolution and can discriminate strains in clinically high-risk populations (2, 14). However, when it is applied to nonselected populations, the limited resolution of ompA sequencing restricts the amount of epidemiological information that can be obtained (6). This is especially true, given that the single serotype E comprises almost half of all urogenital chlamydial infections, and within this serotype, one genotypic variant appears to predominate (3, 4, 6). There is therefore an obvious need to develop better methods for evaluation of the molecular epidemiology of chlamydial infections. Furthermore, it has recently been shown that mutant strains evade systems commonly used for the detection of C. trachomatis (11). At present little is known about the spread of such changed chlamydial strains, but they are known to be prevalent in several regions of Sweden and are probably prevalent elsewhere. In this context, multilocus sequence typing (MLST) is an important tool for the investigation of whether several clones occur simultaneously and whether they have changed over time.
The aim of this study was to develop a high-resolution method for the discrimination of C. trachomatis strains. We explored the differences between the available genomes and tested a set of five candidate target regions for the design of an MLST system. Here we present data showing that the system developed has a high capacity to discriminate between strains of individual major outer membrane protein serovars; i.e., it is capable of identifying high intraserotype variation. The MLST system also demonstrated sequence variation when it was applied to clinical chlamydial specimens of common serotypes. This type of information can be used to gain epidemiological knowledge of C. trachomatis infections, and we present an example of how MLST can be applied in contact tracing. Furthermore, it may provide a better understanding of cell tropism and pathogenesis.
|
|
|---|
All genes from the raw genome sequences were extracted using Glimmer2 software (12). Homologous regions between multiple strains were identified and scanned for maximum nucleotide sequence difference with the help of a database-driven system called GENCOMP (Hans-Henrik Fuxelius and Siv G. E. Andersson, unpublished software). Homologous genes and intergenic regions were aligned with DIALIGN2 software (8) at the nucleotide level. dnapars (DNA parsimony) and dnadist (DNA distance) software from the PHYLIP package (J. Felsenstein, PHYLIP, phylogeny inference package, University of Washington, Seattle, 1993 [distributed by the author]) were applied to the aligned sequences and used to estimate nucleotide sequence divergence levels. Reference to dnadist in the following means that pairwise dnadist values exceeding a cutoff of 0.005 were counted for each cluster to give a scalar value for each cluster showing the highest pairwise variability. Candidate genes and intergenic regions were ranked by their scalar dnapars and dnadist values. Following visual inspection of the alignment, the sequences associated with maximal distance values were selected as candidates for the development of a genotyping system. Selected regions comprised genes coding for known and hypothetical proteins, along with intergenic regions. The regions were named after the gene (hypothetical or defined) dominating them: CT046 (hctB), CT058, CT144, CT172, and CT682 (pbpB).
Strain selection. Reference strains comprising the 16 major serotypes, serotypes A to L3, of C. trachomatis were selected: A/Sa1, B/Iu1226, Ba/AP/2, C/TW3, D/IC-Cal8, D/UW-3/Cx, E/DK-20/ON, F/IC-Cal3, G/UW57/Cx, H/Wash, Ia/Iu-4168, J/UW36, K/UW31/Cx, L1/440, L2/434, and L3/404. A set of 47 isolates of serotypes D, E, G, and K from clinical specimens was collected over 6 months at Malmö University Hospital; and a set of 21 samples of serotypes D and G was collected from February 2004 to March 2005 at Venhälsan Gay Clinic, Karolinska University Hospital, Stockholm, Sweden (5). The samples of chlamydia used to demonstrate the application of MLST in contact tracing were from our earlier study (6).
DNA purification. DNA from Chlamydia cultures and swab samples was purified by using a QIAamp DNA mini kit (QIAGEN, Hilden, Germany). DNA from urine samples was isolated by using a MagAttract DNA Mini M48 kit on a BioRobot M48 workstation (QIAGEN).
PCR amplification. PCR was used to amplify the five different target regions of the genome of C. trachomatis. The primer pairs used amplified regions between 382 and 2,377 bp (Table 1). The reaction mixture contained 0.4 µM each primer, 0.2 mM deoxynucleoside triphosphates, 1.5 mM MgCl2, and 1.3 U Expand high-fidelity polymerase (Roche Applied Science, Mannheim, Germany). The same temperature program was used for most regions: initial denaturation for 2 min at 94°C, followed by 40 cycles of denaturation for 15 s at 94°C, annealing for 30 s at 60°C, and elongation for 1 min at 72°C. The elongation step was increased by 5 s per cycle after 10 cycles. The amplification was terminated with elongation for 7 min at 72°C. Different cycling conditions were required for pbpB, as this region is longer. The annealing temperature was initially 68°C, but it was decreased by 1°C per cycle during the first 10 cycles to 58°C. Elongation was for 80 s for the first 10 cycles, but it was then increased by 5 s per cycle for the remaining 30 cycles. The ompA gene was amplified as described by Lysén et al. (6).
|
View this table: [in a new window] |
TABLE 1. Primers used for PCR
|
|
|
|---|
|
View this table: [in a new window] |
TABLE 2. Genetic variation between and within serotypes
|
![]() View larger version (9K): [in a new window] |
FIG. 1. Schematic view of the five MLST regions. Each region is named after the dominant gene. Lengths are approximate since there are minor deletions/insertions. The hctB region has a repetitive element (gray).
|
![]() View larger version (10K): [in a new window] |
FIG. 2. Variable positions found in the CT058 target region when 16 reference strains of serotypes A to L3 and 47 clinical specimens of serotypes D, E, G, and K were analyzed. Some strains carry multiple variations.
|
|
View this table: [in a new window] |
TABLE 3. Discriminatory capacity of the MLST system
|
![]() View larger version (11K): [in a new window] |
FIG. 3. A contact-tracing chain of males (squares) and females (rounded), with sampling dates. Negative cases are indicated with white boxes. All positive cases were infected by serotype E and had identical ompA sequences. MLST was able to separate them, and the two variants are presented with different allele numbers in two samples. *, no sequence was obtained from the target region.
|
A common ompA variant among MSM is EF017667, which comprised 98% of all serotype G strains in our previous study and which constituted 44% of all cases of chlamydial infection detected. Analysis of 10 such specimens by MLST differentiated them into four genotypes (Table 4). In addition, a single strain of serotype G with only one point mutation in ompA compared with the ompA sequence of variant EF017667 was examined. It was found to be markedly different in the MLST system, with divergent sequence variants in three of the five target regions.
|
View this table: [in a new window] |
TABLE 4. MLST analysis of C. trachomatis specimens of predominant ompA serotypes in MSM
|
|
|
|---|
The C. trachomatis genome contains hypervariable regions called plasticity zones near the origin and terminus of replication. These are characterized by insertion-deletion differences and breakpoints in gene synteny (10). The CT172 locus is located within one of the plasticity zones, and CT144 is immediately next to it. This zone is characterized by a rapid rate of sequence evolution. The adjacent genes pbpB and ompA are centered in another area of the genome that is flanked by the conserved orthologs pkn5 (CT674) and CT695 and that also shows a high degree of sequence variability but no gene order structure differences. The hctB gene is located distantly from these in a shorter zone of six genes. The dnadist variability in these zones ranges from 13.2 (CT144 and CT172) to 17.9 (pbpB and ompA). Positive selection for high sequence variability makes these regions suitable for use for the typing of closely related clinical strains but not for the study of long-term evolutionary histories.
Initial analysis of the five regions separately in 16 reference serotype A to L3 strains indicated fewer genetic variants in those regions than in ompA, with the exception of pbpB, which showed the same number of polymorphisms. However, when the five MLST regions of clinical isolates of the predominant serotype, serotype E, were analyzed, four regions were found to have a similar or a higher number of genetic variants compared with the number for ompA. Although no covariation is expected, it is noteworthy that strains of serotype E have low levels of sequence variation (i.e., the number of nucleotide changes leading to distinct genetic variants is small) in other target regions compared with the numbers in strains of different serotypes. The pbpB target provided the highest intraserotype sequence variation and was more discriminatory than ompA for isolates of all four serotypes tested. This variation was based on the arbitrary selection of samples from 1 year in a laboratory that ran 27,000 annual tests and that detected 2,000 cases of chlamydial infection. Consequently, even higher levels of variation could be expected if sampling were to take place at different points in time and at different geographical sites.
To investigate the usefulness of MLST as a tool for molecular epidemiology in a real-world situation, it was used to test samples from a contact-tracing chain. This showed that MLST was superior to ompA analysis and was able to discriminate cases of chlamydial infection and link them to separate contact chains, whereas the ompA sequences were identical.
Analysis of MSM with chlamydial infections showed that cases with an identical and predominant serotype G variant of ompA were differentiated by MLST into four genotypes. In contrast, for infections of the predominant serotype D variant, all cases were also identical by MLST. Interestingly, this is in agreement with earlier findings of a lower relative frequency of shared ompA mutations in serotype G strains than in serotype D strains (7). Serotype G strains thus appear to be more variable than serotype D strains also in other regions of the genome.
Previous studies of different single genes have failed to find useful targets for a genotyping system of high resolution. This is explained by the genetic isolation of Chlamydia and the exceptionally low degree of horizontal DNA transfer during evolution (1, 13), leading to low sequence variation. When the results for the five candidate regions in our study were combined, the resolution of the MLST system was significantly higher than that obtained with ompA genotyping alone. Thus, within serotypes, between 50% and 82% of the strains tested had a unique MLST sequence pattern. Furthermore, ompA analysis only slightly enhanced the discrimination of C. trachomatis strains when it was added to analysis of the five candidate regions. However, in optimizing an MLST system, both sequence variation and sequence length are important. Omission of the longest target region, pbpB, considerably reduced the sequence information, but the resolution for genotyping was almost maintained. Further analysis based on more sequence information may indicate the optimal combination of target regions for use in MLST. The phylogenetics of these genetically variable regions may also provide information that will improve our understanding of regions important for pathogenesis.
In summary, our findings show that MLST can be used for the discrimination of C. trachomatis strains and can be applied to high-resolution molecular epidemiology. This provides a tool for understanding the spread of chlamydial infections in the community.
Published ahead of print on 28 February 2007. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»