Skip to main content

Bacterial composition and colony structure of the lower respiratory tract in infants and children with recurrent wheezing: a case–control study



The bacterial load of the human lower respiratory tract is at least several times lower than that of the other parts of the body. This study aimed to identify the bacterial composition and colony structure of the lower respiratory tract in infants and children with recurrent wheezing compared with those of children with a bronchial foreign body and clarify whether the length of wheezing in infants can contribute to differences in the lower respiratory tract’s bacterial colony structure.


We collected specimens of alveolar lavage fluid from 48 infants and children who underwent fiberoptic bronchoscopy and were divided into groups: A1 (multiple wheezing: wheezing more than three times in < 1 month), A2 (persistent wheezing: > 1 month), and B (bronchial foreign body; control group). We analyzed the bacterial community structure of alveolar lavage fluid using high-throughput sequencing. The richness and diversity of the microbial communities were assessed by α and β diversity analyses.


A total of 6,644 operational taxonomic units (OTUs) were obtained based on the Illumina Nova sequencing platform and clustered according to those that met the 97% identity threshold, followed by species annotation of the OTU sequences. In the annotation results, 2,608 (39.25%) OTUs were annotated at the genus level. At the genus level, Sphingomonas and Phyllobacterium were significantly higher in group A1 than in group B. There were significantly more Phyllobacterium in group A2 than in group B. Prevotella, Neisseria, and Haemophilus were higher in group B than in groups A1 and A2. The differences in the between-group α and β diversity analyses were statistically significant. The microbial diversity in groups A1 and A2 was significantly less than that in group B, but there was no statistical difference in bacterial community diversity between groups A1 and A2.


Recurrent wheezing in infants and children is more likely due to alterations in the overall bacterial microecology and disruption of host respiration and immune homeostasis than the effects of a single bacterium.


Wheezing is a common respiratory symptom in children, especially in infants and young children, and about one-third of children will have at least one episode of wheezing by the age of 3 years, and 1–2% of children are hospitalized for severe respiratory symptoms [1]. Wheezing resolves with age in some infants and children, while others develop persistent wheezing or asthma [2]. Recurrent wheezing in infants and children is often attributed to airway injury, inflammation, and immune responses due to infection, allergy, and environmental factors. However, few reports have shown whether bacterial flora disorders in the lower respiratory tract are associated with recurrent wheezing in infants and children. The human lower respiratory tract’s bacterial load is at least several times lower than that of the other parts of the body and may be highly plastic and susceptible to external forces [3, 4].

The earliest description of symbiotic microorganisms in the body can be traced back to the study of Enterobacter in the digestive tract of infants by Escherich more than a century ago [5]. Over the decades, the microbiome’s importance has gradually gained recognition, from early descriptions of gut bacteria in response to diet to the detailed mechanisms of host-microbe interactions found in a large body of work documented in recent years, including pathways affecting systemic immunity, cardiovascular health, and neurological development [6,7,8,9]. However, research and understanding of the bacterial community structure of the lung remain scarce as it was long considered a sterile environment, except when infected by respiratory pathogenic microbes. This deeply-rooted notion led to the removal of lungs from the sampling list in the first round of the Human Microbiome Project (HMP) [10]. The second round of the HMP found viruses and fungi in the lung in addition to bacteria, but the health and disease effects of each microorganism are still not well understood [11,12,13,14].

In recent years, recurrent wheezing in infants and preschoolers has become very common. The increase in emergency department visits and hospitalization has societal as well as personal effects [15]. The prognosis of recurrent wheezing in infants and children varies, and the onset and duration of wheezing may sometimes be inevitably associated with asthma [16, 17], which has been proven to be caused by impaired lung function due to bacterial infection. Streptococcus pneumoniae, Moraxella catarrhalis, and Haemophilus influenzae are the microbiota primarily associated with the exacerbation of asthma [18, 19]. However, it remains unclear whether prolonged, persistent wheezing versus multiple wheezing episodes cause different changes in the lower respiratory tract’s colony structure.

The ideal method of specimen collection is to collect lung tissue samples by isolating microbial deoxyribonucleic acid (DNA) from host DNA, but this method cannot be routinely applied to every disease. Although sputum specimens are easy to collect, contact with the oral cavity and upper respiratory tract is unavoidable during the collection process. In contrast, bronchoalveolar lavage fluid is immune to upper airway flora contamination and closely mimics the lower airway’s microbiome but cannot wholly sample the lung’s deepest airways. Bronchial foreign bodies are a typical pediatric emergency, commonly affecting infants and children aged six months to three years. Most children with bronchial foreign bodies have a short duration of illness, and mechanical damage to the airway should not adjust the microbiota’s composition over this short period. Therefore, we chose the results of alveolar lavage fluid from children with bronchial foreign bodies instead of lower airway bacterial colonization from healthy infants and children as a control.

Because persistent and repeated wheezing may cause changes in the colony structure of lower respiratory tract microorganisms, the differences in the genera represented by the lower respiratory tract microbiota in infants and children with recurrent wheezing should be elucidated to provide a basis for future insights into the critical role of lower respiratory tract bacteria in the development of recurrent wheezing in infants and children.

Therefore, the first aim of this study was to describe the composition and diversity of bacterial colonies in the lower respiratory tract of children with recurrent wheezing. The second aim was to determine whether persistent wheezing affects infants' and children's microbial colony structure and diversity. The third objective was to identify the dominant bacterial genera in infants and children with recurrent wheezing, provide new ideas for future insights into the role of bacterial genera in regulating airway immunity, and provide new directions for the prevention of diseases that cause recurrent wheezing.

Materials and methods

Participants and research design

To analyze the relevant bacterial communities, 48 specimens were collected from the alveolar lavage fluid of children with bronchial foreign bodies or recurrent wheezing admitted to Chengdu Women’s and Children’s Center Hospital. The participants were divided into groups A1 (multiple wheezing: repeated wheezing for a short duration each time; n = 13), A2 (persistent wheezing: wheezing for more than one month without previous wheezing; n = 16), and B (bronchial foreign body: control group; n = 19). The criteria for enrolling children in the wheezing group were (1) age 6–36 months; (2) previous wheezing ≥ three times, wheezing duration < 1 month, or wheezing duration > 1 month; (3) bronchoscopy and alveolar lavage; and (4) no history of antibiotic use for 1 week, no heart disease, tuberculosis, bronchial foreign body, bronchopulmonary dysplasia, abnormal immune function, or abnormal imaging tests. Children in group B had a history of foreign body choking and coughing within the past 3 days, and a foreign body was visible on bronchoscopy.

Sample collections and molecular analyses

Pulmonary function tests

Infant tidal pulmonary function (tidal volume, peak time ratio, peak volume ratio, respiratory ratio), fractional exhaled nitric oxide (Fe–NO), and routine blood tests were completed within 24 h after admission.


The guardians of all children signed informed consent forms. Preoperative routine examination and other related tests were performed, and the participants fasted from food and water for 4–6 h. Bronchoscopy was performed with each participant in the supine position under intravenous compound anesthesia with transglottic mask ventilation. The appropriate fiberoptic bronchoscope (OLYMPUS, Tokyo, Japan) was selected according to the child’s age and weight. The bronchoscope was inserted via laryngeal mask airway (LMA) and the child’s epiglottis, tracheal ramus, and each lobe and segment of the bronchus were assessed for lesions.

Extraction of alveolar lavage fluid DNA

The DNA in the sample was extracted by the sodium dodecyl sulfate method, and the concentration and purity of the DNA were detected by agarose gel electrophoresis and diluted. Using the diluted genomic DNA as a template, we used foreign objects with Barcode, selected the sequencing area, and used High-Fidelity Polymerase Chain reaction (PCR) Master Mix with GC Buffer (New England Biolabs Phusion®, US) and high-fidelity enzymes for PCR. The amplified regions included the Archaea 16S V4-V5/Archaea 16S V8, 16S V3-V4/16S V4-V5/16SV5-V7; 18S V9, and ITS2 regions.

Purification of PCR products, library construction, and up-sequencing

The PCR products were electrophoresed using agarose gel; the quantified PCR products were purified using magnetic beads. Using enzyme labeling quantification, equal amounts of PCR products were mixed according to their concentrations, and PCR products were detected by agarose gel electrophoresis. The products were then recovered for the target bands using the gel recovery kit (Qiagen, Germany). The library was constructed using a DNA PCR-free sample preparation kit (TruSeq®, USA). Qubit and quantitative PCR were used to quantify the constructed library, which was then sequenced using NovaSeq 6000 (Novogene Co., Ltd., Beijing, China).

Data processing

The sparse algorithm was used to cluster the valid tag sequences into OTUs, which met the 97% identity threshold by default, for all samples. Species annotation of OTU sequences was performed using the Mothur method with the small subunit ribosomal ribonucleic acid (rRNA) database of SILVA138 (with a threshold value of 0.8–1) to obtain taxonomic information and count each sequence at each taxonomic level (kingdom, phylum, class, order, family). The community composition of each sample was counted at each taxonomic level (kingdom, phylum, class, order, family, and genus). The phylogenetic relationships of all OTUs were obtained by rapid multiple sequence matching using MUSCLE software. Finally, the least amount of data in the samples were used as the criterion for homogenization, and the subsequent α and β diversity analyses were based on the homogenized data.

Statistical analysis and visualization

We drew a light curve to estimate whether the sample’s current sequencing depth could reveal the microbial community’s diversity [20]. QIME software was used to calculate the α diversity index, and differences between the α diversity indices were checked using rank-sum tests [21]. β diversity analysis can detect the similarity of bacterial colony structure between different groups through UniFrac distance, and the non-metric multidimensional scaling (NMDS) method was used for visualization [22,23,24]. The contribution of the differences between the two groups for each species was quantified using the similarity percentage method (SIMPER) [25, 26]. UPGMA (Unweighted Pair-group Method with Arithmetic Mean) is a more commonly used method for cluster analysis; it can be used to study the similarity between different samples and to perform cluster analysis on samples. In addition, linear discriminant analysis of effect size (LEfSe) was used to screen taxa with rich differences between groups and discover biomarkers [27]. Sequence data analysis and data curation were performed using QIIME, R package (v3.2.0), and IBM SPSS Statistics (v21.0). Independent sample t-tests and chi-square tests were used to analyze the statistically significant differences between continuous and categorical variables, respectively. One-way analysis of variance was used for continuous variables among the three groups (*P < 0.05).


Patients’ characteristics

All groups had a higher proportion of male participants, but the between-group sex differences were not statistically significant. The white blood cell count and the percentage of eosinophils in routine blood tests were not significantly different and were within the normal range. We measured fractional exhaled nitric oxide (Fe–NO) and tidal pulmonary function in only 11 cases in group A1 and 13 cases in group A2 (data were missing in the other five cases), in which the mean values of the ratio of tidal expiratory flow to total expiratory time (TPEF/TE) and the ratio of peak expiratory flow to total expiratory flow (VPEF/VE) were lower in group A1 than group A2, suggesting severe obstructive ventilation. However, they were not significantly different between the two groups. The general information regarding the participants is shown in Table 1.

Table 1 Participants’ baseline characteristics

Bacterial abundance and Distribution

The PCR-free library was sequenced by paired-end sequencing using the NovaSeq 6000 (Illumina, USA). By splicing the reads, an average of 87,349 tags were sequenced per sample, and 81,331 valid data points were obtained after quality control, resulting in 61,798 valid data points in quality control and a quality control efficiency of 71.37%. The sequences were clustered into OTUs that met the 97% identity threshold, and 6,644 OTUs were identified. Between 182 and 1326 OTUs (mean = 565) were identified per sample in group A, and 282–1334 OTUs (mean = 580) were identified per sample in group B.

The richness sparsity curves for groups A1, A2, and B confirmed that all the samples had sufficient sequencing depth (Fig. 1A), and the rank clustering curves reflected the samples’ uniform distribution and high species richness (Fig. 1B). Based on the annotation results, the 20 most abundant species at the phylum, order, family, and genus levels were selected for each sample or group, and a cumulative bar graph of the species’ relative abundance was generated (Fig. 1C).

Fig. 1
figure 1

Multiple rarefaction curves for all the alveolar lavage fluid samples. Multiple rarefaction curves of the (A) richness and (B) rank abundance (B). C Columnar cumulative chart of relative abundance of the top 20 species with the greatest abundance at the phylum, class, order, family, and genus levels. D Distribution of dominant bacteria in groups A1, A2, and B at the phylum and genus levels

A total of 11 phyla, 28 classes, 67 orders, 99 families, and 165 genera were detected in the alveolar lavage fluid specimens, which could be annotated to 70.64% at the phylum level, 68.68% at the class level, 64.77% at the order level, 56.17% at the family level, and 39.25% at the genus level. The major bacterial distributions were characterized based on relative taxonomic abundance (Fig. 1D) as Proteobacteria (A1 = 49.1%, A2 = 43.2%, B = 25.5%), Firmicutes (A1 = 30.9%, A2 = 36.9%, B = 34.9%), Bacteroidota (A1 = 9.1%, A2 = 10.1%, B = 25.9%), and Fusobacteriota (A1 = 0.5%, A2 = 0.4%, B = 5.0%), and at the genus level with Streptococcus (A1 = 18.8%, A2 = 16.4%, B = 20.1%) Stenotrophomonas (A1 = 9.2%, A2 = 14.9%, B = 0.1%), Prevotella (A1 = 4.5%, A2 = 1.4%, B = 15.1%), Sphingomonas (A1 = 13.4%, A2 = 6.5%, B = 0.1%), Staphylococcus (A1 = 0.8%, A2 = 4.3%, B = 0.1%), Phyllobacterium (A1 = 11.2%, A2 = 6.3%, B = 0.0%), Neisseria (A1 = 0.1%, A2 = 0.3%, B = 6.1%), Haemophilus (A1 = 0.5%, A2 = 1.1%, B = 9.6%), Brevundimonas (A1 = 4.6%, A2 = 3.7%, B = 0.5%), Porphyromonas (A1 = 0.0%, A2 = 0.0%, B = 3.5%), Fusobacterium (A1 = 0.3%, A2 = 0.01%, B = 4.0%), and Veillonella (A1 = 1.5%, A2 = 2.9%, B = 6.8%)( Table 2).

Table 2 Major bacterial distributions characterized based on relative taxonomic abundance

Genus groups with significant differences at the level of bacterial phylum and genus were assessed by performing t-tests for groups A1, A2, and B. At the phylum level, the number of Bacteroidetes bacteria in groups A1 and A2 were lower than in groupB, and the number of Proteobacteria were higher in groups A1 and A2 than in group B. At the genus level, the numbers of Sphingomonas and Phyllobacterium in group A1 were higher than those in group B. Furthermore, the number of Phyllobacterium in group A2 was higher than that in group B. Prevotella, Neisseria, and Haemophilus were more in group B than in groups A1 and A2.

Bacterial diversity and colony structure analysis

There was no significant difference in the number of OTUs directly observed among the three groups (Fig. 2A). The α diversity in group B was significantly different from the Simpson’s Diversity Index of groups A1 and A2, and the diversity and evenness of group B were higher than those of groups A1 and A2 (Fig. 2B). In the analysis of β-diversity, we found that the Bray–Curtis distance of the bacterial community in group B was different from that in groups A1 and A2. There was no statistically significant difference in the Bray–Curtis distance between group A1 and group A2 (Fig. 2C). NMDS plots showed similar distances between groups A1 and A2 and were distant from group B, suggesting that the colony structures of groups A1 and A2 are highly similar and differ from those of group B (Fig. 2D). The specific information regarding this is shown in Additional file 1: Appendix 1. UPGMA clustering of microbiome classification profiles in all groups with UniFrac distance distances showed relatively significant difference between the control and wheezing groups (Fig. 2E).

Fig. 2
figure 2

Bacterial diversity in alveolar lavage fluid samples of children in groups A1, A2, and B. A The number of observed species. B α diversity shown by Simpson's Diversity Index. C β diversity expressed as the Bray–Curtis distance in the alveolar lavage fluid samples. D NMDS based on the Bray–Curtis distances at the OTU level which met the 97% identity threshold. E UPGMA clustering of the microbiome classification maps in the three groups using the UniFrac distance showed relatively clear separation between groups

Differential microbiota composition

The species abundance data between groups were analyzed using LEfSe analysis to detect the differential species between groups by rank-sum test and to assess the effect size of differential species by linear discriminant analysis (LDA). Histograms of the distribution of LDA values of the differential species and evolutionary branching of the differential species were plotted. We found 47 biomarkers with an LDA score > 4, including Haemophilus, Neisseria, Prevotella, Sphingomonadales, Sphingomonas, Delftia, and Veillonella.

These differentially rich taxa can be considered potential biomarkers (Fig. 3A; see S1 for details) and are shown in LEfSe plots in order of phylum to genus and group importance (the microbial population with similar coloring as the group is more important in that group; Fig. 3B). The contribution of the top ten species in abundance to the variability between each group was quantified by SIMPER (bubble size represents the relative abundance of the species, and contribution is the species’ contribution to the variability between the two groups; Fig. 3C). Anosim’s test showed that the differences between groups B, A1, and A2 were greater than those within the groups (Table 3).

Fig. 3
figure 3

The linear discriminant analysis score (A) and taxonomy tree (B) of the significant between-group differences. The top 10 genera and their abundances in the difference contribution between groups A1, A2, and B

Table 3 Analysis of the differences between Anosim groups


Humans have evolved symbiotically with microorganisms, and host-microbial interactions are gradually being elucidated, with bacterial colonization of the human host impacting a variety of diseases [28]. Microbial colonization can be detected in the naso- and oropharynx of newborns within the first 5 min after birth, indicating that colonization of the upper respiratory tract has begun [29]. Initially, the immune system of newborns will adjust so that the microbiota can colonize the respiratory tract smoothly. The stable symbiotic microbiota will gradually promote the improvement of the neonatal immune system, and the establishment of this microbiota gives the body its ability to fight pathogenic bacteria. The resident microbiota and the infant’s immune system are restricted in many ways and regulated by metabolites and microorganisms. Immunoglobulin A, present in breast milk, can help establish a stable host-microbiota relationship [30]. Bacteria, viruses, and fungi are the primary pathogens responsible for human respiratory infections. However, in 30% of respiratory cases, pathogens remain undetected by conventional methods [31].

Data released in 2017 suggest that the bacterial microbiota differs slightly between bronchial lavage fluid and endotracheal samples [32]. Our analysis of alveolar lavage fluid microbiomes from 48 children showed that the lower respiratory flora of those in the wheezing group (showing both recurrent and persistent wheezing) comprised mainly the phylum Proteobacteria, followed by the thick-walled Firmicutes and the Bacteroidota.

In contrast, the alveolar lavage fluid specimens from children in the foreign body group (control group) were mainly composed of the phylum Firmicutes, followed by Proteobacteria and Bacteroidota. There is a relative lack of data on children’s lower respiratory flora, and studies of oropharyngeal microflora have confirmed that the phyla mentioned above were found to be dominant in both asthmatic and healthy populations [33, 34]. A neonatal-based study showed that H. influenzae, M. catarrhalis, and S. pneumoniae were associated with wheezing outcomes within the first 5 years of life. Although there is no direct evidence for the cellular impact of these bacterial fixations, these associations may be consistent with significant neutrophil elevation in children with severe recurrent wheezing [18, 35].

We compared the wheezing and foreign body groups in terms of microbial genera and observed abundant Streptococcus spp. in each group. The numbers of bacteria of genera Sphingomonas and Phyllobacterium were significantly higher in the recurrent wheezing group than in the control group; Phyllobacterium spp. were also higher in the prolonged wheezing group than in the control group, and Prevotella spp., Neisseria spp., and Haemophilus spp. were higher in the control group than in the wheezing groups.

The low proportions of Moraxella spp., Staphylococcus spp., and Haemophilus spp. in the wheezing group were different from the increased proportions of Staphylococcus spp., Moraxella spp., and Haemophilus spp., which have been previously reported to increase the risk of chronic wheezing in children [35,36,37,38,39,40], and require further investigation. The high proportion of Phyllobacterium spp. has not been reported to be associated with wheezing in infants and children. However, Wen et al. found a higher enrichment of Phyllobacterium foliaris among oropharyngeal microorganisms in children with the influenza A virus; the exact mechanism remains unclear [41].

Reduced numbers of Prevotella spp. were found in both wheezing groups. Prevotella spp. is a common commensal fixer at mucosal sites and is the dominant respiratory genus [42]. The homeostatic role of Prevotella spp. in the healthy lung remains largely unclear. Inflammation and pathology in the chronic obstructive pulmonary disease-like lungs of mice induced by lipopolysaccharide or lipopolysaccharide elastase inhalation were reduced by this genus and mediated the growth of Pseudomonas and Lactobacillus spp. This suggests that the reduction of Prevotella spp. (possibly by creating a microenvironment unsuitable for the growth of the genus, which leads to a decrease in its abundance) in patients with wheezing may not be a pre-disease risk factor.

Some authors have compared the inflammatory properties of Prevotella spp. in the lungs of healthy individuals with those of bacilli associated with pulmonary disease in asthmatic patients (H. influenzae B, H. influenzae, and M. catarrhalis). Prevotella was found to induce similar levels of CD83, CD86, and CD40 activation marker surface expression, but compared with Proteobacteria, it reduced the levels of interleukin (IL)-12p70, IL-23, and IL-10 cytokines in monocyte-derived dendritic cells. The lower inflammatory capacity of Prevotella was also confirmed in mice, where Prevotella reduced the production of macrophage inflammatory protein 2-α (IL-8), tumor necrosis factor-a (TNF-a), and thymic stromal lymphopoietin by lung stromal cells, as well as the production of TNF-a by lung immune cells, suggesting that Prevotella may be well tolerated in the lungs [43,44,45].

A high proportion of Streptococcus spp. was found in all three groups, and it has been shown that increased abundance of Moraxella and Streptococcus spp. in nasopharyngeal swabs from infants with capillary bronchitis is associated with recurrent wheezing in children at 3 years of age [46]. This may suggest that infant wheezing is more likely to result from a disruption of overall bacterial microecology and host respiratory and immune homeostasis rather than the influence of a single bacterium. Bacterial detection in the lower airway may be due to an imbalance between constant exposure to upper respiratory microorganisms through micro-inhalation and clearance by endothelial cilia and coughing in the respiratory tract.

The exact reasons for the differences in the dominant genera are unclear and may lead to molecular changes in the mucin glycoproteins during long-term chronic inflammation; mucin glycans are a significant source of nutrients and adhesion sites for bacteria and fungi in diseased airways, and different microorganisms can utilize specific mucins, which influence airway fixation [47]. The physical and chemical factors in the airway (including, oxygen tension, pH, temperature, and mucus) vary with the airway’s anatomical position [48], which, together with irregular medication use in children with multiple and persistent wheezing, may lead to the establishment of new, stable, and complex microflora in the presence of immune dysregulation.

The microbial diversity of the wheezing group was lower than that of the foreign body group, and the diversity of microbial colonies was reduced. The exact mechanism by which bacteria cause wheezing in infants and children remains unclear. The immune balance between the bacteria and host organism may be disrupted. Bacteria and their products (such as endotoxins) can stimulate the production of cytokines in airway epithelial cells. Some cytokines (such as nuclear factor kappa B) can activate the IL-8 gene and then initiate positive feedback regulation to increase the production of neutrophils and inflammatory factors, such as IL-8, IL-6, and TNF-α. The increase in TNF-α production can increase the reactivity of the airway’s smooth muscles, thereby aggravating airway obstruction and causing wheezing [49,50,51].

Pulmonary function tests can help in judging the severity of diseases with wheezing as a symptom and also in evaluating clinical efficacy and prognosis [52]. The most commonly used lung function test for infants and young children is the tidal breathing lung function test [53]. Both TPEF/TE and VPEF/VE are important parameters reflecting obstructive ventilation disorders, with a normal range of 28–55%. Depending on the degree of obstruction, it is classified as mild obstruction (23–28%), moderate obstruction (15–22%), and severe obstruction (< 15%) [54]. The pulmonary function of the children we observed with wheezing suggested a decrease in severe obstructive ventilation. However, the mean TPEF/TE and VPEF/VE in the multiple wheezing groups were worse than those in the persistent wheezing group, and we speculated that perhaps the persistent wheezing was due to persistent airway spasm caused by some infection, while the more severe pulmonary impairment in the recurrent wheezing group may have caused a greater chance of asthma. However, why the duration of disease differs in patients with similar colony structure is subject to further research. We speculate that this may be related to insufficient and ineffective treatment of infants or children.

The main advantage of our study is that the matched case–control excludes the deviation caused by the influence of age, time, and sex. We also elucidated that the duration of wheezing can affect the structure of the bacterial microbial community. We used alveolar lavage fluid samples obtained via bronchoscopy, which is closer to the proper lower respiratory tract microbial colony composition, and the bronchoscopic foreign body group as the control, which is closer to the community composition of respiratory bacteria in children of the same age.

However, our study had some limitations. We primarily focused on the bacterial group of microorganisms and could not exclude the possibility that viruses and fungi might be involved in recurrent wheezing in children. As in other observational studies, our findings do not necessarily suggest causality. 16S rRNA sequencing allows the annotation between genus- and species-level identification of bacteria but does not provide the resolution of macrogenomic techniques (e.g., birdshot sequencing), especially for closely related species, such as Streptococcus spp. and Prevotella spp. Specific species may need further discussion.


Our study showed that the lower respiratory flora of children with wheezing consisted mainly of the phylum Proteobacteria, followed by the phylum Firmicutes and the phylum Bacteroidota, whereas the lower respiratory flora of healthy children may be mainly phylum Firmicutes, followed by Proteobacteria and Bacteroidota. A decrease in the number of Prevotella spp might be associated with recurrent wheezing in infants and children. It might be involved in the establishment of homeostasis in the pulmonary environment, the mechanism of this is still unclear. When comparing the diversity of bacterial colony structure in the lower respiratory tract of children with persistent wheezing and those with repeated wheezing, we found that there may not be any particular differences in microbiology between them. Wheezing in infants and children is more likely to be due to disruption of the overall bacterial microecology and host respiratory immune balance rather than the effect of a single bacterium.

Overall, our findings suggest that diagnostic methods for microbiota should be further refined and that longitudinal studies related to animal models, such as the alveolar lavage fluid overall fixation in children with recurrent wheezing, could be refined to explore the role of microbial colonization in children with wheezing. Additional studies are needed to confirm whether changes in overall microbial colonization can reverse the clinical symptoms of recurrent wheezing in infants and children.

Availability of data and materials

All data generated or analyzed in this study are included in this published article and its supplementary files.



Deoxyribonucleic acid


Fractional exhaled nitric oxide


Human Microbiome Project




Linear discriminant analysis


Linear discriminant analysis of effect size


Non-metric multidimensional scaling


Operational taxonomic unit


Ribosomal ribonucleic acid


Tumor necrosis factor-a


The ratio of tidal expiratory flow to total expiratory time


The ratio of peak expiratory flow to total expiratory flow


  1. Inoue Y, Shimojo N. Epidemiology of virus-induced wheezing/asthma in children. Front Microbiol. 2013;4:391.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Rodríguez-Martínez CE, Sossa-Briceño MP, Castro-Rodriguez JA. Factors predicting persistence of early wheezing through childhood and adolescence: a systematic review of the literature. J Asthma Allergy. 2017;10:83–98.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Charlson ES, Bittinger K, Haas AR, Fitzgerald AS, Frank I, Yadav A, et al. Topographical continuity of bacterial populations in the healthy human respiratory tract. Am J Respir Crit Care Med. 2011;184:957–63.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Hilty M, Burke C, Pedro H, Cardenas P, Bush A, Bossley C, et al. Disordered microbial communities in asthmatic airways. PLoS ONE. 2010;5:e8578.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  5. Escherich T. The intestinal bacteria of the neonate and breastfed infant. Rev Infect Dis. 1988;10:1220–5.

    Article  CAS  PubMed  Google Scholar 

  6. Torrey JC. The regulation of the intestinal flora of dogs through diet. J Med Res. 1919;39:415–47.

    CAS  PubMed  PubMed Central  Google Scholar 

  7. Belkaid Y, Hand TW. Role of the microbiota in immunity and inflammation. Cell. 2014;157:121–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Wang Z, Klipfell E, Bennett BJ, Koeth R, Levison BS, Dugar B, et al. Gut flora metabolism of phosphatidylcholine promotes cardiovascular disease. Nature. 2011;472:57–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Diaz Heijtz R, Wang S, Anuar F, Qian Y, Björkholm B, Samuelsson A, et al. Normal gut microbiota modulates brain development and behavior. Proc Natl Acad Sci U S A. 2011;108:3047–52.

    Article  PubMed  Google Scholar 

  10. NIH HMP Working Group, Peterson J, Garges S, Giovanni M, McInnes P, et al. The NIH Human Microbiome Project. Genome Res. 2009;19:2317–23.

    Article  CAS  Google Scholar 

  11. Integrative HMP (iHMP) Research Network Consortium. The Integrative Human Microbiome Project: dynamic analysis of microbiome-host omics profiles during periods of human health and disease. Cell Host Microbe. 2014;16:276–89.

    Article  CAS  Google Scholar 

  12. Lozupone C, Cota-Gomez A, Palmer BE, Linderman DJ, Charlson ES, Sodergren E, et al. Widespread colonization of the lung by Tropheryma whipplei in HIV infection. Am J Respir Crit Care Med. 2013;187:1110–7.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Lysholm F, Wetterbom A, Lindau C, Darban H, Bjerkner A, Fahlander K, et al. Characterization of the viral microbiome in patients with severe lower respiratory tract infections, using metagenomic sequencing. PLoS ONE. 2012;7:e30875.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Willger SD, Grim SL, Dolben EL, Shipunova A, Hampton TH, Morrison HG, et al. Characterization and quantification of the fungal microbiome in serial samples from individuals with cystic fibrosis. Microbiome. 2014;2:40.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Soh JE, Kim KM, Kwon JW, Kim HY, Seo JH, Kim HB, et al. Recurrent wheeze and its relationship with lung function and airway inflammation in preschool children: a cross-sectional study in South Korea. BMJ Open. 2017;7:e018010.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Sears MR. Predicting asthma outcomes. J Allergy Clin Immunol. 2015;136:829–36 quiz 837.

    Article  PubMed  Google Scholar 

  17. Owora AH, Becker AB, Chan-Yeung M, Chan ES, Chooniedass R, Ramsey C, et al. Wheeze trajectories are modifiable through early-life intervention and predict asthma in adolescence. Pediatr Allergy Immunol. 2018;29:612–21.

    Article  PubMed  Google Scholar 

  18. Bisgaard H, Hermansen MN, Buchvald F, Loland L, Halkjaer LB, Bønnelykke K, et al. Childhood asthma after bacterial colonization of the airway in neonates. N Engl J Med. 2007;357:1487–95.

    Article  CAS  PubMed  Google Scholar 

  19. Schwerk N, Brinkmann F, Soudah B, Kabesch M, Hansen G. Wheeze in preschool age is associated with pulmonary bacterial infection and resolves after antibiotic therapy. PLoS ONE. 2011;6:e27913.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Lundberg DS, Yourstone S, Mieczkowski P, Jones CD, Dangl JL. Practical innovations for high-throughput amplicon sequencing. Nat Methods. 2013;10:999–1002.

    Article  CAS  PubMed  Google Scholar 

  21. Li B, Zhang X, Guo F, Wu W, Zhang T. Characterization of tetracycline resistant bacterial community in saline activated sludge using batch stress incubation with high-throughput sequencing analysis. Water Res. 2013;47:4207–16.

    Article  CAS  PubMed  Google Scholar 

  22. Kruskal JB. Nonmetric multidimensional scaling: a numerical method. Psychometrika. 1964;29:115–29.

    Article  Google Scholar 

  23. Noval Rivas M, Burton OT, Wise P, Zhang YQ, Hobson SA, Garcia Lloret M, et al. A microbiota signature associated with experimental food allergy promotes allergic sensitization and anaphylaxis. J Allergy Clin Immunol. 2013;131:201–12.

    Article  CAS  PubMed  Google Scholar 

  24. Lozupone C, Knight R. UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol. 2005;71:8228–35.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Clarke KR. Non-parametric multivariate analyses of changes in community structure. Austral Ecol. 1993;18:117–43.

    Article  Google Scholar 

  26. Warton DI, Wright ST, Wang Y. Distance-based multivariate analyses confound location and dispersion effects. Methods Ecol Evol. 2012;3:89–101.

    Article  Google Scholar 

  27. Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, et al. Metagenomic biomarker discovery and explanation. Genome Biol. 2011;12:R60.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Moran CP, Shanahan F. Gut microbiota and obesity: role in aetiology and potential therapeutic target. Best Pract Res Clin Gastroenterol. 2014;28:585–97.

    Article  CAS  PubMed  Google Scholar 

  29. Dominguez-Bello MG, Costello EK, Contreras M, Magris M, Hidalgo G, Fierer N, et al. Delivery mode shapes the acquisition and structure of the initial microbiota across multiple body habitats in newborns. Proc Natl Acad Sci U S A. 2010;107:11971–5.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Marcobal A, Sonnenburg JL. Human milk oligosaccharide consumption by intestinal microbiota. Clin Microbiol Infect. 2012;18(suppl 4):12–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Heikkinen T, Järvinen A. The common cold. Lancet. 2003;361:51–9.

    Article  PubMed  PubMed Central  Google Scholar 

  32. McCauley LM, Webb BJ, Sorensen J, Dean NC. Use of tracheal aspirate culture in newly intubated patients with community-onset pneumonia. Ann Am Thorac Soc. 2016;13:376–81.

    Article  PubMed  Google Scholar 

  33. Qin T, Zhang F, Zhou H, Ren H, Du Y, Liang S, et al. High-level PM2.5/PM10 exposure is associated with alterations in the human pharyngeal microbiota composition. Front Microbiol. 2019;10:54.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Faner R, Sibila O, Agustí A, Bernasconi E, Chalmers JD, Huffnagle GB, et al. The microbiome in respiratory medicine: current challenges and future perspectives. Eur Respir J. 2017;49:28404649.

    Article  Google Scholar 

  35. Teo SM, Mok D, Pham K, Kusel M, Serralha M, Troy N, et al. The infant nasopharyngeal microbiome impacts severity of lower respiratory infection and risk of asthma development. Cell Host Microbe. 2015;17:704–15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Smits HH, Hiemstra PS, Prazeres da Costa C, Ege M, Edwards M, Garn H, et al. Microbes and asthma: opportunities for intervention. J Allergy Clin Immunol. 2016;137:690–7.

    Article  PubMed  Google Scholar 

  37. Fedosenko SV, Ogorodova LM, Karnaushkina MA, Kulikov ES, Deev IA, Kirillova NA. The airways microbial community composition in healthy individuals and bronchial asthma patients. Vestn Ross Akad Med Nauk. 2014;69:71–6.

    Article  Google Scholar 

  38. Depner M, Ege MJ, Cox MJ, Dwyer S, Walker AW, Birzele LT, et al. Bacterial microbiota of the upper respiratory tract and childhood asthma. J Allergy Clin Immunol. 2017;139:826-834.e13.

    Article  PubMed  Google Scholar 

  39. Castro-Nallar E, Shen Y, Freishtat RJ, Pérez-Losada M, Manimaran S, Liu G, et al. Integrating microbial and host transcriptomics to characterize asthma-associated microbial communities. BMC Med Genomics. 2015;8:50.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  40. Hyde ER, Petrosino JF, Piedra PA, Camargo CA Jr, Espinola JA, Mansbach JM. Nasopharyngeal Proteobacteria are associated with viral etiology and acute wheezing in children with severe bronchiolitis. J Allergy Clin Immunol. 2014;133:1220–2.

    Article  PubMed  Google Scholar 

  41. Wen Z, Xie G, Zhou Q, Qiu C, Li J, Hu Q, et al. Distinct nasopharyngeal and oropharyngeal microbiota of children with influenza A virus compared with healthy children. BioMed Res Int. 2018;2018:6362716.

    PubMed  PubMed Central  Google Scholar 

  42. Human Microbiome Project Consortium. Structure, function and diversity of the healthy human microbiome. Nature. 2012;486:207–14.

    Article  CAS  Google Scholar 

  43. Larsen JM, Steen-Jensen DB, Laursen JM, Søndergaard JN, Musavian HS, Butt TM, et al. Divergent pro-inflammatory profile of human dendritic cells in response to commensal and pathogenic bacteria associated with the airway microbiota. PLoS ONE. 2012;7:e31976.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Larsen JM, Musavian HS, Butt TM, Ingvorsen C, Thysen AH, Brix S. Chronic obstructive pulmonary disease and asthma-associated Proteobacteria, but not commensal Prevotella spp., promote toll-like receptor 2-independent lung inflammation and pathology. Immunology. 2015;144:333–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Yadava K, Pattaroni C, Sichelstiel AK, Trompette A, Gollwitzer ES, Salami O, et al. Microbiota promotes chronic pulmonary inflammation by enhancing IL-17a and autoantibodies. Am J Respir Crit Care Med. 2016;193:975–87.

    Article  CAS  PubMed  Google Scholar 

  46. Mansbach JM, Luna PN, Shaw CA, Hasegawa K, Petrosino JF, Piedra PA, et al. Increased Moraxella and Streptococcus species abundance after severe bronchiolitis is associated with recurrent wheezing. J Allergy Clin Immunol. 2020;145:518–527.e8.

    Article  CAS  PubMed  Google Scholar 

  47. Taylor SL, Wesselingh S, Rogers GB. hosts. Host-microbiome interactions in acute and chronic respiratory infections. Cell Microbiol. 2016;18:652–62.

    Article  CAS  PubMed  Google Scholar 

  48. Dickson RP, Erb-Downward JR, Freeman CM, Walker N, Scales BS, Beck JM, et al. Changes in the lung microbiome following lung transplantation include the emergence of two distinct Pseudomonas species with distinct clinical associations. PLoS ONE. 2014;9:e97214.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  49. McNamara PS, Flanagan BF, Hart CA, Smyth RL. Production of chemokines in the lungs of infants with severe respiratory syncytial virus bronchiolitis. J Infect Dis. 2005;191:1225–32.

    Article  CAS  PubMed  Google Scholar 

  50. De Schutter I, Dreesman A, Soetens O, De Waele M, Crokaert F, Verhaegen J, et al. In young children, persistent wheezing is associated with bronchial bacterial infection: a retrospective analysis. BMC Pediatr. 2012;12:83.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Pennings HJ, Kramer K, Bast A, Buurman WA, Wouters EF. Tumour necrosis factor-alpha induces hyperreactivity in tracheal smooth muscle of the guinea-pig in vitro. Eur Respir J. 1998;12:45–9.

    Article  CAS  PubMed  Google Scholar 

  52. Bevacqua BK. Pre-operative pulmonary evaluation in the patient with suspected respiratory disease. Indian J Anaesth. 2015;59:542–9.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Nguyen TT, Hoo AF, Lum S, Wade A, Thia LP, Stocks J. New reference equations to improve interpretation of infant lung function. Pediatr Pulmonol. 2013;48:370–80.

    Article  PubMed  Google Scholar 

  54. Chen J, Liu X, Du W, Srivastava R, Fu J, Zheng M, et al. Pulmonary function testing in pediatric pneumonia patients with wheezing younger than 3 years of age. Glob Pediatr Health. 2019;6(X19840357):2333794.

    Google Scholar 

Download references




This work was supported by the National Health and Family Planning Commission Fund Project of Sichuan province of China (No. 16PJ067).

Author information

Authors and Affiliations



LZ and TA conceived and designed the study. WX and JY collected the data. JY and TA analyzed the data. LZ performed the laboratory analysis. JY wrote the paper. TA, YF, and CX critically revised the paper. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Lei Zhang.

Ethics declarations

Ethics approval and consent to participate

We confirm that all methods were carried out in accordance with the ethical standards as laid down in the 1964 Declaration of Helsinki and its later amendments or comparable ethical standards. We confirm that written informed consent has been obtained from all participants from their parents and/or legal guardians.The study was approved by the ethics committees of Chengdu Women’s and Children’s Central Hospital (Ethnic audit number 2016 [22]).

Consent for publication

Not applicable.

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: Table 1.

Alpha diversity indices of all oropharyngeal samples. Fig. 1. Multiple rarefaction curves of the richness and Rank Abundance. Fig. 2 NMDS based on the Bray-Curtis distances at the OTU level at 97% identity. Fig. 3 All significantly different phyla and genera detected by LefSe.

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

Yao, J., Ai, T., Xia, W. et al. Bacterial composition and colony structure of the lower respiratory tract in infants and children with recurrent wheezing: a case–control study. Ital J Pediatr 48, 120 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Microbiome
  • Bronchoalveolar lavage
  • Children
  • Recurrent wheezing