Previous Article | Next Article ![]()
Journal of Clinical Microbiology, April 2009, p. 914-923, Vol. 47, No. 4
0095-1137/09/$08.00+0 doi:10.1128/JCM.01935-08
Copyright © 2009, American Society for Microbiology. All Rights Reserved.

C. Cazanave,
A. Charron,
H. Renaudin,
C. Bébéar, and
C. M. Bébéar*
Laboratoire de Bactériologie EA 3671, Université Victor Segalen Bordeaux 2, and CHU de Bordeaux, 146 rue Léo Saignat, Bordeaux 33076 Cedex, France
Received 7 October 2008/ Returned for modification 18 December 2008/ Accepted 30 January 2009
|
|
|---|
|
|
|---|
Strain subtyping by molecular methods is a powerful tool for surveillance and outbreak investigation. However, molecular typing is hampered by the fact that M. pneumoniae is a genetically homogeneous species and isolates are poorly differentiated by PCR-restriction fragment length polymorphism (RFLP) analysis of the P1 gene, the most common genotyping method. PCR-RFLP analysis of the gene encoding the P1 protein, a major adhesin that induces a strong humoral immune response during infection, enables the separation of isolates into two types, types 1 and 2 (1, 2, 17). More recent studies used repetitive regions, RepMp2/3 and RepMp4, which can be found in the P1 gene, for molecular typing. Besides the subtype classification, those methods allow the identification of variants of subtypes 1 and 2 (3, 5, 8).
In addition to these molecular typing methods that target only one gene, other methods based on the study of the whole genome have been adapted to M. pneumoniae, such as pulsed-field electrophoresis (1) and multilocus sequence typing (MLST) (4). Pulsed-field electrophoresis is able to differentiate strains into two groups (groups 1 and 2), like P1 gene PCR-RFLP analysis can, and can subdivide group 2 into two subgroups, subgroups 2a and 2b (1). Due to the homogeneity of the M. pneumoniae species, very little polymorphism can be found in M. pneumoniae housekeeping genes. Thus, the attempt to use MLST with housekeeping and structural genes was not useful for molecular typing (4). Among the new genotyping methods, multiple-locus variable-number tandem-repeat (VNTR) analysis (MLVA) has successfully been tested with many bacterial species (20). MLVA is a molecular typing method based on the variation in the copy number of tandemly repeated sequences, called VNTRs, found at different loci on the genome. The variation of the copy number of these tandem repeats (TRs) depends on the isolate tested. Recently, MLVA has proven its suitability for the subtyping of isolates of one animal Mycoplasma species, Mycoplasma mycoides subsp. mycoides SC (12). Furthermore, a multilocus genotyping system has recently been developed for the human species Mycoplasma genitalium, based in part on the analysis of short TR sequences (10).
As current genotyping methods have low discriminatory powers, an alternative method with a high discriminatory ability is needed to differentiate M. pneumoniae isolates according to geographical origin, specimen type, isolation date, and patient age or to differentiate a reinfection from a relapse. The present study aims to develop a new epidemiological tool based on whole-genome analysis, an MLVA scheme, by using VNTR loci selected from the sequence of the sole published genome of M. pneumoniae, that of strain M129, and to apply it to a large collection of isolates.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Characteristics of the 263 M. pneumoniae isolates used in this study
|
Identification of TRs. VNTR markers were identified in the sequenced genome of M. pneumoniae M129 (6) by using the Tandem Repeats Finder program and the Microorganisms Tandem Repeats database, available at http://minisatellites.u-psud.fr. Loci were chosen on the basis of matches of more than 80% between the DNA sequences of the repeat units. The different TR loci were designated by the letters Mpn, followed by a number corresponding to the order of detection and the position on the genome.
To screen for variability in the number of TRs, PCR primers targeting the regions flanking TR loci were designed. As a preliminary step, a total of 17 selected VNTRs were tested with specific primers to amplify DNA from a set of 13 diverse M. pneumoniae strains, including the 2 reference strains, 5 type 1 clinical isolates, 5 type 2 clinical isolates, and 1 isolate that could not be classified by P1 gene PCR-RFLP analysis. These strains were collected over a long time frame and from diverse geographical origins. Each locus was amplified individually and was analyzed by conventional agarose gel electrophoresis. To confirm that length polymorphisms were the result of repeat copy number variations, the PCR products were purified with a Wizard PCR preps DNA purification system (Promega) and were sequenced by using the fluorescence-labeled dideoxynucleotide technology, according to the manufacturer's recommendations (Applied Biosystems). By this approach, five VNTR loci were selected and taken forward for full assessment.
PCR amplification for MLVA. The five loci ultimately selected for MLVA were multiplexed in two solutions, solutions M1 and M2, and the amplifications were performed with a Mastercycler epgradient S apparatus (Eppendorf, Hamburg, Germany) in a final volume of 25 µl. The reaction mixtures contained 1x Qiagen PCR buffer containing 1.5 mM MgCl2; 0.2 mM deoxynucleotide triphosphates; 2 mM and 3 mM MgCl2 for solutions M1 and M2, respectively; 1.25 U of HotStar Taq polymerase (Qiagen, Hilden, Germany); 0.5 µM each primer pairs Mpn1-F-Mpn1-R, Mpn14-F-Mpn14-R, and Mpn16-F-Mpn16-R for solution M1 and primer pairs Mpn13-F-Mpn13-R and Mpn15-F-Mpn15-R for solution M2; and 1 µl of template DNA. Forward or reverse primers were fluorescently labeled at the 5' end with 4,4,7,2',4',5',7'-hexachloro-6-carboxyfluorescein (HEX), 6-carboxyfluorescein (FAM; Eurogentec, France), or NED (2'-chloro-5'-fluoro-7',8'-fused phenyl-1,4-dichloro-6-carboxyfluorescein; Applied Biosystems), depending on the locus to be amplified (Table 2). Both multiplex solutions were run with the same cycling conditions of 95°C for 15 min, followed by 25 cycles of 95°C for 1 min, 60°C for 1 min, and 72°C for 1 min and with a final extension of 72°C for 10 min. Prior to GeneScan analysis, the amplicons were diluted 1:50 in distilled water, and 1 µl of the diluted samples was mixed with 10 µl of Hi-Di formamide (Applied Biosystems) and 0.5 µl of GeneScan ROX 500 size standard (Applied Biosystems). After heat denaturation for 5 min at 95°C and rapid cooling on ice, the fragments were separated on an ABI 3130 genetic analyzer (Applied Biosystems) and were electrophoresed at 15,000 V for 20 min at 60°C. The GeneScan data were analyzed with GeneMapper software (version 4.0; Applied Biosystems) to perform sizing and to calculate the number of repeats in the PCR fragments. Each locus was identified according to color. An allele number string, based on the number of repeats at each locus, was assigned to all isolates. The number of repeats was rounded up to an integer value.
|
View this table: [in a new window] |
TABLE 2. Oligonucleotide primers used for MLVA
|
Stability determination. Five M. pneumoniae isolates, including reference strains M129 and FH, were chosen for use in the stability study. Each strain was passaged 10 times in Hayflick modified broth medium supplemented with glucose to determine the stability of each locus before and after the 10 passages. DNA was extracted and subjected to multiplexing PCR.
|
|
|---|
The stability of the six polymorphic VNTR markers was studied for five strains after 10 passages in Hayflick broth medium. Analysis of the five strains resulted in identical MLVA profiles for all markers except one, Mpn6, which harbored a variation in the amplification size between the 1st and the 10th passages of one strain.
Thus, an MLVA scheme that included five markers (Mpn1, Mpn13, Mpn14, Mpn15, and Mpn16) and that allowed unambiguous type assignment by using agarose gels was proposed. The PCR products ranging in size from 219 to 419 bp. The five loci were distributed around the genome from base pair positions 111920 to 736123. Four of the five VNTRs were located in open reading frames, and three of them encoded hypothetical proteins. Marker Mpn1 was located in the hsdS gene, which encodes the S subunit of a type I restriction-modification enzyme. The characteristics of each VNTR marker are shown in Table 3. The sizes of the unit repeats ranged from 12 bp for Mpn1 to 47 bp for Mpn16. Sequencing of the PCR products of different sizes at each of the five loci from each of the 13 strains confirmed the sizes and sequences of the individual VNTR loci.
|
View this table: [in a new window] |
TABLE 3. Characteristics of VNTR markers
|
The five VNTRs were very efficiently amplified from all 265 strains, and the size variations of the amplicons were exact multiples of the repeats. The results are shown in Table 4. Marker Mpn1 was the most discriminatory VNTR, displaying seven different alleles with repeat copy numbers ranging from 1 to 7, depending on the isolate. Markers Mpn13 and Mpn15 showed three different allele sizes and had repeat units that ranged in number from 2 to 4 and 6 to 8, respectively. Markers Mpn14 and Mpn16 gave two different-sized PCR products. Marker Mpn16 was the most homogeneous marker, with all except two isolates harboring 2 repeat units of 47 bp. This was also reflected by the diversity index of each VNTR, estimated from the HGDI, with a value of 0.784 for the most discriminatory marker (Mpn1) and a value of 0.015 for the less discriminatory one (Mpn16) (Table 3). The overall discriminatory index of the MLVA method was 0.915.
|
View this table: [in a new window] |
TABLE 4. Number of repeat units of the VNTR loci in 265 isolates of M. pneumoniae
|
As a complementary analysis, the MLVA profiles were imported into the Bionumerics software package (Applied Maths), and the genetic relationships of the 265 isolates were deduced by construction of an MST (Fig. 1). This population modeling highlighted the diversity among the isolates tested, but each type was closely related to each other, differing by no more than one locus and emphasizing the lack of distinct clusters. There was no obvious link between the MLVA type and the year when the strains were collected, the patient's age, or resistance to macrolides. Indeed, the 12 macrolide-resistant isolates were dispersed into nine MLVA types (Fig. 1). No association could be established between the 26 unique allele strings and the geographical origin, except for the Tunisian and Spanish isolates. The eight Tunisian isolates were obtained between 2006 and 2008 from different geographical areas and belonged to the same MLVA type, type P, the most important type. The two Spanish strains isolated from two different patients in 1995, both of which were from pericardial fluid or tissue, belonged to MLVA type E (Fig. 1). Furthermore, a clustering was observed on the basis of the P1 subtype, although there were a few exceptions (Fig. 1). P1 variant 1a (strain 4817) belonged to MLVA type J, which included only type 1 isolates, while the five 2a variants (strains F59, Ho, Sta, 79809, and 79692) belonged to MLVA types B, M, and V, which included only type 2 isolates.
![]() View larger version (22K): [in a new window] |
FIG. 1. MST of the MLVA profiles of the 265 M. pneumoniae isolates. Each circle denotes a particular MLVA type, as indicated by a letter and with the number of isolates corresponding to this genotype given in parentheses. The size of the circle is proportional to the number of isolates belonging to the indicated MLVA genotype. The distance between neighboring genotypes is expressed as the number of allelic changes and is outlined by short bold lines for one change. The color of the circles indicates the P1 adhesin type: gray for type 1, white for type 2, and black for characterized variants. Asterisks in the MLVA type indicate the presence of macrolide-resistant isolates.
|
|
|
|---|
The stability of the VNTRs was evaluated in vitro and was confirmed before and after 10 repeated passages of five strains in broth medium. For one patient, two distinct respiratory tract specimens (one sputum specimen and one bronchoalveolar lavage fluid specimen) that were positive for M. pneumoniae were available at the same time. Both isolates were genotyped and were found to harbor the same MLVA type, which also suggested the in vivo stability of the VNTRs chosen (data not shown).
The discriminatory power of this MLVA scheme with the five markers was 0.915 for the collection of 265 isolates. The discriminatory index should ideally be 1.00, but in practice, it should at least be on the order of 0.90 to 0.95 for a typing system to be considered more or less ideal (19). Thus, the global HGDI of this MLVA scheme in our study was very satisfactory. This value must be interpreted in comparison with the discriminatory powers of the other available typing methods. In comparison, the HGDI of the P1 PCR-RFLP typing method for the 265 isolates was only 0.511. The allelic diversity of the VNTR loci varied significantly at each locus, with the Mpn1, Mpn13, Mpn14, and Mpn15 loci being more discriminatory than the Mpn16 locus. The association of this set of markers with variable HGDIs makes this MLVA, in theory, the optimal MLVA scheme for epidemiological studies. Of the five VNTRs tested, only one was located in an intergenic region. Interestingly, marker Mpn1 was located in a coding region containing the hsdS gene, which encodes the S subunit of a type I restriction-modification enzyme involved in protection against foreign DNA (15). The recognition sequence of this enzyme depends on the S subunit, which is known to be variable (16) and which could explain the importance of the allelic variation of the Mpn1 marker.
No clustering was emphasized on the basis of geographical location or specimen origin for any of the isolates except the Spanish and Tunisian isolates. It is interesting to note that two clinical strains, isolated in Spain in 1995 from the pericardial site from two different patients, were of the same MLVA type (13), suggesting that these strains may have a particular virulence mechanism. However, 22 other strains belonged to the same MLVA type as the 2 clinical strains isolated in Spain in 1995. The eight isolates from Tunisia that were studied clustered in the main MLVA type, suggesting a very low level of genetic diversity and the diffusion of a homogeneous population of M. pneumoniae in that country. The 12 macrolide-resistant isolates of M. pneumoniae, which comprised 9 Japanese isolates, 2 French isolates, and 1 Danish isolate, were not related and clustered into nine MLVA types, suggesting the absence of a link between this typing method and macrolide resistance. Nevertheless, this observation could also be explained by the absence of a particular emerging macrolide-resistant clone.
Furthermore, our MLVA scheme was used to study the course of infection in a single patient. Indeed, two clinical isolates obtained 2 months apart from the same patient, isolates B2829 and B2892, showed two different MLVA profiles, as well as two different P1 PCR-RFLP subtypes, suggesting more a M. pneumoniae reinfection than a relapse in this case.
The P1 gene PCR-RFLP typing method (1, 17) divided the clinical isolates only into P1 types 1 and 2. One or the other of the two types tended to predominate among the clinical isolates over time. These changes in the P1 adhesin type may play a role in the development of outbreaks. However, in disagreement with this theory, clear shifts between the M. pneumoniae types did not appear to be correlated with M. pneumoniae epidemic cycles (17). In France, since 1998, both types were present in about the same proportion without a clear dominance of one type over the other one (14). Compared to the findings of the P1 gene PCR-RFLP typing method, which is usually performed to subtype M. pneumoniae, MLVA analysis showed a better power to discriminate among the different isolates, as shown by comparison of the HGDIs between the two methods (see above). Furthermore, while P1 PCR-RFLP typing gave only three patterns (type 1, type 2, and not classified), 26 MLVA types were obtained for the 265 isolates tested. Interestingly, there was a correlation between the typing results obtained by MLVA and PCR-RFLP analysis of the P1 gene. As shown by the MST modeling, M. pneumoniae strains belonging to type 1 or 2 were well separated into two parts, with the exception of a few strains. Beside subtypes 1 and 2, one 1a variant (3) and five 2a variants characterized previously (5) belonged to MLVA types which included only type 1 and type 2 isolates, respectively. Furthermore, one isolate, isolate B3896, was previously found to be different and could not be classified into type 1, type 2, or one of the different subtype variants already described (3, 5, 8, 14). This isolate belonged to MLVA type S, in which only type 2 isolates clustered (Fig. 1), suggesting that isolate B3896 was more related to type 2 than to type 1. These data are in agreement with the sequence analysis of the P1 gene of this isolate, which was more homologous to that of type 2 than to that of type 1 (14).
MST modeling illustrates the relationship of subtypes on the basis of the number of VNTR loci that differ between two MLVA genotypes. Analysis of this population modeling revealed that the 26 MLVA types identified did not differ by more than one locus. Furthermore, no cluster could be individualized, highlighting the homogeneous character of the M. pneumoniae species, despite the high HGDI of the typing method developed in this study. As M. pneumoniae is a relatively homogeneous species, few typing methods have previously been able to detect significant differences between strains (2, 5). One attempt to subtype M. pneumoniae by MLST with housekeeping and structural genes was not sufficiently discriminatory to be used for epidemiological purposes (4).
Although we have studied an important set of clinical M. pneumoniae strains that originated from various countries and that were isolated in different years, unfortunately, no isolates were obtained from an epidemic outbreak of M. pneumoniae disease. Further community outbreaks of acute respiratory tract infections due to M. pneumoniae have been described worldwide (9, 22), and it would be very interesting to evaluate our MLVA scheme in such epidemic situations. The level of discrimination of this typing method should be confirmed by comparing outbreak-related strains to a set of control strains that were isolated from a similar time period and geographical area but that are not epidemiologically related.
In conclusion, this is, to our knowledge, the first attempt to perform the molecular typing of M. pneumoniae by the MLVA method with a consequential number of clinical isolates. Other markers could be identified to improve the MLVA scheme proposed in this study and to find a better correlation between strain information and typing by VNTR.
One of the advantages of MLVA typing is that it is PCR based and does not require the growth of bacteria. Therefore, another important issue will be to adapt the method directly for clinical specimens, since the value of M. pneumoniae culture is far from being satisfactory, is much more fastidious, and is less sensitive than detection by PCR. Furthermore, the data obtained by MLVA could be shared in a web-based database for M. pneumoniae, as has already been done for other bacterial species.
Published ahead of print on 9 February 2009. ![]()
S.D. and C.C. contributed equally to the work. ![]()
|
|
|---|
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»