The Human Salivary Microbiome Is Shaped by Shared Environment Rather than Genetics: Evidence from a Large Family of Closely Related Individuals
ABSTRACT
The human microbiome is affected by multiple factors, including the environment and host genetics. In this study, we analyzed the salivary microbiomes of an extended family of Ashkenazi Jewish individuals living in several cities and investigated associations with both shared household and host genetic similarities. We found that environmental effects dominated over genetic effects. While there was weak evidence of geographical structuring at the level of cities, we observed a large and significant effect of shared household on microbiome composition, supporting the role of the immediate shared environment in dictating the presence or absence of taxa. This effect was also seen when including adults who had grown up in the same household but moved out prior to the time of sampling, suggesting that the establishment of the salivary microbiome earlier in life may affect its long-term composition. We found weak associations between host genetic relatedness and microbiome dissimilarity when using family pedigrees as proxies for genetic similarity. However, this association disappeared when using more-accurate measures of kinship based on genome-wide genetic markers, indicating that the environment rather than host genetics is the dominant factor affecting the composition of the salivary microbiome in closely related individuals. Our results support the concept that there is a consistent core microbiome conserved across global scales but that small-scale effects due to a shared living environment significantly affect microbial community composition.
IMPORTANCE Previous research shows that the salivary microbiomes of relatives are more similar than those of nonrelatives, but it remains difficult to distinguish the effects of relatedness and shared household environment. Furthermore, pedigree measures may not accurately measure host genetic similarity. In this study, we include genetic relatedness based on genome-wide single nucleotide polymorphisms (SNPs) (rather than pedigree measures) and shared environment in the same analysis. We quantify the relative importance of these factors by studying the salivary microbiomes in members of a large extended Ashkenazi Jewish family living in different locations. We find that host genetics plays no significant role and that the dominant factor is the shared environment at the household level. We also find that this effect appears to persist in individuals who have moved out of the parental household, suggesting that aspects of salivary microbiome composition established during upbringing can persist over a time scale of years.
INTRODUCTION
The human microbiome is the name given to the collected communities of bacteria that live on and in the human body. The oral microbiome is one of the most diverse (1) of any human-associated microbial community (2). The oral microbiome is a causative factor in conditions such as dental caries (3), periodontal disease (4), and halitosis (5) and has also been implicated as a reservoir for infection at other body sites (2) and in the pathogenesis of nonoral diseases, such as inflammatory bowel disease (6). Strictly speaking there is no single “oral microbiome,” as its composition is highly heterogeneous at different sites in the mouth (7, 8), but the term is commonly used to encompass all of these sites. Site-specific microbiomes can be observed in the periodontal sulcus, dental plaque, tongue, buccal mucosa, and saliva (9). The salivary microbiome exhibits long-term stability and can be considered an important reservoir that contains microorganisms from all distinct ecological niches of the oral cavity. Characterizing and understanding the factors defining the composition of the salivary microbiome are thus crucial to understanding the oral microbiome (10, 11).
Some factors that are thought to influence the human microbiome include environment, diet, disease status, and host genetics (12). The relative importance of these factors for the oral microbiome is still under debate, with the majority of previous studies focusing on the gut microbiome (7–9), although it seems reasonable to assume some potential interaction between the salivary microbiome and microbial communities in other parts of the human body, including the intestinal tract (10).
There is evidence that genetically related individuals tend to share more gut microbes than unrelated individuals do, whether or not they are living in the same house at the time of sampling (13, 14). However, the levels of covariation are similar in monozygotic and dizygotic twins, suggesting that a shared early environment may be a more important factor than genetics (13, 15). The effect of cohabitation with direct and frequent contact is greatest when considering the skin microbiome, with a less-evident effect on the gut and salivary microbiomes (11).
There is also evidence that genetic variation is linked to microbiome composition across other body sites, including the mouth (12), with a recent genome-wide association study (GWAS) identifying several human loci associated (P < 5 × 10−8) with microbial taxonomies in the gut microbiome (16). However, no study thus far has incorporated both genetic relatedness (as a continuous variable) and shared environment into the same analysis of the salivary microbiome.
Despite high diversity between individuals, the salivary microbiome appears to have little geographical structure at the genus level at the global scale (17). Nevertheless, at smaller geographical scales, it appears that the environment plays a role in the oral microbiome. Song et al. studied 60 household units and found that the bacterial composition of dorsal tongue bacterial samples was more similar between cohabiting family members than for individuals from different households, with partners and mother-child pairs having significantly more similar communities (18). However, this did not include information on genetic relatedness in addition to family relationships. It appears that household-level differences in the salivary microbiome may also apply to genetically unrelated individuals and nonpartners, with a similar pattern observed in analysis of 24 household pairs of genetically unrelated individuals, only half of whom were considered romantic couples at the time of sampling (19).
The establishment of the oral microbiome appears to proceed rapidly in the first few years of life, with a notable increase in diversity from 0 to 3 years (18), especially after the eruption of teeth (20). The plaque microbiome also appears stable within adult individuals over a period of at least 3 months, with a unique “fingerprint” of oligotypes discernible even within a single bacterial genus (21). Another study indicates that the salivary microbiome is relatively stable over a year, despite measurable effects of interventions like flossing (22). Taken together, these findings suggest the intriguing hypothesis that once a particular oral microbiome is established earlier in life, it can potentially persist for months and perhaps even years, particularly if external factors such as diet remain fixed. If this were true, shared upbringing effects would continue to be detected in the salivary microbiome even after individuals are no longer living in the same household (15).
A recently described large Ashkenazi Jewish family (23) offers an opportunity to investigate the effect of both environment and genetics in closely related individuals. The availability of host genetic data for this cohort means that we can calculate similarity between individuals based on single nucleotide polymorphisms (SNPs), rather than using measures of relatedness from pedigrees that do not precisely correspond to shared genetic content (24). We hypothesized that using this more accurate measure of host genetic similarity could lead to different conclusions about the proportion of shared microbiome composition attributable to genetics compared to previous studies. While, like other studies, we lack information on potential confounders such as diet and lifestyle (17), due to shared cultural practices between members of the ultraorthodox Ashkenazi Jewish community (25), we believe that confounding factors are likely to be more controlled for in this cohort than in others. For this reason, this cohort represents a unique opportunity to compare the salivary microbiome within a large number of individuals living in separate locations but nevertheless sharing a similar diet, lifestyle, and genetic background and to investigate the long-term effect of shared upbringing on salivary microbiome composition.
RESULTS
Description of cohort.
The families analyzed in this study have been already described in detail by Levine et al. (23). All individuals sampled were from the ultraorthodox Ashkenazi Jewish community. Family A comprised more than 800 individuals living in at least eight cities in four countries. Family B comprised more than 200 individuals living in at least four cities in three countries. The unrelated controls were sampled from the same community as the two families. In total, data were generated from samples from 133 individuals in family A, 18 individuals in family B, and 27 controls.
Using minimum entropy decomposition (see Materials and Methods), we found 271 phylotypes in the total microbiome data set, all of which were present when considering just family A. Of these 271 phylotypes, 49 were present in >95% of individuals within family A, with the Firmicutes the most abundant phyla (Fig. 1a) as observed in previous oral microbiome studies (15, 26). The most abundant genera were Streptococcus (30.4%), Rothia (18.5%), Neisseria (17.1%), and Prevotella (17.1%). The composition of samples was similar between the two families (families A and B) and the unrelated controls (Fig. 1b). These groupings had a small but significant effect in an analysis of variance (R2 = 0.015; P < 0.01), but this is typical of comparisons between such large groups that may differ in an unknown number of confounded variables (e.g., diet, genetics, lifestyle). We concluded that family A was at the very least a representative sample capturing the majority of the variation present in the wider Ashkenazi Jewish population, if not also individuals who are not Ashkenazi Jewish (for comparison with Human Microbiome Project data, see Fig. S3 in the supplemental material).
FIG 1
This cohort was originally collected for a study of the genetics of Crohn’s disease (23), and 28 individuals within our sample had a diagnosis of the disease at the time of saliva sample acquisition. We found no significant effect of Crohn’s disease on salivary microbiome composition with an exploratory analysis of variance (R2 = 0.009; P = 0.101; n = 148) accounting for other variables. It was therefore not included as a covariate in further analysis.
Host genetic similarity is weakly correlated with salivary microbiome similarity.
We performed an exploratory analysis on individuals in family A with both genetic and microbiome data available (n = 111) (Fig. S4) and found that genetic kinship was weakly but significantly associated with salivary microbiome dissimilarity computed using Bray-Curtis dissimilarity (Fig. S5a; r = 0.065 by Mantel test; P = 0.001). This analysis does not take into account confounding by the shared environment and therefore sets a probable upper bound on the variation that can be attributed to host genetics. An exploratory analysis of microbiome variation across a subfamily within family A (n = 44) showed that individuals from the same household had a more similar microbiome composition as measured by Bray-Curtis dissimilarity (mean ± standard deviation [SD], 0.623 ± 0.088) compared with individuals from different households (0.652 ± 0.084), and this difference was significant (P < 0.001 by two-sided t test). An exploratory visual representation of this variation showed further clustering by household, although there was large overlap between households in two dimensions (Fig. 2). However, such an analysis is insufficient; household is obviously correlated with variation in host genetics (Fig. S6) because parents tend to live with their children. This emphasizes the need for a quantitative approach looking at the effect of both household and genetics simultaneously as well as other potential confounders.
FIG 2
The approach we chose to use was adonis, which performs a permutational analysis of variance in community composition using a sequential sum-of-squares approach (27). We used Bray-Curtis dissimilarities to quantify differences in salivary microbiome composition between individuals. The following sections present our analyses attempting to quantify the effects of shared environment and genetics. The analysis groups were effectively nested as follows: individuals cohabiting with at least one other individual (n = 26; Table 1), individuals who had cohabited with at least one other individual, either at the time of sampling or beforehand (n = 61; Table 2), individuals living in four different cities who were not necessarily cohabiting with another individual (n = 82; Table S1), and individuals with host genetic information available (n = 111; Table 3).
TABLE 1
| Variableb | City only | Household only | City and householdc | |||
|---|---|---|---|---|---|---|
| R 2 | P | R 2 | P | R 2 | P | |
| Sequencing plate | 0.048 | 0.190 | 0.048 | 0.075 | 0.048 | 0.458 |
| Gender | 0.032 | 0.724 | 0.032 | 0.4 | 0.032 | 0.467 |
| Age | 0.069 | 0.017 | 0.069 | 0.004 | 0.069 | 0.013 |
| MDS1 | 0.031 | 0.757 | 0.031 | 0.537 | 0.031 | 0.727 |
| MDS2 | 0.050 | 0.142 | 0.050 | 0.052 | 0.050 | 0.099 |
| MDS3 | 0.030 | 0.807 | 0.030 | 0.585 | 0.030 | 0.862 |
| MDS4 | 0.049 | 0.162 | 0.049 | 0.054 | 0.049 | 0.097 |
| MDS5 | 0.029 | 0.824 | 0.029 | 0.614 | 0.029 | 0.791 |
| City | 0.080 | 0.400 | 0.080 | 0.178 | ||
| Household | 0.300 | 0.001 | 0.220 | 0.001 | ||
| Residuals | 0.582 | 0.362 | 0.362 | |||
| Total | 1.000 | 1.000 | 1.000 | |||
a
Sharing a household was always significant and explained the most variance of any variable (>18%) even in a model that nests permutations within cities.
b
The order of variables in the model is given by their order in the table. MDS1 to MDS5 are the five axes of metric multidimensional scaling (MDS).
c
Permutations stratified by city in this analysis.
TABLE 2
| Variableb | City only | Household only | City and householdc | |||
|---|---|---|---|---|---|---|
| R 2 | P | R 2 | P | R 2 | P | |
| Sequencing plate | 0.029 | 0.018 | 0.029 | 0.012 | 0.029 | 0.013 |
| Gender | 0.018 | 0.258 | 0.018 | 0.219 | 0.018 | 0.257 |
| Age | 0.038 | 0.002 | 0.038 | 0.001 | 0.038 | 0.002 |
| MDS1 | 0.014 | 0.668 | 0.014 | 0.607 | 0.014 | 0.740 |
| MDS2 | 0.017 | 0.362 | 0.017 | 0.305 | 0.017 | 0.440 |
| MDS3 | 0.020 | 0.173 | 0.020 | 0.141 | 0.020 | 0.263 |
| MDS4 | 0.020 | 0.150 | 0.020 | 0.118 | 0.020 | 0.147 |
| MDS5 | 0.012 | 0.783 | 0.012 | 0.744 | 0.012 | 0.943 |
| City | 0.056 | 0.149 | 0.056 | 0.934 | ||
| Household | 0.239 | 0.021 | 0.183 | 0.044 | ||
| Residuals | 0.777 | 0.594 | 0.594 | |||
| Total | 1.000 | 1.000 | 1.000 | |||
a
Sharing a household was always significant and explains the most variance of any variable (>18%) even in a model that nests permutations within cities.
b
The order of variables in the model is given by their order in the table.
c
Permutations stratified by city in this analysis.
TABLE 3
| Variableb | Pedigree (kinship2) | Genome-wide SNPs (LDAK) | ||
|---|---|---|---|---|
| R 2 | P | R 2 | P | |
| Sequencing plate | 0.028 | <0.001 | 0.028 | <0.001 |
| Gender | 0.011 | 0.094 | 0.011 | 0.096 |
| Age | 0.023 | <0.001 | 0.023 | <0.001 |
| MDS1 | 0.010 | 0.174 | 0.011 | 0.119 |
| MDS2 | 0.007 | 0.706 | 0.010 | 0.231 |
| MDS3 | 0.012 | 0.063 | 0.011 | 0.131 |
| MDS4 | 0.016 | 0.009 | 0.011 | 0.111 |
| MDS5 | 0.009 | 0.325 | 0.007 | 0.617 |
| Parental household | 0.215 | <0.001 | 0.217 | <0.001 |
| Residuals | 0.670 | 0.671 | ||
| Total | 1 | 1 | ||
a
Using pedigree information to produce kinship results in a significant association with human genetics via the fourth MDS axis, which is not present using kinships calculated with LDAK based on genome-wide SNPs.
b
The order of variables in the model is given by their order in the table.
Shared household is the dominant factor affecting salivary microbiome composition.
We performed a permutational analysis of variance on the salivary microbiome dissimilarities for 26 individuals within family A, each of whom lived in a household with at least one other individual in the cohort. At the time of sampling, these cohabiting individuals lived in a total of 16 households in four cities (cities I, II, III, and IV). To account for host genetics, we included axes from a metric multidimensional scaling (MDS) of pairwise genetic distances between individuals as explanatory variables (see Materials and Methods and Fig. S6). The magnitude of the effect of a significant variable is given by the amount of variance explained (R2 in tables).
There was no significant effect of any of the MDS axes, suggesting that host genetics in closely related individuals does not significantly affect microbiome composition. We investigated the effect of environment using two levels of geography: city and household (Table 1). A city-only model showed no significant effect of environment (R2 = 0.08; P = 0.4), whereas a household-only model showed a significant effect (R2 = 0.30; P = 0.001). This was reproduced in a model containing both geographical variables, with permutations stratified by city, where household was still a significant effect (R2 = 0.22; P = 0.001), suggesting that differences at the level of household are more important than at larger geographical scales. We confirmed that city-level effects were small by extending our sample to 82 individuals across the four cities who were not necessarily cohabiting with others (48 individuals for city I, 13 for city II, 12 for city III, and 9 for city IV) and found that city still had a small effect, although it was significant (R2 = 0.053; P < 0.01). In this analysis, we also found no significant effect of genetics, but age was significant (R2 = 0.028; P = 0.0101) (Table S1).
Spouses share taxa at the subgenus level.
Restricting the analysis to only married couples within family A (n = 16; eight couples), shared household explained even more of the variance (R2 = 0.591; P = 0.001). Subtle variations in the relative abundance of phylotypes within the same genus between households were observable, even within the same city location. For example, Leptotrichia phylotypes qualitatively varied consistently between spouse pairs, and these patterns were also seen in children living at home (Fig. 3). Minimum entropy decomposition (MED) phylotype X2772 was present in both spouses in household A1.7 and was also present in the two youngest children within that household (≤10 years old). Similarly, within household A2.4, the Leptotrichia phylotypes of the two children who were ≤10 years old were more similar than the Leptotrichia phylotype of an older child. Quantitatively, repeating the permutational analysis of variance based only on the composition of phylotypes within Leptotrichia showed that spousal pair explained 68.4% of variance, although this was not significant (R2 = 0.684; P = 0.068), and spouses had significantly more similar subgenus phylotype composition than nonspouses did (Fig. S7b). Similar patterns with spouses were also visible in other abundant genera (Fig. S7a), with spouses on average having a significantly more similar subgenus phylotype composition than nonspouses did (mean ± standard error [SE] difference in Bray-Curtis dissimilarity for each genus, −0.048 ± 0.013; Fig. S7b).
FIG 3
Household effects persist in individuals who are no longer cohabiting.
There were an additional 35 individuals who had grown up in a household with at least one other individual present but who no longer lived together at the time of sampling. To see whether the effects of household persisted, we repeated analysis of variance with these individuals included along with the cohabiting individuals (n = 61; Table 2). The effect of household remained significant (R2 = 0.183; P = 0.044), and no axes of human genetic variation were significant (P > 0.05). Age had a significant effect (R2 = 0.038; P < 0.01).
Other variables such as age and sequencing plate had smaller effects than household in all our analyses of variance. However, we chose the order of variables as presented in our tables (Tables 1, 2, and 3) to test for the effect of household after controlling for other variables. Because of the sequential sum-of-squares approach used by adonis, this ordering of variables can have an effect with an unbalanced design. To check that this was not biasing our results and therefore our conclusions about the important factors for salivary microbiome composition, we also investigated the effect of randomly permuting the order of variables in our model formula (see supplemental material). The results of this analysis confirmed that household was always significant (false-discovery rate [q] < 0.05, with Benjamini-Hochberg multiple testing correction), as was age.
Relying on pedigree kinships produces a genetic signal.
To test whether our conclusions required using kinships estimated from genome-wide SNP data for individuals or whether pedigree information was sufficient, we also repeated our analyses using pedigree kinships (see Materials and Methods). Using pedigree kinships resulted in a small but significant amount of variation in microbiome composition being attributable to host genetics via MDS axis 4 (R2 = 0.016; P < 0.01; Table 3).
DISCUSSION
We have conducted, to our knowledge, the first simultaneous investigation of the role of environment and host genetics in shaping the human salivary microbiome in a cohort of closely related individuals within a large Ashkenazi Jewish family. We found a weak correlation between host kinship and salivary microbiome dissimilarity before taking shared household into account and an apparent small but significant effect of genetics when using kinships based on the family pedigree as proxies for genetic similarity. However, when using kinship estimates based on genome-wide SNPs between individuals and simultaneously controlling for shared household with a permutational analysis of variance, we find no support for any clear effect of human genetics, suggesting that shared environment has a much larger effect than genetics and is the dominant factor affecting the salivary microbiome. Typically shared household had an order of magnitude greater effect compared with other significant variables. For example, in our analysis where city was also used as an environmental variable, the variance explained was as follows: household (18.3%), age (3.8%), and sequencing plate (2.9%) (Table 2).
We also observed that younger children living in the same household shared subtle variations in phylotype abundance within genera with their parents (Fig. 3). However, despite a persistence of household effects, it would be wrong to conclude that the salivary microbiome is completely fixed once established, as it clearly has aspects that can change over time. For example, shared household explained more variation for spousal pairs (likely due to frequent contact between them) and that phylotypes observed in younger children and their parents were not seen in older children (likely due to less frequent contact between them). Taken together, these observations support the view that human genetics does not play a major role in shaping the salivary microbiome, at least not in individuals of the same ethnicity, compared to the environment and contact with other individuals.
Our results confirm the seemingly paradoxical situation that the salivary microbiome is largely consistent across global geographical scales but can show large variation between households in the same city. Previous studies have also found evidence of small variations in salivary microbiome composition comparing samples across a global scale (17). As noted previously, this variation could be influenced by differences in environmental or cultural factors, in which case controlling for these differences would decrease the amount of geographical variation. All individuals in our study follow a traditional Ashkenazi Jewish lifestyle and subsequently are thought to share a similar diet and lifestyle regardless of geographical location (25), which may reduce the variation attributable to city-level differences.
The establishment of the oral microbiome early in life may lead to the persistence of a similar composition over several years. The microbial composition of sites within the mouth has been previously observed to be persistent within individuals over periods of months (21) to a year (22), and we see similar strain-level variation between spouses and their young children as observed between individuals by Utter et al. (21) (Fig. 3). Our results indicate the persistence of household effects in individuals no longer cohabiting, suggesting that the salivary microbiome composition established early in life via shared upbringing is able to persist for at least several years. It has been observed that monozygotic twins do not have significantly more similar gut microbiomes than dizygotic twins (13). Stahringer et al. observed the same effect in the salivary microbiome and also found that the salivary microbiomes of twins became less similar as they grew older and ceased cohabiting, concluding that “nurture trumps nature” in the salivary microbiome (15). Our findings from a large number of related individuals rather than twins support this view, including the persistence of shared upbringing effects. Shared upbringing appears to be the dominant factor affecting microbiome composition in both the gut and the mouth, rather than genetic similarity. This may have implications for understanding the familial aggregation of diseases such as inflammatory bowel disease, which has been suggested to have an environmental component (28).
The salivary microbiome appears to be far more resilient to perturbation than the gut microbiome is, with a rapid return to baseline composition after a short course of antibiotics (29). While this could be because of the pharmacokinetics of the antibiotics involved, Zaura et al. speculate that this difference may be due to the salivary microbial ecosystem’s higher intrinsic resilience to stress, as the mouth is subject to more frequent perturbation (30). Our work supporting the dominant role of the environment in affecting salivary microbiome composition suggests that another important factor in long-term persistence may be the regular reseeding of the ecosystem with bacteria from the external environment.
The fact that we reached our conclusion about the lack of effect of genetics only after including kinship based on genome-wide SNP markers casts doubt on the reliability of pedigrees for calculating relatedness. There are several possible reasons for a discrepancy between kinship estimates from pedigrees and allele sharing (24). One possibility is errors in the pedigree, most likely due to extrapair paternities, although this explanation can be ruled out in this data set. More importantly, inherent stochasticity in the Mendelian process of inheritance means that although parents always pass on 50% of their genes to their offspring, SNPs are inherited together in blocks (i.e., haplotypes), meaning that the relatedness between two offspring in a family can be substantially different from 50%. Finally, and most importantly for this closely related population, shallow pedigrees cannot fully capture complex inbreeding patterns. Thus, while pedigrees are a good model for host relatedness in microbiome studies of large randomly mating populations, they should be used with caution in closely related large families like this one.
Limitations.
Because all individuals in our main cohort were members of the same extended Ashkenazi Jewish family, the genetic variation in our data set is therefore much lower than between individuals from a wider population. It is conceivable that host genetics between more distantly related individuals may play a significant role in affecting salivary microbiome composition. However, we note that a recent study of the nasopharyngeal microbiome among Hutterite individuals (a founder population in North America) detected associations between host variation and microbial composition with a similar cohort size (31), demonstrating that limited genetic variation can be associated with the composition of other microbiomes; it may simply be that the salivary microbiome is relatively unaffected by such variation.
Furthermore, our study looked at only overall genetic similarity, assessed using community comparison metrics based on taxon abundances. They therefore do not preclude the existence of fine-scale links between particular microbial taxa and individual genetic loci, particularly in immune-sensing genes such as those identified in the gut microbiome by Bonder et al. using a much larger cohort (32), although our study was not designed or have the statistical power to detect such associations.
Additionally, we lack detailed information on diet and lifestyle factors of individuals in this study. However, the shared cultural practices within this ultraorthodox Ashkenazi Jewish family mean that it is not unreasonable to assume that they share similar lifestyles and diet despite living in different locations around the world (25).
The apparent persistence of shared upbringing could be confounded by the fact that individuals may continue living near the household where they grew up. If this were the case, then our observation could instead be due to the persistence of a shared environment beyond the household at a level intermediate between household and city, rather than the persistence of a stable salivary microbiome following environmental change. Finally, our samples represent only a single cross-sectional snapshot in time. More long-term longitudinal studies like the work of Stahringer et al. on twins (15) are necessary to investigate the persistence of the salivary microbiome after its establishment early in life in a variety of relatedness settings.
Conclusion.
In summary, our results incorporating a measure of genetic relatedness using SNPs demonstrate that the overall composition of the human salivary microbiome in a large Ashkenazi Jewish family is largely influenced by shared environment rather than host genetics. An apparent significant effect of host genetics using pedigree-based estimates disappears when using genetic markers instead, which shows that in future microbiome research, the use of pedigree relatedness as a proxy for host genetic similarity should be done with caution. Geographical structuring occurs to a greater extent at the household level within cities than between cities on different continents. Living in the same household is associated with a more similar salivary microbiome, and this effect persists after individuals have left the household. This is consistent with the long-term persistence of the salivary microbiome composition established earlier in life due to shared upbringing.
MATERIALS AND METHODS
Ethics.
Ethical and research governance approval was provided by the National Research Ethics Service London Surrey Borders Committee and the UCL Research Ethics Committee. Written informed consent was provided by all participants.
Cohort.
Our cohort contained data from 133 individuals within the same extended family (family A) living in four disparate cities (I, II, III, and IV) across three continents (see reference 23 for more information). We also had samples available from 18 individuals from a separate smaller family (family B) and 27 unrelated Ashkenazi Jewish controls. All individuals studied were of genetically confirmed Ashkenazi Jewish ancestry (23, 25). When information was not directly available, shared household was inferred according to age; individuals within this community marry and subsequently leave the family home at a median age of 21 years (95% confidence interval, 19 to 26 years) (25). Therefore, we assumed that individuals aged 18 or younger at the time of sampling were living with their parents and individuals aged 25 or older were not.
For analysis of the effects of household, we included only households with two or more individuals so as to remove the possibility that we were measuring only interindividual differences, which can be large in the salivary microbiome (17, 21). Twenty-six individuals were living with at least one other individual at the time of sampling in a total of nine households. An additional 35 individuals who had grown up in a shared household with at least one other individual in the cohort but who were no longer living together were subsequently included in the analysis.
Sampling.
Saliva samples were collected in sterile tubes containing saliva preservative buffer by the method of Quinque et al. (33). For the full protocol, see the supplemental material. Five hundred milliliters of saliva/preservative buffer was used with PurElute bacterial genomic kit (Edge Biosystems, Gaithersburg, MD) for DNA extraction. After bacterial DNA extraction, three spike DNAs were added to all samples at a final concentration of 4 pg/ml, 0.4 pg/ml, and 0.08 pg/ml, respectively.
PCR amplification, purification, and sequencing.
The Mastermix 16S basic PCR kit containing MolTaq 16S DNA polymerase (Molzym GmbH & Co. KG, Bremen, Germany) was used to generate PCR amplicons. PCR amplicons were purified in two rounds using the Agencourt AMPure system (Beckman Coulter, Inc., Beverly, MA) in a Hamilton StarLet (Hamilton Company, Boston, MA) automated liquid handler. DNA quantitation and quality control were performed using the Agilent 2100 Bioanalyzer system (Agilent Technologies, Santa Clara, CA).
We used 785F (F stands for forward) and 1175R (R stands for reverse) 16S rRNA primers (see Text S1 in the supplemental material) that amplified the V5-V7 region of the 16S rRNA gene on the Illumina MiSeq system (Illumina, San Diego, CA).
Quality control.
To assess technical variation across runs, we spiked samples during library preparation with a fixed amount of synthetic DNA (see supplemental material). Three unique spike sequences (length of 350 bases) that could be easily identified for quality control purposes were designed. We found, as expected, that the number of spike sequences and the number of putative 16S sequences (length between 350 and 380 bases) were negatively correlated with each other due to the limited total sequencing depth of the Illumina MiSeq system (Fig. S1a). The variation in reads corresponding to this spike across samples was independent of run. We also resequenced a subset of samples without spikes to verify whether spikes affected our analyses and observed the same qualitative differences (Fig. S1b), implying that the addition of spikes did not have a negative impact on downstream analysis. Paired-end reads were merged with fastq-mergepairs in VSEARCH v1.11.1 (34), discarding reads with an expected error of >1. As the expected length of the V5-V7 region was 369 bases, we discarded sequences with <350 or >380 bases.
Clustering and taxonomic classification.
Sequences were clustered with minimum entropy decomposition (MED) (35). MED requires that the variation in read depth across samples does not differ by several orders of magnitude, so we discarded samples with fewer than 5,000 reads and subsampled to a maximum number of 20,000 sequences, resulting in 6,353,210 sequences. We ran MED v2.1 with default parameters (see Text S1 in the supplemental material), identifying 271 phylotypes in the data set (Table S2). MED offers higher resolution than operational taxonomic unit (OTU) picking methods do, and it has previously been shown to differentiate the composition of the oral microbiomes of individuals over time even within the same genus in a study of plaque samples (21). We verified that using MED phylotypes gave very similar compositional dissimilarities compared to using de novo OTUs clustered at 98.5% sequence similarity (Fig. S2) but allowed slightly increased statistical power in analysis of variance (see supplemental material), consistent with the literature (35). MED phylotypes had taxonomy assigned using RDP (36) against the Human Oral Microbiome Database (HOMD) (37). Comparison to Human Microbiome Project (HMP) samples from various sites in the mouth also indicated that Ashkenazi Jewish individuals do not have a significantly different oral microbiome from those of other populations, with Ashkenazi Jewish saliva samples clustering with nonplaque samples from individuals in the HMP (Fig. S3). However, the use of different primers makes it difficult to reach a robust conclusion on this point.
Inclusion of host genetics.
We investigated the effect of relatedness between individuals on salivary microbiome composition using both genetic kinships (based on SNPs) and pedigree kinships (based on the pedigree). We calculated pedigree kinships with kinship2 (38) and genetic kinships with LDAK v5.94 (39) using genome-wide SNP data from either the Illumina HumanCytoSNPv12 (Illumina, USA) or the Illumina HumanCoreExome-24, as described previously (23). These genetic kinships kg are normalized to have a mean of zero and correspond to genetic similarity between individuals. kg correlates with the pedigree kinship kp, but there can be substantial spread around the expected values due to the random nature of genetic inheritance (Fig. S5b), making kg a more accurate measure of true genetic similarity between individuals (24). We converted these kinships to dissimilarities and then Euclidean distances (supplemental material) which were used in a multidimensional scaling (MDS) ordination (Fig. S6). Following Blekhman et al. (12), we used MDS with five axes as covariates in a permutational analysis of variance of salivary microbiome dissimilarities.
Statistical analysis.
We calculated Bray-Curtis dissimilarities between samples based on relative abundances of phylotypes, excluding samples with fewer than 1,000 reads. Variance explained in Bray-Curtis dissimilarities was calculated using the adonis function from the vegan v2.4.1 package in R, which performs a permutational analysis of variance of distance matrices (27). We used 9,999 permutations, with permutations stratified by geographical sample location where appropriate.
ACKNOWLEDGMENTS
We thank all the individuals who kindly participated in this study.
This project was supported through the following charities and research councils: EFFORT (Eastman Foundation for Oral Research and Training), the Charles Wolfson Charitable Trust, and the Medical Research Council (MR/L000261/1). A.P.L. was supported by the Irwin Joffe Memorial Fellowship. L.S. was supported by the Engineering and Physical Sciences Research Council (EP/F500351/1). A.L.R.R. was supported by the CAPES Foundation of the Ministry of Education of Brazil (0698130).
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Supplemental Material
File (mbo004173481s1.docx)
- Download
- 40.97 KB
File (mbo004173481sf1.tif)
- Download
- 12.09 MB
File (mbo004173481sf2.tif)
- Download
- 11.59 MB
File (mbo004173481sf3.tif)
- Download
- 13.45 MB
File (mbo004173481sf4.tif)
- Download
- 42.53 MB
File (mbo004173481sf5.tif)
- Download
- 15.21 MB
File (mbo004173481sf6.tif)
- Download
- 6.02 MB
File (mbo004173481sf7.tif)
- Download
- 14.34 MB
File (mbo004173481st1.docx)
- Download
- 48.92 KB
File (mbo004173481st2.xlsx)
- Download
- 295.44 KB
REFERENCES
1.
Human Microbiome Project Consortium. 2012. Structure, function and diversity of the healthy human microbiome. Nature486:207–214.
2.
Wade WG. 2013. The oral microbiome in health and disease. Pharmacol Res69:137–143.
3.
Yang F, Zeng X, Ning K, Liu K-L, Lo C-C, Wang W, Chen J, Wang D, Huang R, Chang X, Chain PS, Xie G, Ling J, and Xu J. 2012. Saliva microbiomes distinguish caries-active from healthy human populations. ISME J6:1–10.
4.
Teles R, Teles F, Frias-Lopez J, Paster B, and Haffajee A. 2013. Lessons learned and unlearned in periodontal microbiology. Periodontol 200062:95–162.
5.
Ren W, Zhang Q, Liu X, Zheng S, Ma L, Chen F, Xu T, and Xu B. 2016. Supragingival plaque microbial community analysis of children with halitosis. J Microbiol Biotechnol26:2141–2147.
6.
Lucas López R, Grande Burgos MJ, Gálvez A, and Pérez Pulido R. 2017. The human gastrointestinal tract and oral microbiota in inflammatory bowel disease: a state of the science review. APMIS125:3–10.
7.
Eren AM, Maignien L, Sul WJ, Murphy LG, Grim SL, Morrison HG, and Sogin ML. 2013. Oligotyping: differentiating between closely related microbial taxa using 16S rRNA gene data. Methods Ecol Evol4:1111–1119.
8.
Mark Welch JL, Rossetti BJ, Rieken CW, Dewhirst FE, and Borisy GG. 2016. Biogeography of a human oral microbiome at the micron scale. Proc Natl Acad Sci U S A113:E791–E800.
9.
Krishnan K, Chen T, and Paster BJ. 2017. A practical guide to the oral microbiome and its relation to health and disease. Oral Dis23:276–286.
10.
Takeshita T, Kageyama S, Furuta M, Tsuboi H, Takeuchi K, Shibata Y, Shimazaki Y, Akifusa S, Ninomiya T, Kiyohara Y, and Yamashita Y. 2016. Bacterial diversity in saliva and oral health-related conditions: the Hisayama Study. Sci Rep6:22164.
11.
Belstrøm D, Holmstrup P, Bardow A, Kokaras A, Fiehn N-E, and Paster BJ. 2016. Temporal stability of the salivary microbiota in oral health. PLoS One11:e0147472.
12.
Blekhman R, Goodrich JK, Huang K, Sun Q, Bukowski R, Bell JT, Spector TD, Keinan A, Ley RE, Gevers D, and Clark AG. 2015. Host genetic variation impacts microbiome composition across human body sites. Genome Biol16:191.
13.
Turnbaugh PJ, Hamady M, Yatsunenko T, Cantarel BL, Duncan A, Ley RE, Sogin ML, Jones WJ, Roe BA, Affourtit JP, Egholm M, Henrissat B, Heath AC, Knight R, and Gordon JI. 2009. A core gut microbiome in obese and lean twins. Nature457:480–484.
14.
Yatsunenko T, Rey FE, Manary MJ, Trehan I, Dominguez-Bello MG, Contreras M, Magris M, Hidalgo G, Baldassano RN, Anokhin AP, Heath AC, Warner B, Reeder J, Kuczynski J, Caporaso JG, Lozupone CA, Lauber C, Clemente JC, Knights D, Knight R, and Gordon JI. 2012. Human gut microbiome viewed across age and geography. Nature486:222–227.
15.
Stahringer SS, Clemente JC, Corley RP, Hewitt J, Knights D, Walters WA, Knight R, and Krauter KS. 2012. Nurture trumps nature in a longitudinal survey of salivary bacterial communities in twins from early adolescence to early adulthood. Genome Res22:2146–2152.
16.
Zhernakova A, Kurilshikov A, Bonder MJ, Tigchelaar EF, Schirmer M, Vatanen T, Mujagic Z, Vila AV, Falony G, Vieira-Silva S, Wang J, Imhann F, Brandsma E, Jankipersadsing SA, Joossens M, Cenit MC, Deelen P, Swertz MA, Weersma RK, Feskens EJM, Netea MG, Gevers D, Jonkers D, Franke L, Aulchenko YS, Huttenhower C, Raes J, Hofker MH, Xavier RJ, Wijmenga C, and Fu J. 2016. Population-based metagenomics analysis reveals markers for gut microbiome composition and diversity. Science352:565–569.
17.
Nasidze I, Li J, Quinque D, Tang K, and Stoneking M. 2009. Global diversity in the human salivary microbiome. Genome Res19:636–643.
18.
Song SJ, Lauber C, Costello EK, Lozupone CA, Humphrey G, Berg-Lyons D, Caporaso JG, Knights D, Clemente JC, Nakielny S, Gordon JI, Fierer N, and Knight R. 2013. Cohabiting family members share microbiota with one another and with their dogs. Elife2:e00458.
19.
Abeles SR, Jones MB, Santiago-Rodriguez TM, Ly M, Klitgord N, Yooseph S, Nelson KE, and Pride DT. 2016. Microbial diversity in individuals and their household contacts following typical antibiotic courses. Microbiome4:39.
20.
Cephas KD, Kim J, Mathai RA, Barry KA, Dowd SE, Meline BS, and Swanson KS. 2011. Comparative analysis of salivary bacterial microbiome diversity in edentulous infants and their mothers or primary care givers using pyrosequencing. PLoS One6:e23503.
21.
Utter DR, Mark Welch JL, and Borisy GG. 2016. Individuality, stability, and variability of the plaque microbiome. Front Microbiol7:564.
22.
David LA, Materna AC, Friedman J, Campos-Baptista MI, Blackburn MC, Perrotta A, Erdman SE, and Alm EJ. 2014. Host lifestyle affects human microbiota on daily timescales. Genome Biol15:R89.
23.
Levine AP Pontikos N Schiff ER Jostins L Speed D NIDDK Inflammatory Bowel Disease Genetics Consortium, Lovat LBBarrett JCGrasberger HPlagnol VSegal AW. 2016. Genetic complexity of Crohn’s disease in two large Ashkenazi Jewish families. Gastroenterology151:698–709.
24.
Speed D and Balding DJ. 2015. Relatedness in the post-genomic era: is it still useful?Nat Rev Genet16:33–44.
25.
Levine A. 2015. The genetics of inflammatory bowel disease in extended multiplex Ashkenazi Jewish kindreds. PhD thesis. University College London, London, United Kingdom.
26.
Dewhirst FE, Chen T, Izard J, Paster BJ, Tanner AC, Yu WH, Lakshmanan A, and Wade WG. 2010. The human oral microbiome. J Bacteriol192:5002–5017.
27.
Oksanen J, Kindt R, Legendre P, O’Hara B, Henry M, and Stevens H. 2016. vegan: community ecology package version 2.4.1.
28.
Nunes T, Fiorino G, Danese S, and Sans M. 2011. Familial aggregation in inflammatory bowel disease: is it genes or environment?World J Gastroenterol17:2715–2722.
29.
Zaura E, Brandt BW, Teixeira de Mattos MJ, Buijs MJ, Caspers MPM, Rashid M-U, Weintraub A, Nord CE, Savell A, Hu Y, Coates AR, Hubank M, Spratt DA, Wilson M, Keijser BJF, and Crielaard W. 2015. Same exposure but two radically different responses to antibiotics: resilience of the salivary microbiome versus long-term microbial shifts in feces. mBio6:e01693-15.
30.
Marsh PD, Head DA, and Devine DA. 2015. Dental plaque as a biofilm and a microbial community—implications for treatment. J Oral Biosci57:185–191.
31.
Igartua C, Davenport ER, Gilad Y, Nicolae DL, Pinto J, and Ober C. 2017. Host genetic variation in mucosal immunity pathways influences the upper airway microbiome. Microbiome5:16.
32.
Bonder MJ, Kurilshikov A, Tigchelaar EF, Mujagic Z, Imhann F, Vila AV, Deelen P, Vatanen T, Schirmer M, Smeekens SP, Zhernakova DV, Jankipersadsing SA, Jaeger M, Oosting M, Cenit MC, Masclee AAM, Swertz MA, Li Y, Kumar V, Joosten L, Harmsen H, Weersma RK, Franke L, Hofker MH, Xavier RJ, Jonkers D, Netea MG, Wijmenga C, Fu J, and Zhernakova A. 2016. The effect of host genetics on the gut microbiome. Nat Genet48:1407–1412.
33.
Quinque D, Kittler R, Kayser M, Stoneking M, and Nasidze I. 2006. Evaluation of saliva as a source of human DNA for population and association studies. Anal Biochem353:272–277.
34.
Rognes T, Flouri T, Nichols B, Quince C, and Mahé F. 2016. VSEARCH: a versatile open source tool for metagenomics. PeerJ4:e2584.
35.
Eren AM, Morrison HG, Lescault PJ, Reveillaud J, Vineis JH, and Sogin ML. 2015. Minimum entropy decomposition: unsupervised oligotyping for sensitive partitioning of high-throughput marker gene sequences. ISME J9:968–979.
36.
Wang Q, Garrity GM, Tiedje JM, and Cole JR. 2007. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol73:5261–5267.
37.
Chen T, Yu W-H, Izard J, Baranova OV, Lakshmanan A, and Dewhirst FE. 2010. The Human Oral Microbiome Database: a web accessible resource for investigating oral microbe taxonomic and genomic information. Database (Oxford)2010:baq013.
38.
Sinnwell JP, Therneau TM, and Schaid DJ. 2014. The kinship2 R package for pedigree data. Hum Hered78:91–93.
39.
Speed D, Hemani G, Johnson MR, and Balding DJ. 2012. Improved heritability estimation from genome-wide SNPs. Am J Hum Genet91:1011–1021.
Information & Contributors
Information
Published In
mBio
Volume 8 • Number 5 • 8 November 2017
eLocator: 10.1128/mbio.01237-17
Editor: Claire M. Fraser, University of Maryland, School of Medicine
Copyright
© 2017 Shaw et al. This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license .
History
Received: 18 July 2017
Accepted: 9 August 2017
Published online: 12 September 2017
Keywords
Contributors
Editor
Claire M. Fraser
Editor
University of Maryland, School of Medicine
Metrics & Citations
Metrics
Note: There is a 3- to 4-day delay in article usage, so article usage will not appear immediately after publication.
Citation counts come from the Crossref Cited by service.
Citations
If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.