Skip to main content

Comprehensive virome analysis of the viral spectrum in paediatric patients diagnosed with Mycoplasma pneumoniae pneumonia



Among hospitalized children suffering from community-acquired pneumonia, Mycoplasma pneumoniae (MP) is one of the most common pathogens. MP often exists as a co-infection with bacteria or viruses, which can exacerbate the clinical symptoms. We investigated the pathogen spectrum in MP-positive and MP-negative samples from hospitalized children with respiratory tract infections in Beijing, China.


This study included 1038 samples of nasopharyngeal aspirates obtained between April, 2017 and March, 2018 from hospitalized children under 6 years of age with respiratory tract infections. To explore the impact of MP infection on the composition of the pathogen spectrum, 185 nasopharyngeal aspirates (83 MP-positive/102 MP-negative) were randomly selected for next-generation sequencing and comprehensive metagenomics analysis. Real-time PCR was used to detect and verify common respiratory viruses.


Of the 1038 samples, 454 (43.7%) were infected with MP. In children < 6 years of age, the MP infection rate gradually increased with age, with the highest rate of 74.2% in 5–6-year-olds. The results of metagenomics analysis revealed 11 human, animal and plant virus families, and bacteriophages, including common respiratory viruses, enteroviruses and anelloviruses. The virus family with the highest number of reads in both MP-positive and MP-negative samples was the Pneumoviridae, and the number of reads for human respiratory syncytial virus (HRSV) in MP-positive samples was higher than that in MP-negative samples. Among the 83 MP-positive samples, 47 (56.63%) were co-infected with viruses, the most common of which was influenza virus (IFV). The durations of hospitalization and fever were higher in patients with MP co-infection than MP single infection, but the difference was not statistically significant.


The viral family with the highest number of reads in both groups was Pneumoviridae, and the number of reads matched to HRSV in MP-positive samples was much higher than MP-negative samples. Co-infection of MP and IFV infection were the most cases.


Mycoplasma pneumoniae (MP) is a pathogenic microorganism that contains DNA and RNA and lacks a cell wall. It is a common pathogen that causes respiratory tract infections (RTIs). Mycoplasma pneumoniae pneumonia (MPP) is a major cause of community-acquired pneumonia in hospitalized children [1]. It is estimated that MPP is responsible for approximately 3–10% of respiratory infections in children [2], and approximately 5–12% of the hospitalized children with MPP were admitted to the intensive care unit [3,4,5]. The clinical presentation of MPP is usually tracheobronchitis, pneumonia, nonspecific symptoms of upper RTI (such as headache, sore throat, rhinitis, otitis media) and extrapulmonary manifestations (such as encephalitis, urticaria, allergic purpura) [6]. MPP is generally considered a mild and self-limiting disease, but in some cases can lead to sever pneumonia, which can be life-threatening in children [7]. It is reported that MP can damage human airway epithelial cells and cilia, and can affect the function of the mucus-ciliary clearance system and host immune function [8]. Additionally, recent research showed that indirect (immune-mediated) mechanisms are mainly implicated in Mycoplasma pneumoniae-related extra-pulmonary diseases (MpEPDs), such as Henoch-Schonlein purpura, Myositis, Kawasaki disease, Nephritis [9]. A cohort study had found that incident cases of early-onset and late-onset asthma might be related to MPP [10]. MP occurs worldwide and is easily transmitted by direct respiratory droplet contact. MPP can cause infections in people of any age, and is one of the leading causes of community-acquired pneumonia (CAP), especially in children [5, 11].

In addition to MP, respiratory viruses also cause RTIs in children. Human respiratory syncytial virus (HRSV), human adenovirus (HAdV) and influenza virus (IFV) are the most common respiratory viruses infecting infants and children under 5 years of age [12]. In children with lower RTIs, approximately 31–66% of cases are caused by co-infection of respiratory viruses and bacteria [13]. In recent years, it has been reported that MP can co-infect a host along with bacteria and viruses (such as Streptococcus pneumoniae, HRSV) and the co-infection rate can reach 52% [14]. Compared with mono-infection, MP co-infection cause serious clinical symptoms and the clinical manifestations are more diverse [7]. The relationship between MP and viruses is yet to be fully elucidated.

Metagenomic next generation sequencing (NGS) can simultaneously detect multiple viruses, as well as new and highly-differentiated viruses, providing comprehensive detection and quantitative analysis of all microorganisms present in clinical samples [15]. In this study, to analyze the pathogenic spectrum of MP-positive and MP-negative samples in hospitalized children with RTIs, nasopharyngeal aspirate (NPA) samples were collected and sequenced using NGS technology. By analyzing the NGS sequences of the two groups, we investigated the impact of MP infection on the composition of the respiratory virus spectrum, and explore whether MP infection affect the diversity and complexity of virus co-infection.


Participants and sample collection

A total of 1038 NPA samples were collected from patients under 6 years of age diagnosed with a RTIs in the pediatric inpatient unit of Beijing Friendship Hospital between April, 2017 and March, 2018. In this study, we described the epidemiological characteristics of MP infection, compared the pathogen spectrum and clinical manifestation between MP-positive group and MP-negative group. The protocol was approved by the Ethical Committee (IVDC2017-021, approval date: Mar 20, 2017). Informed consents were obtained by the legal guardians. In total, 185 samples of NPA were randomly selected using the RAND function of Excel for NGS, including 83 MP-positive samples and 102 MP-negative samples. MPP was diagnosed according to the following criteria: clinical symptoms of a cough and fever; wet rales on auscultation or infiltrative inflammatory manifestations on a chest radiograph; ≥ four-fold antibody titers increase of paired sera; a positive IgM antibody test; PCR-positive for MP; or isolation of MP from cultures of respiratory specimens [16, 17]. In the present study, Serum MP antibody titers were detected with SERODIA-MYCO II kit (Fujirebio, Japan), and titer ≥ 1:80 was considered positive.

NPAs samples were collected by professional personnel and stored in virus transport solution in an ice box during transport to the laboratory. Samples were then stored at ˗80°℃ for later use.

Sample pretreatment and viral nucleic acid extraction

NPA samples in solution (800 µl) were centrifuged at 3000 rpm for 20 min, then 20,000 rpm for 15 min. The supernatant was filtered through a 0.22 μm filter to remove eukaryotic and bacterial cell-sized particles. The filtered sample was subjected to nuclease treatment at 37°℃ for 2 h to eliminate the free and unprotected nucleic acids, including benzonase nuclease (Novagen, Germany), TURBO™DNase (Thermo Fisher Scientific, USA) and RNaseI (Thermo Fisher Scientific). Viral nucleic acids were extracted using the QIAamp MinElute Virus Spin Kit (Qiagen, Germany) in accordance with the manufacturer’s instructions.

cDNA library construction and high throughput sequencing

First-strand cDNA synthesis was performed using the SuperScript III First-Strand Synthesis System (Invitrogen, USA). Subsequently, double-stranded cDNA (ds-cDNA) was synthesized using 3ʹ–5ʹ exo-Klenow fragment (Biolabs, New England). The obtained ds-cDNA was amplified using the QuantiTect Whole Transcriptome Kit (Qiagen, Germany). The amplification products were sent to Beijing Liuhe Huada Gene Technology Company for NGS. The Illumina HiSeq2000 platform was used for NGS. A total of 8G clean data were obtained for each sample.

Pathogen verification

Real-time PCR was performed on the NPA samples to detect common respiratory virus, including HAdV, HRSV, IFV (A/B/C), human parainfluenza virus (HPIV1–4), human metapneumovirus (HMPV), human polyomavirus (HPyV 3/4), human rhinovirus (HRV), human bocavirus (HBoV), and human coronavirus (HCoV 229E/OC43/NL63/ HKU1). Samples tested positive with a cycle threshold (CT) value < 37 (primers and probes are shown in T Additional file 1: Table S2). DNA viruses were detected with the TaqMan™ Gene Expression Master Mix kit (Thermo Fisher Scientific, USA) and RNA viruses were detected with the AgPath-ID™ One-Step RT-PCR kit (Ambion, USA) in accordance with the corresponding manufacture’s protocols.

Bioinformatics analysis

The obtained sequences were processed and analyzed as follows. Prinseq-lite software was used to filter out low-quality data (QC cut-offs was listed in Additional file 1: Table S1). Then Bowtie2 alignment was performed to remove the sequences of the host, and the remaining reads were assembled to obtain contig sequences using MIRA with default parameters [15]. Blastn and Blastx were used to align the obtained reads at the nucleotide and amino acid levels against the virus reference database (RefSeq) download from Genbank, respectively, to determine the virus species and coverage. The E value was set to 1 × 10− 5 for blastn and blastx to maintain high sensitivity and low false-positive rate.

Statistical analysis

SAS 9.4 software was used for statistical analysis of the data. We compared the basic information of children hospitalized with and without MP. Additionally, the clinical symptoms of MP and virus co-infection versus single MP infection were also assessed. Categorical data, such as age and gender, were tested using Pearson’s chi-square test. However, continuous variables including pneumonia diagnosis cases, white blood cell count, lymphocyte percentage, hospital stay, fever duration and max temperature were analyzed using independent-samples two-tailed t-test. Since the C-reactive protein data is non-normal, it was analyzed using Wilcoxon’s rank sum test. P < 0.05 was considered to indicate a statistically significant difference.


Epidemiological characteristics of MP-positive and MP-negative samples

From April, 2017 to March, 2018, a total of 1,038 NPA samples from hospitalized children, < 6 years of age, with RTIs were collected. The age distribution was 1 day to 6 years old, and the average age was 29 months. Of the total cases, 454 were positive for MP and 584 were negative. The MP infection rate was 43.7%. 953 patients were diagnosed with pneumonia in 1038 cases, and 96.7% (439/454) of MP-positive patients were diagnosed with pneumonia. As the age of children increased, the MP infection rate gradually increased, with the highest rate of 74.2% (46/62) detected among those aged 5–6 years (Fig. 1). Among the 454 MP-positive cases, most of the children were older than 2 years (75.5%, 343/454). Among the 584 MP-negative children, 359 cases (61.5%, 359/584) were under 2 years. The difference in age distribution between the two groups was statistically significant (P < 0.01). In the MP-positive group, males accounted for 48.7% and females accounted for 51.3% (sex ratio: 1:1.05). In the MP-negative group, 58.9% cases were from males, 41.1% were from females (sex ratio: 1.43:1) The difference in sex ratio between the two groups was statistically significant (P < 0.01) (Table 1). The laboratory information including white blood cell (WBC) count, lymphocyte percentage (LY%) and C-reactive protein (CRP) were also collected, among which, LY% and CRP were different between MP-positive and MP-negative groups, as shown in Table 1. Seasonal distribution results showed that the rate of MP infection in September was the highest (62.6%, 57/91), but there was no obvious seasonal epidemic (Fig. 2). Additionally, common complications in all patients were secondary thrombocythemia (562/1038), abnormal electrocardiography (227/1038), electrolyte disturbances (225/1038), lymphadenopathy (100/1038) and urticarial (12/1038). The severe complications were respiratory failure (52/1038) and heart failure (15/1038).

Fig. 1
figure 1

Age distribution of RTI hospitalized children infected with MP.

Table 1 Basic information of 1038 hospitalized children
Fig. 2
figure 2

Seasonal distribution of RTI hospitalized children infected with MP.

NGS data

After sequencing 185 randomly-selected NPA samples, the obtained data were processed to removed irrelevant sequences such as host and bacterial sequences, and then compared with the NCBI virus database. The total number of reads for samples from the MP-positive group ranged from 6538 to 666,756 (median 19,811), with a mean read length of 72–130 bp, and the number of contigs ranged from 230 to 8,745 (median 703). The total number of reads for samples from the MP-negative group ranged from 4281 to 438,405 (median 15,006), with a mean read length of 57–194 bp, and the number of contigs ranged from 183 to 3341 (median 531). These reads matched to diverse viruses including human-related viruses, insect viruses, mammalian viruses, fungal viruses, phages and other viruses.

Analysis of the pathogen spectrum

In this study, we separately analyzed torque teno virus (TTV) and other common human-related viruses. Viral spectrum analysis was performed on the common human-related viruses after removing TTV (Additional file 1: Table S3).

A variety of common virus-related sequences were detected in the respiratory tract samples of the MP-positive group. The total number of reads was 637,248. These virus sequences were from the following 11 viral families (Fig. 3): Pneumoviridae, Picornaviridae, Polyomaviridae, Orthomyxoviridae, Adenoviridae, Paramyxoviridae, Coronaviridae, Parvoviridae, Caliciviridae, Astroviridae and Papillomaviridae. Among them, reads from Pneumoviridae, Picornaviridae and Polyomaviridae were more abundant, with relative abundances of 56.98%, 20.01% and 10.54%, respectively. Furthermore, the virus-related sequence reads were mainly consisted of HRSV (56.98%), HRV (19.78%), HPyV (10.54%), IFV (4.93%), HAdV (3.94%) and HPIV (1.42%). Other viruses were also detected, but the number of related reads was relatively small, such as for norovirus, mamastrovirus, human rubulavirus and human papillomavirus.

Fig. 3
figure 3

Distribution of different viral families in MP positive and negative groups. The reads of 11 viral families were sequenced by NGS. A MP-positive group; B MP-negative group

In the MP-negative group, a variety of common virus-related sequences were also detected. The total number of reads was 807,989. These virus sequences were from the following 11 viral families: Pneumoviridae, Picornaviridae, Polyomaviridae, Parvoviridae, Paramyxoviridae, Adenoviridae, Orthomyxoviridae, Coronaviridae, Caliciviridae, Astroviridae and Papillomaviridae (Fig. 3). Different from the MP-positive group, reads from Picornaviridae, Polyomaviridae, Parvoviridae and Paramyxoviridae were more abundant, with relative abundances of 22.27% (P < 0.001), 20.46% (P < 0.001), 11.95% (P < 0.001) and 10.47% (P < 0.001), respectively. Further analysis of the viruses showed that HRV (22.04%), HRSV (24.07%), HPyV (20.46%), HBoV (11.94%), human respirovirus (10.45%), HAdV (6.27%) and IFV (3.65%) accounted for the majority of the sequencing data from MP-negative samples. Other viruses were also detected, but the number of related reads was relatively small, such as for human rubulavirus, human papillomavirus and foot-and-mouth disease virus.

The composition of the pathogen spectrum was compared between the two groups. The virus families with the highest representation among the two groups were the same, both Pneumoviridae and Picornaviridae. The proportion of Pneumoviridae-related virus reads among the total number of reads in the MP-positive group was 56.98%, which was much higher than that in the MP-negative group (24.10%) (P < 0.001). The proportion of Polyomaviridae, Paramyxoviridae and Parvoviridae virus-related sequences in the MP-negative group was higher than that in the MP-positive group (P < 0.001). Regarding specific viruses, the proportion of HRSV-related reads in the MP-positive samples was 57.0%, which was much higher than that in the MP-negative samples (24.1%) (P < 0.001).

Common respiratory viruses

The positive viral infection samples detected by NGS was confirmed by real-time PCR and a number of respiratory viruses were detected including HRSV, HMPV, HRV, HPyV, IFV, HAdV, HCoV and HBoV (Fig. 4 and Table 2). Among 83 cases in the MP-positive group, 47 (56.62%) were mixed infections, of which 38 (45.78%) were single virus infections, 7 (8.43%) were double virus infections, and 2 (2.41%) were triple virus infections. MP and IFV co-infection had the highest detection rate (24.10%, 20/83), followed by MP and HRSV co-infection (12.05%, 10/83). Among 102 MP-negative cases (Fig. 4), 66 (64.71%) were infected with virus. HRSV was the most commonly detected virus (20.59%, 21/102), followed by IFV (15.69%, 16/102) and HRV (15.69%, 16/102).

Fig. 4
figure 4

Common respiratory virus co-infection of all samples sequenced by NGS and real-time PCR. A 83 MP-positive cases; B 102 MP-negative cases. HAdV (human adenoviru), HRSV (Human respiratory syncytial virus), IFV (influenza virus A/B/C), HPIV (human parainfluenza virus 1–4), HMPV (human metapneumovirus), HPyV (human polyomavirus 3/4), HRV (human rhinovirus), HBoV (human bocavirus), and HCoV (human coronavirus, 229E/OC43/NL63/HKU1). Samples tested positive with a cycle threshold (CT) value < 37

Table 2 MP co-infection with other respiratory viruses detected by real-time PCR

Comparison of the clinical symptoms for MP co-infection with single MP infection

To compare the clinical symptoms of MP virus co-infection with single MP infection, the clinical information for patients from the two groups was compared (Table 3). All NPAs were detected by real-time PCR to verify co-infection with common respiratory virus. The results showed that the mean number of days of hospitalization were 9.64 ± 3.87 and 8.72 ± 2.02 days for MP co-infection and MP single infection patients, respectively, but this difference was not significant (P > 0.05). The fever duration for patients with MP co-infection (9.74 ± 3.61 days) was similar to that of single infection patients (9.44 ± 3.08 days) (P > 0.05). Serum Creatine Kinase MB Isoenzyme (CK-MB) levels and the white blood cell count were approximate between MP co-infection and mono-infection patients, no significant differences were observed (P > 0.05).

Table 3 Comparison of clinical symptoms of hospitalized children infected with MP

Anelloviruses and phage

In this study, a large number of sequence reads (34.7%) matched to the Anelloviridae family, mainly TTV, torque teno mini virus (TTMV), torque teno midi virus (TTMDV) and small anellovirus (SAV).

In addition to the common respiratory viruses and anelloviruses, a large number of phage-related sequences were detected, such as Bacillus phage, Escherichia phage and Enterobacter phage. Since bacteriophage are viruses that infect bacteria, the presence of phage indirectly reflects bacterial infection.


Respiratory infections are considered a major health threat to infants and young children [18]. MP is one of the most commonly detected bacteria responsible for mild to severe lower RTIs among older children [19]. Some respiratory viruses can also cause pneumonia, such as IFV, HRSV, HAdV, HBoV and HRV [12]. A previous study showed that in cases of MPP, co-infection with viruses and bacteria is common and the co-infection rate can reach 52% [14]. Compared with mono-infection, the symptoms associated with co-infection can be relatively more serious [20]. This study collected NPA samples between April, 2017 and March, 2018 from hospitalized children under 6 years of age with RTIs in Beijing, China. NGS and real-time PCR were used to detect the causative agents of the RTIs in the children. We compared the virus spectrum and identified the common infecting respiratory viruses between the MP-positive and MP-negative groups, found that both Pneumoviridae and Picornaviridae are the largest proportion among the two groups. 113 of the 185 NPA samples were co-infected with respiratory viruses. By contrast with MP isolated infection children, the clinical symptoms of children with MP combined with virus infection had no statistically difference.

Accurately and promptly detecting respiratory viruses is a key step prior to initiating antiviral treatment [21]. There are noticeable differences in the detection rates of viral pathogens as a result of the method used for detection. NGS is used for detecting pathogens in samples containing a mixture of different species, without the need for sequence-specific amplification [22]. The advantages of NGS include the unbiased, sensitive and simultaneous detection of multiple viral pathogens. In our study, we employed NGS to obtain comprehensive and reliable data. The results were verified by real-time PCR, and 113/185 samples (61.1%) were positive for common respiratory viruses. These data were consistent with those of a previous study that reported a 63% detection rate among samples using a PCR method [23].

RTIs present a major threat to the health of infants and young children [18], and MP is one of the most common pathogens causing mild to severe lower RTIs in children [19]. It was reported that 70.9% of MP infections occurred in children from 1 to 6 years of age in South Korea [24]. In China, 80% of MP infections occurred in children under 7 years of age [25]. In the current study of children aged under 6 years, the proportion of MP infection increased with age, peaking among children aged 5–6 years. A study in the United States showed that 73.3% of children hospitalized with MP infection were over 5 years of age [26]. This may be because younger children are less exposed to the general population, thereby reducing the risk of MP infection. In addition, the results of the current study showed that MP infection could occur throughout the year, with no obvious epidemic seasonality. However, the MP infection rate in fall was slightly higher than that in the other seasons, which was consistent with previous reports that MP infection is more common in the summer or early fall, but may occur at any time of the year [1].

Previous studies have shown that some pediatric patients with RTIs are infected with multiple respiratory pathogens at the same time. MP co-infection with bacteria and viruses is common. For example, MP and IFV, MP and HPIV, and MP and Legionella pneumophila all have high co-detection rates [2]. In our study, the ratio of co-infection with multiple pathogens in MP-positive cases was 56.63%, which was higher than previous reports. In MP-positive cases, the co-infection rate with IFV was highest among all cases (24.1%) and this was consistent with the findings of Zheng, et al. [27] and Kalenahalli, et al. [28]. It has been reported that viral infection can make patients susceptible to MP infection, and that this interaction may lead to the aggravation of asthma [29]. It may be that MP infection can cause cellular and humoral immunity dysfunction and induce immunosuppression in the body, damaging the patient’s systemic and local defense system, and thereby promoting virus infection [30]. These findings suggest that clinicians should consider the co-infection of virus and MP when making clinical diagnosis, and initiate treatment accordingly.

Compared with single infection with MP, the clinical symptoms in co-infected children are relatively more severe and more diverse [20]. MP and virus co-infection can prolong the clearance time of pathogenic bacteria, disrupt the host’s immune response, and prolong the inflammatory response time [8]. The results of this study showed that compared with children with MP single infection, children with MP combined with viral infection had longer hospital stays and longer fever duration, but these differences were not statistically significant. Chiu et al. compared the clinical manifestations of MP single infection, MP combined with Streptococcus pneumoniae infection, and MP combined with virus infection. They found that the duration of fever in children with MP and S. pneumoniae co-infection was longer than that of children with MP single infection, but similar to our study, there was no significant difference [31]. However, in another study, the duration of fever in children with MP combined with HAdV infection was significantly longer than for MP single infection [8]. These differences between studies may reflect the different pathogens analyzed, and the impact of particular virus on MP infection outcomes and the clinical manifestations of MPP require further investigation.

TTV was detected in almost every sample in this study. TTV is a human virus with a circular negative-strand DNA genome of approximately 3.8 kb in length [32, 33]. Although this virus has a high prevalence among the population, its interaction with the host and the etiological relationship with specific diseases are not yet fully understood. There is a significant correlation between TTV load and airflow limitation in the peripheral airways, as well as the severity of bronchiectasis and decreased lung function [34]. Therefore, TTV may cause influenza-like symptoms. However, whether TTV is the cause of influenza-like symptoms, and the reason why TTV resides in certain individuals, even healthy individuals, requires further research.

Several studies have shown that children with MP-related extra-pulmonary diseases higher IgE levels exhibit more severe clinical symptoms and complications. Therefore, it is speculated that IgE might be a biomarker for complications after MP infection [35, 36]. Respiratory virus infection is a further clinical complication for children with MP, and may increase the severity of disease [37]. Although clinical presentation and laboratory tests were not statistically significant between children with single MP infection and children with MP and virus coinfection in this study, the impact of particular virus on MP infection outcomes and the clinical manifestations of MPP require further investigation to find virus-related biomarkers. Therefore, it is recommended that pediatricians consider the co-infection of MP and respiratory viruses in the diagnosis and treatment of respiratory diseases. Medical workers should adhere to strict hand washing procedures to avoid more serious complications in children infected with MP.

Our study has several limitations. First, some important laboratory information such as neutrophils, thrombocytes count, hemoglobin level, lactate dehydrogenase, aspartate transaminase and alanine transaminase are not included and do not provide full definition of study population characteristics. Second, the sample size was limited and all the samples were collected from one sentinel hospital in urban areas of Beijing, our findings may not be representative of viral spectrum in children diagnosed with MPP in China. Finally, we only focused on co-infection with multiple viruses, the co-infection of MP with a specific virus may provide greater insight into the impact of certain virus co-infections on MPP.


NGS analysis revealed 11 virus families in both MP-positive and MP-negative samples, the highest number of reads was the Pneumoviridae and HRSV in MP-positive samples, which was much higher than MP-negative samples. Among MP-positive cases, IFV co-infection were the most cases. By contrast with MP isolated infection children, the clinical symptoms of children with MP combined with virus infection had no statistically difference. Our study provided a theoretical basis for the development of effective prevention and treatment strategies for MPP.

Data availability

The data set supporting the conclusions of this paper is included in the article.



Community-acquired pneumonia


Mycoplasma pneumoniae


Nasopharyngeal aspirates


Mycoplasma pneumoniae pneumonia


Human respiratory syncytial virus


Human adenovirus


Influenza virus


Metagenomic next generation sequencing


Respiratory tract infections


Human parainfluenza virus


Human metapneumovirus


Human polyomavirus


Human rhinovirus


Human bocavirus


Human coronavirus


Torque teno virus


Cycle threshold


Next-generation sequencing


Polymerase chain reaction


Legionella pneumophila


Influenza-like symptoms.


  1. Waites KB, Xiao L, Liu Y, Balish MF, Atkinson TP. Mycoplasma pneumoniae from the respiratory tract and beyond. Clin Microbiol Rev. 2017;30(3):747–809.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Liu J, Wang M, Zhao Z, Lin X, Zhang P, Yue Q, et al. Viral and bacterial coinfection among hospitalized children with respiratory tract infections. Am J Infect Control. 2020;48(10):1231–6.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Lee KL, Lee CM, Yang TL, Yen TY, Chang LY, Chen JM, et al. Severe Mycoplasma pneumoniae pneumonia requiring intensive care in children, 2010–2019. J Formos Med Assoc. 2021;120(1 Pt 1):281–91.

    Article  PubMed  Google Scholar 

  4. Khoury T, Sviri S, Rmeileh AA, Nubani A, Abutbul A, Hoss S, et al. Increased rates of intensive care unit admission in patients with Mycoplasma pneumoniae: a retrospective study. Clin Microbiol Infect. 2016;22(8):711–4.

    Article  CAS  PubMed  Google Scholar 

  5. Kutty PK, Jain S, Taylor TH, Bramley AM, Diaz MH, Ampofo K, et al. Mycoplasma pneumoniae among children hospitalized with community-acquired pneumonia. Clin Infect Dis. 2019;68(1):5–12.

    CAS  PubMed  Google Scholar 

  6. Meyer SP, Unger WW, Nadal D, Berger C, Vink C, van Rossum AM. Infection with and carriage of mycoplasma pneumoniae in children. Front Microbiol. 2016;7:329.

    Google Scholar 

  7. Sondergaard MJ, Friis MB, Hansen DS, Jorgensen IM. Clinical manifestations in infants and children with mycoplasma pneumoniae infection. PLoS ONE. 2018;13(4):e195288.

    Article  Google Scholar 

  8. Gao J, Xu L, Xu B, Xie Z, Shen K. Human adenovirus coinfection aggravates the severity of Mycoplasma pneumoniae pneumonia in children. BMC Infect Dis. 2020;20(1):420.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Poddighe D. Extra-pulmonary diseases related to Mycoplasma pneumoniae in children: recent insights into the pathogenesis. Curr Opin Rheumatol. 2018;30(4):380–7.

    Article  PubMed  Google Scholar 

  10. Yeh JJ, Wang YC, Hsu WH, Kao CH. Incident asthma and Mycoplasma pneumoniae: a nationwide cohort study. J Allergy Clin Immunol. 2016;137(4):1017-23.e6.

    Article  PubMed  Google Scholar 

  11. Krafft C, Christy C. Mycoplasma pneumonia in children and adolescents. Pediatr Rev. 2020;41(1):12–9.

    Article  PubMed  Google Scholar 

  12. Oumei H, Xuefeng W, Jianping L, Kunling S, Rong M, Zhenze C, et al. Etiology of community-acquired pneumonia in 1500 hospitalized children. J Med Virol. 2018;90(3):421–8.

    Article  PubMed  Google Scholar 

  13. Honkinen M, Lahti E, Osterback R, Ruuskanen O, Waris M. Viruses and bacteria in sputum samples of children with community-acquired pneumonia. Clin Microbiol Infect. 2012;18(3):300–7.

    Article  CAS  PubMed  Google Scholar 

  14. Michelow IC, Olsen K, Lozano J, Rollins NK, Duffy LB, Ziegler T, et al. Epidemiology and clinical characteristics of community-acquired pneumonia in hospitalized children. Pediatrics. 2004;113(4):701–7.

    Article  PubMed  Google Scholar 

  15. Wang C, Zhou S, Xue W, Shen L, Huang W, Zhang Y, et al. Comprehensive virome analysis reveals the complexity and diversity of the viral spectrum in pediatric patients diagnosed with severe and mild hand-foot-and-mouth disease. Virology. 2018;518:116–25.

    Article  CAS  PubMed  Google Scholar 

  16. Chen Z, Shang Y, Zhao S, et al. Expert consensus on diagnosis and treatment of mycoplasma pneumoniae pneumonia in children (2015). Chin J Appl Clin Pediatr. 2015;30(17):1304–8.

    Google Scholar 

  17. Meyer Sauteur PM, Unger WWJ, van Rossum AMC, Berger C. The art and science of diagnosing Mycoplasma pneumoniae Infection. Pediatr Infect Dis J. 2018;37(11):1192–5.

    Article  PubMed  Google Scholar 

  18. El SM, Fouda EM, Ibrahim HM, Fathy MM, Husseiny AA, Khater WS, et al. Microbial etiology of community-acquired pneumonia among infants and children admitted to the pediatric hospital, Ain Shams university. Eur J Microbiol Immunol (Bp). 2016;6(3):206–14.

    Article  Google Scholar 

  19. Kayser FH. Changes in the spectrum of organisms causing respiratory tract infections: a review. Postgrad Med J. 1992;68(Suppl 3):17–23.

    Google Scholar 

  20. Zheng X, Lee S, Selvarangan R, Qin X, Tang YW, Stiles J, et al. Macrolide-resistant Mycoplasma pneumoniae, United States. Emerg Infect Dis. 2015;21(8):1470–2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Rao S, Nyquist AC. Respiratory viruses and their impact in healthcare. Curr Opin Infect Dis. 2014;27(4):342–7.

    Article  PubMed  Google Scholar 

  22. Leung CM, Li D, Xin Y, Law WC, Zhang Y, Ting HF, et al. MegaPath: sensitive and rapid pathogen detection using metagenomic NGS data. BMC Genom. 2020;21(Suppl 6):500.

    Article  CAS  Google Scholar 

  23. Richter J, Panayiotou C, Tryfonos C, Koptides D, Koliou M, Kalogirou N, et al. Aetiology of acute respiratory tract infections in hospitalised children in cyprus. PLoS ONE. 2016;11(1):e147041.

    Article  Google Scholar 

  24. Kim EK, Youn YS, Rhim JW, Shin MS, Kang JH, Lee KY. Epidemiological comparison of three Mycoplasma pneumoniae pneumonia epidemics in a single hospital over 10 years. Korean J Pediatr. 2015;58(5):172–7.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Tian DD, Jiang R, Chen XJ, Ye Q. Meteorological factors on the incidence of MP and RSV pneumonia in children. PLoS ONE. 2017;12(3):e173409.

    Article  Google Scholar 

  26. Diaz MH, Winchell JM. The evolution of advanced molecular diagnostics for the detection and characterization of Mycoplasma pneumoniae. Front Microbiol. 2016;7:232.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Zheng YX, Chen J, Kong DC, Pan H, Zhou YQ, Chen ML, et al. Pathogenic characteristics of hospitalized severe acute respiratory infections in Shanghai, China, 2015–2017. Zhonghua Liu Xing Bing Xue Za Zhi. 2019;40(8):911–6.

    CAS  PubMed  Google Scholar 

  28. Kalenahalli KJ, Kumar NA, Chowdary KV, Sumana MS. Fatal swine influenza A H1N1 and Mycoplasma pneumoniae coinfection in a child. Tuberk Toraks. 2016;64(3):246–9.

    Article  PubMed  Google Scholar 

  29. Duenas ME, Jaramillo CA, Correa E, Torres-Duque CA, Garcia C, Gonzalez M, et al. Virus and Mycoplasma pneumoniae prevalence in a selected pediatric population with acute asthma exacerbation. J Asthma. 2016;53(3):253–60.

    Article  Google Scholar 

  30. He J, Liu M, Ye Z, Tan T, Liu X, You X, et al. Insights into the pathogenesis of Mycoplasma pneumoniae (Review). Mol Med Rep. 2016;14(5):4030–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Chiu CY, Chen CJ, Wong KS, Tsai MH, Chiu CH, Huang YC. Impact of bacterial and viral coinfection on Mycoplasmal pneumonia in childhood community-acquired pneumonia. J Microbiol Immunol Infect. 2015;48(1):51–6.

    Article  PubMed  Google Scholar 

  32. Hino S. TTV, a new human virus with single stranded circular DNA genome. Rev Med Virol. 2002;12(3):151–8.

    Article  CAS  PubMed  Google Scholar 

  33. Peng YH, Nishizawa T, Takahashi M, Ishikawa T, Yoshikawa A, Okamoto H. Analysis of the entire genomes of thirteen TT virus variants classifiable into the fourth and fifth genetic groups, isolated from viremic infants. Arch Virol. 2002;147(1):21–41.

    Article  CAS  PubMed  Google Scholar 

  34. Pifferi M, Maggi F, Caramella D, De Marco E, Andreoli E, Meschi S, et al. High torquetenovirus loads are correlated with bronchiectasis and peripheral airflow limitation in children. Pediatr Infect Dis J. 2006;25(9):804–8.

    Article  PubMed  Google Scholar 

  35. Zhou L, Li Y, Xu Z, Peng X, Gong X, Yang L. Increased total serum immunoglobulin E is likely to cause complications of Mycoplasma pneumoniae pneumonia in children. Front Cell Infect Microbiol. 2021;11:783635.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Poddighe D, Marseglia GL. Is There any relationship between extra-pulmonary manifestations of Mycoplasma Pneumoniae infection and atopy/respiratory allergy in children? Pediatr Rep. 2016;8(1):6395.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Shinozaki T, Sasahara K, Iwami E, Kuroda A, Matsuzaki T, Nakajima T, et al. A case of influenza B and Mycoplasma pneumoniae coinfection in an adult. Case Rep Infect Dis. 2018;2018:3529358.

    PubMed  PubMed Central  Google Scholar 

Download references


We would like to thank the Department of Pediatrics, Beijing Friendship Hospital (Capital Medical University) for providing the NPA specimens. We would also like to thank Liwen Bianji, PhD, from Edanz Group (, for editing the English text of this manuscript.


This work was supported by the Key Technologies R&D Program of the National Ministry of Science, China (2018ZX10305-410-003-002, 2018ZX10713-002).

Author information

Authors and Affiliations



QG, CW and LSZ conceived and designed the experiments. QG, LHY, CW, FLM and AJC performed the experiments. QG, LLL, YMH and SSC analyzed the data; QG, LLL, CW and JJT contributed reagents/materials/analysis tools; QG and LSZ wrote the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Lishu Zheng.

Ethics declarations

Ethics approval and consent to participate

The project was approved by the Ethical Committee of National Institute for Viral Disease Control and Prevention, China CDC, and the committee’s reference number is 021, IVDC2017.

Consent for publication

Written informed consent for specimen collection, testing and publication was obtained from the individuals or their parents.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

. The tables of the filtering parameters setting using Prinseq-lite software, primers and probes used to detect co-infection respiratory viruses, type analysis of common virus-related sequence viral in MP positive and MP negative (TTV has been removed). Table S1: The filtering parameters setting using Prinseq-lite software. Table S2: Primers and probes used to detect co-infection respiratory viruses. Table S3: Type analysis of common virus-related sequence viral in MP positive and MP negative (TTV has been removed).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Guo, Q., Li, L., Wang, C. et al. Comprehensive virome analysis of the viral spectrum in paediatric patients diagnosed with Mycoplasma pneumoniae pneumonia. Virol J 19, 181 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Next-generation sequencing
  • Pathogen spectrum
  • Metagenomics analysis
  • Respiratory virus
  • Mycoplasma pneumoniae