Figures
Abstract
A duplex real-time PCR assay was designed for simultaneous detection and genotyping of Mycoplasma pneumoniae (M. pneumoniae). The detection/typing performance of this duplex PCR method, targeting specific genes for M. pneumoniae type 1 (mpn 459) and type 2 (mpna 5864), was compared to that of the previously published MpP1 real-time PCR assay and the genotyping method for the adhesin P1 gene (mpn 141). A total of 1,344 throat swab specimens collected from patients in Beijing, China were tested for M. pneumoniae by bacterial culture, MpP1 real-time PCR assay, and our duplex PCR assay, and positive detection rates of 26.9%, 34.4%, and 33.7%, respectively, were obtained. The duplex PCR method demonstrated high sensitivity and accuracy for detecting and genotyping M. pneumoniae, and significant differences in genotyping ability were observed when compared to the conventional P1 gene-based method. M. pneumoniae type 1 was the predominate genotype from 2008 to 2012 in Beijing, and a shift from type 1 to type 2 began to occur in 2013. To our knowledge, this is the first reported incidence of a type shift phenomenon of M. pneumoniae clinical isolates in China. These genotyping results provide important information for understanding recent changes in epidemiological characteristics of M. pneumoniae in Beijing.
Citation: Zhao F, Liu L, Tao X, He L, Meng F, Zhang J (2015) Culture-Independent Detection and Genotyping of Mycoplasma pneumoniae in Clinical Specimens from Beijing, China. PLoS ONE 10(10): e0141702. https://doi.org/10.1371/journal.pone.0141702
Editor: Qijing Zhang, Iowa State University, UNITED STATES
Received: August 26, 2015; Accepted: October 12, 2015; Published: October 28, 2015
Copyright: © 2015 Zhao et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited
Data Availability: All relevant data are within the paper.
Funding: This work was supported by the National Natural Science Foundation of China (Grant No. 81401637) and the National Key Scientific Instrument and Equipment Development Projects (Grant No. 2012YQ18011709).
Competing interests: The authors have declared that no competing interests exist.
Introduction
Mycoplasma pneumoniae (M. pneumoniae), a small prokaryote devoid of a cell wall, is one of the most common etiological agents of human respiratory tract infections. About 10–40% of all cases of community acquired pneumonia can be attributed to M. pneumoniae, which occurs in epidemic peaks at intervals of 3–7 years [1–3]. Occasionally, severe cases with extrapulmonary involvement can result in hospitalization and death due to neurological disease, such as encephalitis [4, 5]. The high rate of morbidity and the occasional mortality reinforces the need for timely diagnosis. While the culturing of M. pneumoniae is time-consuming (2–8 weeks) and unreliable, serological detections require the availability of paired sera exhibiting IgG antibodies, which is age- and time-dependent [6]. Molecular approaches, especially real-time PCR methods, are currently used for rapid, sensitive, and specific detection of M. pneumoniae in clinical specimens [7–10].
In addition to detection, typing of clinical isolates is practically important for understanding the epidemiology of M. pneumoniae infections and for analysis of endemic outbreaks [11, 12]. In general, M. pneumoniae is regarded as a genetically highly uniform species. Previous investigations have revealed 2 types of clinical isolates (type 1 and type 2) that significantly differ in their sequences of the P1gene, mpn 141 [13]. This P1 gene, which codes for an important 170-kDa adhesin and antigenic factor, is the most frequently used typing target for M. pneumoniae. Typing schemes based on conventional PCR [14], restriction fragment length polymorphism [15], real-time PCR with high-resolution melt-analysis [16], and sequencing [17] have been developed to differentiate the P1 genotypes of M. pneumoniae. About 8% of the M. pneumoniae genome consists of multiple copies of four RepMP elements [18–20], two of which (RepMP 2/3 and RepMP 4) are within mpn 141. Besides the two major M. pneumoniae types, in-depth analysis of the RepMp2/3 and RepMp4 repeat sequences within the P1 gene has revealed one sequence variant of type 1 (i.e., variant 1 or V1) [21] and four sequence variants of type 2 (i.e., V2a, b, c, and d) [17, 22–24]. This variation can be generated by homologous recombination of the RepMP elements located both inside and outside the P1 gene [21, 25]. Preservation of these repetitive sequences during presumed genome minimalization reinforces their importance, as they form pools of sequences for homologous recombinations that yield antigenic variations of M. pneumoniae proteins. It is therefore highly likely that novel variants of the P1 gene will continue to appear with homologous recombination, thereby increasing the ambiguity of the traditional P1 typing approach. Furthermore, some studies have indicated that typing targets other than the P1 gene may exist. Musatovova et al. [20] revealed that the sequence divergence involving the RepMP 1-containing genes (mpn 130, mpn 137, and mpn 138), is strictly two types specific. Spuesens et al. [26] and Catrein et al. [27] also confirmed two highly conserved groups of M. pneumoniae strains based on sequence divergence of the RepMP5-containing gene, mpn 142.
In recent years, attempts have been made to find additional proteomic or genomic markers for typing of M. pneumoniae. Pereyre et al. [28] and our group [29] classified M. pneumoniae strains into two distinct groups, based on polypeptide levels, by matrix-assisted laser desorption ionization–time of flight mass spectrometry (MALDI-TOF MS). Genomic analyses [30, 31] have also demonstrated that M. pneumoniae strains can be classified into two distinct types and further classified into sequence variants of each type. This notion was supported by the analysis of phylogenetic trees generated from all RepMP2/3 and RepMP4 sequences. In each of these trees, clear separation was observed between the type 1 and type 2 strains [25]. In 2015, two reports of the whole-genome analysis of 38 M. pneumoniae were published [32, 33]. Whole-genome sequences were obtained by sequencing 20 clinically isolated strains, including type 1, type 2, V2a, and V2c strains, at our laboratory (unpublished data). We compared our genome sequences with the following M. pneumoniae whole-genome sequences reported on the National Center for Biotechnology Information (NCBI) website: M129 (GenBank: U00089.2), FH (GenBank: CP002077.1), resequenced FH (GenBank: CP010546.1), 309 (GenBank: AP012303.1), M29 (GenBank: CP008895.1), M129-B7 (GenBank: CP003913.2), PO1 (GenBank: CP010551.1), 19294 (GenBank: CP010539.1), 39443 (GenBank: CP010540.1), 85138 (GenBank: CP010545.1), 85084 (GenBank: CP010544.1), 54524 (CP010543.1), 54089 (CP010542.1), 51494 (CP010541.1), MAC (CP010550.1), M2592 (GenBank: CP010549.1), M2192 (GenBank: CP010548.1), M1139 (GenBank: CP010547.1), and PI1428 (GenBank: CP010538.1). We found that the genome sequences exhibited distinct characteristics consistent with two biological types of M. pneumoniae. Taken together, these previous reports confirm the existence of two highly conserved types of M. pneumoniae strains, each with respective variants based on minor sequence variations in the P1 gene.
Nevertheless, broad application of the aforementioned conventional bacterial typing schemes is hampered by the time-consuming and unreliable cultivation of M. pneumoniae from clinical specimens. Thus, these lengthy typing methods have low sensitivity compared to real-time PCR performed directly in clinical specimens without cultivation of M. pneumoniae. To simplify the detection/typing procedure and increase sensitivity, our aims were to develop a culture-independent, duplex real-time PCR method for detecting and typing M. pneumoniae simultaneously in clinical specimens and to evaluate the detection/typing performance of this method by comparison to that of conventional methods.
Materials and Methods
Ethics statement
The study was approved by the Research Ethics Committee of the National Institute for Communicable Disease Control and Prevention (ICDC). No specimens were collected for this study. We used 1344 preserved throat swab specimens that were obtained in our laboratory during previous studies. The consent procedure for obtaining the throat samples was approved by the ethics board committee of Beijing Chao-Yang Hospital, Beijing Children’s Hospital, and Beijing Centers for Disease Control and Prevention, separately. Because throat swab is noninvasive, verbal consent was obtained from all patients or their guardians (for minor patients).
Clinical specimens and strains
From August 2008 to December 2014, a total 1344 throat swab specimens from 388 pediatric (< 14 years of age) and 956 adolescent and adult patients were collected from Beijing Chao-Yang Hospital, Beijing Children’s Hospital, and Beijing Centers for Disease Control and Prevention. All clinical specimens were obtained from respiratory tract infection patients according to clinical symptoms. The following American Type Culture Collection (ATCC) reference and clinically isolated strains were used to assess the sensitivity, specificity, lowest detectable limit (LDL), and typing ability of the duplex real-time PCR method: M. pneumoniae (ATCC 29342, 15531, 39505, 29343, 29085, 15377, 15492, 15293, and 49894; ATCC, Manassas, VA), Mycoplasma salivarium (ATCC 23064), Mycoplasma orale (ATCC 23714), Mycoplasma faucium (ATCC 25293), Mycoplasma genitalium (ATCC 33530), Ureaplasma urealyticum (ATCC 27618), Mycoplasma fermentans (ATCC 19989), Mycoplasma hominis (ATCC 23114), Mycoplasma penetrans (ATCC 55252), Mycoplasma hyorhinis (ATCC 17981), Mycoplasma pirum (ATCC 25960), Escherichia coli (ATCC 11229), Streptococcus pneumoniae (ATCC 49619), Staphylococcus aureus (ATCC 29213), Staphylococcus epidermidis (ATCC 35984), Pseudomonas aeruginosa (clinical isolate), Chlamydophila pneumoniae (clinical isolate), Mycobacterium tuberculosis (clinical isolate), Legionella pneumophila (clinical isolate), Haemophilus influenzae (clinical isolate), Neisseria meningitidis (clinical isolate), and 56 other M. pneumoniae clinical isolates.
Duplex real-time PCR assay
By comparing the genome sequences of 20 M. pneumoniae clinical isolates with those of the reference strains, we observed remarkable differences in large fragment insertions or deletions, mainly between genes mpn 459 (corresponding to the type 1 strain, M129) and mpna 5864 (corresponding to the type 2 strain, 309). The type 1 strains had an approximately 3-kb insertion encoding three genes, mpn 457–459, whereas the type 2 strains had an approximately 5-kb insertion encoding five genes, mpna 5861–5865. Therefore, the primers and probes, designed with Primer Express software (version 3.0; Life Technologies-Applied Biosystems, Grand Island, NY, USA), were based on the conserved sequences of the 3kb and 5kb insertion fragments within the type 1 and type 2 strain, respectively. These primers and probes are listed in Table 1. Each PCR mixture was prepared in a total volume of 25 μl and contained the following per reaction: 12.5 μl Platinum Quantitative PCR SuperMix-UDG (Life Technologies-Invitrogen), 1.5 μl MgCl2 (50 mM), 0.5 μM final concentration of each primer, 0.2 μM final concentration of each probe, 1.25 U Platinum Taq DNA polymerase (5 U/μl; Life Technologies-Invitrogen), 1 μl PCR nucleotide mix (10 mM), 5 μl nucleic acid extracted from each specimen, and nuclease-free water to achieve a 25 μl final volume. Real-time PCR for each target was performed in the CFX96 Real-time PCR Detection System (Bio-Rad, Hercules, CA, USA) under the following conditions: predenaturation at 95°C for 2 min, followed by 45 cycles at 95°C for 15 s and 56°C for 15 s. The data were analyzed with the CFX Manager Software (version 2.1; Bio-Rad).
LDL, sensitivity, and specificity of duplex real-time PCR assay
DNA was extracted from the reference strains and clinical isolates mentioned above with the QIAamp DNA Blood Mini Kit (QIAGEN, Hilden, Germany). The standard curves and LDL for each M. pneumoniae type were determined by analyzing serial 10-fold dilutions of Tris-buffered DNA extracted from M. pneumoniae type 1 (M129 strain [ATCC29342], 0.162 fg/μl to 0.162 ng/μl) and M. pneumoniae type 2 (FH strain [ATCC 15531] 0.183 fg/μl to 0.183 ng/μl). M. pneumoniae type 1 DNA was detected with FAM fluorescence, and type 2 DNA was detected with VIC fluorescence. Each dilution of DNA was assayed in triplicate.
DNA mixtures for duplex real-time PCR detection analysis
The detection ability was determined by analyzing the following four mixtures comprised of different compositions of ATCC 29342 and ATCC 15531 DNA: mixture A (1.62 pg/μl ATCC 29342 and 1.83 fg/μl ATCC 15531), mixture B (1.62 fg/μl ATCC 29342 and 1.83 pg/μl ATCC 15531), mixture C (1.62 pg/μl ATCC 29342 and 1.83 pg/μl ATCC 15531), and mixture D (1.62 fg/μl ATCC 29342 and 1.83 fg/μl ATCC 15531). Each DNA mixture was assayed in triplicate.
M. pneumoniae strains for comparative typing analysis by culture-based and duplex real-time PCR methods
DNA extracted from 9 M. pneumoniae ATCC strains and 56 clinical isolates were amplified with the duplex real-time PCR method for typing analysis. All 65 M. pneumoniae strains were genotyped by full-length sequencing of the P1 gene [23] with the amplification primers SeqP1-F (5’-ATGCACCAAACCAAAAAAACTGCCT-3’) and SeqP1-R (5’-CTAAGCGGGTTTTTTAGGTGGTTGC-3’). The same M. pneumoniae strains were also genotyped by MALDI-TOF MS, based on their peptide mass fingerprints, as described in our previous report [29].
Analysis of M. pneumoniae detection and typing in clinical specimens
A total of 1344 throat swab specimens were cultured with selective Mycoplasma Broth (Thermo Scientific-Oxoid Limited, Basingstoke Hampshire, UK). Each culture-positive specimen was subcultured and purified by the filtration-cloning technique for M. pneumoniae clinical isolates. Genomic DNA of each isolate was then extracted with the QIAamp DNA Blood Mini Kit, and subsequently genotyped with the P1 gene-based PCR method as previously reported [12, 32]. The duplex real-time PCR assay and the previously reported MpP1 real-time PCR assay [34] were used simultaneously to detect M. pneumoniae DNA in the throat swab specimen extracts.
Statistical analysis
Quantitative data are presented as the mean ± standard deviation (SD). The results were calculated in the Statistical Analysis System software package (version 9.3; Cary, NC, USA). Differences in number of positive results between each two of the three detection/typing methods were tested with the paired chi-square test. Differences in M. pneumoniae typing differentiation ability and in number of positive typing results between the P1 gene-based PCR method and the duplex PCR method were tested with the chi-square test and paired chi-square test, respectively, with a P value < 0.01 considered statistically significant.
Results
LDL, sensitivity, and specificity of duplex real-time PCR assay
Based on the crossing threshold (CT) values for the serial 10-fold dilutions of DNA, the LDLs of M. pneumoniae type 1 and type 2 were approximately 8.1 fg and 9.3 fg DNA, respectively (Table 2). Both these amounts are equivalent to approximately 3 colony forming units (CFU). The standard curves (CT vs. log CFU) for M. pneumoniae type 1 and type 2, generated by the serial 10-fold dilutions, yielded coefficient of determination (r2) values of 0.996 and 0.997, respectively. The duplex real-time PCR assay amplified DNA from all M. pneumoniae species tested (9 ATCC strains and 56 clinical isolates). This assay did not amplify DNA from any of the other Mycoplasma species, other respiratory bacteria (10 mollicute species and 10 common respiratory bacteria, listed in Materials and Methods), or human DNA, even when tested at concentrations (at ng level) markedly higher than the LDL.
Duplex real-time PCR detection analysis with DNA mixtures
Mixture C, with higher (pg) concentrations of M. pneumoniae type 1 and type 2 DNA, detected positive for both types, but the CT for each fluorescent (FAM and VIC) signal was delayed about 1–2 cycles compared to that for each signal of amplified M. pneumoniae DNA detected individually at the same concentration. Mixture D, with lower (fg) concentrations of type 1 and type 2 DNA, detected positive for both types only once. In mixtures A and B, with different concentrations of type 1 and type 2 DNA, only the type with the higher concentration could be detected (Table 3).
Comparative M. pneumoniae typing analysis by culture-based and duplex real-time PCR methods
Duplex real-time PCR typing of 9 M. pneumoniae ATCC strains and 56 clinical isolates revealed 44 type 1 strains and 21 type 2 strains. These typing results were completely consistent with those produced by MALDI-TOF MS and P1 gene-based typing. The 21 type 2 strains revealed by the P1 gene-based method comprised of three traditional type 2, two V2a, and 16 V2c strains, based on the minor variations in the P1 gene. No variant 1 strain was detected (Table 4).
M. pneumoniae detection and typing analysis in clinical specimens
Of the 1344 throat swab specimens, the M. pneumoniae positive detection rates of the culture method, MpP1 real-time PCR assay, and our duplex real-time PCR assay were 26.9% (362/1344), 34.4% (462/1344), and 33.7% (453/1344), respectively. The positive detection rates of the MpP1 and duplex real-time PCR assays were not significantly different (paired χ2, P = 0.0389). Differences were observed between those of the culture method and each of the PCR methods (both paired χ2, P < 0.01). A total of 356 specimens tested positive by all three methods, 92 specimens tested positive by the two real-time PCR methods, and only 4 specimens tested positive by the culture method alone (Table 5). Of the 362 M. pneumoniae isolates cultured from the throat swab specimens, 98 were genotyped by full-length sequencing of the P1 gene, and the other 264 isolates were genotyped by the P1 gene-based PCR method [12, 35]. Based on genotyping analysis by the culture-based methods, 85.4% (309/362) of the isolates were traditional type 1 (no variant 1 detected), and 14.6% (53/362) were type 2 (2 isolates were variant 2a, and 51 isolates were variant 2c).
Compared to the aforementioned culture-based typing methods (362/1344), our duplex real-time PCR (453/1344) analysis of the throat swabs identified a greater number of M. pneumoniae genotypes (paired χ2, P < 0.01). Among the 453 clinical specimens that tested positive by duplex real-time PCR, both M. pneumoniae types were identified in 3 specimens (0.7%). MpP1 real-time PCR identified only one of the 3 specimens as positive, and each of these 3 specimens was negative by culture and P1 gene-based typing. Based on genotyping analysis by the duplex real-time PCR assay, 84.4% (385/456 identified genotypes) were type 1, and 15.6% (71/456 identified genotypes) were type 2. A total of 97 culture negative specimens without a P1 gene-based typing result were successfully genotyped with the duplex PCR assay, and only 6 culture-positive specimens genotyped with the P1 gene-based typing method were not detected positive and typed by the duplex real-time PCR assay (Table 5).
From 2008 to 2014, the yearly percentages of M. pneumoniae type 2 genotypes identified by the P1 gene-based typing method were 11.5% (7/61), 20.0% (4/20), 3.1% (2/64), 9.9% (13/131), 15.2% (5/33), 42.3% (11/26), and 40.7% (11/27), respectively. During the same period, the yearly percentages of M. pneumoniae type 2 genotypes identified by the duplex real-time PCR assay were 10.8% (8/74), 12.9% (4/31), 5.1% (4/79), 10.9% (18/165), 19.5% (8/41), 45.5% (15/33), and 42.4% (14/33), respectively. The total percentages of M. pneumoniae type 2 genotypes identified by the P1 gene-based typing and duplex real-time PCR methods during this period were 14.6% (53/362) and 15.6% (71/456), respectively. These genotype percentages were not different from those generated by the culture-based typing methods mentioned above (χ2, P = 0.7128). Between 2008 and 2012, the percentage of identified M. pneumoniae type 2 genotypes in clinical specimens from Beijing did not exceed 20%. However, the percentage of type 2 genotypes rapidly increased to over 40% in2013 and 2014 (Table 6).
Discussion
In this study, we demonstrated that our duplex real-time PCR assay can simultaneously detect and type M. pneumoniae directly from extracted DNA in clinical specimens, without cultivation of the pathogen. The design of our primers/probes was based on the conserved sequences of the type 1-specific 3kb and type 2-specific 5kb insertion fragments. In the NCBI database, the M. pneumoniae FH strain (GenBank: CP002077.1) is the only type 2 strain that does not have the 5-kb insertion fragment. Unexpectedly, the resequenced FH strains (GenBank: CP010546.1) had the 5-kb insertion at position mpna 5861–5865. Furthermore, the ATCC 15531 (FH) strain that we acquired and used in our analysis still had the 5-kb insertion. This finding was also reported in a previous study by Kenri et al. [36]. Whether the ATCC FH strain in our laboratory and the genome-sequenced FH strain (GenBank: CP002077.1) in the NCBI database are in fact different is unclear.
Results from the paired chi-square analyses of the positive detection rates generated by the culture method and the two real-time PCR methods indicate that the detectability of our duplex PCR assay was similar to that of the MpP1 PCR assay and higher than that of the culture method. Only 6 of the 362 culture-positive throat swabs were detected negative by the duplex real-time PCR assay. After further purification of the M. pneumoniae DNA from these specimens, we repeated the assay and found that all 6 specimens tested positive and were correctly genotyped. It is likely that the DNA extracts from these clinical specimens initially contained unknown substances that inhibited PCR amplification. Our assessment of the duplex real-time PCR typing ability revealed that the ATCC strain and clinical isolate typing results from the duplex real-time PCR method were in perfect agreement with those from the culture-based methods. Furthermore, duplex real-time PCR analysis (453/1344) of the clinical specimens identified a greater number of M. pneumoniae genotypes compared to the culture-based typing methods (362/1344), and there were no significant differences in the percentages of each M. pneumoniae type between the duplex PCR and the culture-based methods. Taken together, these findings demonstrate that this novel duplex real-time PCR assay has high accuracy for genotyping as well as detecting M. pneumoniae in clinical specimens.
Coinfection with both types of M. pneumoniae in humans has never been previously reported. In this study, the duplex real-time PCR assay yielded a coinfection rate of 0.7%. However, this rate is likely inaccurate, because analysis of the M. pneumoniae type 1 and type 2 DNA mixtures revealed that the fluorescent signal corresponding to the DNA type of lower concentration could be quenched by that corresponding to the DNA type of higher concentration. Furthermore, rare strains may exist that contain a single P1 gene, as well as both the 3- and 5-kb insertions. Although our findings do not prove their existence, the possibility of such rare strains cannot be excluded. Unfortunately, M. pneumoniae strains could not be isolated from 3 specimens. Therefore, the clinical significance of coinfection with both types and the possible existence of a rare strain of M. pneumoniae remain unclear and require further study.
The annual rates of M. pneumoniae genotypes in clinical specimens from Beijing revealed that type 1 was the predominate genotype from 2008 to 2012. This trend continued for about 5 years, and then a shift toward type 2 began in 2013 and continued through 2014. In previous years, type 1 was always the predominate genotype in China [37–42]. Therefore to our knowledge, this is the first report of M. pneumoniae type shifting in China. Of the 53 type 2 strains detected by the culture-based methods, the majority of them were V2c. Unfortunately, the type 2 variants could not detected by duplex real-time PCR. Previous reports have indicated that a certain strain of M. pneumoniae could confer a type-specific immunity [43]. Although type-specific protection against M. pneumoniae in response to previous infections in the human population has been reported [1, 44], the relationship between the M. pneumoniae type shift phenomenon and pneumonia outbreaks needs to be further verified with more extensive surveillance data. This type shift phenomenon has already been reported in Japan and France [12, 15, 45]. Given that, as reported by Kenri et al. [12], a type shift phenomenon occurs every 8–10 years, and a shift from one type to another requires 2–3 years, we can deduce that a shift from M. pneumoniae type 1 to type 2 in Japan would have occurred around 2012–2013. A recent study by Kubota et al. [46] found that type 1 was still the predominate genotype among M. pneumoniae strains isolated in Japan between 2011 and 2012. Another report by Ishiguro et al. [47] revealed that the percentage of M. pneumoniae type 2 isolates in Japan increased from December 1, 2012 to July 31, 2014. To a certain extent, these results support our deduction. Given that Beijing and Japan are located in the same latitude of the East Asia region where there is frequent population migration, it is possible that an M. pneumoniae type shift phenomenon occurred in these two regions simultaneously. Based on our own findings and inference from the Kenri et al. report, we hypothesize that an M. pneumoniae type shift phenomenon began in Beijing in 2013 and will last for a few years. Currently, M. pneumoniae in our population is still in the process of type shifting. Type 2c (V2c) will most likely become the predominate genotype in our population over the next few years. Although M. pneumoniae type 2 variants have been detected in many regions [38, 47, 48], they comprised a very small percentage of the M. pneumoniae isolates and never became the predominant genotype. In this study, we found that M. pneumoniae in our population started shifting from type 1 to 2c at about the same time that this type 2 variant was becoming increasingly prevalent in Japan [47]. Based on these findings, M. pneumoniae type 2c will likely become the first predominant type 2 variant in the world, perhaps as a result of selection pressure on the human population.
Interestingly, a recent publication failed to find a clear increasing or decreasing trend in P1 type 1, P1 type 2, or their proportion from 2003 to 2012 in Germany [49]. Rather than supporting the hypothesis of type shifts before or during epidemic phases, this result favors the influence of a decreasing level of specific antibodies within the human population as the reason for the start of epidemics. Exactly why these results were different is not clear and requires further study. In any event, we will continue to monitor the status of V2c in Beijing over the next decade. Although the relationship between the type shift phenomenon and pneumonia outbreaks needs to be further verified, increased attention on M. pneumoniae surveillance in Beijing is warranted.
Conclusion
Compared to the conventional, culture-based M. pneumoniae typing methods, our duplex real-time PCR assay has distinct advantages, including time saving and higher typing sensitivity. Although this assay can rapidly differentiate the two M. pneumoniae types, it is unable to differentiate the variant strains. Based on the findings from this study, this duplex PCR assay is most useful as a rapid surveillance tool for M. pneumoniae rather than a refined scientific research method. The sensitivity and specificity of the assay should be validated by testing more known M. pneumoniae isolates in future studies. In conclusion, this novel duplex real-time PCR assay can detect and type M. pneumoniae in clinical specimens with high sensitivity and accuracy. Our genotyping results provide important information for understanding recent type shifts that can change the epidemiological characteristics of M. pneumoniae in Beijing.
Author Contributions
Conceived and designed the experiments: FZ JZ. Performed the experiments: FZ LL XT LH FM. Analyzed the data: FZ JZ. Contributed reagents/materials/analysis tools: FZ. Wrote the paper: FZ JZ.
References
- 1. Waites KB, Talkington DF. Mycoplasma pneumoniae and its role as a human pathogen. Clin Microbiol Rev. 2004. 17: 697–728. pmid:15489344
- 2. Morozumi M, Takahashi T, Ubukata K. Macrolide-resistant Mycoplasma pneumoniae: characteristics of isolates and clinical aspects of community-acquired pneumonia. J Infect Chemother. 2010. 16:78–86. pmid:20094751
- 3. McIntosh K. Community-acquired pneumonia in children. N Engl JMed. 2002. 346:429–437.
- 4. Tsiodras S, Kelesidis I, Kelesidis T, Stamboulis E, Giamarellou H. Central nervous system manifestations of Mycoplasma pneumoniae infections. J Infect. 2005. 51:343–354. pmid:16181677
- 5. Sanchez-Vargas FM, Gomez-Duarte OG. Mycoplasma pneumoniae-an emerging extra-pulmonary pathogen. Clin Microbiol Infect. 2008. 14: 105–117. pmid:17949442
- 6. Daxboeck F, Krause R, Wenisch C. Laboratory diagnosis of Mycoplasma pneumoniae infection. Clin Microbiol Infect. 2003. 9:263–273. pmid:12667235
- 7. Daxboeck F, Khanakah G, Bauer C, Stadler M, Hofmann H, Stanek G. Detection of Mycoplasma pneumoniae in serum specimens from patients with mycoplasma pneumonia by PCR. IntJMedMicrobiol. 2005. 295:279–285.
- 8. Winchell JM, Thurman KA, Mitchell SL, Thacker WL, Fields BS. Evaluation of three real-time PCR assays for detection of Mycoplasma pneumoniae in an outbreak investigation. JClinMicrobiol. 2008. 46:3116–3118.
- 9. Ursi D, Dirven K, Loens K, Ieven M, Goossens H. Detection of Mycoplasma pneumoniae in respiratory samples by real-time PCR using an inhibition control. J Microbiol Methods. 2003. 55:149–153. pmid:14500006
- 10. Dumke R, Schurwanz N, Lenz M, Schuppler M, Lück C, Jacobs E. Sensitive detection of Mycoplasma pneumoniae in human respiratory tract samples by optimized real-time PCR approach. J Clin Microbiol. 2007. 45:2726–2730. pmid:17537933
- 11. Jacobs E, Vonski M, Oberle K, Opitz O, Pietsch K. Are outbreaks and sporadic respiratory infections by Mycoplasma pneumoniae due to two distinct subtypes? Eur J Clin Microbiol InfectDis. 1996. 15:38–44.
- 12. Kenri T, Okazaki N, Yamazaki T, Narita M, Izumikawa K, Matsuoka M, et al. Genotyping analysis of Mycoplasma pneumoniae clinical strains in Japan between 1995 and 2005: type shift phenomenon of M. pneumoniae clinical strains. J Med Microbiol. 2008. 57:469–475. pmid:18349367
- 13. Su CJ, Dallo SF, Chavoya A, Baseman JB. Sequence divergency of thecytadhesin gene of Mycoplasma pneumoniae. Infect Immun. 1990. 58:2669–2674. pmid:1973413
- 14. Kong F, Gordon S, Gilbert GL. Rapid-cycle PCR for detection and typing of Mycoplasma pneumoniae in clinical specimens. JClinMicrobiol. 2000. 38:4256–4259.
- 15. Sasaki T, Kenri T, Okazaki N, Iseki M, Yamashita R, Shintani M, et al. Epidemiological study of Mycoplasma pneumoniae infections in Japan based on PCR-restriction fragment length polymorphism of the P1 cytadhesin gene. J Clin Microbiol. 1996. 34:447–449. pmid:8789036
- 16. Schwartz SB, Mitchell SL, Thurman KA, Wolff BJ, Winchell JM. Identification of P1 variants of Mycoplasma pneumoniae by use of high-resolution melt analysis. J Clin Microbiol. 2009. 47:4117–4120. pmid:19828737
- 17. Dumke R, Lück PC, Noppen C, Schaefer C, von Baum H, Marre R, et al. Culture-independent molecular subtyping of Mycoplasma pneumoniae in clinicalsamples. J Clin Microbiol. 2006. 44:2567–2570. pmid:16825381
- 18. Himmelreich R, Hilbert H, Plagens H, Pirkl E, Li BC, Herrmann R. Complete sequence analysis of the genome of the bacterium Mycoplasma pneumoniae. Nucleic Acids Res. 1996. 24:4420–4449. pmid:8948633
- 19. Dandekar T, Huynen M, Regula JT, Ueberle B, Zimmermann CU, Andrade MA, et al. Re-annotating the Mycoplasma pneumoniae genome sequence: adding value, function and reading frames. Nucleic Acids Res. 2000. 28:3278–3288. pmid:10954595
- 20. Musatovova O, Kannan TR, Baseman JB. Mycoplasma pneumoniae large DNA repetitive elements RepMP1 show type specific organization among strains. 2012. PLoS One. 7:e47625. pmid:23091634
- 21. Dorigo-Zetsma JW, Wilbrink B, Dankert J, Zaat SA. Mycoplasma pneumoniae P1 type 1- and type 2-specific sequences within the P1cytadhesin gene of individual strains. Infect Immun. 2001.69:5612–5618. pmid:11500436
- 22. Kenri T, Taniguchi R, Sasaki Y, Okazaki N, Narita M, Izumikawa K, et al. Identification of a new variable sequence in the P1 cytadhesin gene of Mycoplasma pneumoniae: evidence for the generation of antigenic variation by DNA recombination between repetitive sequences. Infect Immun. 1999. 67:4557–4562. pmid:10456900
- 23. Zhao F, Cao B, Li J, Song S, Tao X, Yin Y, et al. Sequence analysis of the p1 adhesin gene of Mycoplasma pneumoniae in clinical isolates collected in Beijing in 2008 to 2009. 2011. J Clin Microbiol. 49:3000–3003. pmid:21697320
- 24. Xiao J, Liu Y, Wang M, Jiang C, You X, Zhu C. Detection of Mycoplasma pneumoniae P1 subtype variations by denaturing gradient gel electrophoresis. Diagn Microbiol Infect Dis. 2014. 78:24–28. pmid:24268834
- 25. Spuesens EB, Oduber M, Hoogenboezem T, Sluijter M, Hartwig NG, van Rossum AM, Vink C. Sequence variations in RepMP2/3 and RepMP4 elements reveal intragenomic homologous DNA recombination events in Mycoplasma pneumoniae. Microbiology. 2009. 55:2182–2196.
- 26. Spuesens EB, van de Kreeke N, Estevão S, Hoogenboezem T, Sluijter M, Hartwig NG, et al. Variation in a surface-exposed region of the Mycoplasma pneumoniae P40 protein as a consequence of homologous DNA recombination between RepMP5 elements. Microbiology. 2011. 157:473–483. pmid:20966086
- 27. Catrein I, Dumke R, Weiner J, Jacobs E, Herrmann R. Cross-complementation between the products of the genes P1 and ORF6 of Mycoplasma pneumoniae subtypes 1 and 2. Microbiology. 2004. 150:3989–4000. pmid:15583152
- 28. Pereyre S, Tardy F, Renaudin H, Cauvin E, Del Prá Netto Machado L, Tricot A, et al. Identification and subtyping of clinically relevant human and ruminant mycoplasmas by use of matrix-assisted laser desorption ionization-time of flight mass spectrometry. J Clin Microbiol. 2013. 51:3314–3323. pmid:23903545
- 29. Xiao D, Zhao F, Zhang H, Meng F, Zhang J. Novel strategy for typing Mycoplasma pneumoniae isolates by use of matrix-assisted laser desorption ionization-time of flight mass spectrometry coupled with ClinProTools. J Clin Microbiol. 2014. 52:3038–3043. pmid:24920781
- 30. Musatovova O, Kannan TR, Baseman JB. Genomic analysis reveals Mycoplasma pneumoniae repetitive element 1-mediated recombination in a clinical isolate. Infect Immun. 2008. 76:1639–1648. pmid:18212079
- 31. Dumke R, Catrein I, Pirkil E, Herrmann R, Jacobs E. Subtyping of Mycoplasma pneumoniae isolates based on extended genome sequencing and on expression profiles. Int J Med Microbiol. 2003. 292: 513–525. pmid:12635934
- 32. Xiao L, Ptacek T, Osborne JD, Crabb DM, Simmons WL, Lefkowitz EJ, et al. Comparative genome analysis of Mycoplasma pneumoniae. 2015. BMC Genomics. 16: 610–615. pmid:26275904
- 33. Lluch-Senar M, Cozzuto L, Cano J, Delgado J, Llórens-Rico V, Pereyre S, et al. Comparative "-omics" in Mycoplasma pneumoniae Clinical Isolates Reveals Key Virulence Factors. 2015. PLoS One.10:e0137354. pmid:26335586
- 34. Zhao F, Cao B, He LH, Yin YD, Tao XX, Song SF, et al. Evaluation of a new real-time PCR assay for detection of Mycoplasma pneumoniae in clinical specimens. Biomed Environ Sci. 2012. 25:77–81. pmid:22424630
- 35. Yan X, Zhao F, Zhang J. Development and application of a new p1-based genotyping method for Mycoplasma pneumoniae. Wei Sheng Wu Xue Bao. 2012. 52:262–267. pmid:22587007
- 36. Kenri T, Horino A, Matsui M, Sasaki Y, Suzuki S, Narita M, et al. Complete genome sequence of Mycoplasma pneumoniae type 2a strain 309, isolated in Japan. J Bacteriol. 2012. 194:1253–1254. pmid:22328753
- 37. Zhou Z, Li X, Chen X, Luo F, Pan C, Zheng X, et al. Macrolide-resistant Mycoplasma pneumoniae in adults in Zhejiang, China. Antimicrob Agents Chemother. 2015. 59:1048–1051. pmid:25451048
- 38. Xue G, Wang Q, Yan C, Jeoffreys N, Wang L, Li S, et al. Molecular characterizations of PCR-positive Mycoplasma pneumoniae specimens collected from Australia and China. J Clin Microbiol. 2014. 52:1478–1482. pmid:24574282
- 39. Tian XJ, Dong YQ, Dong XP, Li JY, Li D, Jiang Y, et al. P1 gene of Mycoplasma pneumoniae in clinical isolates collected in Beijing in 2010 and relationship between genotyping and macrolide resistance. Chin Med J (Engl). 2013. 126:3944–3948.
- 40. Zhao F, Liu G, Cao B, Wu J, Gu Y, He L, et al. Multiple-locus variable-number tandem-repeat analysis of 201 Mycoplasma pneumoniae isolates from Beijing, China, from 2008 to 2011. J Clin Microbiol. 2013. 51:636–639. pmid:23224090
- 41. Zhang L, Chen ZM, Shen Z, Li SX, Xu D, Tang LF. Genotyping of Mycoplasma pneumoniae in nasopharyngeal aspirate from children with pneumonia in Zhejiang Province. Zhonghua Er Ke Za Zhi. 2011. 49:750–754. pmid:22321180
- 42. Qu J, Yu X, Liu Y, Yin Y, Gu L, Cao B, et al. Specific multilocus variable-number tandem-repeat analysis genotypes of Mycoplasma pneumoniae are associated with diseases severity and macrolide susceptibility. PLoS One. 2013. 8:e82174. pmid:24367502
- 43. Dumke R, Catrein I, Herrmann R, Jacobs E. Preference, adaptation and survival of Mycoplasma pneumoniae subtypes in an animal model. Int J Med Microbiol. 2004. 294:149–155. pmid:15493825
- 44. Razin S, Jacobs E. Mycoplasma adhesion. J Gen Microbiol. 1992. 138:407–422. pmid:1593256
- 45. Pereyre S, Charron A, Renaudin H, Bébéar C, Bébéar CM. First report of macrolide-resistant strains and description of a novel nucleotide sequence variation in the P1 adhesin gene in Mycoplasma pneumoniae clinical strains isolated in France over 12 years. J Clin Microbiol. 2007. 45:3534–3539. pmid:17881549
- 46. Kubota H, Okuno R, Hatakeyama K, Uchitani Y, Sadamasu K, Hirai A, et al. Molecular typing of Mycoplasma pneumoniae isolated from pediatric patients in Tokyo, Japan. Jpn J Infect Dis. 2015. 68:76–78. pmid:25702792
- 47. Ishiguro N, Koseki N, Kaiho M, Kikuta H, Togashi T, Oba K, et al. Regional differences in rates of macrolide-resistant Mycoplasma pneumoniae in Hokkaido, Japan. Jpn J Infect Dis. 2015. [Epub ahead of print]
- 48. Spuesens EB, Hoogenboezem T, Sluijter M, Hartwig NG, van Rossum AM, Vink C. Macrolide resistance determination and molecular typing of Mycoplasma pneumoniae by pyrosequencing. J Microbiol Methods. 2010. 82:214–222. pmid:20547188
- 49. Jacobs E, Ehrhardt I, Dumke R. New insights in the outbreak pattern of Mycoplasma pneumoniae. 2015. Int J Med Microbiol. S1438-4221(15)00084-3.