Skip to main content

The developing airway and gut microbiota in early life is influenced by age of older siblings

Abstract

Background

Growing up with siblings has been linked to numerous health outcomes and is also an important determinant for the developing microbiota. Nonetheless, research into the role of having siblings on the developing microbiota has mainly been incidental.

Results

Here, we investigate the specific effects of having siblings on the developing airway and gut microbiota using a total of 4497 hypopharyngeal and fecal samples taken from 686 children in the COPSAC2010 cohort, starting at 1 week of age and continuing until 6 years of age. Sibship was evaluated longitudinally and used for stratification. Microbiota composition was assessed using 16S rRNA gene amplicon sequencing of the variable V4 region.

We found siblings in the home to be one of the most important determinants of the developing microbiota in both the airway and gut, with significant differences in alpha diversity, beta diversity, and relative abundances of the most abundant taxa, with the specific associations being particularly apparent during the first year of life. The age gap to the closest older sibling was more important than the number of older siblings. The signature of having siblings in the gut microbiota at 1 year was associated with protection against asthma at 6 years of age, while no associations were found for allergy.

Conclusions

Having siblings is one of the most important factors influencing a child’s developing microbiota, and the specific effects may explain previously established associations between siblings and asthma and infectious diseases. As such, siblings should be considered in all studies involving the developing microbiota, with emphasis on the age gap to the closest older sibling rather than the number of siblings.

Video abstract

Key messages

  • Siblings are a strong determinant of the developing airway and gut microbiota, possibly serving as a mediator for the previously observed health effects of growing up with older siblings.

  • Age gap to the closest sibling is more important than the number of siblings.

Background

From birth, humans are continuously exposed to a multitude of microorganisms, many of which take up residence on the surfaces of our bodies, thus establishing the human microbiota [1]. Especially in the first years of life, these microbial communities are shaped and matured by environmental factors as well as host factors. Simultaneously, the microbiota affects our health, e.g., through complex interactions with the immune system and the metabolism [2].

Major factors shaping the early microbiota include mode of delivery, breastfeeding, and antibiotic exposure [3]. All of these factors have in turn been linked to numerous health outcomes during childhood [4,5,6]. Additionally, they are all modifiable to a certain extent, further justifying their intense scrutiny. Likewise, siblings have been linked to numerous health outcomes. These include decreased risk of atopic diseases [7] and inflammatory bowel diseases [8], as well as increased risk of infectious diseases [9] and even some malignancies associated with infectious diseases [10, 11].

The infant’s microbiota is in early life affected by its closest microbiota reservoirs. Here, the mother is of large importance [12], but various other exposures are related to the family, including family size and living environment [13]. Among these, the presence of older siblings in the home seems to be very important. However, previous studies on siblings and the microbiota have either been limited in size [14,15,16] and only reported associations of siblings as one among many other exposures [3, 17,18,19,20,21,22,23,24] and no studies report long-term associations. Some of these studies have described the influence of the number of older siblings, while the age gap to the older child has not been explored. Furthermore, most of these studies are evaluating only differences in the gut microbiota composition after exposure to siblings, while only few evaluate the airway composition [3, 17]. Curiously, in the development of the airway microbiota in early life, siblings have been shown to have an even larger impact than all of the above-mentioned factors (delivery mode, breastfeeding, and antibiotic exposure) [3, 17].

We therefore aimed to investigate the specific effects of siblings on the developing microbiota in both the airway and the gut, using data from 16S rRNA gene amplicon sequencing of hypopharyngeal and fecal samples from the 700 children in the Copenhagen Prospective Studies on Asthma in Childhood 2010 (COPSAC2010) cohort, collected between 1 week and 6 years of age, and furthermore, to explore the sibling signature in the microbiota in relation to asthma, allergic rhinitis, and sensitization at age of 6 years. While the presence or absence of siblings can hardly be considered a modifiable risk factor in the traditional sense, further knowledge on these microbial effects may facilitate intervention elsewhere in the causal chain and is essential to interpret the underlying mechanisms of microbiota–disease associations.

Methods

Study population

The COPSAC2010 cohort is an ongoing population-based mother–child cohort with the main objective to study and reduce the burden of disease from asthma, eczema, and allergy through deep phenotyping [25]. A total of 738 women were recruited at pregnancy week 24 from 2008 to 2010. A total of 700 of their children were included at 1 week of age and were followed prospectively by the in-house pediatricians with scheduled visits at 1 week; 1, 3, 6, 12, and 18 months; and 2, 2 1/2, 3, 4, 5, and 6 years of age. Data on both younger and older siblings were collected or updated prospectively at each scheduled visit registering sex and date of birth.

Sample collection

Hypopharyngeal aspirates were obtained by a study physician at ages 1 week, 1 month, and 3 months using a soft suction catheter passed through the nose [17]. Aspirates were diluted in 1 mL sterile 0.9% NaCl and transported to the microbiological laboratory at Statens Serum Institut (Copenhagen, Denmark), where they were stored at −80 °C until DNA extraction.

Fecal samples were collected at ages 1 week, 1 month, 1 year, 4 years, and 6 years, either at the research clinic or by the parents at home using detailed instructions [26]. Each sample was mixed on arrival to the laboratory with 1 mL of 10% vol/vol glycerol broth and likewise stored at −80 °C.

Sequencing and bioinformatics

We examined the microbiota using 16S rRNA gene amplicon sequencing, as previously described [26, 27]. In brief, DNA was extracted from both sample types using the PowerMag® Soil DNA Isolation Kit (MO-BIO Laboratories), after which the hypervariable V4 region of the 16S rRNA gene was amplified, using the modified broad range primers 515F and 806R. Paired-end sequencing was performed on the Illumina MiSeq System (Illumina) (please see Supplemental Materials for more details). Sequencing adapters were removed with Cutadapt v1.15 [28]. Reads were analyzed using QIIME 2 v2018.2.0 [29] and denoised using DADA2 [30] to infer the ASVs present and their relative abundances across the samples. Forward and reverse reads were trimmed at the 5′ end till 8 bp; other quality parameters used dada2 default values. Taxonomy was assigned using a pre-trained Naïve Bayes classifier (Silva database, release 132, 99% ASV) [31]. All 16S sequences have been deposited in the Sequence Read Archive repository with the accession numbers PRJNA340273 (hypopharyngeal aspirates) and PRJNA417357 (fecal samples).

Covariates

Information on all covariates was obtained through interviews with the primary caregivers during scheduled visits to the COPSAC clinic. Information on antibiotic use was validated using national registries [32].

Clinical endpoints

The clinical endpoints, asthma, allergic rhinitis, and allergic sensitization were evaluated cross-sectionally at age 6 years, while the number of lower respiratory tract infections was summarized across the second and third years of life.

Asthma was diagnosed prospectively and required all following points: (1) recurrent episodes of troublesome lung symptoms (5 episodes within 6 months, each lasting at least 3 consecutive days, 4 weeks of continuous symptoms, or an acute severe exacerbation requiring oral prednisone or hospital admission); (2) symptomatology typical of asthma, including exercise-induced symptoms, prolonged nocturnal cough, and persistent cough outside the common cold; (3) need for intermittent rescue use of inhaled β2-agonist; and (4) response to a 3-month course of inhaled corticosteroids and relapse upon ended treatment [25].

Allergic rhinitis was based on allergic sensitization and clinical interviews of the parents on the history of significant nasal congestion, sneezing, and/or runny nose outside periods with the common cold and in the relevant period of the sensitized aeroallergen.

Allergic sensitization was defined as any positive skin prick test (SPT) ≥3 mm (ALK-Abello, Horsholm, Denmark) or specific IgE (sIgE) in serum ≥0.35 kUa/L against common allergens (birch, grass, mugwort, horse, dog, cat, D. pteronyssinus, D. farinae, Aspergillus fumigatus, Cladosporium herbarum, Penicillium notatum, Alternaria alternata) and/or food allergens (milk, egg, wheat flour, soybean, peanut, cod, shrimp, rye flour, pork, and beef) (ImmunoCAP; Thermo Fisher Scientific, Allerod, Denmark). Children classified as “not sensitized” were both SPT and sIgE negative for all tested allergens.

Lower respiratory infection was defined as a diagnosis of pneumonia or bronchiolitis. Pneumonia was doctor diagnosed and based on clinical appearance, whereas bronchiolitis was defined as cough, tachypnea, chest retractions, and auscultating widespread crepitation and/or rhonchi in a child below 1 year of age, independent of identified pathogens, X-ray or laboratory findings [33].

Statistics

All data analyses were performed in the statistical software “R” version 3.6.1 [34], with notable add-on packages “tidyverse” v1.3.1 [35] (for general data handling), ggplot2 v3.3.5, ggpubr v0.4.0, and cowplot v1.1.1 (for plotting), “phyloseq” v1.38 [36] (to handle the microbiota data), vegan v2.5-7 (for beta diversity and marginal adonis calculation), rabuplot (https://github.com/jstokholm/rabuplot for relative abundance plotting and test), caret v6.0-91, MLeval v0.3 (for classification using LDA), and Tw2 [37] (https://github.com/alekseyenko/Tw2/ for multivariate Welsh-type t-test). Alpha diversity was measured using the Shannon diversity index and ASV richness and tested using Wilcoxon rank-sum or Kruskal–Wallis test, and Dunn’s test with false discovery rate (FDR) control for multiple testing. Beta diversity was computed using weighted UniFrac distances from the raw ASV counts after excluding samples below 2000 reads. The beta diversity is calculated for each compartment but across all time points. Differences in beta diversity were visualized with principal coordinates analysis (PCoA) plots, in which the PCoA ordinations were calculated based on all samples from that sample site. Beta diversity was tested for inference using permutational multivariate analysis of variance (PERMANOVA, Adonis from the “R” package “vegan” [38]) with 9999 permutations. The assumption of variance homogeneity between groups in the PERMANOVA test was evaluated by the betadisper function in vegan [38]. For cases with statistically heterogenous variance, robustness of the results was evaluated by post hoc analysis using a Welsh-type t-test for distance data [37]. Differences in relative abundances were tested using the Wilcoxon rank-sum or Kruskal–Wallis test for the top ten most abundant genera and Benjamini–Hochberg FDR.

A microbiota-derived sibling score was created by training a linear discriminant analysis model to classify the status of older siblings (yes/no) from the corresponding microbiota profiles represented by the first 10 multidimensional scaling (MDS) components from a weighted UniFrac compression of the community distribution. To avoid this score being overfitted towards class separation, training was done on these 10 components, and the cross-validated class probabilities were used as a sibling score, where a score close to 1 resembles the microbiota of children with older siblings.

In the study of longitudinal carryover effects between time points of the microbiota in relation to siblings, the sibling score was used as a univariate measure of the microbiota-related sibling signature: At each age, the current sibling score was regressed on to all prior sibling scores recording the cumulative explained variance as well as current sibling status, using consecutive type II analysis of variance (ANOVA) inference.

The associations between the sibling microbiota scores and asthma, sensitization, and allergic rhinitis at age 6 years were investigated univariately by odds ratios (OR) from logistic regression models with a standardized logit version of the sibling probability score, such that the OR reflects a change in the population of one standard deviation. These analyses were conducted both as unadjusted, as well as adjusted for and stratified by sibling status to reflect the net microbiota relation.

Results

Seven hundred children were initially enrolled in the COPSAC2010 cohort. Ten twins were excluded post hoc for the purpose of this study. Of the remaining 690 children, 686 (99.4%) had at least one sample successfully sequenced for assessment of the airway and/or gut microbiota at a depth of 2000 reads or more. One fecal sample taken at 1 year of age was excluded post hoc, as it represented the only sample from a child with younger siblings at that sample time. All other samples (nairway = 1832, ngut = 2665) were included in the further analyses. The number of children sampled at each time point, as well as sibling status and other covariates, is shown in Supplemental Table 1.

A detailed overview of sibship at each sample time is shown in Supplemental Table 2. At 1 week, 58.2% (of the 638 children with either an airway or a gut sample) had one or more older siblings. This changed minimally until the 1-year sample time, with minor fluctuations exclusively due to slightly different sample populations with concordant microbiome data at each sample time. After the 1-year sample time, numerous younger siblings were introduced: At 6 years, only 11.1% had no siblings, whereas 45.2% had exclusively older siblings, 33.2% had exclusively younger siblings, and 10.5% had both.

The airway microbiota and older siblings

The overall composition of the airway microbiota in the COPSAC2010 cohort, including temporal changes, has previously been described using operational taxonomic units (OTUs) [17]. For the current study, ASVs were used instead. In brief, a total of 7876 ASVs were identified in the airway, most frequently belonging to the phyla Firmicutes (mean relative abundance: 60.9%), Proteobacteria (30.3%), Actinobacteria (5.5%), and Bacteroidetes (1.7%). Alpha diversity increased significantly between 1 week, 1 month, and 3 months of age (Shannon diversity index (SDI): median [IQR] 1.07 [0.61–1.59], 1.36 [0.89–1.77], and 1.52 [1.11–1.99]; Wilcoxon p < 0.001 for both time spans). Likewise, there was a significant separation in beta diversity between time points (PERMANOVA, p < 0.001, r2 = 0.095).

Figure 1 shows alpha diversity in the airway across all sample times, stratified by siblings. Older siblings were associated with a significantly higher alpha diversity at 1 week (SDI: older = 1.1 [0.72–1.65] vs none = 0.93 [0.51–1.51], Wilcoxon p = 0.02). At 1 month, there was no significant difference in alpha diversity (older = 1.33 [0.9–1.74] vs none = 1.46 [0.87–1.80], p = 0.24), whereas older siblings were associated with a significantly lower alpha diversity at 3 months (older = 1.47 [1.09–1.91] vs none = 1.65 [1.16–2.07], p = 0.0057).

Fig. 1
figure 1

Alpha diversity (Shannon diversity index), stratified by sample site, time, and siblings. Boxplots demonstrate medians and IQR. p-values determined by Dunn’s test with FDR control for multiple testing. The overall statistical difference at 4 and 6 years is p = 0.014 and p = 0.09, respectively (determined by the Kruskal–Wallis test). Richness results and Shannon diversity based on rarefied data are shown in Supplemental Fig. 1 and Supplemental Table 3, respectively

Older siblings were also associated with a significant shift in airway beta diversity at all sample times, as shown in Fig. 2 (PERMANOVA, all p < 0.001).

Fig. 2
figure 2

Beta diversity stratified by sample site, time, and siblings. PcoA plots of weighted UniFrac distances showing the first two axes. Ellipses demonstrate means ± 1 SD. p-value and r2 determined by PERMANOVA (9999 permutations) for siblings (yes vs no) after adjustment for covariates

The presence or absence of siblings accounted for an increasing explained variance in airway beta diversity with each sample time (1 week: r2=0.019; 1 month: r2=0.031; 3 months: r2=0.049). When including all covariates in the same model, no other covariate had a higher impact than siblings on airway microbiota composition at any sample time, with the closest being maternal education (1 week: r2=0.015), exclusive breastfeeding (1 month: r2=0.005), and use of antibiotics in 3 months prior to sampling (1 week: r2=0.005), as shown in Fig. 3.

Fig. 3
figure 3

Microbial variance of weighted UniFrac distances explained by 15 different covariates, including siblings, stratified by sample site and time. Horizontal bars and corresponding labels show the variance (r2) explained by each covariate as well as the p-value, determined by PERMANOVA (9999 permutations). All factors are mutually adjusted apart from the two breastfeeding variables, which are not adjusted for each other. Breastfeeding is defined as yes/no at sampling for the 1-week, 1-month, and 3-month time points and for total breastfeeding additionally the 1-year time point, while the time points thereafter reflect the breastfeeding duration. Distribution of covariates across time is shown in Supplemental Table 1. Test for variance homogeneity and post hoc robustness test are shown in Supplemental Tables 4 and 5

Differential abundance analyses at the phylum level (Supplemental Fig. 2) showed a significantly lower abundance of Firmicutes at all sample times in the airways of children with older siblings (Wilcoxon, all p < 0.001) and a corresponding significantly higher abundance of Proteobacteria at all sample times (all p < 0.001). Actinobacteria were significantly more abundant at 1 month (p = 0.002) and 3 months (p < 0.001) in children with older siblings (p = 0.002). No significant difference was found for Actinobacteria at 1 week.

At the genus level (Fig. 4), the presence of older siblings was associated with significant changes in the relative abundances of each of the 10 most abundant genera in the airways at one or more sample times. Both Moraxella (Wilcoxon, all p < 0.001) and Neisseria (1 week: p = 0.036; 1 month; p < 0.001, 3 months; p = 0.001) were significantly more abundant at all time points in children with older siblings, whereas Staphylococcus was less abundant at all time points (1 month: p = 0.003; 1 week and 3 months; p < 0.001).

Fig. 4
figure 4

Relative abundances of the 10 most abundant genera (highest mean relative abundance) in each sample site, stratified by sample site, time, and siblings. p-values determined by the Wilcoxon test and q-values from FDR. A pseudocount (+1e−06) was added to all abundances for the log-scale presentation. The black dots indicate median values. Children with younger siblings or both younger and older siblings (only applicable at 4 and 6 years of age) are included in Supplemental Fig. 3. Summary statistics are shown in Supplemental Table 8 and Supplemental Table 10

The gut microbiota and older siblings

The overall composition of the gut microbiota in the COPSAC2010 cohort, including temporal changes, has previously been described [26]. In brief, a total of 12,644 ASVs were identified in the gut, most frequently belonging to the phyla Firmicutes (mean relative abundance: 33.0%), Bacteroidetes (31.8%), Proteobacteria (19.7%), Actinobacteria (13.3%), and Verrucomicrobia (1.9%). Alpha diversity in the gut decreased initially from 1 week to 1 month, but increased with each subsequent sample time, albeit most markedly during the interval from 1 month to 1 year (SDI median [IQR] 1 week = 1.51 [1.08–1.96]; 1 month = 1.4 [1.01–1.8]; 1 year = 2.38 [1.94–2.77]; 4 years = 3.52 [3.21–3.82]; 6 years = 3.67 [3.31–3.89], Wilcoxon p ≤ 0.002 for all intervals). Likewise, age at sampling accounted for a large part of the variance in beta diversity (PERMANOVA, p < 0.001, r2 = 0.255).

When stratifying alpha diversity in the gut by siblings (Fig. 1), children with older siblings had a significantly higher alpha diversity than children with no siblings at the 3 intermediary sample times, though attenuating at 4 years: 1 month (SDI 1.47 [1.03–1.83] vs 1.34 [0.92–1.72], p = 0.03), 1 year (2.45 [2.05–2.86] vs 2.27 [1.81–2.61], p < 0.001), and 4 years (3.57 [3.3–3.87] vs 3.49 [3.21–3.74], p = 0.02).

The presence or absence of siblings was also associated with a significant difference in gut beta diversity, albeit only at the first 3 sample times (PERMANOVA, all p≤ 0.001), as shown in Fig. 2. Siblings accounted for a stable explained variance in gut beta diversity during this period (1 week: r2 = 0.008; 1 month: r2 = 0.009; 1 year: r2 = 0.012), with no significant difference at 4 and 6 years (Fig. 3). At 1 week, the variance explained by siblings was overshadowed only by the effect of delivery mode (r2 = 0.052), whereas at 1 month, both delivery mode (r2 = 0.018) and exclusive breastfeeding (r2 = 0.013) accounted for a higher explained variance than siblings. However, at 1 year, siblings were the single most important determinant of gut beta diversity, followed by delivery mode (r2 = 0.004).

Differential abundance analyses at the phylum level in the gut microbiota (Supplemental Fig. 2) showed significant associations between siblings and the abundance of each of the 4 most abundant phyla at one or more sample times. Comparable patterns were found at 1 week and 1 month, with older siblings being associated with lower abundances of Proteobacteria (p = 0.027 and p < 0.001) and higher abundances of Actinobacteria (both p < 0.001) at both sample times. At 1 year of age, older siblings were associated with lower abundance of both Bacteroidetes (p = 0.028) and Firmicutes (p = 0.007). At 4 and 6 years of age, there was no significant difference in abundances when comparing older siblings to none.

At the genus level (Fig. 4), having older siblings was significantly associated with the abundance of 8 out of the 10 most abundant genera in the gut at one or more sample times. Again, comparable patterns were found at 1 week and 1 month, with older siblings notably being associated with higher abundances of Bifidobacterium (Wilcoxon, both p < 0.001) and lower abundances of Clostridium (p = 0.009 and p < 0.001) at both sample times. The highest number of associations was found at 1 year of age, where children with older siblings had significantly lower abundances of Escherichia/Shigella (p = 0.029), other Enterobacteriaceae (p = 0.002), and Veillonella (p <0.001), while having significantly higher abundances of Prevotella (p = 0.011). Notably, the increased abundance of Prevotella in children with older siblings was even more pronounced at 4 years (p < 0.001) and persisted at 6 years (p = 0.037).

Younger siblings

At the 4- and 6-year sample times, 136 (30.1%) and 146 (33.2%) had younger siblings, while 38 (8.4%) and 46 (10.5%) had both younger and older siblings. No significant differences were found in gut alpha diversity at either sample time when comparing children with younger siblings to children with no siblings (Fig. 1). Having both younger and older siblings at 4 years of age was associated with an increased alpha diversity when compared to none (SDI 3.67 [3.22–3.9] vs 3.35 [3.03–3.66], Wilcoxon p = 0.018). However, a similar increase in alpha diversity was also found in children with exclusively older children at that sample time, as previously mentioned (SDI 3.57 [3.3–3.87]).

Gut beta diversity showed no significant shifts at either 4 or 6 years when comparing younger children to older, both, and none (Fig. 2).

In differential abundance analyses at the phylum level (Supplemental Fig. 3) and the genus level (Supplemental Fig. 4), we found no consistently significant associations with younger siblings at 4 and 6 years of age.

Age gap, number of older siblings, and sex of closest older sibling

To investigate different aspects of the sibling effect on the developing microbiota, we stratified alpha diversity and beta diversity by either age gap to the closest older sibling (“More than 4 years,” “2 to 4 years,” and “Less than 2 years”), the number of older siblings (“1” and “2 or more”), and sex of the closest older sibling. All subsequent statistical tests were done for differences across all strata except “None,” which have nonetheless been included in plots for visual reference. At 4 and 6 years of age, children with younger siblings or both were excluded from these analyses. The distribution of strata across sample times is shown in Supplemental Table 2.

When stratifying alpha diversity by age gap to the closest older sibling (Supplemental Fig. 5), we found a significant difference in the airways at 3 months (Kruskal–Wallis, p = 0.029) and gut at 1 year (p = 0.023), with a stepwise increase in sibling effect with decreasing age gap. The same stepwise pattern was evident in airway beta diversity (Supplemental Fig. 6), where significant differences were found across strata at all 3 sample times (PERMANOVA, 1 week, p = 0.0011; 1 month, p < 0.001; 3 months, p < 0.001). In the gut, age gap to the closest older sibling was associated with a shift in beta diversity at 1 year (p = 0.0028), where the general effects of siblings were previously found to be most pronounced (Supplemental Fig. 7).

In contrast, an increasing number of older siblings was only weakly associated with an increased effect on airway alpha diversity at 1 week (Wilcoxon, p = 0.022). No other significant associations were found for alpha diversity or beta diversity at other sample times or sample sites. When stratified by sex of the closest older sibling, no significant associations were found on any of the tested outcomes.

To further investigate the relationship between age gap, number of older siblings, and microbial impacts, a microbial sibling score was defined as the probability of having older siblings given the microbiota and calculated for each time point, by learning the microbial profiles most associated with having an older sibling: A cross-validated prediction of older siblings (yes/no) was performed using linear discriminant analysis on the first 10 MDS components from a weighted UniFrac matrix for each compartment and sample time individually (Supplemental Fig. 8). The strongest results were obtained for airways at 3 months of age and gut 1 year resembling the beta diversity results (Fig. 2) (for results using other beta diversity metrics, see Supplemental Fig. 9).

The score was then plotted against the age gap to the closest older sibling and the number of older siblings. In the airway (Fig. 5), a sigmoid relationship was observed between sibling score and age gap to the closest older sibling, showing a steep drop in sibling score with age gaps larger than 4 to 6 years. In the gut (Supplemental Fig. 10), the curve appeared logarithmic and only dropped with age gaps over 10 years, where the number of observations grew increasingly sparse.

Fig. 5
figure 5

Airway sibling score plotted as a function of age gap to the closest older sibling (top) and the number of older siblings (bottom), stratified by sample time. Boxplots denote upper and lower quartiles. The middle of each boxplot denotes the mean sibling score, with the mean sibling score of children with no older siblings being extrapolated horizontally (dashed line). Age gap truncated at 10 years; number of older siblings truncated at 3. The curves were fitted by local (alpha = 0.75) quadratic polynomial regression

A numerical total sibling burden was defined from the age gap to all older siblings and used to assess the influence of having several older siblings and the age gap to these. For the microbiota of the gut, a slight marginal increase in correlation is observed, pointing to additional effects of the age gap to other older siblings as well (Supplemental Fig. 11).

As suggested by the stratified analysis (Fig. 5 lower panel), an increasing number of siblings was associated with little to no variation in sibling score, regardless of compartment and sample time. However, the higher number of older siblings was positively associated with a larger age gap to the closest older sibling (Supplemental Fig. 12, Spearman’s rank correlation: rho = 0.21, p < 0.001). Thus, the potential effect of an increasing number of siblings could be canceled out by a larger, opposite effect of an increasing age gap.

Longitudinal effects

The same sibling score was used to test for microbial carryover effects between sample times to evaluate when the siblings, as an environmental exposure, are most prominent, and how strong the potential microbial carryover effect between time points is: Linear models were constructed to predict current sibling scores, based on all previous sibling scores and current sibling status to demonstrate how much of the microbial sibling signature is due to having siblings at present and how much is carried over from previous time points. Estimates of explained variance from these models are visualized in Fig. 6.

Fig. 6
figure 6

Effect of older siblings at sample time compared to carryover from previous sample times. Sibling score at each sample time was modeled as a function of current sibling status as well as previous sibling scores, using consecutive type II ANOVA inference

In the airway, a remarkably large carryover effect was found alongside an equally large direct effect at both 1 month and 3 months of age (all p < 0.001, Supplemental Table 11). However, the carryover effect was limited to temporally adjacent sample times, with no significant carryover from 1 week of age to 3 months of age.

In the gut, significant carryover was also limited to temporally adjacent sample times: from 1 week to 1 month (p < 0.001) and from 1 month to 1 year (p = 0.011). At the age of 1 month, the sibling effect is solely a carryover effect from 1 week, while at 1 year, there is an additional direct effect.

The sibling microbiota and asthma, allergic rhinitis, and sensitization

The microbiota sibling scores from each compartment at each age were used as predictors for asthma, sensitization, and allergic rhinitis status at age 6 years. Of the children with both microbiota samples as well as 6-year disease outcomes, there were approximately 7.0% with asthma, 6.4% with allergic rhinitis, and 25.7% with sensitization. For asthma, the gut microbiota sibling score at 1 year was found significantly associated with reduced risk (OR = 0.97 CI: 0.95–0.99, relative to a score change of 1 SD), where the effect size and statistical strength were unaffected by both adjusting for and stratifying by siblings. None of the sibling scores from the other compartments and time points was associated with asthma. For allergic rhinitis and sensitization, none of the associations was significant (see Supplemental Table 12 for details). No associations were found for lower respiratory tract infections after adjustment (Supplemental Table 13).

Discussion

Primary findings

We found siblings to be one of the most important determinants of the developing microbiota in both the airways and the gut, particularly during the first year of life, with a significant impact on alpha diversity, beta diversity, and relative abundance of the most abundant taxa. These sibling associations were highly dependent on a small age gap to the closest older sibling, with the number of older siblings having comparatively little influence. The sibling signature in the gut microbiota at age 1 year was associated with a reduced risk of asthma at 6 years.

Strengths and limitations

The major strength of the present study is the size and reach: 4497 samples taken from 2 different compartments in 686 children, starting as early as 1 week of age and continuing until 6 years of age. As the effect of siblings on the airway microbiota was continuously increasing in magnitude until the last airway sample time at 3 months of age, later airway samples could have potentially revealed an even greater effect. Likewise, more frequent sampling of both compartments between the ages of 1 and 4 years could potentially provide greater insight into the effects of younger siblings than we were able to glean with the present data.

Information on all covariates was collected and validated in a carefully monitored prospective birth cohort with high follow-up, minimizing recall bias. Despite this, direct comparisons between the effect size of sibship and other covariates should be interpreted with caution. Such comparisons are entirely dependent on how data on each covariate is handled subsequent to collection, and particularly the use of antibiotics has the potential to be evaluated in a multitude of ways. Furthermore, the present study population exhibited a high degree of homogeneity in daycare use, and partly also the duration of breastfeeding. This was to be expected from a Danish study population but will inevitably have limited the detection of microbial variance due to these factors.

16S rRNA gene amplicon sequencing of the V4 region is limited in taxonomic resolution [39]. While this limits our understanding of which specific taxa the siblings introduce, it will not limit our overall findings on a changed microbiome from having siblings. Furthermore, the carryover effect as described in Fig. 6 may not necessarily reflect retainment of a specific microbial community, but may just as well reflect a persistent environmental exposure to siblings.

Interpretation

To the best of our knowledge, the present study is the most extensive investigation of the impact of siblings on the developing microbiota thus far, with the majority of previous studies being limited in size [14,15,16] or only tangentially addressing the topic [3, 17,18,19,20,21,22,23,24].

Studies on the gut microbiota are still by far the most common, and our results are in general agreement with those previously published in this area, where siblings have been associated with increased microbial diversity in the first years of life [14, 22], increased abundance of Bifidobacterium, and decreased abundances of Enterobacteriaceae [23] and Clostridium [21]. We found a consistent long-term increase in the genus Prevotella in the gut after exposure to older siblings beginning from 1 year and lasting to age 6 years. Prevotella is associated with carbohydrate-based diets and fermentation of dietary fibers to short-chain fatty acids [40]. Even at 4 years of age, having older siblings was continuously associated with increased alpha diversity in the gut (beginning at 1 month of age), though increased stratification at 4 and 6 years due to the arrival of younger siblings resulted in a steady loss of statistical power. For the same reason, disentangling the respective effects of younger and older siblings also proved difficult, with analyses on alpha diversity nonetheless suggesting an additive rather than opposing effect of younger and older siblings. We speculate that older siblings are direct seeding sources of microbiota inducing higher diversity and as well as transfer of specific bacteria. This transfer may decrease with the age of the older child, as hygiene measures probably improve with age. Also, we speculate that parental behavior such as hygiene and dietary practices could differ between the first-born child and later children. Studies on the airway microbiota are comparatively sparse, though also in general agreement with our results: Older siblings have thus previously been reported to be associated with lower alpha diversity at 3 months [16] and increased abundances of particularly Moraxella [16, 24] and Pasteurellaceae, including the genus Haemophilus [3], perfectly in line with our findings and may be explained by the close contact between the older siblings and the studied child, which is also a plausible explanation for a higher risk of infections from the presence of older children [9].

Our results are supportive of the “Hygiene hypothesis” [41] and its successor the “Old Friends hypothesis” [42], where the microbial impact of siblings is hypothesized to reduce the risk of atopic diseases. We only found associations between the sibling signature in the gut microbiota at 1 year and asthma at age 6 years, while no associations were found at other ages or for allergic sensitization and allergic rhinitis. The 1-year gut microbiota is however especially interesting, as we have previously found robust associations between the gut maturity here and later asthma [26] as well as an increased asthma risk in children born by cesarean section if their microbiota did not normalize by age 1 year [43]. Importantly, in these studies, older siblings in the home were associated both with higher gut maturity and with normalization of a cesarean section perturbed gut microbiota by age 1 year. Likewise, the increased abundance of Moraxella and Haemophilus in particular is likely mediators of the increased risk of respiratory tract infections associated with having siblings [44].

Sporadic evidence of differential effects of brothers and sisters on the risk of atopic diseases has been around for as long as the hygiene hypothesis itself [45]. This is not entirely implausible, as sex has also been implicated in differences in the gut microbiome [20]. Furthermore, behavioral differences between older sisters and brothers [46] might be hypothesized to result in a differential impact on the developing microbiota of their younger sibling. Despite the above, we found no evidence for such sex-dependent differences.

A stepwise correlation between the number of older siblings and their impact on gut microbiota diversity has previously been shown [14]. However, here, we demonstrate the importance of the age gap to the closest older sibling and how it has a comparatively much larger modulating effect. One explanation may be that shorter age gap results in a closer relationship, as closer relationships have previously been associated with higher gut microbiota similarity in an adult study population comprising both sibling pairs and spouses [47]. Also, children have a higher rate of respiratory tract infections and gastrointestinal infections during their first few years of life [48], as well as different approaches to hygiene, possibly allowing them to exert a larger microbial impact than their older counterparts.

An additional explanation is related to the maturation of the human microbiota: While this maturation is a continuous process throughout life, the most striking development of the gut microbiota occurs within the first few years of life [49]. During this early maturation, it is plausible that the microbiota would be more likely to be influenced by other microbiotas on a not-too-distant level of maturity. This hypothesis is supported by previous findings that parents share significantly more tongue and gut communities with their own children at older ages (3–18 years) than with other children, while the same is not the case for parents and their younger infants (<3 years) [50]. This is also supported by our findings that older siblings aged 4 years and above may have a larger impact on the airway microbiota at later sample times, whereas the 1-week time point was mainly affected by siblings aged 3 years and below.

Regardless of the mechanism behind them, the implications of the above results are not to be underestimated and may extend beyond the hygiene hypothesis. Indeed, care should also be exercised when using siblings as healthy controls in microbiota studies, as the results may be confounded by sibship [51, 52]. Finally, we suggest that the frequent stratification by number of siblings in microbiota-related studies should be reconsidered in favor of stratification by age gap to the closest sibling.

Conclusions

We found that the presence of older siblings is associated with increased bacterial diversity in the gut and decreased bacterial diversity in the airways particularly early in childhood. These results are in line with previous reports of siblings being associated with microbial changes leading to decreased risk of asthma. We urge colleagues to consider the age gap to the closest sibling rather than just the number of siblings in future studies involving the developing microbiota.

Availability of data and materials

The 16S sequences have been deposited in the Sequence Read Archive repository with the accession numbers PRJNA340273 (hypopharyngeal aspirates) (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA340273/) and PRJNA417357 (fecal samples) (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA417357). All other relevant data are available from the corresponding authors upon reasonable requests, after signing a data access agreement.

Abbreviations

ANOVA:

Analysis of variance

ASV:

Amplicon sequence variant

COPSAC2010 :

Copenhagen Prospective Studies on Asthma in Childhood 2010

PERMANOVA:

Permutational multivariate analysis of variance

rRNA:

Ribosomal RNA

SDI:

Shannon diversity index

PCoA :

Principal coordinates analysis

WUF:

Weighted UniFrac

References

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

  2. Thorsen J, Rasmussen MA, Waage J, Mortensen M, Brejnrod A, Bønnelykke K, et al. Infant airway microbiota and topical immune perturbations in the origins of childhood asthma. Nat Commun. 2019;10:1–8.

    Article  CAS  Google Scholar 

  3. Bosch AATM, de Steenhuijsen Piters WAA, van Houten MA, Chu MLJN, 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:1582–90.

    Article  PubMed  Google Scholar 

  4. Sevelsted A, Stokholm J, Bisgaard H. Risk of asthma from cesarean delivery depends on membrane rupture. J Pediatr. 2016;171:38–42.e1-4.

    Article  PubMed  Google Scholar 

  5. Pedersen TM, Stokholm J, Thorsen J, Mora-Jensen A-RC, Bisgaard H. Antibiotics in pregnancy increase children’s risk of otitis media and ventilation tubes. J Pediatr. 2017.

  6. Giwercman C, Halkjaer LB, Jensen SM, Bønnelykke K, Lauritzen L, Bisgaard H. Increased risk of eczema but reduced risk of early wheezy disorder from exclusive breast-feeding in high-risk infants. J Allergy Clin Immunol. 2010;125:866–71.

    Article  PubMed  Google Scholar 

  7. Karmaus W, Botezan C. Does a higher number of siblings protect against the development of allergy and asthma? A review. J Epidemiol Commun Health. 2002;56:209–17.

    Article  CAS  Google Scholar 

  8. Cholapranee A, Ananthakrishnan AN. Environmental hygiene and risk of inflammatory bowel diseases: a systematic review and meta-analysis. Inflamm Bowel Dis. 2016;22:2191–9.

    Article  PubMed  Google Scholar 

  9. Sibship structure and risk of infectious mononucleosis: a population-based cohort study. International Journal of Epidemiology. Available from: https://academic.oup.com/ije/article/43/5/1607/2949582. Cited 2020 Jun 24.

  10. Liu Z, Fang F, Chang ET, Adami H-O, Ye W. Sibship size, birth order and risk of nasopharyngeal carcinoma and infectious mononucleosis: a nationwide study in Sweden. Int J Epidemiol. 2016;45:825–34.

    Article  PubMed  Google Scholar 

  11. Westergaard T, Melbye M, Pedersen JB, Frisch M, Olsen JH, Andersen PK. Birth order, sibship size and risk of Hodgkin’s disease in children and young adults: a population-based study of 31 million person-years. Int J Cancer. 1997;72:977–81.

    Article  CAS  PubMed  Google Scholar 

  12. Ferretti P, Pasolli E, Tett A, Asnicar F, Gorfer V, Fedi S, et al. Mother-to-infant microbial transmission from different body sites shapes the developing infant gut microbiome. Cell Host Microbe. 2018;24:133–145.e5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Lehtimäki J, Thorsen J, Rasmussen MA, Hjelmsø M, Shah S, Mortensen MS, et al. Urbanized microbiota in infants, immune constitution, and later risk of atopic diseases. J Allergy Clin Immunol. 2021;148:234–43.

    Article  PubMed  CAS  Google Scholar 

  14. Laursen MF, Zachariassen G, Bahl MI, Bergström A, Høst A, Michaelsen KF, et al. Having older siblings is associated with gut microbiota development during early childhood. BMC Microbiol. 2015;15:154.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  15. Azad MB, Konya T, Maughan H, Guttman DS, Field CJ, Sears MR, et al. Infant gut microbiota and the hygiene hypothesis of allergic disease: impact of household pets and siblings on microbiota composition and diversity. Allergy Asthma Clin Immunol. 2013;9:15.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Hasegawa K, Linnemann RW, Mansbach JM, Ajami NJ, Espinola JA, Fiechtner LG, et al. Household siblings and nasal and fecal microbiota in infants. Pediatr Int. 2017;59:473–81.

    Article  CAS  PubMed  Google Scholar 

  17. Mortensen MS, Brejnrod AD, Roggenbuck M, Abu Al-Soud W, Balle C, Krogfelt KA, et al. The developing hypopharyngeal microbiota in early life. Microbiome. 2016;4:70.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Gut microbiota and development of atopic eczema in 3 European birth cohorts. J Allergy Clin Immunol. 2007;120:343–50.

  19. Penders J, Thijs C, Vink C, Stelma FF, Snijders B, Kummeling I, et al. Factors influencing the composition of the intestinal microbiota in early infancy. Pediatrics. 2006;118:511–21.

    Article  PubMed  Google Scholar 

  20. Martin R, Makino H, Yavuz AC, Ben-Amor K, Roelofs M, Ishikawa E, et al. Early-life events, including mode of delivery and type of feeding, siblings and gender, shape the developing gut microbiota. PLoS One. 2016;11:e0158498.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  21. Penders J, Gerhold K, Stobberingh EE, Thijs C, Zimmermann K, Lau S, et al. Establishment of the intestinal microbiota and its role for atopic dermatitis in early childhood. J Allergy Clin Immunol. 2013;132:601–607.e8.

    Article  PubMed  Google Scholar 

  22. Stewart CJ, Ajami NJ, O’Brien JL, Hutchinson DS, Smith DP, Wong MC, et al. Temporal development of the gut microbiome in early childhood from the TEDDY study. Nature. 2018;562:583.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Yap GC, Chee KK, Hong P-Y, Lay C, Satria CD, Sumadiono, et al. Evaluation of stool microbiota signatures in two cohorts of Asian (Singapore and Indonesia) newborns at risk of atopy. BMC Microbiol. 2011;11:193.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Huy Ta LD, Yap GC, Tay CJX, Lim ASM, Huang CH, Chu CW, et al. Establishment of the nasal microbiota in the first 18 months of life – correlation with early onset rhinitis and wheezing. J Allergy Clin Immunol. 2018;142:86–95.

    Article  PubMed Central  Google Scholar 

  25. Bisgaard H, Vissing NH, Carson CG, Bischoff AL, Følsgaard NV, Kreiner-Møller E, et al. Deep phenotyping of the unselected COPSAC2010 birth cohort study. Clin Exp Allergy. 2013;43:1384–94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Stokholm J, Blaser MJ, Thorsen J, Rasmussen MA, Waage J, Vinding RK, et al. Maturation of the gut microbiome and risk of asthma in childhood. Nat Commun. 2018;9:141.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  27. Gupta S, Hjelmsø MH, Lehtimäki J, Li X, Mortensen MS, Russel J, et al. Environmental shaping of the bacterial and fungal community in infant bed dust and correlations with the airway microbiota. Microbiome. 2020;8:115.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet.journal. 2011;17:10–2.

    Article  Google Scholar 

  29. Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet C, Al-Ghalith GA, et al. QIIME 2: reproducible, interactive, scalable, and extensible microbiome data science. PeerJ Inc.; 2018 Oct. Report No.: e27295v1. Available from: https://peerj.com/preprints/27295.

  30. Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: high-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13:581–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41:D590–6.

    Article  CAS  PubMed  Google Scholar 

  32. Stokholm J, Schjørring S, Pedersen L, Bischoff AL, Følsgaard N, Carson CG, et al. Prevalence and predictors of antibiotic administration during pregnancy and birth. PLoS One. 2013;8:e82932.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  33. Vissing NH, Chawes BLK, Bisgaard H. Increased risk of pneumonia and bronchiolitis after bacterial colonization of the airways as neonates. Am J Respir Crit Care Med. 2013;188:1246–52.

    Article  PubMed  Google Scholar 

  34. Team RC. R: a language and environment for statistical computing [Computer software]. Vienna: R Foundation for Statistical Computing; 2016.

    Google Scholar 

  35. Wickham H, Averick M, Bryan J, Chang W, McGowan L, François R, et al. Welcome to the Tidyverse. JOSS. 2019;4:1686.

    Article  Google Scholar 

  36. McMurdie PJ, Holmes S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One. 2013;8:e61217.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Alekseyenko AV. Multivariate Welch t-test on distances. Bioinformatics. 2016;32:3552–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  38. Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, McGlinn D, et al. vegan: community ecology package. 2019; Available from: https://CRAN.R-project.org/package=vegan.

  39. Bukin YS, Galachyants YP, Morozov IV, Bukin SV, Zakharenko AS, Zemskaya TI. The effect of 16S rRNA region choice on bacterial community metabarcoding results. Sci Data. 2019;6:190007.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Roager HM, Christensen LH. Personal diet-microbiota interactions and weight loss. Proc Nutr Soc. 2022;1–12.

  41. Strachan DP. Hay fever, hygiene, and household size. BMJ. 1989;299(6710):1259–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Rook GAW, Brunet LR. Microbes, immunoregulation, and the gut. Gut. 2005;54:317–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Stokholm J, Thorsen J, Blaser MJ, Rasmussen MA, Hjelmsø M, Shah S, et al. Delivery mode and gut microbial changes correlate with an increased risk of childhood asthma. Sci Transl Med. 2020;12.

  44. Christensen ED, Thorsen J, Stokholm J, Pedersen TM, Brix S, Krogfelt KA, et al. Early life bacterial airway colonization, local immune mediator response and risk of otitis media. J Med Microbiol. 2020.

  45. Strachan DP. Family size, infection and atopy: the first decade of the “hygiene hypothesis”. Thorax. 2000.

  46. Stoneman Z, Brody GH, MacKinnon CE. Same-sex and cross-sex siblings: activity choices, roles, behavior, and gender stereotypes. Sex Roles. 1986;15:495–511.

    Article  Google Scholar 

  47. KA D-MF, Tang Z-Z, Kemis JH, Kerby RL, Chen G, Palloni A, et al. Close social relationships correlate with human gut microbiota composition. Sci Rep. 2019;9:703.

    Article  CAS  Google Scholar 

  48. Vissing NH, Chawes BL, Rasmussen MA, Bisgaard H. Epidemiology and risk factors of infection in early childhood. Pediatrics. 2018;141:e20170933.

    Article  PubMed  Google Scholar 

  49. Derrien M, Alvarez A-S, de Vos WM. The gut microbiota in the first decade of life. Trends Microbiol. 2019;27:997–1010.

    Article  CAS  PubMed  Google Scholar 

  50. Song SJ, Lauber C, Costello EK, Lozupone CA, Humphrey G, Berg-Lyons D, et al. Cohabiting family members share microbiota with one another and with their dogs. eLife. 2013;2 Available from: http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3628085/. Cited 2015 Apr 14.

  51. Son JS, Zheng LJ, Rowehl LM, Tian X, Zhang Y, Zhu W, et al. Comparison of fecal microbiota in children with autism spectrum disorders and neurotypical siblings in the Simons Simplex Collection. PLOS ONE. Public Library of. Science. 2015;10:e0137725.

    Google Scholar 

  52. Duytschaever G, Huys G, Bekaert M, Boulanger L, De Boeck K, Vandamme P. Cross-sectional and longitudinal comparisons of the predominant fecal microbiota compositions of a group of pediatric patients with cystic fibrosis and their healthy siblings. Appl Environ Microbiol. 2011;77:8015–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We would like to express our deepest gratitude to the children and families of the COPSAC2010 cohort study for all their support and commitment. We would also like to acknowledge and appreciate the unique efforts of the entire COPSAC research team.

Funding

All funding received by COPSAC is listed on www.copsac.com. The Lundbeck Foundation (Grant no R16-A1694), The Ministry of Health (Grant no 903516), Danish Council for Strategic Research (Grant no 0603-00280B), and The Capital Region Research Foundation have provided core support to the COPSAC research center. The funding agencies did not have any role in the design and conduct of the study; collection, management, and interpretation of the data; or preparation, review, or approval of the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

The guarantors of the study are MR and JS, from conception and design to the conduct of the study and acquisition of data, data analysis, and interpretation of data. SJS has coordinated all microbiome characterizations from DNA isolation and sequencing to bioinformatic analyses. MHH and EDC have written the first draft of the manuscript. JS, JT, SH, TR, CEP, UT, JR, SG, BC, KB, and MR have provided important intellectual input and contributed considerably to the analyses and interpretation of the data. All authors guarantee that the accuracy and integrity of any part of the work have been appropriately investigated and resolved and all have approved the final version of the manuscript.

Authors’ information

Not applicable.

Corresponding authors

Correspondence to Søren Johannes Sørensen, Morten Arendt Rasmussen, Hans Bisgaard or Jakob Stokholm.

Ethics declarations

Ethics approval and consent to participate

The study was conducted in accordance with the guiding principles of the Declaration of Helsinki and was approved by the Local Ethics Committee (H-B-2008-093) and the Danish Data Protection Agency (2015-41-3696). Both parents gave written informed consent before enrolment.

Consent for publication

Both parents gave written informed consent before enrolment.

Competing interests

EDC, MHH, JS, BC, KB, UT, and MAR all have younger siblings. JS, TR, KB, UT, SJS, and JR have older siblings. The authors declare 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: Supplemental Table 1

: Distribution of all covariates shown in Fig. 3, stratified by sample site and time. For day care and breastfeeding, asterisks denote the encoding used in analyses at each sample time. Supplemental Table 2: Detailed overview of sibship characteristics at each sample time, including sibship categories used in stratified analyses. Supplemental Table 3: Relative difference in Shannon diversity index (siblings versus no siblings) Tested by Wilcoxon rank sum test. Supplemental Table 4: Test of variance homogeneity of weighted UniFrac betadiversity between groups of covariates. The statistics reported is the ratio of the average distance to the centroid between the group with largest over the smallest variance. Supplemental Table 5: Comparison of inference results based on adonis (assumimg variance homogeneity) and Welch t-test for distances (allowing heterogenous variance between groups). Included is only comparisons where the variance homogeneity assumption is problematic (see Supplemental Table 4), for binary covariates, and as crude unadjusted tests. Supplemental Table 6: Age at scheduled visits as function of sibling status at visit timepoint. Statistical inferences were calculated by t-test. Possible age-confounding is present at the 1-week visit. We tested this for both airway and gut results. We observe an effect of siblings (R2airways = 0.0185 and R2gut = 0.0086), and when accounting for the actual age this drops relatively ~5% to (air = 0.0174 and R2gut = 0.0083). However, all signals remain significant. In the airways the actual age accounts for R2 = 0.0036 (p = 0.078), while in the gut this value is R2 = 0.00742 (p = 0.005). Siblings are in general associated with a faster microbiome maturation, and as children with siblings attend the 1-week visit earlier, the confounding effect is in the opposite direction of the sibling effect at this visit. Supplemental Table 7: Summary statistics of the 10 most abundant phyla in the airways, stratified by sample time and siblings. [truncated due to space limitations]. Supplemental Table 8: Summary statistics of the 20 most abundant genera in the airways, stratified by sample time and siblings. [truncated due to space limitations]. Supplemental Table 9: Summary statistics of the 10 most abundant phyla in the gut, stratified by sample time and siblings. [truncated due to space limitations]. Supplemental Table 10: Summary statistics of the 20 most abundant genera in the gut, stratified by sample time and siblings. [truncated due to space limitations]. Supplemental Table 11: Estimates of the direct effect of siblings at current sample time vs. carryover of effects at earlier sample times. Sibling score at each sample time, modeled as a function of current sibling status as well as previous sibling scores, using consecutive type II anova inference. Supplemental Table 12: Association between allergic rhinitis, asthma and sensitization at age six years and sibling microbiome score from gut (type=Fecal) and airways (type = Trach) at various ages (visit). Model indicates type of analysis: unadjusted (all crude), adjusted for current siblings (all adjusted), stratified by siblings (siblings and no siblings respectively). Supplemental Table 13: Association between number of lower respiratory tract infections between zero and three years of age and sibling microbiome score from gut (type=Fecal) and airways (type = Trach) at various ages (time) analyzed by quasi-poisson regression. Model indicates type of analysis: unadjusted (all crude), adjusted for current siblings (all adjusted), stratified by siblings (siblings and no siblings respectively). Supplemental Figure 1: Alpha diversity (Richness), stratified by sample site, time, and siblings. Boxplots demonstrate medians and IQR. P-values determined by Dunn’s test with FDR control for multiple testing. Supplemental Figure 2: Relative abundance of the 6 most abundant phyla (highest mean relative abundance) in each sample site, stratified by sample site, time, and siblings. P-values determined by Wilcoxon test. A pseudocount (+1e-06) was added to all abundances for the log-scale presentation. The black dots indicate median values. Children with younger siblings or both younger and older siblings (only applicable at 4 and 6 years of age) are included in Supplemental Figure 2. Summary statistics shown in Supplemental Table 7 and Supplemental Table 9. Supplemental Figure 3: Relative abundance of the 8 most abundant phyla (highest mean relative abundance) in the gut at 4 and 6 years of age, stratified by time and siblings. P-values determined by Kruskal-Wallis test. A pseudocount (+1e-06) was added to all abundances for the log-scale presentation. The black dots indicate median values. Summary statistics shown in supplemental table 9. Supplemental Figure 4: Relative abundance of the 10 most abundant genera (highest mean relative abundance) in the gut at 4 and 6 years of age, stratified by time and siblings. P-values determined by Kruskal-Wallis test. A pseudocount (+1e-06) was added to all abundances for the log-scale presentation. The black dots indicate median values. Supplemental Figure 5: Alpha diversity (Shannon diversity index), stratified by sample site, time, and either age gap to closest older sibling, number of older siblings or sex of closest older sibling. Boxplots demonstrate medians and IQR. P-values determined by Wilcoxon or Kruskal Wallis test on all categories except ‘None’. Supplemental Figure 6: Airway beta diversity stratified by sample time, and either age gap to closest older sibling, number of older siblings or sex of closest older sibling. PCoA plots of weighted UniFrac distances showing the first two axes. Ellipses demonstrate means ± 1 SD. P-value and r2 determined by PERMANOVA (9,999 permutations) on all categories except ‘None’. Supplemental Figure 7: Gut beta diversity stratified by sample time, and either age gap to closest older sibling, number of older siblings or sex of closest older sibling. PCoA plots of weighted UniFrac distances showing the first two axes. Ellipses demonstrate means ± 1 SD. P-value and r2 determined by PERMANOVA (9,999 permutations) on all categories except ‘None’. Supplemental Figure 8: Upper: Coverage of variation from the global NMDS component representation on the individual timepoints. Lower: 10 fold Cross-validated AUC for classifying siblings based on increasing number of NMDS components using linear discriminant analysis. Supplemental Figure 9: 10 fold Cross-validated AUC for classifying siblings based on increasing number of NMDS components from Bray Curtis (bray), Jaccard (jaccard), Unifrac (uf) and Weighted unifrac (wuf) ordination, using linear discriminant analysis. Supplemental Figure 10: Gut sibling score plotted as a function of age gap to closest older sibling (top) and number of older siblings (bottom), stratified by sample time. Boxplots denote upper and lower quartile. The middle of each boxplot denotes the mean sibling score, with the mean sibling score of children with no older siblings being extrapolated horizontally (dashed line). Age gap truncated at 10 years; number of older siblings truncated at 3. Supplemental Figure 11: Correlation (reflected as R squared values from univariate linear models for alpha div and sibling score and adonis for weighted UniFrac betadiversity) as a function of sibling burden based on closest youngest sibling (closest – blue) and all older siblings (all older – red). Supplemental Figure 12: Age gap to closest older sibling stratified by number of older siblings showing a positive association (Spearman’s rank correlation, rho = 0.21, p < 0.001). Censored at 3 older siblings (n3 older siblings = 17, n4 older siblings = 10, n5 older siblings = 2).

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 http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) 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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Christensen, E.D., Hjelmsø, M.H., Thorsen, J. et al. The developing airway and gut microbiota in early life is influenced by age of older siblings. Microbiome 10, 106 (2022). https://doi.org/10.1186/s40168-022-01305-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40168-022-01305-z

Keywords