Viral respiratory infections and the oropharyngeal bacterial microbiota in acutely wheezing children
Authors:
Leah Cuthbertson aff001; Stephen W. C. Oo aff002; Michael J. Cox aff001; Siew-Kim Khoo aff002; Des W. Cox aff002; Glenys Chidlow aff005; Kimberley Franks aff002; Franciska Prastanti aff002; Meredith L. Borland aff002; James E. Gern aff008; David W. Smith aff002; Joelene A. Bizzintino aff002; Ingrid A. Laing aff002; Peter N. Le Souëf aff002; Miriam F. Moffatt aff001; William O. C. Cookson aff001
Authors place of work:
National Heart and Lung Institute, Imperial College, London, England, United Kingdom
aff001; Division of Paediatrics, Faculty of Health and Medical Sciences, University of Western Australia, Perth, Australia
aff002; Respiratory Department, Perth Children’s Hospital, Perth, Western Australia
aff003; Telethon Kids Institute, Perth, Australia
aff004; Department of Microbiology, PathWest Laboratory Medicine WA, QEII Medical Centre, Perth, Australia
aff005; Emergency Department, Perth Children’s Hospital, Perth, Australia
aff006; Division of Emergency Medicine, Faculty of Health and Medical Sciences, University of Western Australia, Perth, Australia
aff007; Department of Pediatrics, University of Wisconsin-Madison, Madison, Wisconsin, United States of America
aff008; Medical School, Faculty of Health and Medical Sciences, University of Western Australia, Perth, Australia
aff009; Royal Brompton and Harefield NHS Foundation Trust, London, England, United Kingdom
aff010
Published in the journal:
PLoS ONE 14(10)
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pone.0223990
Summary
Acute viral wheeze in children is a major cause of hospitalisation and a major risk factor for the development of asthma. However, the role of the respiratory tract microbiome in the development of acute wheeze is unclear. To investigate whether severe wheezing episodes in children are associated with bacterial dysbiosis in the respiratory tract, oropharyngeal swabs were collected from 109 children with acute wheezing attending the only tertiary paediatric hospital in Perth, Australia. The bacterial community from these samples was explored using next generation sequencing and compared to samples from 75 non-wheezing controls. No significant difference in bacterial diversity was observed between samples from those with wheeze and healthy controls. Within the wheezing group, attendance at kindergarten or preschool was however, associated with increased bacterial diversity. Rhinovirus (RV) infection did not have a significant effect on bacterial community composition. A significant difference in bacterial richness was observed between children with RV-A and RV-C infection, however this is likely due to the differences in age group between the patient cohorts. The bacterial community within the oropharynx was found to be diverse and heterogeneous. Age and attendance at day care or kindergarten were important factors in driving bacterial diversity. However, wheeze and viral infection were not found to significantly relate to the bacterial community. Bacterial airway microbiome is highly variable in early life and its role in wheeze remains less clear than viral influences.
Keywords:
Respiratory infections – Viral transmission and infection – children – Microbiome – Sequence databases – Ribosomal RNA – Bronchiolitis – Rhinovirus infection
Introduction
Respiratory infections are a significant cause of morbidity in young children, with over half of hospitalisations in children under 5 years of age being related to respiratory disease [1]. Evidence suggests that environmental and infectious exposures in early life are a major risk factor for the development of disease [2, 3], including asthma [4–6].
Acute wheeze is the most common reason for children to present to hospital with between 80–90% of these cases being attributed to viral infection [7, 8]. While a number of respiratory viruses have been implicated in wheeze the most common are respiratory syncytial virus (RSV) and rhinovirus (RV)[9]. In infants aged less than one year RSV infections predominate, but after this RV is also viewed as an important risk factor for acute wheeze [10].
The influence of the respiratory tract bacterial community in the context of wheeze is still being established. Recent studies in this area have implicated organisms such as Moraxella catarrhalis, Haemophilus influenzae and Staphylococcus aureus in increased rates of wheeze [11, 12]. It has been suggested this may be due to immune modulation by these organisms, in infants resulting in increased risk of asthma development in later life [13].
Using 16s rRNA gene sequencing of respiratory samples from children presenting to hospital with acute wheezing, this study aimed to examine whether the bacterial community in the airways of children with acute respiratory wheeze was altered compared with that of non-wheezing children. Changes in the bacterial community were also explored to determine if acute RV infection or species had a significant effect on the airway microbiota.
Methods
Children between 0–16 years were recruited as part of the MAVRIC (Mechanisms of Acute Viral Respiratory Infection in Children) study on presentation to Princess Margaret Hospital (PMH) for Children, Perth, Western Australia between January 2004 and January 2014. Cases and controls were recruited through all seasons, and a questionnaire was administered to all subjects to determine symptoms of any current illness, including coryzal symptoms, and risk factors, including birth history, postnatal and in utero cigarette smoke exposure, daycare attendance, atopy and allergy history, diagnoses of acute illness (children were excluded if they had any chronic illness other than asthma), recent antibiotic use, systemic steroids.
Recurrence data was collected on cases included in the study from birth. Hospital presentation records were used to determine frequency of respiratory presentations both prior to and following presentation. Five patterns were determined described as “few”, “persistent”, “multiple A”, “multiple B”, and “atypical”. (See supplemental methods for further definitions).
Cases recruited had acute wheezing illness with no other co-morbid conditions besides asthma, eczema, or atopy. Controls with no pre-existing chronic disease including chronic respiratory illness were recruited from four sources: siblings and relatives of cases, PMH patients (presenting with minor injury/fractured limbs), volunteers from the local community or day care facilities. A proportion of controls had symptoms of mild acute respiratory infection but no wheeze. Blood, oropharyngeal (OP) swabs and nasal samples (wash or blow) were collected from each participant.
Several cases were followed up within 9 months of recruitment and viral, OP swab and blood samples were repeated at this time.
This study was approved by the PMH for Children, Perth, Western Australia Human ethics committee (Reference: 1761EP). At least one parent or guardian provided informed written consent for children, prior to participation in the study.
Atopy
Skin prick tests were performed to 11 allergens (cow’s milk, whole egg, cat pelt, dog dander, rye grass, mixed grasses, D. pteronyssinus, D. farinae, cockroach and Aspergillus fumigatus, Alternaria tenuis), a wheal size ≥3mm or self-reported allergic reaction to an allergen was used determined atopy.
Bloods cell counts and cathelicidin measurements
Blood cell counts were determined by PathWest Laboratory Medicine WA (PathWest) at PMH hospital pathology lab.
Cathelicidin (LL-37) was measured in plasma using a commercially available ELISA kit (Hycult Biotech).
Viral detection
A nasal blow or wash was collected at recruitment and placed immediately on ice prior to storage at -80°C. Nasal specimens were typed for RV species using methods previously described [14] using modified primers [15]. Briefly, a semi-nested RT-PCR with primers that amplify the 260-bp variable region of the 5’ untranslated region of the RV genome was completed. RV-positive samples were sequenced and assigned an RV strain type and species following sequence alignment with sequences of the 101 classic serotypes and the 53 newly assigned genotypes, using ClustalX software (University College Dublin, Dublin, Ireland) [16].
A second aliquot from each nasal specimen was tested for RSV, influenza A and B, parainfluenza 1–4, RSV, and human metapneumovirus (hMPV) at PathWest using routine diagnostic methods including PCR [17], direct or indirect fluorescent antibody testing, or immunofluorescence after cell culture as previously described [18]. Whenever possible, samples were also tested for enterovirus, coronavirus, and bocavirus.
Bacterial detection
Sterile dry rayon swabs (Copan) were used to sample the soft palate of the oropharynx. Specific care was made not to contaminate samples with any other part of the mouth. Any swab that made contact with the tongue or cheek were disposed of and a new swab was performed, for full details see supplementary materials. Samples were immediately put on ice and stored at -80°C.
DNA was extracted from samples utilizing an MPbio FastDNA SPIN Kit for soil as per manufacturers’ instructions and frozen at -80°C. Samples were transported on dry ice to Imperial College London, for bacterial analysis.
Quantitative PCR
Total bacterial burden was measured using a SYBR green quantitative PCR assay using the primers 520F, 5’- AYTGGGYDTAAAGNG and 820R, 5’-TACNVGGGTATCTAATCC, targeting the V4 region of the 16S rRNA gene as described in Cuthbertson et al 2017[19]. All reactions were performed in triplicate and included standards and non-template controls on the ViiA 7 Real-time PCR system (Life Technologies, Paisley, UK) using SYBR Fast qPCR Master mix (KAPA Biosystems, Wilmington, MA, USA). Standards were generated from near full length cloned 16S rRNA gene of Vibrio natregens. Plasmids quantified using Quantit picogreen dsDNA Assay kit (Promega, Madison, USA), samples were then serially diluted 10 fold to form standards ranging from 1 x 108–1 x 104.
16S rRNA sequencing
Community analysis was carried out using 16S rRNA gene sequencing. Custom dual barcoded fusion primers were used to target the previously quantified region of the 16S rRNA gene as previously described [19]. Each sequencing run contained a PCR negative control and a Mock community, consisting of 34 16S rRNA gene clones of known bacterial species in equal proportions. Sequencing was carried out using the Illumina MiSeq platform using the Illumina V2 2x250bp cycle kit. Sequences were submitted to the European nucleotide database, project number PRJEB32061.
Sequencing analysis
Downstream sequencing analysis was carried out using Quantitative Insights in Microbial Ecology (QIIME) Version 1.9.0 due to the dual barcoded indexs used in this study. All sequences were trimmed to 200bp and joined with a minimum of 150bp overlap, a maximum of 10% mismatch was stipulated. Sequences were then demultiplexed and any phiX reads were removed, prior to OTU picking using open reference UCLUST OTU picking [20], clustering at 97% similarity using the Silva reference database (www.arb-silva.de), reads with less that 60% id were discarded and 10% of sequences that failed to id were included for de novo clustering. Representative sequences were picked from the most abundant read in the cluster. PYNAST [21] was used to align representative sequences before running the nearest alignment space termination (NAST) algorithm [22]. ChimeraSlayer (http://microbiomeutil.sourceforge.net/) was used to identify and remove any chimeric sequences. The Ribosomal Database Project (RDP) naive Bayesian classifier was used to apply taxonomic identification using Silva 115 NR database. Finally, an OTU biom table was created for further downstream analysis.
Exact sequence variants (ESVs) may in some circumstances improve identification of microbial taxa, but genomic sequencing of the airway microbiota is at an early stage. In order to avoid over-splitting of taxa and false inflation of diversity, we have taken the conservative approach of using OTUs in our analyses rather than amplicon sequence variants.
Statistical analysis
All further analysis was carried out using R version 3.3.2 [23]. Pre-processing and primary analysis was carried out in Phyloseq [24]. Contamination was removed using Decontam [25]. A minimum threshold of 2,000 reads was applied to all, samples with less than 2,000 reads were removed from further analysis. All remaining samples were then rarefied to the minimum number of reads present in the data subset.
Non-parametric Wilcoxon sign ranked tests were used to test significant differences between means. Pearson correlations were used to test the relationship between continuous clinical variables and diversity measures. Adonis permutational ANOVA was used to investigate changes in community composition while Random forest analysis was used to identify OTU’s associated with disease.
Results
After removal of sequencing controls and those with less than 2000 reads, a total of 201 samples were taken forward for analysis, see Table 1 (further information in S1 Table). This included 109 samples from children with acute wheeze (Fig 1) and paired stable follow-up samples from 17 of these children. Control samples from 75 children without symptoms of wheeze were collected, this cohort included children from day-care and siblings of those with acute wheeze (Fig 1).
Comparison of acute wheeze and controls
Samples from 109 children with acute wheeze and 75 control children were compared (Table 1 and Fig 1). No significant difference in alpha diversity was observed between the 2 groups (Richness, P = 0.363; Shannon-Weiner, P = 0.98; Inverse Simpsons, P = 0.654). Adonis permutational ANOVA revealed a significant difference in Bray-Curtis dissimilarity (R2 = 0.016, P = 0.003), Fig 2A. However, the variation explained was only 1%, suggesting a relationship driven by a small number of samples (see S1 Fig).
A random forest model was used to predict wheeze based on OTUs in the bacterial community. The predictive power of the model for cases and controls was found to be poor (out of box (OOB) error rate = 34.24%, p = 0.37), therefore the model results were considered unreliable and was not used in further investigation. However, indicator species analysis revealed a significant association of a single Veillonella OTU with the acute group. DeSeq2 analysis indicated the same Veillonella OTU associated with cases while several Haemophilus OTUs were found to be associated with controls (S2 Fig)
Clinical variables were investigated to see if they were able to better explain bacterial community variation within this population. Wilcoxon rank sum test was used to investigate 33 binary variables across all samples (S2 Table). After Bonferroni correction for multiple testing, only the diagnosis of bronchiolitis was found to show a significant difference in bacterial richness (p = 0.006), Shannon-Weiner (p = 0.017) and inverse Simpsons (p = 0.062). Significant differences in Shannon-Weiner and inverse Simpsons were observed in those that required oxygen and those that had ever attended kindergarten. It was noted, however, that none of the controls were positive for bronchiolitis or required oxygen. Therefore, this relationship was further tested in the acute wheeze population. Kruskal-Wallis with Dunns test was used to investigate categorical variables, however no significant results were observed. We did not find any differences in the microbiota that were attributable to recorded breastfeeding. No data was available to explore the effect of delivery method in this study.
Pearson’s correlations were used to investigate relationships between 20 continuous variables and alpha diversity measures (S3 Table). Only bacterial biomass was found to have a weak significant correlation with bacterial richness. Significant changes in blood counts were not found to hold up to multiple testing and expected differences in the antimicrobial peptide cathlecidin were not found.
Adonis permutational ANOVA with 99,999 iterations was used to investigate changes in community composition with 58 variables. Only bacterial biomass, bronchiolitis diagnosis and regular attendance at day care were found to be significant (S4 Table).
Bronchiolitis diagnosis was found to show significant reduction in alpha diversity. Due to the control samples coming from a slightly older population true age matched case controls were unable to be compared (wheeze, median = 0.195 (min = 0.08 –max = 1.87), control, 1.38 (0.6–1.91)). However, when the microbiota in children with bronchiolitis were compared to those closest in age there was a significant reduction in alpha diversity (richness; p = 0.01, Shannon-Weiner, p = 0.003, inverse Simpsons, p = 0.005) that was not associated with change in bacterial biomass (p = 0.874). Changes in community composition using Bray-Curtis dissimilarity were found to explain 9.8% of the variation (p = 0.012), S3 Fig.
Acute wheeze
To explore variation within the group suffering from acute wheeze, this group was examined as a subset from the whole data set and clinical variables were explored (Fig 1). In the acute wheeze group, significant decrease in all alpha diversity measures (differences in bacterial diversity between subjects) were seen with bronchiolitis. Fourteen patients with doctor diagnosed bronchiolitis were included in the study, all of these patients were under the age of 2 years. Of those diagnosed, all patients also had wheeze. Nine of the patients were positive for RSV at the time of sampling, 5 of these were also positive for RV (Fig 1).
A significant increase was observed with attendance at kindergarten (S5 Table). Attendance at preschool showed a significant increase in bacterial richness (number of OTUs) only. Only patient age was found to show significant positive correlation with bacterial richness (S6 Table).
In the children with acute wheeze, bacterial biomass was found to explain 8.9% of the variation in beta diversity (p < 0.001). Bronchiolitis (r2 = 0.042, p < 0.001) and attendance at kindergarten (r2 = 0.052, p < 0.001) were also found to be significant in the acute wheeze dataset (S7 Table).
Within the 109 children with acute wheeze, 71 were below the age of 5. In the acute wheeze group, the under 5s bacterial biomass remained a significant driver of bacterial community composition (r2 = 0.06, p < 0.001). However, no other variables were significant in the younger cohort.
Kruskal-Wallis was used to investigate wheeze recurrence patterns for the five patterns of wheeze observed in wheeze cases (S4 Fig). No significant difference in alpha bacterial diversity was observed within or between any of the above recurrence groups.
Rhinovirus (RV)
Within the participants suffering from acute wheeze, 71 patients had RV infection. RV strains from two patients were unable to be typed. Twenty-six patients were infected with RV-A, 3 with RV-B and 40 with RV-C. The sample numbers only allowed comparison of the bacterial community from patients with RV-A with those with RV-C. No significant difference in alpha diversity was observed. In the under 5s group, 15 patients had RV-A, 1 had RV-B and 27 had RV-C. A significant difference in richness (W = 123.5, p-value = 0.039) was observed between those with RV-A and RV-C, however no differences in diversity were observed. No significant difference in community composition was observed using adonis (R2 = 0.039, p = 0.328).
Random forest modelling was unable to identify OTUs associated with viral species (p = 0.504).
Follow-up samples
Stable follow-up samples were collected from 17 individuals after an acute wheezing episode. No significant difference in bacterial biomass (V = 44, p-value = 0.132) or alpha diversity (richness, p = 0.289, Shannon, p = 0.145, Simpson, p = 0.109) was found when comparing paired samples with Wilcoxon signed rank test.
Between and within sample beta diversity was compared (Fig 2B). Significant differences between paired samples were observed with 5.8% of the variation explained by time of sampling (p = 0.013), however 57.6% of the variation was explained by individual (p = 0.013).
Discussion
In this cohort of children, we found substantial diversity and heterogeneity in OP microbiota composition regardless of wheeze status. There were no significant microbiome differences between acute wheeze cases and non-wheezing controls on alpha diversity measures. Beta diversity was examined using a Bray-Curtis dissimilarity adonis model and while significant difference (p = 0.003) was found, only 1.6% of the variance is explained by differences between cases and controls. Further interrogation with Bray-Curtis dissimilarity hierarchical clustering (S1 Fig) reveals that the correlation is likely driven by a small number of samples and no pattern was observed between those that had wheeze and those that did not (Fig 2). Overall, despite cases having an acute wheezing illness serious enough to cause presentation to a children’s hospital emergency department, there were no clear differences in the bacterial community between cases and controls.
In the small number of cases with matched acute and convalescent samples taken up to 9 months later, we were similarly unable to find differences in alpha diversity between the cases and follow-up samples. Significant differences in beta diversity were observed between cases and follow-up samples, however the variation explained was 5.8% compared to within subject changes which explained 57.6% of the variation. The heterogeneous nature of the microbial community between patients makes investigation into microbial changes a challenge and supports the need for large-scale longitudinal investigation into community composition in health and disease.
Random forest analysis has been used in a number of recent studies to investigate predictive OTUs in microbiome analysis [26, 27]. In this study, random forest was used to attempt to identify OTU predictors for wheeze. We found that the variation in the microbial community observed between these groups was high and this led to low confidence in the predictive power of the model. Random forest analysis of the bacterial community found bacterial OTUs to be poor predictors of wheeze. This analysis was also performed to determine if any OTUs were significantly associated with RV infection or between RV species. Neither model was found to have strong predictive power and estimated accuracy was low.
Diagnosis with bronchiolitis showed a significant difference in bacterial richness compared with the entire cohort, however this comprised only a small subset of patients (n = 13). When those with a diagnosis of bronchiolitis were compared to other cases that wheezed or age matched control subjects, significant differences in alpha and beta diversity were observed however this may have been driven by subjects with bronchiolitis being significantly younger than other groups of subjects. Bacterial diversity was found to be significantly reduced in these subjects with many patients dominated with Streptococcus sp. Further investigation into the streptococcal species, which is unable to be elucidated from 16S rRNA sequencing, would be important in further investigations. Longitudinal studies in young children to determine if low bacterial diversity in the respiratory tract is a risk factor for or the result of bronchiolitis would be important for future investigations.
In contrast to our findings, there are a number of studies that show significant differences in microbiome in relation to viral infections compared to healthy subjects [11, 28–30]. These studies show haemophilus, streptococcus, moraxella associated with acute respiratory tract infections and staphylococcus, alloiococcus and corynebacterium in healthy samples. A number of factors account for the differences we find; firstly, they analyse nasopharyngeal (NP) samples compared to the OP samples analysed in our study. A number of studies have shown that these areas show significantly different bacterial communities, and that OP swabs are a better proxy for the lower respiratory tract microbiome. There are also age differences between other cohorts and ours.
Lastly analytical methods to look at the microbiome differ between studies. Analytical methods for microbiota are slowly being standardised, and much of the analytic techniques are borrowed from ecological analysis to try and accurately portray biodiversity and differences. In order to find differences, earlier studies have defaulted to simplistic groupings based on dominant organisms. However, we have observed that many samples are diverse and without a singular dominant organism. For example, Man et al, comparing those with lower respiratory tract infection and healthy controls, found no difference in alpha diversity, and using beta diversity measures found, similar to us, very small but significant differences between the groups (r2 = 0.0031, p = <0.0001) [28]. To find significant differences in species, previous studies grouped samples based on perceived dominance and not necessarily by a dominance measure or beta diversity based on hierarchical clustering.
The greater variation in age within our population compared to other studies may have influenced our ability to find significant differences. Bacterial biomass was found to explain much of the variation in beta diversity within the acute population. Bacterial biomass was also found to show a weak positive correlation with age. Age was found to be a key variable in the analysis of this dataset and this relationship has been demonstrated in previous studies [31]. Age is interrelated with developmental changes in blood cell counts and immune function [32]. With these considerations, we considered changes in the microbiome while controlling for age.
A strength of our study is that we extracted blood from patients to examine differences in blood counts and antimicrobial peptides. However, variation in age in the study resulted in small and non-significant trends in blood cell counts. Cathelicidin was also examined as a candidate innate immunity marker as it has both anti-viral and anti-bacterial activity and lower levels have been shown to be associated with worsened severity in bronchiolitis [33]. However, we were unable to find a significant change in the bacterial community with changes in cathelicidin.
Seasonal influences in the same study location (Perth, Australia) have shown to have limited effects on microbiome in comprehensive longitudinal studies at least in infants [29]. Risk factors that are associated with hospitalisation as a result of respiratory infections in early life include maternal smoking during pregnancy, season of birth, delivery mode and gestational age[34]. When considered in this cohort of mostly hospitalised wheezing children (96%), these factors were not found to be significantly associated with bacterial or viral infection suggesting more work is required to fully understand any potential relationship.
At least one other study exists comparing OP microbiome samples between healthy adult patients and severe asthmatics showing no significant diversity differences [35]. An abundance of moraxella is found in prior studies in NP. Hilty et al demonstrate that Moraxella (proteobacteria) exists in oropharynx and may account for differences [36]. OP samples were specifically chosen as they showed greater representation to lower airway samples at least in stable patients [36, 37]. Whether acute OP swabs reflect lower airway samples in acute wheeze or acute viral episodes is still unclear but, interestingly, Man et al did show that intensive care patients had reasonable similarity between NP samples and lower respiratory tracheal samples (Bray Curtis similarity p = 0.61), but did find key differences, specifically that staphylococcus, Corynebacterium, and Dolosigranulum sp. were almost exclusively present in NP samples and absent from endotracheal aspirates [28].
The bacterial community in the paediatric population is diverse and heterogeneous. The wide range of clinical factors tested did not fully explain the wide variation in the bacterial community in these subjects. Age had a significant influence on both the bacterial community and blood cell counts in this study and the wide age range within the study population helped explain some of the variation observed.
The simplest interpretation of our results is that acute wheezing illnesses are driven by viral infections and that these infections have little influence on the bacterial community during the acute phase of the illness. However, the heterogeneous nature of the subjects made it difficult to test for significant associations between clinical variables and the oropharyngeal bacterial community. Prospective longitudinal investigation of children pre, during and post viral infection may help identify if the bacterial community is either protective or a risk factor for viral infection and respiratory wheeze.
Supporting information
S1 Table [docx]
Participant demographics continued.
S2 Table [docx]
Comparison of categorical clinical variables with Wilcoxon rank sum test in the complete cohort to alpha diversity measures, richness, Shannon-Weiner and inverse Simpsons.
S3 Table [docx]
Pearson’s correlations of continuous clinical variables and alpha diversity measures, richness, Shannon-Weiner and inverse Simpsons, in the complete cohort.
S4 Table [docx]
Results of Bray-Curtis adonis permutational ANOVA examining clinical variables with 99,999 iterations.
S5 Table [docx]
Comparison of categorical clinical variables in those with acute wheeze to alpha diversity measures, richness, Shannon-Weiner and inverse Simpsons.
S6 Table [docx]
Pearson’s correlations of continuous clinical variables in those with acute wheeze and alpha diversity measures.
S7 Table [docx]
Results of Bray-Curtis adonis permutational ANOVA examining clinical variables from individuals with acute wheeze with 99,999 iterations.
S1 Fig [tiff]
Stacked bar plot comparing OP samples from children with acute wheeze to healthy controls.
S2 Fig [pdf]
Volcano plot showing differential abundance of OTUs based on results from DeSeq2 analysis.
S3 Fig [tiff]
Stacked bar plot comparing OP samples from children with bronchiolitis and acute wheeze to healthy controls.
S4 Fig [tiff]
Boxplots showing diversity measures associated with recurrence data from acute wheeze cases.
S1 File [docx]
Supplementary methods.
Zdroje
1. Goto T, Tsugawa Y, Mansbach JM, Camargo CA Jr, Hasegawa K. Trends in Infectious Disease Hospitalizations in US Children, 2000 to 2012. Pediatr Infect Dis J. 2016;35(6):e158–63. Epub 2016/03/12. doi: 10.1097/INF.0000000000001134 26967815; PubMed Central PMCID: PMC4912127.
2. Bisgaard H, Hermansen MN, Buchvald F, Loland L, Halkjaer LB, Bonnelykke K, et al. Childhood asthma after bacterial colonization of the airway in neonates. N Engl J Med. 2007;357(15):1487–95. Epub 2007/10/12. doi: 10.1056/NEJMoa052632 17928596.
3. Ege MJ, Mayer M, Normand AC, Genuneit J, Cookson WO, Braun-Fahrlander C, et al. Exposure to environmental microorganisms and childhood asthma. N Engl J Med. 2011;364(8):701–9. Epub 2011/02/25. doi: 10.1056/NEJMoa1007302 21345099.
4. Castro-Rodriguez JA, Forno E, Rodriguez-Martinez CE, Celedon JC. Risk and Protective Factors for Childhood Asthma: What Is the Evidence? J Allergy Clin Immunol Pract. 2016;4(6):1111–22. Epub 2016/06/12. doi: 10.1016/j.jaip.2016.05.003 27286779; PubMed Central PMCID: PMC5107168.
5. Stern DA, Morgan WJ, Halonen M, Wright AL, Martinez FD. Wheezing and bronchial hyper-responsiveness in early childhood as predictors of newly diagnosed asthma in early adulthood: a longitudinal birth-cohort study. Lancet. 2008;372(9643):1058–64. Epub 2008/09/23. doi: 10.1016/S0140-6736(08)61447-6 18805334; PubMed Central PMCID: PMC2831297.
6. Mirzakhani H, Carey VJ, Zeiger R, Bacharier LB, O'Connor GT, Schatz MX, et al. Impact of parental asthma, prenatal maternal asthma control, and vitamin D status on risk of asthma and recurrent wheeze in 3-year-old children. Clin Exp Allergy. 2018. Epub 2018/11/22. doi: 10.1111/cea.13320 30461089.
7. Bizzintino J, Lee WM, Laing IA, Vang F, Pappas T, Zhang G, et al. Association between human rhinovirus C and severity of acute asthma in children. Eur Respir J. 2011;37(5):1037–42. Epub 2010/08/10. doi: 10.1183/09031936.00092410 20693244; PubMed Central PMCID: PMC3024467.
8. Bacharier LB. Viral-induced wheezing episodes in preschool children: approaches to therapy. Curr Opin Pulm Med. 2010;16(1):31–5. Epub 2009/10/10. doi: 10.1097/MCP.0b013e32833303e6 19816178.
9. Rossi GA, Colin AA. Infantile respiratory syncytial virus and human rhinovirus infections: respective role in inception and persistence of wheezing. Eur Respir J. 2015;45(3):774–89. Epub 2014/11/02. doi: 10.1183/09031936.00062714 25359340.
10. Miller EK, Bugna J, Libster R, Shepherd BE, Scalzo PM, Acosta PL, et al. Human rhinoviruses in severe respiratory disease in very low birth weight infants. Pediatrics. 2012;129(1):e60–7. Epub 2011/12/28. doi: 10.1542/peds.2011-0583 22201153; PubMed Central PMCID: PMC3255465.
11. Biesbroek G, Tsivtsivadze E, Sanders EA, Montijn R, Veenhoven RH, Keijser BJ, et al. Early respiratory microbiota composition determines bacterial succession patterns and respiratory health in children. Am J Respir Crit Care Med. 2014;190(11):1283–92. Epub 2014/10/21. doi: 10.1164/rccm.201407-1240OC 25329446.
12. Bosch A, de Steenhuijsen Piters WAA, van Houten MA, Chu M, Biesbroek G, Kool J, et al. Maturation of the Infant Respiratory Microbiota, Environmental Drivers, and Health Consequences. A Prospective Cohort Study. Am J Respir Crit Care Med. 2017;196(12):1582–90. Epub 2017/07/01. doi: 10.1164/rccm.201703-0554OC 28665684.
13. Folsgaard NV, Schjorring S, Chawes BL, Rasmussen MA, Krogfelt KA, Brix S, et al. Pathogenic bacteria colonizing the airways in asymptomatic neonates stimulates topical inflammatory mediator release. Am J Respir Crit Care Med. 2013;187(6):589–95. Epub 2013/02/02. doi: 10.1164/rccm.201207-1297OC 23370914.
14. Lee WM, Grindle K, Pappas T, Marshall DJ, Moser MJ, Beaty EL, et al. High-throughput, sensitive, and accurate multiplex PCR-microsphere flow cytometry system for large-scale comprehensive detection of respiratory viruses. J Clin Microbiol. 2007;45(8):2626–34. Epub 2007/06/01. doi: 10.1128/JCM.02501-06 17537928; PubMed Central PMCID: PMC1951217.
15. Bochkov YA, Grindle K, Vang F, Evans MD, Gern JE. Improved molecular typing assay for rhinovirus species A, B, and C. J Clin Microbiol. 2014;52(7):2461–71. Epub 2014/05/03. doi: 10.1128/JCM.00075-14 24789198; PubMed Central PMCID: PMC4097758.
16. Lee WM, Kiesner C, Pappas T, Lee I, Grindle K, Jartti T, et al. A diverse group of previously unrecognized human rhinoviruses are common causes of respiratory illnesses in infants. Plos One. 2007;2(10):e966. Epub 2007/10/04. doi: 10.1371/journal.pone.0000966 17912345; PubMed Central PMCID: PMC1989136.
17. Chidlow GR, Harnett GB, Shellam GR, Smith DW. An economical tandem multiplex real-time PCR technique for the detection of a comprehensive range of respiratory pathogens. Viruses. 2009;1(1):42–56. Epub 2009/06/01. doi: 10.3390/v1010042 21994537; PubMed Central PMCID: PMC3185464.
18. Subrata LS, Bizzintino J, Mamessier E, Bosco A, McKenna KL, Wikstrom ME, et al. Interactions between innate antiviral and atopic immunoinflammatory pathways precipitate and sustain asthma exacerbations in children. J Immunol. 2009;183(4):2793–800. Epub 2009/07/22. doi: 10.4049/jimmunol.0900695 19620293.
19. Cuthbertson L, Craven V, Bingle L, Cookson W, Everard ML, Moffatt MF. The impact of persistent bacterial bronchitis on the pulmonary microbiome of children. Plos One. 2017;12(12):e0190075. Epub 2017/12/28. doi: 10.1371/journal.pone.0190075 29281698; PubMed Central PMCID: PMC5744971.
20. Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26(19):2460–1. Epub 2010/08/17. doi: 10.1093/bioinformatics/btq461 20709691.
21. Caporaso JG, Bittinger K, Bushman FD, DeSantis TZ, Andersen GL, Knight R. PyNAST: a flexible tool for aligning sequences to a template alignment. Bioinformatics. 2010;26(2):266–7. Epub 2009/11/17. doi: 10.1093/bioinformatics/btp636 19914921; PubMed Central PMCID: PMC2804299.
22. DeSantis TZ Jr, Hugenholtz P, Keller K, Brodie EL, Larsen N, Piceno YM, et al. NAST: a multiple sequence alignment server for comparative analysis of 16S rRNA genes. Nucleic Acids Res. 2006;34(Web Server issue):W394-9. Epub 2006/07/18. doi: 10.1093/nar/gkl244 16845035; PubMed Central PMCID: PMC1538769.
23. Team RDC. R: A Language and Environment for Statistical Computing. R Found Stat Comput. 2015.
24. McMurdie PJ, Holmes S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. Plos One. 2013;8(4):e61217. Epub 2013/05/01. doi: 10.1371/journal.pone.0061217 23630581; PubMed Central PMCID: PMC3632530.
25. Davis NM, Proctor DM, Holmes SP, Relman DA, Callahan BJ. Simple statistical identification and removal of contaminant sequences in marker-gene and metagenomics data. Microbiome. 2018;6(1):226. Epub 2018/12/19. doi: 10.1186/s40168-018-0605-2 30558668; PubMed Central PMCID: PMC6298009.
26. Belk A, Xu ZZ, Carter DO, Lynne A, Bucheli S, Knight R, et al. Microbiome Data Accurately Predicts the Postmortem Interval Using Random Forest Regression Models. Genes (Basel). 2018;9(2). Epub 2018/02/22. doi: 10.3390/genes9020104 29462950; PubMed Central PMCID: PMC5852600.
27. Dahl C, Stanislawski M, Iszatt N, Mandal S, Lozupone C, Clemente JC, et al. Gut microbiome of mothers delivering prematurely shows reduced diversity and lower relative abundance of Bifidobacterium and Streptococcus. Plos One. 2017;12(10):e0184336. Epub 2017/10/27. doi: 10.1371/journal.pone.0184336 29069100; PubMed Central PMCID: PMC5656300.
28. Man WH, van Houten MA, Merelle ME, Vlieger AM, Chu M, Jansen NJG, et al. Bacterial and viral respiratory tract microbiota and host characteristics in children with lower respiratory tract infections: a matched case-control study. Lancet Respir Med. 2019. Epub 2019/03/20. doi: 10.1016/S2213-2600(18)30449-1 30885620.
29. 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(5):704–15. Epub 2015/04/14. doi: 10.1016/j.chom.2015.03.008 25865368; PubMed Central PMCID: PMC4433433.
30. Bogaert D, Keijser B, Huse S, Rossen J, Veenhoven R, van Gils E, et al. Variability and diversity of nasopharyngeal microbiota in children: a metagenomic analysis. Plos One. 2011;6(2):e17035. Epub 2011/03/10. doi: 10.1371/journal.pone.0017035 21386965; PubMed Central PMCID: PMC3046172.
31. Mika M, Mack I, Korten I, Qi W, Aebi S, Frey U, et al. Dynamics of the nasal microbiota in infancy: a prospective cohort study. J Allergy Clin Immunol. 2015;135(4):905–12 e11. Epub 2015/02/01. doi: 10.1016/j.jaci.2014.12.1909 25636948.
32. Valiathan R, Ashman M, Asthana D. Effects of Ageing on the Immune System: Infants to Elderly. Scand J Immunol. 2016;83(4):255–66. Epub 2016/01/26. doi: 10.1111/sji.12413 26808160.
33. Mansbach JM, Piedra PA, Borregaard N, Martineau AR, Neuman MI, Espinola JA, et al. Serum cathelicidin level is associated with viral etiology and severity of bronchiolitis. J Allergy Clin Immunol. 2012;130(4):1007–8 e1. Epub 2012/09/05. doi: 10.1016/j.jaci.2012.07.044 22944482; PubMed Central PMCID: PMC3462235.
34. Moore HC, Hall GL, de Klerk N. Infant respiratory infections and later respiratory hospitalisation in childhood. Eur Respir J. 2015;46(5):1334–41. Epub 2015/08/22. doi: 10.1183/13993003.00587-2015 26293501.
35. Brusselle G, Canvin J, Weiss S, Sun SX, Buhl R. Stratification of eosinophilic asthma patients treated with reslizumab and GINA Step 4 or 5 therapy. ERJ Open Res. 2017;3(3). Epub 2017/08/29. doi: 10.1183/23120541.00004–2017 28845430; PubMed Central PMCID: PMC5570512 openres.ersjournals.com.
36. Hilty M, Burke C, Pedro H, Cardenas P, Bush A, Bossley C, et al. Disordered microbial communities in asthmatic airways. Plos One. 2010;5(1):e8578. Epub 2010/01/07. doi: 10.1371/journal.pone.0008578 20052417; PubMed Central PMCID: PMC2798952.
37. 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(8):957–63. Epub 2011/06/18. doi: 10.1164/rccm.201104-0655OC 21680950; PubMed Central PMCID: PMC3208663.
Článek vyšel v časopise
PLOS One
2019 Číslo 10
- Řekněte mi, co jíte, a já vám řeknu, jaká je úroveň vaší kognice
- Původní český test ALBA pro odhalení mírných kognitivních deficitů získal certifikaci MZ ČR
- Duševní onemocnění se u teenagerů přenáší v rámci školní třídy
- Od kostiček k lidské kůži. Jak si postavit 3D biotiskárnu?
- Prof. Libor Grubhoffer: Nová varianta mpox je příležitostí pokročit v poznání původce varioly
Nejčtenější v tomto čísle
- Correction: Low dose naltrexone: Effects on medication in rheumatoid and seropositive arthritis. A nationwide register-based controlled quasi-experimental before-after study
- Combining CDK4/6 inhibitors ribociclib and palbociclib with cytotoxic agents does not enhance cytotoxicity
- Risk factors associated with IgA vasculitis with nephritis (Henoch–Schönlein purpura nephritis) progressing to unfavorable outcomes: A meta-analysis
- Prevalence of pectus excavatum (PE), pectus carinatum (PC), tracheal hypoplasia, thoracic spine deformities and lateral heart displacement in thoracic radiographs of screw-tailed brachycephalic dogs