Cohabiting family members share microbiota with one another and with their dogs

  1. Se Jin Song
  2. Christian Lauber
  3. Elizabeth K Costello
  4. Catherine A Lozupone
  5. Gregory Humphrey
  6. Donna Berg-Lyons
  7. J Gregory Caporaso
  8. Dan Knights
  9. Jose C Clemente
  10. Sara Nakielny
  11. Jeffrey I Gordon
  12. Noah Fierer
  13. Rob Knight  Is a corresponding author
  1. University of Colorado, Boulder, United States
  2. Stanford University School of Medicine, United States
  3. Northern Arizona University, United States
  4. Institute for Genomics and Systems Biology, United States
  5. University of Minnesota, United States
  6. BioTechnology Institute, University of Minnesota, United States
  7. Howard Hughes Medical Institute, University of California, San Francisco, United States
  8. Washington University School of Medicine, United States
  9. Howard Hughes Medical Institute, University of Colorado, Boulder, United States
  10. Biofrontiers Institute, University of Colorado, Boulder, United States
5 figures and 9 tables

Figures

Community similarity within and between families across body sites, and taxa contributing to these differences.

Panels (AD) show average unweighted UniFrac distances between family members (blue) and between members of different families (red). ‘Child’ refers to all offspring aged 3–18 years who cohabit with the parents. ‘Infants’ were considered to be individuals aged 0–12 months. Palm/Paw refers to the right palm in the human comparisons and the back left paw in the dog comparison. Although there are distinguishable differences between the left and right palm communities within and across individuals (Fierer et al., 2008), the same analysis using the left palms showed a similar pattern (Table 2) and neither composition nor diversity were different enough between palms or among the four dog paws to affect the overall patterns. Mean ± 95% CI and R values (ANOSIM) are shown. *p<0.05 and **p<0.001 based on 10,000 permutations. Panel (E) shows the families of bacteria that exhibit the greatest differences in the number of phylotypes (OTUs) shared within and between adult partners on the right palm. Bars represent the average number of shared phylotypes for a given bacterial family within partners from the same family (blue) and between partners of different families (red). Mean ± 95% CI shown. *p<0.05 after Bonferroni correction (Wilcoxon test).

https://doi.org/10.7554/eLife.00458.006
Approach towards or departure from the ‘adult’ state in each body site with age.

(A) Each point represents the average distance (unweighted UniFrac in red; weighted UniFrac in blue) between each participant and all other participants in the ‘adult’ age bracket. Here we define baseline ‘adult’ as 30–45 years in age (the results are not sensitive to this threshold). R2 values (linear regression model) are shown. *p<0.01, **p<0.001. (B) Phylogenetic diversity (PD) of the communities on each body site is plotted for all of the offspring in the study (aged 0–18 years).

https://doi.org/10.7554/eLife.00458.008
Community similarity and phylotype sharing between dogs-owners and their dogs.

The left panel shows the average unweighted UniFrac distance between adult dog-owners and their dogs (blue), between dog-owners and other (not their own) dogs (red), and between adults who do not own dogs and dogs (green). The right panel shows the number of phylotypes shared for the same categories. Comparisons are labeled on the y-axis such that the first body site listed corresponds to the dog and the second site corresponds to the human. Mean ± 95% CI shown. The presence of asterisks lacking brackets indicates that all pairwise comparisons within that group are significant. Generally, dog-owners tend to share more similar communities and more phylotypes with their own dogs than with other dogs. *p<0.05, **p<0.001 after Bonferroni correction (Wilcoxon test).

https://doi.org/10.7554/eLife.00458.011
Alpha diversity and shared phylotypes in couples with and without dogs and children.

The left panels show rarefaction curves for skin communities of couples (including seniors) who have dogs (top, in red), those without dogs (top, in blue), couples (excluding seniors) with infants/children (bottom, in red), and those without infants/children (bottom, in blue). Mean ± 95% CI shown. The right panels show the average number of phylotypes shared among individuals from the same categories shown in the left panels. Mean ± 95% CI shown. *p<0.05, **p<0.001 after Bonferroni correction (Wilcoxon test).

https://doi.org/10.7554/eLife.00458.012
Variation within and between the communities of skin, oral, and fecal samples from humans and dogs.

Panel (A) shows a PCoA plot of all the body habitats, using unweighted UniFrac distances of human and dog samples, rarefied at 5000 sequences/sample. Panels (BD) show select body habitats from the full plot. Panel (E) shows a summary of the taxa shaded by relative abundance at the phylum level broken down by specific body habitat; the seven most abundant taxa are shown in the legend.

https://doi.org/10.7554/eLife.00458.016

Tables

Table 1

Summary of the number, age classification, and gender of humans surveyed in each family and the number and types of animals in families with pets

https://doi.org/10.7554/eLife.00458.003
Adults (Sex)Infants (Sex)Adolescents (Sex)Seniors (Sex)Dogs (breed)Other pets
2 (M, F)1 (F)1 (F)Cat
2 (M, F)1 (Unknown)
2 (M, F)Cats
2 (M, F)1 (F)
2 (M, F)2 (Unknown)
2 (M, F)1 (F)2 (Unknown)
2 (M, F)1 (F)1 (M)
2 (M, F)Cats, guinea pigs
2 (M, F)
2 (M, F)2 (Unknown)
2 (M, F)3 (Unknown)
2 (M, F)1 (Unknown)
2 (M, F)
2 (M, F)1 (F)Cat, fish
2 (M, F)1 (M)Cat
2 (M, F)
2 (M, F)
2 (M, F)Fish
2 (M, F)1 (F)1 (Jack Russell Terrier)Cat
2 (M, F)1 (Australian Cattle mix)
2 (M, F)2 (M, F)Cat, tarantula
2 (M, F)1 (Labrador/Golden mix)Cat
2 (M, F)1 (Springer Spaniel)Cat
2 (M, F)1 (M)1 (M)Cat
2 (M, F)Cats
2 (M, F)Reptiles, amphibians
2 (M, F)Cat
2 (M, F)Cats
2 (M, F)Fish
2 (M, F)1 (F)2 (M, M)
2 (M, F)2 (M, F)Cat
2 (M, F)1 (F)
3 (M, F, M)1 (M)Cat
2 (M, F)1 (Border Collie)
2 (M, F)2 (Kelpie, Standard Poodle)
2 (M, F)1 (F)1 (Border Collie mix)
2 (M, F)1 (M)2 (Boxer, Boxer)
2 (M, F)2 (Labrador, Labrador)
2 (M, F)
2 (M, F)2 (English Setter, Labrador)
2 (M, F)2 (Labrador, Australian Shepherd/Spaniel mix)Chickens
2 (M, F)2 (M, M)Cats, rabbit, reptiles
2 (M, F)1 (F)
2 (M, F)1 (F)1 (German Shepherd mix)
2 (M, F)1 (Bernese Mountain)
2 (M, F)2 (M, F)
2 (M, F)1 (M)1 (M)Cat
2 (M, F)3 (F, M, F)1 (German Shepherd/Malamute mix)
2 (M, F)1 (F)1 (F)
2 (M, F)2 (Australian Shepherd, Australian Shepherd)
2 (M, F)1 (M)2 (Border Collie/German Shepherd mix, Labrador mix)
2 (M, F)1 (M)Cat
1 (F)1 (M)Cat
2 (M, F)1 (F)2 (Siberian Husky, Greater Swiss Mountain)
2 (M, F)
2 (M, F)Cat
2 (M, F)1 (Unknown)*
2 (M, F)
2 (M, F)1 (Unknown)*
2 (M, F)
10612261536*
  1. *

    These dogs were not sampled and thus not included in the total number of dogs.

  2. Each row is one family and the last row contains the total for each column. Infants were considered to be individuals aged 0–12 months, children/adolescents as 1–17 years, adults as 18–59 years and seniors as ≥60 years

Table 2

Summary of the optimal linear mixed model explaining microbial phylogenetic diversity of body sites in relation to the main factors

https://doi.org/10.7554/eLife.00458.004
Age groupBody siteTermTypeEstimateSE%Variability
AllPalms (L&R)NoDogFixed−0.190.06
FamilyRandom38.98
AgeRandom17.14
PlateRandom4.34
FamSizeRandom0.38
LaneRandom0
BSRandom0
Residual39.16
ForeheadFamilyRandom24.85
AgeRandom15.41
LaneRandom8.50
FamSizeRandom7.81
PlateRandom6.33
Residual37.10
FecalAgeRandom45.64
PlateRandom7.15
FamSizeRandom3.08
FamilyRandom2.24
LaneRandom2.97 × 10-9
Residual41.88
OralAgeRandom35.92
FamilyRandom14.54
LaneRandom3.64
PlateRandom4.60 × 10-10
FamSizeRandom0
Residual45.90
AdultsPalms (L&R)NoDogsFixed−0.220.068
MaleFixed−0.110.031
FamilyRandom46.90
PlateRandom2.77
FamSizRandom5.17 × 10-11
LaneRandom2.16 × 10-12
BSRandom0
Residual50.33
ForeheadNoDogsFixed−0.270.087
FamilyRandom34.91
FamSizeRandom9.08
PlateRandom0.80
LaneRandom0
Residual55.21
FecalFamilyRandom20.20
LaneRandom9.77
FamSizeRandom0
PlateRandom0
Residual70.03
OralFamilyRandom26.16
FamSizeRandom14.59
LaneRandom3.43
PlateRandom3.42
Residual52.40
  1. The model takes into account the variability between age groups (Age), families (Family), family sizes (FamSize), sequencing lanes (Lane), and primer plates (Plate). Variability between left and right palms (BS) is also controlled for in the palm model. The table gives parameter estimates and standard errors for the significant terms in the model and the percentage of explained variability for each of the random effects ordered from highest to lowest.

Table 3

Summary of a permutational analysis of variance (PERMANOVA) assessing the effect of family membership on unweighted UniFrac distances between families

https://doi.org/10.7554/eLife.00458.005
Body siteSource of variationDfMSPseudo-FPV
Palms (L&R)Family membership690.32742.07020.0010.21
Residual2100.158150.40
ForeheadFamily membership700.237581.46030.0010.18
Residual910.16270.40
FecalFamily membership660.2121.27590.0010.14
Residual920.166150.41
OralFamily membership680.125441.68030.0010.15
Residual960.074650.27
  1. V are the estimates of the components of variance for the factors. Statistically significant effects are in bold.

Table 4

Summary of ANOSIM analyses of the differences between within-family and between-family community comparisons using unweighted (unwtd) and weighted (wtd) UniFrac

https://doi.org/10.7554/eLife.00458.007
ComparisonBody siteFam (N)Ind (N)R (unwtd)pR (wtd)p
FamiliesForehead601510.49<0.00010.0100.39
Right palm601410.62<0.00010.180.0008
Left palm561220.61<0.00010.30<0.0001
Fecal591510.19<0.00010.0680.056
Oral591550.28<0.00010.32<0.0001
SpousesForehead601140.28<0.00010.0910.0048
Right palm591050.35<0.00010.110.0042
Left palm55910.31<0.00010.130.002
Fecal591140.21<0.00010.0770.0032
Oral591170.18<0.00010.18<0.0001
Father–ChildForehead14310.170.0250.0260.36
Right palm14310.210.00850.110.096
Fecal14330.170.0530.0660.20
Oral14330.070.0690.230.0051
Mother–ChildForehead14310.250.00240.0180.41
Right palm14300.200.0130.00190.51
Fecal14330.100.110.100.095
Oral14330.210.00470.260.0012
Father-InfantForehead12240.250.00340.0300.99
Right Palm12210.160.0160.00630.96
Fecal12230.0570.96−0.100.97
Oral12240.0610.61−0.00860.70
Mother–InfantForehead12230.230.0025−0.0250.99
Right palm12210.34<0.00010.0120.87
Fecal12230.0660.94−0.0360.86
Oral12240.0920.230.0740.42
DogsForehead12220.710.00020.56<0.0001
Back left paw12230.83<0.00010.68<0.0001
Fecal12250.250.00790.200.039
Oral12250.250.0170.300.023
  1. ‘Child’ refers to all offspring aged 3–18 years who cohabit with the parents. ‘Infant’ is considered to be an individual aged 0–12 months. The number of families and individuals used in each analysis is shown. Statistically significant comparisons (p<0.05) are bolded. Those both statistically significant and of relatively large magnitude (R>0.25) are bolded and italicized.

Table 5

Summary of a permutational distance-based linear model testing the effect of age (Age), gender (Sex), pet ownership (Dog or Cat), and family size (FamSize) on unweighted and weighted UniFrac distance (measures of community dissimilarity)

https://doi.org/10.7554/eLife.00458.009
AgeGroupBody siteBest modelAICcPsuedo-Fpr2
UnweightedAllPalms (L&R)(+)Age−453.495.91820.0012.2
(+)Dog−455.944.48370.0011.19
(+)FamSize−454.952.03780.0010.62
Forehead(+)Age−264.312.75790.0011.69
Fecal(+)Age−269.344.37750.0012.71
Oral(+)Age−386.582.47740.0021.5
AdultsPalms (L&R)(+)Dog−284.034.3030.0012.43
Forehead(+)Dog−166.12.38780.0012.33
FecalNoneNANANANA
OralNoneNANANANA
InfantsPalms (L&R)NoneNANANANA
ForeheadNoneNANANANA
FecalNoneNANANANA
OralNoneNANANANA
SeniorsPalms (L&R)NoneNANANANA
ForeheadNoneNANANANA
FecalNoneNANANANA
OralNoneNANANANA
WeightedAllPalms (L&R)(+)Age−727.6623.5970.0019.05
(+)Sex−730.845.2180.0011.57
(+)FamSize−733.394.57950.0020.39
Forehead(+)Age−418.219.4630.00115.11
(+)Sex−423.787.69390.0011.29
(+)FamSize−424.893.16750.0280.22
Fecal(+)Age−353.082.40440.0611.51
Oral(+)FamSize−416.323.40730.0192.05
AdultsPalms (L&R)(+)Sex−462.244.50580.0025.36
(+)FamSize−464.113.91490.0031.09
(+)Dog−465.653.58580.0030.22
Forehead(+)Sex−295.944.36280.0084.18
FecalNoneNANANANA
OralNoneNANANANA
  1. This analysis was performed for the entire data set (All) as well as separately for age groups. The terms from the best model for each body site are shown, along with the percent of total variation explained (r2). Those terms of statistically significant and largest effect are bolded for each body site.

Table 6

Summary of taxon abundances (%) present on the forehead for each age group

https://doi.org/10.7554/eLife.00458.010
TaxonInfantsChildren/AdolescentsAdultsSeniors
Actinobacteria
Propionibacteriaceae* (8.78 × 10-13)6.3125131
Corynebacteriaceae0.72.04.29.3
Micrococcaceae* (0.0072)2.02.71.22.5
Bacteroidetes
Prevotellaceae* (7.02 × 10-9)7.75.61.41.5
Porphyromonadaceae* (8.53 × 10-12)2.83.10.40.8
Flavobacteriaceae1.21.41.12.0
Bacteroidaceae0.51.41.11.9
Firmicutes
Streptococcaceae* (5.32 × 10-37)47275.58.2
Staphylococcaceae* (0.0088)2.12.6124.3
Carnobacteriaceae* (5.99 × 10-24)4.83.40.50.7
Veillonellaceae* (1.25 × 10-13)4.61.90.71.3
Alphaproteobacteria
Sphingomonadaceae0.91.41.30.9
Betaproteobacteria
Neisseriaceae2.44.12.36.9
Comamonadaceae0.41.21.52.9
Gammaproteobacteria
Pasteurellaceae* (1.64 × 10-6)5.37.01.32.4
Moraxellaceae0.72.21.71.8
  1. *

    A significant effect of age (p<0.05 after Bonferroni correction; exact p-values are shown in parentheses). Infants were considered to be individuals aged 0–12 months, children/adolescents as 1–17 years, adults as 18–59 years and seniors as ≥60 years.

  2. Family level abundances of >1% were subjected to ANOVA analysis in QIIME.

Table 7

Summary of a linear mixed effects analysis on the response of bacterial diversity across the body sites using the full data set and then filtered to include just adults

https://doi.org/10.7554/eLife.00458.013
Age groupBody siteFixed effectsRandom effectsModelΔAICPr(Chi)
AllPalms (L&R)Dog + Cat + Sex + KidBS + Ag + Fa + FS + Pl + LaFull0NA
(−)Dog2.910.027
(−)Cat−5.601
(−)Sex1.730.053
(−)Kid−5.621
ForeheadDog + Cat + Sex + KidAg + Fa + FS + Pl + LaFull0NA
(−)Dog0.570.11
(−)Cat−5.251
(−)Sex0.970.085
(−)Kid−4.221
FecalDog + Cat + Sex + KidAg + Fa + FS + Pl + LaFull0NA
(−)Dog−0.0490.16
(−)Cat−0.150.17
(−)Sex0.330.13
(−)Kid−0.130.17
OralDog + Cat + Sex + KidAg + Fa + FS + Pl + LaFull0NA
(−)Dog1.220.072
(−)Cat−1.950.82
(−)Sex−1.880.73
(−)Kid−1.780.64
AdultsPalms (L&R)Dog + Cat + Sex + KidBS + Fa + FS + Pl + LaFull0NA
(−)Dog4.090.014
(−)Cat−5.291
(−)Sex4.110.013
(−)Kid−5.691
ForeheadDog + Cat + Sex + KidFa + FS + Pl + LaFull0NA
(−)Dog2.520.033
(−)Cat−4.991
(−)Sex−0.370.20
(−)Kid−4.381
FecalDog + Cat + Sex + KidFa + FS + Pl + LaFull0NA
(−)Dog0.300.13
(−)Cat0.240.13
(−)Sex−1.440.45
(−)Kid−0.420.21
OralDog + Cat + Sex + KidFa + FS + Pl + LaFull0NA
(−)Dog−0.290.19
(−)Cat−6.501
(−)Sex−7.221
(−)Kid−5.411
  1. The factors tested are co-habitation of dogs (Dog), cats (Cat), children (Kid) and host gender (Sex). The base model takes into account the variability between age groups (Ag), families (Fa), family sizes (FS), sequencing lanes (La), and primer plates (Pl). Variability between left and right palms (BS) is also controlled for in the palm model. The change in model fit resulting from exclusion of each fixed effect based on Akaike's Information Criterion (AIC) is shown. Statistically significant values are in bold text.

Table 8

Summary of taxon abundances (%) on the external body sites of dogs

https://doi.org/10.7554/eLife.00458.014
TaxonBack left pawBack right pawFront left pawPaws averagedForehead
Actinobacteria
Corynebacteriaceae1.81.30.71.21.0
Microbacteriaceae4.04.24.34.21.8
Micrococcaceae2.22.62.42.41.7
Nocardioidaceae3.33.43.23.31.5
Propionibacteriaceae3.03.93.13.34.5
Bacteroidetes
Bacteroidaceae2.11.92.62.23.2
Porphyromonadaceae2.81.61.72.05.7
Prevotellaceae1.11.83.02.02.4
Flavobacteriaceae2.62.52.52.53.3
Flexibacteraceae1.61.71.71.71.5
Firmicutes
Staphylococcaceae1.42.01.01.51.2
Streptococcaceae1.21.61.51.42.8
Lachnospiraceae1.10.91.71.31.2
Veillonellaceae0.30.51.80.90.7
Fusobacteria
Fusobacteriaceae2.31.71.82.03.8
Alphaproteobacteria
Bradyrhizobiaceae0.60.80.50.61.0
Hyphomicrobiaceae1.11.21.11.20.7
Methylobacteriaceae0.40.70.30.51.0
Sphingomonadaceae5.36.65.25.76.1
Betaproteobacteria
Comamonadaceae1.91.91.51.82.6
Oxalobacteraceae1.51.81.51.61.3
Neisseriaceae1.31.01.41.22.8
Gammaproteobacteria
Enterobacteriaceae2.25.15.94.42.5
Oceanospirillaceae2.31.03.12.10.3
Pasteurellaceae2.42.52.32.46.9
Moraxellaceae1.91.41.61.62.2
Pseudomonadaceae7.14.85.55.83.4
  1. Family level abundances >1% are shown. The front right paw was sampled but failed to amplify and is therefore not shown.

Table 9

Summary of taxon abundances (%) present on the palms for each age group

https://doi.org/10.7554/eLife.00458.015
InfantsChildren/ AdolescentsAdultsSeniors
Actinobacteria
Corynebacteriaceae0.7 (0.6)3 (2.3)4.2 (3.9)4.3 (4.0)
Micrococcaceae4.7 (3.9)3.9 (3.9)2.8 (2.6)2.7 (2.8)
Propionibacteriaceae* (0.00016)3.1 (2.2)11 (11)27 (27)20 (15.9)
Bacteroidetes
Bacteroidaceae0.5 (0.6)1.6 (1.4)1.6 (3.3)3.4 (8.7)
Flavobacteriaceae1 (1.4)1.9 (1.4)2.7 (2.0)3.8 (4.1)
Porphyromonadaceae2 (2.1)1.9 (1.5)1.4 (0.7)0.6 (1.0)
Prevotellaceae5.6 (4.2)4.2 (5.7)3.1 (3.0)2.1 (1.7)
Firmicutes
Carnobacteriaceae* (1.15 × 10-11)6.4 (5.8)5 (3.5)1.7 (1.9)1 (0.7)
ClostridialesFamilyXI.IncertaeSedis0.2 (NA)1.5 (NA)1 (NA)1.4 (NA)
Lachnospiraceae* (8.57 × 10-5)0.3 (0.5)0.9 (1.3)1 (1.5)3.1 (7.2)
Lactobacillaceae0.1 (NA)0.2 (NA)1.5 (NA)4.2 (NA)
(Ruminococcaceae*) (1.17 × 10-5)NA (0.2 )NA (0.7)NA (0.8)NA (3.5)
Staphylococcaceae3.2 (2.1)5.1 (7.2)6.7 (7.3)2.8 (2.2)
Streptococcaceae* (3.91 × 10-10)49 (54)27 (26)15 (16)13 (9.1)
Veillonellaceae* (1.76 × 10-7)5.5 (4.5)2.2 (2.1)1.7 (2.0)1.9 (2.1)
Fusobacteria
Fusobacteriaceae1.6 (1.6)1.8 (1.2)1.4 (1.2)1 (1.0)
Betaproteobacteria
Comamonadaceae* (9.90 × 10-5)0.7 (0.3)0.7 (1.0)1.5 (1.7)3.6 (3.7)
Neisseriaceae2.6 (2.4)3.5 (2.3)1.6 (1.1)2 (1.6)
Gammaproteobacteria
Moraxellaceae1.5 (0.7)1.2 (1.7)3.2 (2.7)3.3 (2.7)
Pasteurellaceae* (0.017)2.6 (3.3)4.4 (2.9)1.8 (1.6)1.2 (1.1)
Pseudomonadaceae* (0.049)0.2 (0.2)0.6 (0.7)1.1 (1.0)2.4 (3.0)
(Enterobacteriaceae*) (0.027)NA (0.6)NA (0.8)NA (0.8)NA (2.7)
  1. *

    A significant effect of age (p<0.05 after Bonferroni correction; exact p-values are shown in parentheses).

  2. Shown only for the right palm (left palm showed similar trends). Infants were considered to be individuals aged 0–12 months, children/adolescents as 1–17 years, adults as 18–59 years and seniors as ≥60 years. Abundances for the left palm are shown in parentheses. Family level abundances of greater than 1% were subjected to ANOVA analysis in QIIME. Taxa present at >1% on the left palm but <1% on the right are shown in parentheses.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Se Jin Song
  2. Christian Lauber
  3. Elizabeth K Costello
  4. Catherine A Lozupone
  5. Gregory Humphrey
  6. Donna Berg-Lyons
  7. J Gregory Caporaso
  8. Dan Knights
  9. Jose C Clemente
  10. Sara Nakielny
  11. Jeffrey I Gordon
  12. Noah Fierer
  13. Rob Knight
(2013)
Cohabiting family members share microbiota with one another and with their dogs
eLife 2:e00458.
https://doi.org/10.7554/eLife.00458