key: cord-0715888-kdx0kier
authors: Mansbach, Jonathan M.; Luna, Pamela N.; Shaw, Chad A.; Hasegawa, Kohei; Petrosino, Joseph F.; Piedra, Pedro A.; Sullivan, Ashley F.; Espinola, Janice A.; Stewart, Christopher J.; Camargo, Carlos A.
title: Increased Moraxella and Streptococcus species abundance after severe bronchiolitis is associated with recurrent wheezing
date: 2019-11-15
journal: J Allergy Clin Immunol
DOI: 10.1016/j.jaci.2019.10.034
sha: 3bcd82d8c1991aa3b1448ed237df0c427d672685
doc_id: 715888
cord_uid: kdx0kier

BACKGROUND: The role of the airway microbiome in the development of recurrent wheezing and asthma remains uncertain, particularly in the high-risk group of infants hospitalized for bronchiolitis. OBJECTIVE: We sought to examine the relation of the nasal microbiota at bronchiolitis-related hospitalization and 3 later points to the risk of recurrent wheezing by age 3 years. METHODS: In 17 US centers researchers collected clinical data and nasal swabs from infants hospitalized for bronchiolitis. Trained parents collected nasal swabs 3 weeks after hospitalization and, when healthy, during the summer and 1 year after hospitalization. We applied 16S rRNA gene sequencing to all nasal swabs. We used joint modeling to examine the relation of longitudinal nasal microbiota abundances to the risk of recurrent wheezing. RESULTS: Among 842 infants hospitalized for bronchiolitis, there was 88% follow-up at 3 years, and 31% had recurrent wheezing. The median age at enrollment was 3.2 months (interquartile range, 1.7-5.8 months). In joint modeling analyses adjusting for 16 covariates, including viral cause, a 10% increase in relative abundance of Moraxella or Streptococcus species 3 weeks after day 1 of hospitalization was associated with an increased risk of recurrent wheezing (hazard ratio [HR] of 1.38 and 95% high-density interval [HDI] of 1.11-1.85 and HR of 1.76 and 95% HDI of 1.13-3.19, respectively). Increased Streptococcus species abundance the summer after hospitalization was also associated with a greater risk of recurrent wheezing (HR, 1.76; 95% HDI, 1.15-3.27). CONCLUSIONS: Enrichment of Moraxella or Streptococcus species after bronchiolitis hospitalization was associated with recurrent wheezing by age 3 years, possibly providing new avenues to ameliorate the long-term respiratory outcomes of infants with severe bronchiolitis.

Increased Moraxella and Streptococcus species abundance after severe bronchiolitis is associated with recurrent wheezing Bronchiolitis is an acute respiratory tract infection and the leading cause of hospitalization for US infants, accounting for approximately 130,000 hospitalizations annually. 1 In addition to this acute severity of illness, these hospitalized infants are 3 to 4 times more likely than healthy control subjects to have recurrent wheezing and childhood asthma. 2 Unfortunately, the infants who will have this long-term respiratory morbidity after a bronchiolitis hospitalization remain incompletely defined.

To date, most longitudinal bronchiolitis cohort studies have used 1 of the 2 most common viral causes of bronchiolitis, respiratory syncytial virus (RSV) and rhinovirus, as the primary exposure. 2 However, respiratory tract viruses infect infants in airways colonized with highly functional bacteria. 3, 4 Emerging evidence suggests that the bacterial composition of the airway is associated not only with acute severity outcomes 5 but also longterm respiratory outcomes. 6 More than a decade ago, Bisgaard et al, 7 using conventional culture, reported that hypopharyngeal colonization by Moraxella, Haemophilus, or Streptococcus species in asymptomatic 1-month-old infants was associated with recurrent wheezing and childhood asthma. More recently, cross-sectional, culture-independent studies have shown that when compared with healthy control subjects, children and adults with asthma have airway microbiota compositions enriched with Proteobacteria (eg, Moraxella and Haemophilus species). 4, [8] [9] [10] [11] Furthermore, longitudinal microbiota studies in the first 2 years of life have found that Moraxella, Haemophilus, and Streptococcus species dominate nasopharyngeal communities during acute respiratory tract infection 6, 12 and in sensitized children are associated with chronic wheeze at age 5 years. 6 Abbreviations used HDI: High-density interval HR: Hazard ratio MARC-35: 35th Multicenter Airway Research Collaboration NPA: Nasopharyngeal aspirate RSV: Respiratory syncytial virus

To date, however, no study has examined in a multicenter cohort of infants hospitalized for bronchiolitis the association between the longitudinal composition of the nasal microbiota and risk of recurrent wheezing of childhood. Indeed, this is the first longitudinal microbiota analysis of the present cohort. 5, 13 We hypothesized that a posthospitalization increase in Moraxella species relative abundance is associated with a greater risk of recurrent wheezing by age 3 years.

As previously described, 5 site teams at 17 hospitals in 14 US states enrolled infants (age <1 year) hospitalized for bronchiolitis into the 35th Multicenter Airway Research Collaboration (MARC-35; http://www.emnet-usa.org/ Marc_35/M35.htm). Attending physicians diagnosed an infant's bronchiolitis based on the American Academy of Pediatrics definition. 14 Enrollment occurred over 3 consecutive winter seasons (ie, November-April) from 2011-2014. The institutional review board at the 17 participating hospitals approved this study, and all families signed an informed consent form before participation.

Research teams at each of the sites not only extracted clinical data from the emergency department and inpatient charts but also conducted structured interviews with parents/legal guardians during hospitalization for demographic, historical, and environmental information.

All site researchers collected nasopharyngeal aspirates (NPAs) within 24 hours of hospitalization by using the same equipment (eg, sample traps and suction catheters from Medline Industries [Mundelein, Ill] ) and a standardized protocol, as previously described. 15 Nasal swab collection Within 24 hours of hospitalization, site teams collected an index nasal swab from both anterior nares by using a single nylon pediatric FLOQSwab (Copan, Brescia, Italy). 15 This swab was collected either before the NPA collection or more than 2 hours after the NPA collection. During hospitalization, site teams taught parents how to collect the nasal swab sample. When possible, during the index sample collection, the site team swabbed 1 nare and had the parent practice their technique by swabbing the other nare. Three weeks after the date of hospitalization, using the same sample collection methods, trained parents collected the ''clearance swab.'' The summer (June, July, or August) after the hospitalization, when the children were healthy, the parents collected the ''summer swab.'' The year after the hospitalization, children were randomly assigned to a specific season for their healthy swab collection. Therefore one fourth of the parents were assigned to collect the ''seasonal swab'' over summer (June, July, and August), fall (September, October, and November), winter (December, January, and February), or spring (March, April, and May). All collected swabs were processed, as previously described. 15 16S rRNA gene sequencing As previously described, 15 we sequenced the 16S rRNA gene V4 region of the nasal swab bacteria on the Illumina MiSeq platform (Illumina, San Diego, Calif). We used microbiota data with sufficient sequence depth (ie, > _1000 reads per sample) for all analyses.

Quantitative RT-PCR assays By using the NPA samples from hospitalization, quantitative RT-PCR assays were conducted at Baylor College of Medicine, as previously described, for 17 viruses: RSV (types A and B), rhinovirus, influenza (types A and B and 2009 novel H1N1), human metapneumovirus, human bocavirus, parainfluenza viruses (types 1, 2, and 3), adenovirus, enterovirus, and coronaviruses (NL-63, OC-43, HKU1, and 229E). 16 

Based on biannual parent telephone interviews, the primary outcome for the present analysis was recurrent wheezing by age 3 years, as defined by the 2007 National Institutes of Health's Expert Panel Report 3: having at least 2 corticosteroid-requiring exacerbations in 6 months or having at least 4 wheezing episodes in 1 year that last at least 1 day and affect sleep. 17 There is no required or specified time frame between episodes of wheezing. Because wheezing in early childhood can be transient, 18 we extended this outcome to children who had recurrent wheezing by age 3 years and also had asthma at age 4 years. We defined asthma at age 4 years using an epidemiologic definition: having both a physician's diagnosis of asthma and either asthma medication use (eg, albuterol inhaler, and inhaled corticosteroids) or asthma-related symptoms (eg, wheezing and nocturnal cough) between the ages of 3 and 4 years. 19 

Analyses only included subjects with an index microbiota sample of sufficient quality. Microbiota samples collected after the onset of recurrent wheezing were excluded to maintain consistency between cross-sectional and longitudinal time-to-event analyses. Covariates for patients' characteristics and clinical presentation collected at the index visit were compared between children with and without recurrent wheezing by using x 2 or Wilcoxon rank sum tests, as appropriate.

The unrarefied 16S rRNA gene sequence read counts were combined at the genus level. We defined the top 10 genera as those with the greatest mean relative abundance across all samples. Kruskal-Wallis rank sum tests assessed the significance of differences in age, a-diversity (Shannon index), and relative abundances of the top genera across all 4 time points, with P values for the top genera comparisons determined by using the Bonferroni correction for multiple testing.

To examine the relation of genus abundance to recurrent wheezing, we first performed cross-sectional analyses on the data. We calculated the overall mean microbial abundances for each genus of interest using the relative abundance across all time points. Based on microbial abundances at each of the 4 time points, we categorized the nasal samples into 2 groups: those greater than and those less than the mean. x 2 Tests determined the significance of the association between these groups and recurrent wheezing. Multiple comparison correction was not used in these analyses because they served as exploratory analyses for the joint modeling.

Next, we performed longitudinal analyses by constructing joint models for the longitudinal microbial data to quantify their relationship with the onset of recurrent wheezing. 20 Specifically, we created a longitudinal mixed-effects generalized linear model to predict relative genus abundances over time using the 8 most abundant taxa and clinical covariates as the fixed effects, time as the random slope, and subject as the random intercept. These relative abundances and their interaction by time point were both scaled such that a 10% change in abundance was equivalent to a unit change and incorporated as time-dependent covariates into a Cox proportional hazards model with the outcome recurrent wheezing by age 3 years. We used Bayesian methods through a modified version of the rstanarm R package, assuming scaled default prior distributions. This Bayesian approach produces conservative model estimates without the multiple comparison problems found in frequentist approaches. 21 Significance was determined by using the 95% high-density interval (HDI) of the posterior distribution of the log hazard ratio (HR) within a region of practical equivalence around zero.

To confirm joint model results, we generated Kaplan-Meier curves stratified by groups of participants identified by the joint model as being at greater risk for recurrent wheezing and used log-rank tests to determine statistical significance. We also conducted a joint model sensitivity analysis restricting the outcome to children who had recurrent wheezing by age 3 years and also had asthma at age 4 years.

For each of the high-risk groups identified by using the joint model, we generated locally weighted scatterplot smoothing (ie, loess) curves with 95% CIs for the relative abundances for the time after day 1 of hospitalization. The process of creating smoothed curves for each of the risk groups was repeated to include 8 of the top 10 genera; we removed 2 genera that showed imperceptible changes over time. To further improve clarity, we used a rug plot of sample times (instead of 95% CIs) to indicate the distribution of the observations.

Overlap between participant memberships in the risk groups identified by using the joint model was assessed with a Venn diagram. After removing group overlap, patients' characteristics and clinical presentation were compared between the risk groups by using Fisher exact and Kruskal-Wallis tests, as appropriate.

The analysis used R version 3.3.3 software. All P values were 2-tailed, with P values of less than .05 were considered statistically significant.

Of 921 infants in the longitudinal cohort, 842 had sufficient sequence depth (ie, > _1000 reads per sample) and were included in the present analysis. The analytic and nonanalytic cohorts differed in antibiotic use before hospitalization (P 5 .02), with a greater percentage of infants in the nonanalytic cohort having used antibiotics (see Table E1 in this article's Online Repository at www.jacionline.org). The analytic cohort had a median age at hospitalization of 3.2 months (interquartile range, 1.7-5.8 months), 60% were male, and 80% had no history of wheezing. We characterized the bacterial microbiota using 2086 nasal swabs (842 index, 599 clearance, 379 summer, and 266 seasonal swabs) collected during and after hospitalization for bronchiolitis. We had 88% follow-up at age 3 years. Overall, 265 (31%) infants had recurrent wheezing by age 3 years. Table I summarizes characteristics among children who did and did not have recurrent wheezing. In general, infants who had recurrent wheezing were older on study entry and had a mother with asthma, a personal history of eczema, and RSV-negative bronchiolitis. (14) 45 (17) .35

All data are presented as numbers (percentages), unless otherwise indicated. P values were computed by using x 2 tests (categorical variables) and Wilcoxon rank sum tests (continuous variables). IQR, Interquartile range. *Detected through NPA samples collected at the index visit. Defined as detection of any positive values for serum allergen-specific IgE at the index visit. àDefined as admission to intensive care unit and/or use of mechanical ventilation (continuous positive airway pressure and/or intubation during inpatient stay, regardless of location) at any time during the index hospitalization.

The median age at the time of nasal swab collection, Shannon index scores, and abundances of top 10 nasal microbiota genera by swab collection time point are shown in Table E2 in this article's Online Repository at www.jacionline.org. Nasal swab collections continued until a median age of 16 months. The a-diversity was lowest at the clearance time point and greatest at the seasonal time point. Among the genera of interest, Moraxella, Haemophilus, and Streptococcus species were all in the top 10 genera by relative abundance.

Relative abundance of Moraxella, Haemophilus, and Streptococcus species by time point and risk of recurrent wheezing by age 3 years

In cross-sectional analyses, on hospitalization, we found that the relative abundances of risk bacteria (ie, Moraxella, Haemophilus, and Streptococcus species) were not associated with recurrent wheezing by age 3 years. By contrast, at the clearance time point (ie, 3 weeks after day 1 of hospitalization), infants with a relative abundance of Moraxella species greater than the mean had significantly greater risk of recurrent wheezing by age 3 years than infants with relative abundances of less than the mean (P 5 .04; Fig 1, A) . In parallel analyses the association of Haemophilus species abundance with recurrent wheezing did not reach significance at any of the time points (see Fig E1 in this article's Online Repository at www.jacionline.org). However, infants with Streptococcus species abundance of greater than the mean at either the clearance or summer time points had significantly greater risk of recurrent wheezing by age 3 years when compared with infants who had less than mean Streptococcus species abundance (P 5 .01 and P 5 .03, respectively; Fig 1, B) .

Longitudinal Moraxella, Haemophilus, and Streptococcus species abundance after hospitalization and hazard of having recurrent wheezing by age 3 years

The aforementioned cross-sectional results were also shown in our longitudinal analysis by using joint modeling. In the joint model increased Moraxella species relative abundance at the clearance time point was associated with development of recurrent wheezing by age 3 years after adjusting for 16 confounders, such as maternal asthma, birth mode, history of eczema, sensitization, viral cause, and intensive care use (HR, 1.38; 95% HDI, 1.11-1.85; Fig 2) . Thus for every 10% increase in the abundance of Moraxella species at the clearance time point, there was a 38% increase in the risk of recurrent wheezing by age 3 years. Similar to our cross-sectional analyses, Haemophilus species abundance did not reach significance in the joint model analyses (see Fig E1) . By contrast, we found that increased Streptococcus species abundance at the clearance time point and the summer after hospitalization were both associated with recurrent wheezing by age 3 years (HR of 1.76 and 95% HDI of 1.13-3.19 and HR of 1.76 and 95% HDI of 1.15-3.27, respectively; Fig 2) .

In sensitivity analyses we created Kaplan-Meier curves stratified by the 3 risk groups (ie, increased Moraxella species abundance at the clearance time point, increased Streptococcus species abundance at the clearance time point, and increased Streptococcus species abundance at the summer time point) identified in our cross-sectional and joint modeling analyses (Fig 3) . In these survival analyses we again demonstrated the increased hazard of recurrent wheezing by age 3 years for infants with high Moraxella or Streptococcus species abundance at the clearance time point (both P 5 .01) or high Streptococcus species abundance at the summer time point (P < .01). Longitudinal Moraxella and Streptococcus species abundance after hospitalization among children with recurrent wheezing by age 3 years who also had asthma at age 4 years

To address the potential heterogeneity of recurrent wheezing, we examined children who not only had recurrent wheezing by age 3 years but also had asthma at age 4 years. Supportive of the recurrent wheezing results, we found in joint models that children with greater Moraxella and Streptococcus species abundance at the clearance time point had an increased risk of recurrent wheezing accompanied by asthma at age 4 years (HR of 1.75 and 95% HDI of 1.13-3.18 and HR of 2.84 and 95% HD of 1.07-9.03, respectively; Fig 4) . Children with increased Streptococcus species abundance at the summer time point had a borderline significant increased risk of recurrent wheeze with asthma at age 4 years (HR, 2.51; 95% HDI, 0.96-8.43; Fig 4) .

The 3 nasal microbiota risk groups, increased Moraxella species abundance at the clearance time point, increased Streptococcus species abundance at the clearance time point, and increased Streptococcus species abundance at the summer time point, have minimal overlap (see Fig E2 in this article's Online Repository at www.jacionline.org). Comparing the characteristics and clinical presentation of infants within each of the 3 risk groups after excluding infants in more than 1 risk group, we found that children with increased Moraxella or Streptococcus species abundance at the clearance time point were wheezing on hospitalization more than children with increased Streptococcus species abundance at the summer time point. Otherwise, the 3 risk groups had similar demographics and clinical presentations (see Table E3 in this article's Online Repository at www.jacionline.org). We also found that the increased abundances found in each of the 3 risk groups were not present at the subsequent sample time point Examining whether the transient nature of the increased abundances could be due to antibiotic administration, we found that most patients did not receive antibiotics either before or during the hospitalization in accordance with the current bronchiolitis practice guidelines. 14 Moreover, there were no significant differences in antibiotic administration between patients with taxa abundances greater than or less than the mean (see Table E4 in this article's Online Repository at www.jacionline. org).

In this prospective multicenter cohort of 842 infants hospitalized for bronchiolitis, we identified in cross-sectional, survival, and joint modeling analyses that having increased relative abundance of nasal Moraxella or Streptococcus species 3 weeks after day 1 of hospitalization or increased nasal Streptococcus species abundance the summer after hospitalization was associated with recurrent wheezing by age 3 years. Our data corroborate the relation of infant Moraxella and Streptococcus species abundance to childhood wheezing outcomes 6, 12 and extend this prior research by focusing not only on a high-risk cohort of infants hospitalized for bronchiolitis but also by analyzing hospitalization and posthospitalization samples using robust longitudinal analytic techniques (ie, joint modeling). Although replication is needed and causality has not been established, we might have identified, for infants hospitalized with bronchiolitis, a short window of opportunity (ie, the weeks after hospitalization) to influence the secondary succession of the nasal microbial ecosystem and potentially the risk of recurrent wheezing by age 3 years 22, 23 and recurrent wheezing accompanied by asthma at age 4 years. Over the past 60 years, multiple cohort studies have demonstrated that up to 60% of infants hospitalized with bronchiolitis will have recurrent wheezing within 3 years of hospitalization. 2 Furthermore, compared with healthy control subjects, children hospitalized with bronchiolitis are also at much greater risk of childhood asthma 2 and even longer-term respiratory morbidity. 24 Many of these bronchiolitis cohorts focused on viral infections (eg, RSV). However, viruses infect infants within a complex microbial ecosystem (ie, microbiome), 25 and research results over the last decade have demonstrated that airway bacteria, including nasal microbiota, are also associated with wheezing outcomes in children. 6, 7, 12, 26 Several studies of airway microbiota have implicated Moraxella, Haemophilus, and/or Streptococcus species as risk bacteria for wheezing and asthma by using both cross-sectional 4, [8] [9] [10] [11] 26 and prospective 6,7,12 study designs. Teo et al 12 found in the Australian Childhood Asthma Study (n 5 234) that asymptomatic nasopharyngeal colonization with a high abundance of Streptococcus species around age 2 months was associated with wheeze at age 5 years. 12 In the same cohort greater frequency of detecting nasopharyngeal Moraxella, Haemophilus, or Streptococcus species at times of health between the ages of 6 months and 2 years was associated in sensitized children (n 5 73) with greater odds of chronic wheeze at age 5 years and in nonsensitized children (n 5 64) with any wheeze in the first 3 years of life (ie, transient wheeze). 6 For the first time, we have examined the longitudinal nasal abundances of these 3 risk bacteria (ie, Moraxella, Haemophilus, and Streptococcus species) in a large and diverse multicenter cohort of US infants hospitalized for bronchiolitis. The relative abundances of these risk bacteria on hospitalization were not associated with recurrent wheezing by age 3 years. However, in the weeks after severe viral infection, enrichment of an infant's nasal ecosystem with either Moraxella or Streptococcus species increased the risk of recurrent wheezing by age 3 years and . Kaplan-Meier curves were generated for each of the 3 risk groups for the outcome of recurrent wheezing by age 3 years. Curves were stratified by membership in the risk group. The table shows the number of subjects still at risk for recurrent wheeze based on event occurrences and study dropout. Infants can be included in more than 1 risk group. Log-rank test results indicated that differences in hazard functions of stratified curves were significant for all risk groups. = J ALLERGY CLIN IMMUNOL VOLUME 145, NUMBER 2 recurrent wheezing accompanied by asthma at age 4 years. Similar to the Childhood Asthma Study, 12 we also found a separate group of children who at a median age of 9 months (ie, around 6 months after hospitalization), although healthy, had increased Streptococcus species abundance and a greater risk of recurrent wheezing by age 3 years. However, it remains unclear why these infants have a later increase in Streptococcus species abundance and, as seen in Fig 3, a later onset of recurrent wheezing. The lack of a significant association with Haemophilus species might have been due to inadequate power, a weaker association, or both. Importantly, we confirmed all of these findings using multiple analytic techniques. Moreover, our joint models demonstrated that the present microbiota findings are independent of multiple other potential confounding factors associated with childhood asthma, including viral cause and IgE sensitization in infancy.

These findings suggest that interventions that alter the nasal microbiome after hospitalization for bronchiolitis might be beneficial. Indeed, Beigelman et al 23 found in a small (n 5 39) randomized placebo-controlled trial that azithromycin for 14 days in infants hospitalized with RSV-related bronchiolitis reduced nasal lavage Moraxella species abundance at 2 weeks. 22 Furthermore, in concordance with our findings, decreased Moraxella species abundance at 14 days, no matter the group assignment, was associated with reduced odds of recurrent wheezing (ie, > _3 wheezing episodes) over the subsequent 12 months. 22 This potential approach to bronchiolitis treatment needs further study and thoughtful consideration about possible unintended consequences of widespread antibiotic use. Indeed, bronchiolitis and its long-term respiratory outcomes might be best understood as a complex adaptive system. 27 Taken together, the potential health effects of the airway microbiome are tantalizing. However, we note that early Moraxella species colonization in healthy infants was associated with lower rates of respiratory tract infection. 28 Thus understanding the potential protective role of Moraxella species in healthy populations, compared with its consistent identification as a risk bacterium (not only before and during acute respiratory tract infections 6, 11 but also in asthmatic patients 8 ), will be critical to advancing our understanding of the role of this genus (and ultimately the strain of Moraxella species) in respiratory conditions. It is also possible that the posthospitalization increases in Moraxella or Streptococcus species abundance might simply be markers of children who are already predisposed to recurrent wheezing. Alternatively, the association between posthospitalization increases in Moraxella or Streptococcus species abundance and recurrent wheezing might be causal (ie, microbiota might influence host immune response to viral infection). 23, 29, 30 If true, influencing microbiome reassembly after hospitalization for bronchiolitis either with antibiotics 22, 23 or, preferably, with intranasal administration of protective species 31 might be an approach to help prevent asthma in this high-risk group. This study has potential limitations. First, although bronchiolitis is the leading cause of hospitalization in infants, we suggest caution when generalizing our results beyond the approximately 130,000 US infants hospitalized with bronchiolitis each year. 32 Second, bronchiolitis is a condition that affects the entire airway from the nose to the lungs. 32 Although our samples from the nose might not represent pathobiology in the lung, nasal and bronchial cells have similar responses to stimulation and significant mediator release correlations. 33 Third, there was no healthy control group in MARC-35, but our study objective was not to examine the association between the microbiome and the development of bronchiolitis within the general population but rather to determine the relation of the nasal microbiota to recurrent wheezing within the high-risk population of infants hospitalized with bronchiolitis.

Fourth, although we have 88% follow-up at age 3 years, the number of nasal swabs included in the analyses decreased over time because of parents forgetting to collect swabs and the exclusion of samples collected after the onset of recurrent wheezing, as required by the joint model.

Fifth, the present results are unable to untangle how host environmental conditions, including bacterial-bacterial interactions at the time of infection, might enhance the survival of certain species. 34 However, our group is currently examining the role of airway microbe-host relations in the development of childhood asthma.

Sixth, these analyses will need to be repeated when the children in the cohort turn 6 years of age and are formally evaluated for asthma.

Finally, the present analysis accounted for neither potential genetic differences within our sample 32 nor different strains of Moraxella and Streptococcus species or microbiome function. We are currently pursuing these important research questions.

In summary, we have identified in a large, prospective, multicenter cohort of diverse US infants hospitalized with bronchiolitis that increased abundance of nasal Moraxella or Streptococcus species in the weeks and months after hospitalization was related to recurrent wheezing by age 3 years. These results not only confirm that Moraxella and Streptococcus species are risk bacteria but also identify a potential window of opportunity for clinicians to ameliorate the chronic respiratory outcomes for this high-risk group of children. Indeed, these findings might serve as a starting point to begin investigations that would change management for a common pediatric condition that has no proven effective intervention beyond supportive care. 5 102) ; B, increased Streptococcus species abundance at the clearance time point (n 5 48); and C, increased Streptococcus species abundance at the summer time point (n 5 58). Infants were classified into risk groups using the mean cutoffs for Moraxella (mean 5 0.08) and Streptococcus (mean 5 0.04) species abundance at the relevant time points and might have been included in more than 1 risk group. The rug plot along the x-axis denotes the distribution of observation times for subjects in the specified risk group. Acinetobacter and Bacillus species abundances did not have noticeable variation over time and were excluded to improve readability. (18) .61

All data are presented as numbers (percentages), unless otherwise indicated. P values were computed by using x 2 tests (categorical variables), Wilcoxon rank sum tests (continuous variables), and Kruskal-Wallis rank sum tests (relative abundances). IQR, Interquartile range. *One index swab was lost during shipment. Detected through NPA samples collected at the index visit. àDefined as detection of any positive values for serum allergen-specific IgE at the index visit for both food allergens and aeroallergens. §Defined as admission to the intensive care unit and/or use of mechanical ventilation (continuous positive airway pressure and/or intubation during inpatient stay, regardless of location) at any time during the index hospitalization. All data are presented as numbers (percentages), unless otherwise indicated. Subjects in more than 1 risk group are excluded. P values were computed by using Fisher exact tests (categorical variables) and Wilcoxon rank sum tests (continuous variables). IQR, Interquartile range; RV, rhinovirus. *Detected through NPA samples collected at the index visit. Defined as detection of any positive values for serum allergen-specific IgE at the index visit for both food allergens and aeroallergens. àDefined as admission to the intensive care unit and/or use of mechanical ventilation (continuous positive airway pressure and/or intubation during inpatient stay, regardless of location) at any time during the index hospitalization. (7) 20 (7) All data are presented as numbers (percentages). P values were computed by using x 2 tests.

Trends in bronchiolitis hospitalizations in the United States

Infectious pathogens and bronchiolitis outcomes

Comparative analyses of the bacterial microbiota of the human nostril and oropharynx

Disordered microbial communities in asthmatic airways

Association of nasopharyngeal microbiota profiles with bronchiolitis severity in infants hospitalised for bronchiolitis

Airway microbiota dynamics uncover a critical window for interplay of pathogenic bacteria and allergy in childhood respiratory disease

Childhood asthma after bacterial colonization of the airway in neonates

Bacterial microbiota of the upper respiratory tract and childhood asthma

Airway microbiota and bronchial hyperresponsiveness in patients with suboptimally controlled asthma

Asthma-associated differences in microbial composition of induced sputum

Nasopharyngeal microbiome in respiratory syncytial virus resembles profile associated with increased childhood asthma risk

The infant nasopharyngeal microbiome impacts severity of lower respiratory infection and risk of asthma development

Haemophilus-dominant nasopharyngeal microbiota is associated with delayed clearance of respiratory syncytial virus in infants hospitalized for bronchiolitis

Clinical practice guideline: the diagnosis, management, and prevention of bronchiolitis

The association between anterior nares and nasopharyngeal microbiota in infants hospitalized for bronchiolitis

Prospective multicenter study of viral etiology and hospital length of stay in children with severe bronchiolitis

EPR-3): guidelines for the diagnosis and management of asthma-summary report

Asthma and wheezing in the first six years of life. The Group Health Medical Associates

Cord-blood 25-hydroxyvitamin D levels and risk of respiratory infection, wheezing, and asthma

Joint modelling of time-to-event and multivariate longitudinal outcomes: recent developments and issues

Why we (usually) don't have to worry about multiple comparisons

Azithromycin therapy during respiratory syncytial virus bronchiolitis: upper airway microbiome alterations and subsequent recurrent wheeze

Randomized trial to evaluate azithromycin's effects on serum and upper airway IL-8 levels and recurrent wheezing in infants with respiratory syncytial virus bronchiolitis

Do childhood respiratory infections continue to influence adult respiratory morbidity?

Respiratory syncytial virus and rhinovirus severe bronchiolitis are associated with distinct nasopharyngeal microbiota

The nasal microbiome in asthma

Towards an ecology of the lung: new conceptual models of pulmonary microbiology and pneumonia pathogenesis

Early respiratory microbiota composition determines bacterial succession patterns and respiratory health in children

Nasopharyngeal microbiota, host transcriptome, and disease severity in children with respiratory syncytial virus infection

Pathogenic bacteria colonizing the airways in asymptomatic neonates stimulates topical inflammatory mediator release

Intranasal administration of probiotic Lactobacillus rhamnosus GG prevents birch pollen-induced allergic asthma in a murine model

Viral bronchiolitis in children

Nasal epithelial cells as surrogates for bronchial epithelial cells in airway inflammation studies

Airway response to respiratory syncytial virus has incidental antibacterial effects

Principal investigators at the 17 participating sites in MARC-35

MD Children's Hospital of Los Angeles

PhD Children's Hospital of Philadelphia

MD Children's Mercy Hospital & Clinics

Children's Hospital and Medical Center

MD Phoenix Children's Hospital

We thank the MARC-35 hospitals and research personnel for their ongoing dedication to bronchiolitis and asthma research (see Table E5 in this article's Online Repository at www.jacionline.org). We also thank Alkis Togias, MD, for his ongoing support.

Clinical implications: In a prospective multicenter cohort of infants hospitalized with bronchiolitis, we found the posthospitalization reassembly of the nasal microbial ecosystem was associated with recurrent wheezing by age 3 years.