Skip to main content

Global analysis of saliva as a source of bacterial genes for insights into human population structure and migration studies

Abstract

Background

The genetic diversity of the human microbiome holds great potential for shedding light on the history of our ancestors. Helicobacter pylori is the most prominent example as its analysis allowed a fine-scale resolution of past migration patterns including some that could not be distinguished using human genetic markers. However studies of H. pylori require stomach biopsies, which severely limits the number of samples that can be analysed. By focussing on the house-keeping gene gdh (coding for the glucose-6-phosphate dehydrogenase), on the virulence gene gtf (coding for the glucosyltransferase) of mitis-streptococci and on the 16S-23S rRNA internal transcribed spacer (ITS) region of the Fusobacterium nucleatum/periodonticum-group we here tested the hypothesis that bacterial genes from human saliva have the potential for distinguishing human populations.

Results

Analysis of 10 individuals from each of seven geographic regions, encompassing Africa, Asia and Europe, revealed that the genes gdh and ITS exhibited the highest number of polymorphic sites (59% and 79%, respectively) and most OTUs (defined at 99% identity) were unique to a given country. In contrast, the gene gtf had the lowest number of polymorphic sites (21%), and most OTUs were shared among countries. Most of the variation in the gdh and ITS genes was explained by the high clonal diversity within individuals (around 80%) followed by inter-individual variation of around 20%, leaving the geographic region as providing virtually no source of sequence variation. Conversely, for gtf the variation within individuals accounted for 32%, between individuals for 57% and among geographic regions for 11%. This geographic signature persisted upon extension of the analysis to four additional locations from the American continent. Pearson correlation analysis, pairwise Fst-cluster analysis as well as UniFrac analyses consistently supported a tree structure in which the European countries clustered tightly together and branched with American countries and South Africa, to the exclusion of Asian countries and the Congo.

Conclusion

This study shows that saliva harbours protein-coding bacterial genes that are geographically structured, and which could potentially be used for addressing previously unresolved human migration events.

Background

The human body and its commensal microbiota represent a highly integrated “superorganism”, which is the result of a long-term co-evolution. It is generally recognized that the structure and diversity of the human microbiome is of medical importance not only at the species level but also at the intra-species level given that bacterial strains or clonal variants thereof can differ greatly in virulence [1],[2]. While clonal types can vary among individuals they may be quite similar among family members [2]. The fact that human history is characterized by numerous migration events with the microbiome being a permanent companion suggests that the global distribution of bacterial clones is likely to be biogeographically structured. The concept of a shared history and co-migration between humans and microbes has been particularly exemplified in the case of Helicobacter pylori, whose genetic diversity is geographically and ethnically well structured [3],[4]. For instance, genetic variation in H. pylori exhibits a signature of a recent African origin followed by clonal spread across Europe and Asia, in excellent agreement with inferences from human DNA markers [3]. In addition to providing independent evidence for previously-inferred human migrations and relationships, studies of H. pylori have helped resolve controversies concerning particular issues. For example, a recent study of H. pylori found strong evidence that Polynesians trace their ancestry specifically to Taiwan [5], a view favored by some studies of the human genetic evidence [6],[7] but not others [8],[9]. Moreover, a study of H. pylori was able to provide insights into the history of two religious communities in India that could not be distinguished based on human DNA markers [10], presumably because of the recent time scale involved. Despite the value and the background knowledge that already exists for H. pylori, this organism has two major drawbacks. Studies of H. pylori require stomach biopsy material, which severely limits the populations that can be sampled. In addition, H. pylori is present only in about 50% of the human population [11], and this prevalence is constantly decreasing [12].

Thus, the identification of potentially informative and easily-accessible species for studies of human population history would be highly valuable. The claim for convenient sampling and for a high number of bacterial specimens instantly draws attention to the oral cavity and especially to saliva, which is easily collected non-invasively, without causing pain, discomfort or embarrassment to the individual. In addition, the salivary microbiome constitutes a rich source of hundreds of very diverse bacterial species reflecting the numerous oral niches, such as the tongue, tonsils, throat, hard and soft palate, buccal surfaces and gingivae [13]. Furthermore, members of the oral microbiome are stably preserved for millennia in calcified dental plaque [14],[15], which means that the establishment of suitable oral bacteria or their genes out of saliva may also support the targeted analysis of historic and prehistoric dental calculus of our ancestors. Lastly, saliva is increasingly preferred for sampling as a source of human DNA for epidemiologic and population genetic studies [16].

Whether or not the genetic diversity of oral species may allow distinction of human populations strongly depends on two requirements: first is the degree by which the bacteria are transmitted vertically (i.e. from parents to their children), and such evidence is provided by a number of studies [17]–[20]; second is the degree by which variations arising from mutation are confined to the descendants of the bacterial cell in which they occurred (such a so-called clonal structure would require that genetic exchange is sufficiently rare so that it does not affect population structure). In fact, evidence of a clonal structure and associations between oral bacterial strains and their human host populations exist for a number of species, e.g. streptococci [21], Fusobacterium nucleatum[22], Aggregatibacter actinomycetemcomitans[23], and Haemophilus influenzae[24].

However as with H. pylori, epidemiological data from the above oral species stem generally from clinical isolates grown in routine laboratories for diagnostic purposes. This means that each patient or individual is represented by a single bacterial isolate (i.e. single clone type) as usually routine microbial diagnosis does not account for multiple clones of an identified bacterial species. Given that bacterial strains are often split into numerous clones within the same individual [2], it is legitimate to ask how well a single clone of a bacterial strain represents its host. In addition it is not known (or has not been reported) whether diagnostic laboratories from different areas or countries employ different sampling and cultivation techniques, which may lead to the preferential isolation of certain clonal types over others. Such a possible cultivation bias may contribute to some extent to the genetic variations that has been observed across geographic regions. Therefore, in order to systematically elucidate whether oral bacteria have a biogeographic structure, we here analyzed bacterial genes from the saliva of 120 healthy individuals that were sampled from 12 world-wide regions covering the Northern and South American continent, Europe, Asia and Africa [25]. Based on a non-culture strategy (i.e. PCR amplification directly from salivary DNA-extracts followed by cloning and sequencing) our aim was to analyze and compare gene variation of multiple clones per individual. Instead of a combined analysis of a set of 6 to 8 house keeping genes, as is usually performed with multi-locus sequence typing (MLST) for characterizing bacterial isolates [26], we examined genetic variation within single genes, thereby accounting for the multiple clonal types that exist within individuals. For this purpose we selected three different gene types that are likely to undergo different mutation rates, and consequently may exhibit different discriminatory power among bacterial strains. With a focus on the mitis-streptococci group we used a partial stretch of the “housekeeping gene” gdh, encoding for glucose-6-phosphate dehydrogenase [27],[28], and of the virulence gene gtf, encoding the glucosyltransferase involved in dental biofilm formation [29]. In addition, we analyzed the 16S-23S rRNA internal transcribed spacer region (ITS) of Fusobacterium nucleatum/periodonticum[30],[31]. The above species were initially selected because they are ubiquitous and easily detectable in human saliva, and distinguishable down to sub-species and clonal-cluster level [2],[22]. Furthermore, these species are either pioneers of biofilm formation (i.e. mitis-streptococci), or important bridging organisms at more advanced biofilm stages (i.e. F. nucleatum) which means that there might exist a more or less tight evolutionary relationship between these commensals and their host. While the gdh gene is generally under stabilizing selection for conserving important metabolic functions [32], the ITS may evolve largely neutrally, as mutations at most nucleotides sites lead to only low functional constraints of the transcribed molecule [33]. The overall evolutionary pressure that acts upon the virulence gene gtf may range from stabilizing to positive selection depending on the individual site (e.g. its catalytic domain has recently been reported to be under stabilizing selection [34]). However we restricted the analysis to a rather non-conserved region [35],[36], which may also evolve largely under neutrality. In this pilot study we show that sequence data sets generated from oral bacterial genes differ in their potential to display human population relationships, with gtf exhibiting a geographic signature comparable to the genetic structure of human populations from different continents.

Methods

DNA samples, PCR amplification, cloning and sequencing

DNA-extracts from saliva samples were available from a previous study, in which sample details are described [25]. Briefly, saliva samples encompassed 10 individuals from each of 12 locations and were collected with informed consent. The criteria for selection of individuals were that individuals should not be from the same household, should not exhibit or report any oral diseases, and they should not have been away from the sampled location within the past two months [25]. The 12 geographic regions were: Argentina (AR), Bolivia (BO), California (CA), China (CH), Congo (CO), Georgia (GE), Germany (DE), Louisiana (LO), Philippines (PH), Poland (PO), South Africa (SO), and Turkey (TU), Additional file 1: Figure S1.

The three target genes (gdh, gtf, and ITS) were PCR-amplified directly from DNA-extracts with primers and thermal profiles specified in Table 1. The PCR assay for gdh (encoding the glucose-6-phosphate dehydrogenase gene) covered the three species, Streptococcus oralis, and two closely related species S. mitis, and S. infantis[2]. The PCR assay for gtf (encoding the glucosyltransferase gene) was specific to the gtfR of S. oralis[35], and the PCR assay for the 16S-23S rRNA internal transcribed spacer region (ITS) covered Fusobacterium nucleatum and the closely related F. periodonticum. For amplification of the ITS the forward primer (targeting a specific region of the 16S rRNA gene) was designed by the “probe design” tool of the ARB software package [37]. The reverse primer (targeting a specific region of the 23S rRNA gene) was designed manually based on a conserved region of publicly available 23S rRNA gene sequences. The sensitivity and specificity of both primers was verified separately through “test probe” using SILVA [38].

Table 1 Primer description

PCR amplification products were visualized by electrophoresis on a 1% agarose gel supplemented with ethidium bromide. PCR products were subsequently ligated into the pCR 4-TOPO cloning vector (Invitrogen) and cloned into Escherichia coli One Shot chemocompetent cells (Invitrogen). Cloning and transformation reactions were performed as described by the supplier. Ten to 15 clones were picked at random from each library and bi-directionally sequenced with the plasmid primers M13.

Data analysis

Raw sequences were quality checked manually using the program VectorNT. Sequence alignment was performed with MAFFT [39]. For the analysis of the ITS the co-amplified 16S rRNA and the 23S rRNA gene region were initially used for establishment of a proper alignment and subsequently removed for further analysis. Determination of the start and end points of the ITS region was performed as described previously [30]. ITS-sequences containing tRNAs were not used for further analysis. Sequence clustering into unique sequence types and definition of operational taxonomic units (OTUs) based on a 99%-identity level was performed with mothur [40]). Coverage of OTUs per geographic location was assessed based on the nonparametric estimator of total richness using Chao1 [41]. Pair-wise Fst values as well as the analysis of molecular variance (AMOVA) were calculated in Arlequin (version 3.0) under default settings [42]. Resulting distance matrices were then exported into MEGA 5 [43], and used for Cluster analysis (UPGMA or neighbor-joining). Based on the OTU-distribution the Pearson distance correlations were calculated and resulting distance matrices used for cluster analysis with MEGA 5. Sequence information was also used for phylogenetic tree reconstruction (UPGMA or neighbor-joining) based on MEGA 5 and interpretation of phylogenetic trees (Cluster analysis and principal coordinate analysis) performed with FastUniFrac [44]. Haplotype network analysis was performed with Network, version 4.612 [45]. Correspondence analysis was performed with the “ca”-package in R in analogy to a recent study in which the genetic diversity of the human polyomavirus JC was used for reconstructing human migrations [46]. Sequences determined in this study have been deposited to Genbank and are available under the acc.-no. KM074060 - KM075227 (gtf), KM075228 - KM075859 (gdh), and KM075860 - KM076596 (ITS).

This study was approved by the Ethics Commission of the University of Aachen, Medical Faculty, reference number EK167/09.

Results

Saliva samples from ten individuals each from seven countries (i.e. CH, CO, DE, GE PH, PO, TU) were initially selected for analysis. For all three target genes PCR amplicons of correct fragment size were obtained in the majority of cases (i.e. 65/70 for gdh, 64/70 for gtf, and 68/70 for ITS). Target specificity was verified through comparison with the GenBank database and based on tree reconstructions using publicly available reference sequences (data not shown). Generally 9 to 10 individuals tested positive in each country, except for Georgia which showed the lowest number of positive cases (e.g. only 6/10 for gtf, Table 2). Approximately 700 sequences were analyzed from generated clone libraries for each of the three genes. An overview of the number of individuals, sequences, OTUs and coverage of diversity based on the Chao 1 Richness Estimator is provided in Table 2. The highest number of unique sequence types and OTUs defined at the 99%-identity level were obtained for the gdh, followed by ITS and gtf. At the OTU level, coverage of genetic diversity per country ranged from 28% to 66% (gdh), from 65% to 100% (gtf), and from 43% to 63% (ITS), respectively.

Table 2 Statistics for the analyzed individuals and sequences

The number of polymorphic sites for each gene was 363, 68, and 200 for gdh, gtf, and ITS, respectively. Compared to the lengths of the analyzed gene region, the proportion of polymorphic sites thus was 59% for gdh, 21% for gtf, and 79% for the ITS.

Generally multiple gene variants were found in each individual, with the majority of people exhibiting five to ten different OTUs in the case of gdh and ITS (Figure 1). This was different for the gtf as most individuals were characterized by only one or two different OTUs. A maximum of four OTUs was found in just four individuals (Figure 1).

Figure 1
figure 1

Comparison of the number of different OTUs per individual among the genes gdh , gtf , and ITS. Results are depicted from the individuals of the seven countries CH, CO, DE, GE PH, PO, TU. CH = China, CO = Congo, DE = Germany, GE = Georgia, PH = Philippines, PO = Poland, TU = Turkey.

Classification of the OTUs according to their occurrence in a single or in multiple countries (i.e. co-occurrence in any two countries, any three countries and so on) resulted in very different curve progressions for the three genes. The left panel of Figure 2 displays the type of OTU (i.e. prevalence across countries), while the right panel indicates the number of sequences for each OTU-type. The majority of OTUs were found to be restricted to one particular country in the case of gdh and ITS (95% and 67%, respectively). For gdh those country-specific OTUs were represented by 544 sequences (86% of all sequences), and for ITS by 164 sequences (22% of all sequences). The remaining 5% of gdh-OTUs were distributed over a maximum of three countries. The curve of the ITS-OTUs was characterized by a less steep gradient, with 41 OTUs (20%) shared between any two countries, and 26 OTUs (13%) shared among multiple countries. Four ITS-OTUs were distributed over all seven countries, and were represented by the majority of sequences (i.e. roughly one third of all sequences). The distribution and sequence representation for gtf was different (Figure 2). Here, only 16 OTUs (52%) were unique to one country and they were represented by only 67 sequences (9% of all sequences). Six OTUs (19%) were shared between any two countries, accounting for 11% of all sequences. 22% of OTUs were present in three, four or five countries and accounted together for 26% of all sequences. There was one gtf-OTU present in 6 countries and one gtf-OTU present in all seven countries. Both of these OTUs were represented by roughly equal numbers of sequences which together made up 53% of all sequences in the gtf-dataset (Figure 2).

Figure 2
figure 2

Distribution of OTUs (left panel) and corresponding number of sequences (right panel) across countries. Results are depicted from the individuals of the seven countries CH, CO, DE, GE PH, PO, TU. CH = China, CO = Congo, DE = Germany, GE = Georgia, PH = Philippines, PO = Poland, TU = Turkey.

The sources of genetic variation as calculated via AMOVA were quite distinctive among the three genes (Table 3). Most variation in gdh and ITS were due to sequence diversity within individuals, followed by sequence variations among individuals with virtually little sequence variation due to geographic origin. Conversely, sequence variation among individuals made up the highest source of variation in the case of gtf. Most importantly, geographic origin accounted for 10.7% of the variance.

Table 3 AMOVA values calculated at three hierarchical levels for seven countries

As the gdh- and ITS-dataset revealed no geographic structure further analysis focused on the gtf, for which five additional locations (i.e. AR, BO, CA, LO, and SO, again 10 individuals each) were completed and jointly analyzed with the former data, resulting in the full data set comprising 12 world-wide regions. This dataset comprised 1,168 gtf-sequences, including 332 unique sequences which could be further classified into 57 distinct OTUs (defined at >99% identity). Recalculation of the AMOVA based on these OTUs revealed a partitioning of genetic variation into 34.5% within individuals, 57.3% among individuals per country and 8.2% among countries.

The global distribution pattern revealed two OTUs (OTU001 and OTU002) as conspicuously dominating in prevalence and abundance, whereas the remaining 55 OTUs were notably less abundant and either scattered across two to six countries or found only in a single country (Figure 3). OTU001 was present in all but one country and accounted for 34% of sequences (Figure 3A) with highest prevalence in DE (9 individuals) and lowest in CH and CO (one individual each), Figure 3B. Likewise, OTU002 was ubiquitously present and accounted for 25% of sequences. Highest prevalence was seen in PH (9 individuals) and lowest in PO (1 individual). The abundance of all other OTUs strongly decreased and so did their prevalence across and within countries. For instance, the abundance of OTU003 and OUT004 accounted for only 5.0% and 4.5% of all sequences, yet they were distributed over 5 and 6 countries, respectively. 36 OTUs were found only in a single country, most of which were represented with very few clone sequences. Exceptions to this were OTU008 (CA), OTU011 (CO), OTU015 (PH), and OTU017 (CH), confined to a single country but represented by a relative large number of clone sequences (i.e. 11 to 24).

Figure 3
figure 3

Abundance distribution of gtf -OTUs across all 12 geographic locations, A) # of sequences, B) # of individuals positive for that OTU. AR = Argentina, BO = Bolivia, CA = California, CH = China, CO = Congo, DE = Germany, GE = Georgia, LO = Louisiana, PH = Philippines, PO = Poland, SO = South Africa, TU = Turkey.

Haplotype Network analysis based on those 28 OTUs, represented by at least three clone sequences, revealed some plausible associations in that OTUs shared by multiple countries were generally linked with OTUs shared by largely the same set of countries, e.g. the pairs OTU001/OTU018, OTU002/OTU004, OTU003/OTU027, OTU004/OTU019 (Figure 4). The two major OTUs (OTU001 and OTU002) both of which were dominated by the American and European countries, were not directly connected to each other but instead were linked to a number of minor OTUs. Most of these minor OTUs lacked sequences from American individuals (e.g. OTU006, 007, 011, 013, and 014). Except for OTU011, all OTUs unique to a given country were only connected via one or two nodes with other OTUs.

Figure 4
figure 4

Network analysis based of 28 gtf -OTUs obtained from 12 geographic locations. AR = Argentina, BO = Bolivia, CA = California, CH = China, CO = Congo, DE = Germany, GE = Georgia, LO = Louisiana, PH = Philippines, PO = Poland, SO = South Africa, TU = Turkey.

Phylogenetic tree reconstruction of those 28 OTUs, represented by at least three clone sequences revealed two major branches, which could be subdivided into 9 sub-clusters (A to I, defined at a >97% identity level; Figure 5A). The two most abundant OTUs (OTU001 and OTU002) group separately (i.e. in sub-clusters A and B), and were each affiliated with an additional four OTUs. While BO and LO were almost entirely composed of sub-cluster A and B, these two sub-clusters made up at least 55% of all other countries (except for GE which was devoid of sub-cluster A; Figure 5B). The remaining sub-clusters varied in abundance among countries. For instance, CH and GE were characterized by a high abundance of sub-cluster C, which in turn was very low in BO, CA, SO, and completely absent in LO. Likewise, sub-cluster D was strongly dominant in CA, moderately present in CO, GE, PH, and PO but otherwise absent in all other countries. Sub-cluster E showed the highest abundance in CO and was otherwise confined in low abundance to CH, GE, PO and SO (Figure 5B). Sub-cluster F, composed of a single OTU (i.e. OTU0015, which was identical to OTU002, but with a 54-nucleotide deletion that maintained the reading frame) was exclusively identified in one individual from PH. In order to rule out a possible PCR-artifact, this sample was re-amplified by PCR and new clones analyzed. Again, the identical sequence with a 54-nucleotide deletion was obtained in the new clone library, confirming the authenticity of this sequence variant. Authenticity of this 54-nucleotide deletion was further supported by the fact that it was observed in another (rare) sequence variant in an individual from SO (not displayed in Figure 5A). Sub-cluster G was only found in GE and TU, sub-cluster H only in AR and TU, and sub-cluster I only in SO (Figure 5B).

Figure 5
figure 5

Gtf -subclusters. (A) Phylogenetic UPGMA tree of predominant gtf-variants defined as OTUs at 99% identitiy level with the sub-clusters A to I defined at 97% identity level. Note, that OTU015 (sub-cluster F) differs from OTU002 by a deletion of 54-nucleotides, which is not visible in the treeing analysis. For further information, see main text. The scale bar represents 2% nucleotide differences; (B) Relative abundance of gtf-sub-clusters (A to I) across all 12 geographic locations. AR = Argentina, BO = Bolivia, CA = California, CH = China, CO = Congo, DE = Germany, GE = Georgia, LO = Louisiana, PH = Philippines, PO = Poland, SO = South Africa, TU = Turkey.

We next mapped the 12 countries based on the distribution pattern of all 57 gtf-OTUs using the Pearson Distance analysis. Two major branches emerged, one displaying the European and American countries along with SO, the other displaying the Asian countries along with GE and CO (Figure 6A). While AR and CA were characterized by distinct lineages within the former branch, the countries BO, LO and SO grouped close together next to the European cluster. In contrast, the Asian countries and CO in the second branch were all separated by very long branches. The overall topology of this cluster tree correlated with the OTU-types displayed in Figure 6B. Countries in the first branch were less characterized by rare OTUs but instead were strongly dominated by OTUs present in at least 11 countries (red square, Figure 6B), except for AR and CA (which explains their somewhat unique positions). The pattern for the countries from the second branch was different in that the ratio of rare and common OTUs was more balanced, with a comparably high number of OTUs confined to a single country (blue square, Figure 6B). Two further patterns emerge in Figure 6B: first, the countries from the American continent (AR, BO, CA, and LO) are comparably depauperate in OTUs, which are distributed over four to six countries; second, the entire gtf-gene dataset was devoid of any OTU distributed across seven, eight, nine or ten countries (Figure 6B).

Figure 6
figure 6

GTF -OTU-distribution. (A) Dendrogram displaying the assignment of 12 geographic regions using the Pearson Distance correlation analysis based on gtf-OTU-distribution; (B) heatplot displaying OTU-types and their distribution across 12 geographic locations. OTU-types are defined according to their presence in only one or multiple countries (white, 0% presence, black, 64% presence). By way of example, the upper six countries are characterized by dominance of OTUs, which are distributed over 11 countries (red square), while the lower six countries are comparably high in OTUs confined to a single geographic location (blue square). AR = Argentina, BO = Bolivia, CA = California, CH = China, CO = Congo, DE = Germany, GE = Georgia, LO = Louisiana, PH = Philippines, PO = Poland, SO = South Africa, TU = Turkey.

Countries were also mapped based on pair-wise Fst-values calculated from nucleic acid and deduced amino acid sequences (Figure 7A and B). Both cluster trees largely resembled each other and mirrored the tree that was obtained with the Pearson distance calculation, except for CA, which grouped together with the Asian/CO-branch when Fst-values were used. The distances among countries were markedly lower in the protein-based tree for the countries BO, DE, LO, PO, SO, TU, and CO, GE, and PH, respectively (Figure 7B).

Figure 7
figure 7

(Upper panel) Neighbor-joining trees based on pairwise Fst-Values using DNA (A) or protein (B) sequence; (Lower panel) Neighbor-joining trees based on UniFrac analysis with DNA (C) or protein (D) sequence. AR = Argentina, BO = Bolivia, CA = California, CH = China, CO = Congo, DE = Germany, GE = Georgia, LO = Louisiana, PH = Philippines, PO = Poland, SO = South Africa, TU = Turkey.

The tree topology based on pair-wise Fst-values was largely confirmed with UniFrac-based clustering of countries calculated based on phylogenetic trees using nucleic-acid and deduced amino-acid sequences (Figure 7C and D). However, in the UniFrac analysis the countries were generally separated by longer branches.

The overall influence of the OTU definition at >99% identity level for the assignment of geographic locations was verified by re-defining OTUs using a 97% and 98% threshold, respectively. The resulting cluster-trees based on pair-wise Fst-values were highly similar to each other and to the trees shown in Figure 7, except for the position of CH which now grouped in the upper branch (Additional file 2: Figure S2).

Correspondence analysis confirmed the tight clustering of the European and American countries, with a moderate proximity of PH to those countries (Figure 8). CH is a clear outlier due to several OTUs, as are CO and GE.

Figure 8
figure 8

Relationship between OTUs and countries calculated with correspondence analysis, based on the data presented in Figure3. AR = Argentina, BO = Bolivia, CA = California, CH = China, CO = Congo, DE = Germany, GE = Georgia, LO = Louisiana, PH = Philippines, PO = Poland, SO = South Africa, TU = Turkey.

We also used UniFrac-distances calculated from phylogenetic trees for displaying the relatedness of individuals based on the variation of the gtf. As can be seen in Figure 9A individuals from any given country did not exclusively group together, however certain tendencies were apparent. For instance, individuals from PH grouped comparably close together and were particularly distinct from GE and most individuals of CH. Individuals from GE grouped also tight together and showed some overlap with CH but no overlap with SO. In this plot the most scattered individuals were those from CO. Individuals from America and Europe were also largely scattered and therefore omitted in Figure 9A for transparency. However, re-assignment of corresponding countries at the continent level showed some separation between individuals from Europe and America (Figure 9B).

Figure 9
figure 9

Two-dimensional principle coordinate analysis based on UniFrac distances. Each colored dot is an individual of the following countries (A) CH, CO, GE, PH, and SO; (B) EU = European countries DE, PO, and TU; AM = American countries/states AR, BO, CA, LO; Dashed circles indicate the tentative grouping of individuals from the same geographic location; AR = Argentina, BO = Bolivia, CA = California, CH = China, CO = Congo, DE = Germany, GE = Georgia, LO = Louisiana, PH = Philippines, PO = Poland, SO = South Africa, TU = Turkey.

Discussion

General considerations

Analyses of the human microbiome have strongly progressed in recent years. Cultivation-independent approaches (particularly metagenomics) have provided insights into the tremendous diversity of microbial species, many of which were previously unknown [47],[48] From a medical perspective it seems indisputable that the type of bacterial species, their composition and interactions affect in many (yet to be elucidated) ways our health or predisposition for various diseases [49],[50]. However, a major hurdle towards an advanced comprehension of the human microbiome is its huge intra- and interpersonal variation as a reflection of age, gender, immune system, physiologic status, diet and hygiene, along with many other cultural or lifestyle factors. This complexity has led to an emerging view for the need for interdisciplinary research linking human microbial ecology with various subfields of anthropology [51],[52]. For instance, comparably little is known so far about microbiome variation among individuals from different geographic or ethnic origins. If parts of the microbiome are structured due to geography or ethnicity this lends itself to a more “applied” research question, namely: can we deduce human population relationships and migration history based on the genetic variations within their microbiomes? This is a worthwhile question given that there is still controversy over how humans have migrated in prehistoric times [53],[54], and given that studies on the genetic architecture of human microbes (in particular on H. pylori) have helped resolve some of these issues [4],[5],[10]. With the consideration of a body habitat more conveniently accessible than the stomach, Nasidze et al. [25] recently provided a comprehensive 16S rRNA gene based survey of the salivary microbiota by comparing healthy individuals from 12 geographic regions world wide. A major conclusion of this study was that variations among individuals are persistently/consistently high independently from the geographic distance of any two compared individuals. Hence, no geographic structure of the salivary microbiome at the genus and species level could be ascertained.

However, the fact that many bacterial species are generally represented by numerous clonal variants indistinguishable by the 16S rRNA gene prompted us to re-investigate the above samples with a focus on protein-coding genes and the ITS-region, which are less conserved than the 16S rRNA gene. Our aim was to find out whether or not clonal variants of oral bacterial species correlate with the geographic region/country of their hosts, and we selected the above samples, because of the apparent lack of a geographic structure at the 16S rRNA gene level. In addition, these samples suited our purpose as they were collected independently from general health or dental status and lifestyle of the individuals, such as gender, age, diet or hygiene etc. We believe this to be important, because if a bacterial gene is applicable at all for inferring human population relationships its geographic signature (if existing) must be robust enough to overcome at least partially individual-specific traits or behavior. Only such a bacterial gene would make large-scale population analyses with thousands of individuals feasible, because with increasing number of individuals (involving family members, such as grandparents, parents and children) the consistent controlling for the above factors becomes more and more difficult if not entirely impossible. While individuals of the currently tested sample set were characterized by multiple clone variants, we did not observe a correlation with geography for the genes gdh and ITS. Conversely, the gene gtf exhibited a geographic signature of around 8% for the 12 countries analyzed, which is a notable value because it compares to the human genome structure from different continents [55]. Hence of the three genes tested one exhibited a promising robustness over individual-specific factors, such as diet and hygiene. However, it should be noted that if we had controlled for these individual-specific factors, we may have found significant geographic signature values for the other two genes analyzed.

The multiple clones of bacterial species that can be found in one individual may reflect neutral variants or may be the result of adaptation to the numerous ecological niches in the upper respiratory tract. However, at a given time point and at one particular oral site only few or even only one clonal type might dominate, while others are below the level detectable by traditional culture techniques [2]. Consequently, saliva, which does not contain its own microflora but rather constitutes a representation of other oral sites [13], may also harbor a limited number of highly abundant clonal variants. With restriction to ten clones per individual we aimed at those predominating sequence types rather than assessing the entire clonal diversity. This is because with the long-term objective of large-scale population studies it would be much more practical to retrieve meaningful information with abundant gene variants rather than having to perform an in-depth sequencing for every individual. In addition, if we consider that intra-individual clonal types are more closely related to each other than those found between individuals, the dominating sequence type per individual should be a fair representative for distinguishing clones obtained from geographically different regions or different ethnic groups. In any case, our approach is still more sensitive than a cultivation-based approach, since most clonal types escape cultivation efforts anyway due to their low abundances [2].

Comparison of gdh, gtf, and ITS

Lack of a geographic structure in the gdh and ITS genes contrasts the fact that more polymorphic sites and ultimately much more sequence variants were available for these genes compared to gtf. For instance, the number of gdh- and ITS-OTUs per country exceeded the number of gtf-OTUs by a factor of four to five (Table 2), which may then result in a higher amount of information for proper distinction of geographic locations. This assumption is supported by Figure 2, where both gdh and ITS suggest that each country is unique owing to the large proportions of OTUs only present in a single country. Nonetheless, despite the higher number of country-specific OTUs, when looking at polymorphic sites the AMOVA revealed that most of the genetic variation of gdh and ITS is largely explained by the high within-individual variation (Table 3), which is also supported by the higher richness in OTUs per individual (Figure 1). Conversely, the proportion of shared OTUs among two or more countries is strikingly higher for gtf (Figure 2), and consequently, it would be fair to assume that the gtf gene tentatively reflects similarities rather than differences among the countries (resulting in a low geographic signal). However, despite the relative high proportion of common OTUs the precise combination (distribution pattern) of shared OTUs (Figure 3) along with a comparably low within-individual genetic variation (Table 3) has apparently led to the preservation of some degree of geographic structure of the gtf gene during its diversification and adaptation to the host.

An explanation for this could be that the glucosyltransferases (GTFs) play critical and complex roles in the development of dental plaque. For instance, GTFs of cariogenic S. mutans not only adsorb to enamel synthesizing glucans but also adsorb to surfaces of other oral microorganisms that do not encode for gtf-genes (e.g. Actinomyces spp.), thereby converting them to glucan producers [56],[57]. Although mainly described for S. mutans it is plausible to assume that the GTFs of other oral streptococci are engaged in similar ecological interactions with the oral microbiota. This means, that the adsorption of GTFs to a multitude of oral bacterial species may only work with the existence of universal binding sites, which has to come along with some degree of conservation of the gtf-gene for proper biofilm formation in the human oral cavity. Recently reported evidence of negative selection acting upon the gtf genes of S. mutans supports this assumption [34].

Another possible explanation for the differential characteristics of the three gene datasets could be that with each gene a different spectrum of species was covered. That is, all gtf-variants analyzed in this study could be assigned to a single species, namely S. oralis, while the gdh-variants covered three species of the mitis-streptococci-group (i.e. S. oralis, S. mitis, and S. infantis). In addition, with ITS we covered F. nucleatum and its close relative F. periodonticum. If closely related bacterial species co-evolved differently with the host their joint analysis would probably blur a geographic signature. However, we think that no such scenario has affected our data interpretation, as we performed the analysis also with a separated, species-specific dataset. To this end publicly available reference sequences from bacterial culture isolates were used for joint phylogenetic tree analysis with our sequence data. This enabled an unambiguous assignment of the gdh-OTUs and ITS-OTUs to a given species (data not shown) and those “Species-specific” OTUs could then be analyzed as single datasets. Based on AMOVA the geographic signature (that is, the “among countries genetic variation”) remained close to zero for the split datasets of gdh and ITS. Although the overall statistical power was reduced in this case (due to the lower number of sequences analyzed per species) we believe that a possible geographic signature of the gdh gene and the ITS may only be revealed when looking at the entire clonal diversity per individual. This, however, would be rather impractical for the reasons stated above. Furthermore, since the gtf gene already “scored best” based on the abundant clonal types, analyzing the full breadth of gtf gene diversity appears to be more promising and work is currently underway characterizing the gtf gene by massive in-depth parallel sequencing.

The geographic signature of gtf

Interestingly, the two most abundant clonal variants of gtf (OTU001 and OTU002) are distributed world-wide, yet in no country with a prevalence of 100 percent (Figure 3). In fact, the prevalence fluctuated from one up to nine positive individuals depending on country. Clearly, these clonal variants may be in fact present in every single individual, but were simply by chance not picked up in the clones that were sequenced. However, lack of detection in all samples from GE, nine individuals each from CH and CO (in the case of OTU001) and in nine individuals from PO (in the case of OTU002) suggests that these clonal variants (even if present) are at least markedly less abundant in individuals from the corresponding countries. This implies that they likely might have been overlooked also by traditional culture techniques in those individuals. Hence, while apparently successfully adapted to the human host as reflected by the world-wide ubiquity, some individual factors may preclude the clonal variants OTU001 and OTU002 to colonize every single person in high abundance. Successful colonization of the human host appears to drop dramatically for all subsequent clonal variants (i.e. OTU003-OTU057, Figure 3). These OTUs are scattered among few countries and only with few individuals positive therein. This indicates on the one hand that these clonal variants are rare in prevalence and on the other hand it indicates that more individuals need to be analyzed in order to see which OTUs are truly unique or shared among countries. However, the haplotype network shows that many shared OTUs are connected with other OTUs also shared by the same or a similar range of countries (Figure 4). Such a distribution pattern would probably not emerge, if the geographic regions were massively undersampled. Interestingly, the haplotype network also mirrors some plausible geographical arrangements in that OTUs predominantly shared among individuals from CH or American individuals tentatively are located at opposite sites of the net (e.g. OTU009, OTU010, OTU012, and OTU017 versus OTU004, OTU008, OTU018, and OTU28).

Even when just looking at OTU distribution (that is without sequence information) the Pearson correlation analysis displayed some plausible/fundamental relationships among countries (Figure 6). European countries form one tight cluster, and group with those historically influenced by Europe (i.e. America and SO) and this branch is clearly separated from CO, GE and Asian countries. Figure 6B clearly explains why, namely the first branch is much more strongly characterized by the pre-dominance of the two ubiquitous OTUs (OTU001 and OTU002, with a prevalence in 11 and 12 countries, respectively) while the second branch reflects an elevated number of rare OTUs. Those relationships among countries based on the distribution pattern of OTUs are also largely confirmed based on precise sequence information using pair-wise Fst-values or UniFrac based interpretations of phylogenetic trees, albeit with some minor alterations in the tree topologies (Figure 7). In particular, UniFrac analysis indicates that underlying phylogenetic trees are not heavily blurred by homologous recombination events which most likely occurred at a local rather than at a global scale, thereby probably even contributing to an increased geographic signature. The fact that the topology of the countries in the cluster-trees remained stable even when OTUs were defined at an identity level of 97% (i.e. subclusters A to I, Figure 5, and Additional file 2: Figure S2) indicates the overall robustness of the sequence data set. Even at this relatively high threshold level, the proportions of the resulting nine subclusters (A to I) are quite distinctive among countries with some variants restricted to one or few countries (Figure 5). The relative low abundance of the country-specific subclusters however again indicates that more than ten individuals per country are needed for a proper comparison and differentiation of countries. Another interesting insight stems from the correspondence analysis (Figure 8), which shows, that CH, CO, and GE are distinct from the European/American countries based on several OTUs. These results are consistent with the OTU abundances (Figure 3). Moreover, the OTUs that characterize these outlying countries are not closely-related to one another (Figure 4, Figure 5), indicating that each outlier is characterized by a set of disparate OTUs, rather than by clonal expansion of closely-related OTUs.

In the current study, we did not test the sequence data for potential chimeric artifacts formed during PCR, because such artifacts are indistinguishable from true sequence variants as a result of homologous recombination. However, it is plausible to assume that such PCR artifacts (if they exist) lead to OTUs that are at best abundant in one particular sample or individual, but not spread as identical OTUs among individuals. This means that those OTUs that occur only in one individual, as displayed in Figure 3, are, if any, the most likely candidates to represent chimeric molecules. Therefore, we performed AMOVA and calculated a cluster-tree based on pair-wise Fst-values without those putative chimeric OTUs. The proportion of the total variance due to differences between geographic regions remained at about 8%, and the Fst-cluster tree was highly identical to those shown in Figure 7 (Additional file 3: Figure S3). Hence, our overall interpretation of the data is not distorted by the possible presence of chimeric molecules. Instead, such artifacts would be expected to increase within-individual genetic diversity, thereby reducing the geographic signal. In addition, the haplotype network (Figure 4) has no reticulations, indicating that there is no detectable recombination in the dataset.

In order to relate our work with data from H. pylori we compared the pairwise Fst-values obtained for the 12 countries (Additional file 4: Table S1) with those previously described for human populations from countries across Europe, Middle East and Asia [58], based on H. pylori. The Fst-values in our study range from 0.002 to 0.38 (median 0.11, n = 66) whereas in the case of H. pylori such values only range from 0.002 to 0.13 (median 0.04, n = 210) [58]. Although mainly different countries were analyzed for H. pylori, the notably higher range of Fst-values observed in our study can be viewed as a reasonable estimation of the overall informational potential of the gtf gene. In addition, three countries in the H. pylori study above coincide with our study, namely Germany, Turkey and the Philippines allowing a direct comparison. And in fact, Fst values based on gtf are about four times higher than those based on H. pylori for the pairs (DE-PH and PH-TU) but, interestingly, two times lower for the pair DE-TU (Table 4). Although more data are needed for valid conclusions the particularly high values observed in conjunction with PH indicates quite distinctive oral bacteria in the individuals from the Philippines.

Table 4 Pairwise Fst-values based on the gtf gene (in brackets) and H. pylori 1

Conclusions

We have shown here that protein-coding bacterial genes from human saliva (such as the gtf gene from S. oralis) exhibit a geographic signature, which could potentially be used for addressing previously unresolved human migration events. Clearly, in the current study only a relatively short fragment (around 330 bp) of the gtf gene was analyzed, which means that a sequence dataset solely based thereon hardly provides sufficient information for proper distinction of human populations. This may, in part explain, why individuals from a country did only group together tentatively (Figure 9). A higher resolution with oral bacterial genes would be achieved through the accumulation of sequence information as it is done with H. pylori, for which usually about 4,500 bp, (seven housekeeping genes and two virulence genes) are routinely analyzed [10]. More information could potentially be obtained through joint analysis of other gtf genes and homologous genes from multiple streptococci. Even non-related genes from other oral bacteria (with a geographic signature) could be used in concert, as long as the mere OTU-distribution pattern (e.g. analyzed with the Pearson correlation) allows meaningful inferences. In other words, once tested successfully as single genes, a set of several oral bacterial genes (related or even non-related) could be assembled and the distribution pattern jointly compared among human individuals. This would be in rough analogy to the concatenation of different genes from microbial isolates in order to compare them via MLST. Owing to the functional richness of the oral microbiome [59] a sufficiently high number of biogeographically relevant genes may eventually be identified for the above purpose. The current study provides proof-of-principle and points towards investigations to establish specific saliva bacteria and their genes as a useful alternative to H. pylori for investigating human population structure and migrations.

Additional files

References

  1. Hohwy J, Reinholdt J, Kilian M: Population dynamics of Streptococcus mitis in its natural habitat. Infect Immun. 2001, 69: 6055-6063. 10.1128/IAI.69.10.6055-6063.2001.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  2. Bek-Thomsen M, Tettelin H, Hance I, Nelson KE, Kilian M: Population diversity and dynamics of Streptococcus mitis, Streptococcus oralis, and Streptococcus infantis in the upper respiratory tracts of adults, determined by a nonculture strategy. Infect Immun. 2008, 76: 1889-1896. 10.1128/IAI.01511-07.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  3. Linz B, Balloux F, Moodley Y, Manica A, Liu H, Roumagnac P, Falush D, Stamer C, Prugnolle F, van der Merwe SW, Yamaoka Y, Graham DY, Perez-Trallero E, Wadstrom T, Suerbaum S, Achtman M: An African origin for the intimate association between humans and Helicobacter pylori. Nature. 2007, 445: 915-918. 10.1038/nature05562.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Falush D, Wirth T, Linz B, Pritchard JK, Stephens M, Kidd M, Blaser MJ, Graham DY, Vacher S, Perez-Perez GI, Yamaoka Y, Megraud F, Otto K, Reichard U, Katzowitsch E, Wang X, Achtman M, Suerbaum S: Traces of human migrations in Helicobacter pylori populations. Science. 2003, 299: 1582-1585. 10.1126/science.1080857.

    Article  PubMed  CAS  Google Scholar 

  5. Moodley Y, Linz B, Yamaoka Y, Windsor HM, Breurec S, Wu JY, Maady A, Bernhoft S, Thiberge JM, Phuanukoonnon S, Jobb G, Siba P, Graham DY, Marshall BJ, Achtman M: The peopling of the Pacific from a bacterial perspective. Science. 2009, 323: 527-530. 10.1126/science.1166083.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  6. Melton T, Peterson R, Redd AJ, Saha N, Sofro AS, Martinson J, Stoneking M: Polynesian genetic affinities with Southeast Asian populations as identified by mtDNA analysis. Am J Hum Genet. 1995, 57: 403-414. 10.1002/ajmg.1320570308.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  7. Trejaut JA, Kivisild T, Loo JH, Lee CL, He CL, Hsu CJ, Lee ZY, Lin M: Traces of archaic mitochondrial lineages persist in Austronesian-speaking Formosan populations. PLoS Biol. 2005, 3: e247-10.1371/journal.pbio.0030247.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Oppenheimer SJ, Richards M: Polynesian origins. Slow boat to Melanesia?. Nature. 2001, 410: 166-167. 10.1038/35065520.

    Article  PubMed  CAS  Google Scholar 

  9. Richards M, Oppenheimer S, Sykes B: mtDNA suggests Polynesian origins in Eastern Indonesia. Am J Hum Genet. 1998, 63: 1234-1236. 10.1086/302043.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  10. Wirth T, Wang X, Linz B, Novick RP, Lum JK, Blaser M, Morelli G, Falush D, Achtman M: Distinguishing human ethnic groups by means of sequences from Helicobacter pylori: lessons from Ladakh. Proc Natl Acad Sci U S A. 2004, 101: 4746-4751. 10.1073/pnas.0306629101.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  11. Suerbaum S, Josenhans C: Helicobacter pylori evolution and phenotypic diversification in a changing host. Nat Rev Microbiol. 2007, 5: 441-452. 10.1038/nrmicro1658.

    Article  PubMed  CAS  Google Scholar 

  12. Blaser MJ, Falkow S: What are the consequences of the disappearing human microbiota?. Nat Rev Microbiol. 2009, 7: 887-894. 10.1038/nrmicro2245.

    Article  PubMed  CAS  Google Scholar 

  13. Segata N, Haake SK, Mannon P, Lemon KP, Waldron L, Gevers D, Huttenhower C, Izard J: Composition of the adult digestive tract bacterial microbiome based on seven mouth surfaces, tonsils, throat and stool samples. Genome Biol. 2012, 13: R42-10.1186/gb-2012-13-6-r42.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  14. Warinner C, Rodrigues JF, Vyas R, Trachsel C, Shved N, Grossmann J, Radini A, Hancock Y, Tito RY, Fiddyment S, Speller C, Hendy J, Charlton S, Luder HU, Salazar-Garcia DC, Eppler E, Seiler R, Hansen LH, Castruita JA, Barkow-Oesterreicher S, Teoh Y, Kelstrup CD, Olsen JV, Nanni P, Kawai T, Willerslev E, von Mering C, Lewis CM, Collins MJ, Gilbert MT, et al: Pathogens and host immunity in the ancient human oral cavity. Nat Genet. 2014, 46: 336-344. 10.1038/ng.2906.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  15. Adler CJ, Dobney K, Weyrich LS, Kaidonis J, Walker AW, Haak W, Bradshaw CJ, Townsend G, Soltysiak A, Alt KW, Parkhill J, Cooper A: Sequencing ancient calcified dental plaque shows changes in oral microbiota with dietary shifts of the Neolithic and Industrial revolutions. Nat Genet. 2013, 45: 450-455. 10.1038/ng.2536.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  16. Quinque D, Kittler R, Kayser M, Stoneking M, Nasidze I: Evaluation of saliva as a source of human DNA for population and association studies. Anal Biochem. 2006, 353: 272-277. 10.1016/j.ab.2006.03.021.

    Article  PubMed  CAS  Google Scholar 

  17. Li Y, Ismail AI, Ge Y, Tellez M, Sohn W: Similarity of bacterial populations in saliva from African-American mother-child dyads. J Clin Microbiol. 2007, 45: 3082-3085. 10.1128/JCM.00771-07.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  18. Li Y, Caufield PW: The fidelity of initial acquisition of mutans streptococci by infants from their mothers. J Dent Res. 1995, 74: 681-685. 10.1177/00220345950740020901.

    Article  PubMed  CAS  Google Scholar 

  19. Van Winkelhoff AJ, Boutaga K: Transmission of periodontal bacteria and models of infection. J Clin Periodontol. 2005, 32 (Suppl 6): 16-27. 10.1111/j.1600-051X.2005.00805.x.

    Article  PubMed  CAS  Google Scholar 

  20. Caufield PW: Tracking human migration patterns through the oral bacterial flora. Clin Microbiol Infect. 2009, 15 (Suppl 1): 37-39. 10.1111/j.1469-0691.2008.02680.x.

    Article  PubMed  Google Scholar 

  21. Caufield PW, Ratanapridakul K, Allen DN, Cutter GR: Plasmid-containing strains of Streptococcus mutans cluster within family and racial cohorts: implications for natural transmission. Infect Immun. 1988, 56: 3216-3220.

    PubMed  CAS  PubMed Central  Google Scholar 

  22. Haraldsson G, Holbrook WP, Kononen E: Clonal persistence of oral Fusobacterium nucleatum in infancy. J Dent Res. 2004, 83: 500-504. 10.1177/154405910408300613.

    Article  PubMed  CAS  Google Scholar 

  23. Haubek D, Poulsen K, Westergaard J, Dahlen G, Kilian M: Highly toxic clone of Actinobacillus actinomycetemcomitans in geographically widespread cases of juvenile periodontitis in adolescents of African origin. J Clin Microbiol. 1996, 34: 1576-1578.

    PubMed  CAS  PubMed Central  Google Scholar 

  24. Musser JM, Kroll JS, Granoff DM, Moxon ER, Brodeur BR, Campos J, Dabernat H, Frederiksen W, Hamel J, Hammond G, Arne Høiby EA, Jonsdottir KE, Kabeer M, Kallings I, Khan WN, Kilian M, Knowles K, Koornhof HJ, Law B, Li KI, Montgomery J, Pattison PE, Piffaretti JC, Takala AK, Thong ML, Wall RA, Ward JI, Selander RK: Global genetic structure and molecular epidemiology of encapsulated Haemophilus influenzae. Rev Infect Dis. 1990, 12: 75-111. 10.1093/clinids/12.1.75.

    Article  PubMed  CAS  Google Scholar 

  25. Nasidze I, Li J, Quinque D, Tang K, Stoneking M: Global diversity in the human salivary microbiome. Genome Res. 2009, 19: 636-643. 10.1101/gr.084616.108.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  26. Urwin R, Maiden MC: Multi-locus sequence typing: a tool for global epidemiology. Trends Microbiol. 2003, 11: 479-487. 10.1016/j.tim.2003.08.006.

    Article  PubMed  CAS  Google Scholar 

  27. Do T, Jolley KA, Maiden MC, Gilbert SC, Clark D, Wade WG, Beighton D: Population structure of Streptococcus oralis. Microbiology. 2009, 155: 2593-2602. 10.1099/mic.0.027284-0.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  28. Ip M, Chi F, Chau SS, Hui M, Tang J, Chan PK: Use of the housekeeping genes, gdh (zwf) and gki, in multilocus sequence typing to differentiate Streptococcus pneumoniae from Streptococcus mitis and Streptococcus oralis. Diagn Microbiol Infect Dis. 2006, 56: 321-324. 10.1016/j.diagmicrobio.2006.04.013.

    Article  PubMed  CAS  Google Scholar 

  29. Monchois V, Willemot RM, Monsan P: Glucansucrases: mechanism of action and structure-function relationships. FEMS Microbiol Rev. 1999, 23: 131-151.

    Article  PubMed  CAS  Google Scholar 

  30. Gonzales-Marin C, Spratt DA, Allaker RP: Maternal oral origin of Fusobacterium nucleatum in adverse pregnancy outcomes as determined using the 16S-23S rRNA gene intergenic transcribed spacer region. J Med Microbiol. 2013, 62: 133-144. 10.1099/jmm.0.049452-0.

    Article  PubMed  CAS  Google Scholar 

  31. Conrads G, Claros MC, Citron DM, Tyrrell KL, Merriam V, Goldstein EJ: 16S-23S rDNA internal transcribed spacer sequences for analysis of the phylogenetic relationships among species of the genus Fusobacterium. Int J Syst Evol Microbiol. 2002, 52: 493-499.

    Article  PubMed  CAS  Google Scholar 

  32. Kiratisin P, Li L, Murray PR, Fischer SH: Use of housekeeping gene sequencing for species identification of viridans streptococci. Diagn Microbiol Infect Dis. 2005, 51: 297-301. 10.1016/j.diagmicrobio.2004.12.001.

    Article  PubMed  CAS  Google Scholar 

  33. Alvarez I, Wendel JF: Ribosomal ITS sequences and plant phylogenetic inference. Mol Phylogenet Evol. 2003, 29: 417-434. 10.1016/S1055-7903(03)00208-2.

    Article  PubMed  CAS  Google Scholar 

  34. Argimon S, Alekseyenko AV, DeSalle R, Caufield PW: Phylogenetic analysis of glucosyltransferases and implications for the coevolution of mutans streptococci with their mammalian hosts. PLoS One. 2013, 8: e56305-10.1371/journal.pone.0056305.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  35. Hoshino T, Kawaguchi M, Shimizu N, Hoshino N, Ooshima T, Fujiwara T: PCR detection and identification of oral streptococci in saliva samples using gtf genes. Diagn Microbiol Infect Dis. 2004, 48: 195-199. 10.1016/j.diagmicrobio.2003.10.002.

    Article  PubMed  CAS  Google Scholar 

  36. Fujiwara T, Hoshino T, Ooshima T, Sobue S, Hamada S: Purification, characterization, and molecular analysis of the gene encoding glucosyltransferase from Streptococcus oralis. Infect Immun. 2000, 68: 2475-2483. 10.1128/IAI.68.5.2475-2483.2000.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  37. Ludwig W, Strunk O, Westram R, Richter L, Meier H, Yadhukumar Buchner A, Lai T, Steppi S, Jobb G, Forster W, Brettske I, Gerber S, Ginhart AW, Gross O, Grumann S, Hermann S, Jost R, Konig A, Liss T, Lussmann R, May M, Nonhoff B, Reichel B, Strehlow R, Stamatakis A, Stuckmann N, Vilbig A, Lenke M, Ludwig T: ARB: a software environment for sequence data. Nucleic Acids Res. 2004, 32: 1363-1371. 10.1093/nar/gkh293.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  38. Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, Peplies J, Glockner FO: The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013, 41: D590-D596. 10.1093/nar/gks1219.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  39. Katoh K, Toh H: Recent developments in the MAFFT multiple sequence alignment program. Brief Bioinform. 2008, 9: 286-298. 10.1093/bib/bbn013.

    Article  PubMed  CAS  Google Scholar 

  40. Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, Lesniewski RA, Oakley BB, Parks DH, Robinson CJ, Sahl JW, Stres B, Thallinger GG, Van Horn DJ, Weber CF: Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009, 75: 7537-7541. 10.1128/AEM.01541-09.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  41. Chao A: Estimating the population size for capture-recapture data with unequal catchability. Biometrics. 1987, 43: 783-791. 10.2307/2531532.

    Article  PubMed  CAS  Google Scholar 

  42. Excoffier L, Laval G, Schneider S: Arlequin (version 3.0): an integrated software package for population genetics data analysis. Evol Bioinform Online. 2005, 1: 47-50.

    CAS  PubMed Central  Google Scholar 

  43. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  44. Lozupone C, Hamady M, Knight R: UniFrac–an online tool for comparing microbial community diversity in a phylogenetic context. BMC Bioinformatics. 2006, 7: 371-10.1186/1471-2105-7-371.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Bandelt HJ, Forster P, Rohl A: Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999, 16: 37-48. 10.1093/oxfordjournals.molbev.a026036.

    Article  PubMed  CAS  Google Scholar 

  46. Pavesi A: Utility of JC polyomavirus in tracing the pattern of human migrations dating to prehistoric times. J Gen Virol. 2005, 86: 1315-1326. 10.1099/vir.0.80650-0.

    Article  PubMed  CAS  Google Scholar 

  47. Koren O, Knights D, Gonzalez A, Waldron L, Segata N, Knight R, Huttenhower C, Ley RE: A guide to enterotypes across the human body: meta-analysis of microbial community structures in human microbiome datasets. PLoS Comput Biol. 2013, 9: e1002863-10.1371/journal.pcbi.1002863.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  48. Gill SR, Pop M, Deboy RT, Eckburg PB, Turnbaugh PJ, Samuel BS, Gordon JI, Relman DA, Fraser-Liggett CM, Nelson KE: Metagenomic analysis of the human distal gut microbiome. Science. 2006, 312: 1355-1359. 10.1126/science.1124234.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  49. Turnbaugh PJ, Ley RE, Hamady M, Fraser-Liggett CM, Knight R, Gordon JI: The human microbiome project. Nature. 2007, 449: 804-810. 10.1038/nature06244.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  50. Schloissnig S, Arumugam M, Sunagawa S, Mitreva M, Tap J, Zhu A, Waller A, Mende DR, Kultima JR, Martin J, Kota K, Sunyaev SR, Weinstock GM, Bork P: Genomic variation landscape of the human gut microbiome. Nature. 2013, 493: 45-50. 10.1038/nature11711.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Benezra A, DeStefano J, Gordon JI: Anthropology of microbes. Proc Natl Acad Sci U S A. 2012, 109: 6378-6381. 10.1073/pnas.1200515109.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  52. Fortenberry JD: The uses of race and ethnicity in human microbiome research. Trends Microbiol. 2013, 21: 165-166. 10.1016/j.tim.2013.01.001.

    Article  PubMed  CAS  Google Scholar 

  53. Mellars P: Going east: new genetic and archaeological perspectives on the modern human colonization of Eurasia. Science. 2006, 313: 796-800. 10.1126/science.1128402.

    Article  PubMed  CAS  Google Scholar 

  54. Stoneking M, Harvarti K: Early Old World migrations of Homo sapiens: human biology. The Encyclopedia of Global Human Migration. Edited by: Ness I. 2013, Blackwell, New York

    Google Scholar 

  55. Romualdi C, Balding D, Nasidze IS, Risch G, Robichaux M, Sherry ST, Stoneking M, Batzer MA, Barbujani G: Patterns of human diversity, within and among continents, inferred from biallelic DNA polymorphisms. Genome Res. 2002, 12: 602-612. 10.1101/gr.214902.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  56. Rolla G, Ciardi JE, Schultz SA: Adsorption of glucosyltransferase to saliva coated hydroxyapatite. Possible mechanism for sucrose dependent bacterial colonization of teeth. Scand J Dent Res. 1983, 91: 112-117.

    PubMed  CAS  Google Scholar 

  57. Bowen WH, Koo H: Biology of Streptococcus mutans-derived glucosyltransferases: role in extracellular matrix formation of cariogenic biofilms. Caries Res. 2011, 45: 69-86. 10.1159/000324598.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  58. Breurec S, Guillard B, Hem S, Brisse S, Dieye FB, Huerre M, Oung C, Raymond J, Tan TS, Thiberge JM, Vong S, Monchy D, Linz B: Evolutionary history of Helicobacter pylori sequences reflect past human migrations in Southeast Asia. PLoS One. 2011, 6: e22058-10.1371/journal.pone.0022058.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  59. Belda-Ferre P, Alcaraz LD, Cabrera-Rubio R, Romero H, Simon-Soro A, Pignatelli M, Mira A: The oral metagenome in health and disease. Isme J. 2012, 6: 46-56. 10.1038/ismej.2011.85.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

Download references

Acknowledgement

This work was supported by grants of the German Research Foundation (DFG) to HPH. MS and JL acknowledge support from the Max Planck Society.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Hans-Peter Horz.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

HPH, MS and GC designed the study. AS, HS, and OK carried out the laboratory work. KH, JL, HPH, and MS analyzed the data. HPH, KH, GC and MS wrote the manuscript. All authors read and approved the manuscript.

Electronic supplementary material

12862_2014_190_MOESM1_ESM.pdf

Additional file 1: Figure S1.: Map of the 12 sampling regions. AR = Argentina, BO = Bolivia, CA = California, CH = China, CO = Congo, DE = Germany, GE = Georgia, LO = Louisiana, PH = Philippines, PO = Poland, SO = South Africa, TU = Turkey. (PDF 257 KB)

12862_2014_190_MOESM2_ESM.pdf

Additional file 2: Figure S2.: Neighbor-joining trees based on Pairwise-Fst-values of the gtf gene variants defined as OTUs at 98% identity level (A) and at 97% identitiy level (B). AR = Argentina, BO = Bolivia, CA = California, CH = China, CO = Congo, DE = Germany, GE = Georgia, LO = Louisiana, PH = Philippines, PO = Poland, SO = South Africa, TU = Turkey. (PDF 23 KB)

12862_2014_190_MOESM3_ESM.pdf

Additional file 3: Figure S3.: Neighbor-joining trees based on Pairwise-Fst-values of the gtf gene variants defined as OTUs at 99% identity level, under omission of all OTUs that were unique to a single individual. This calculation was performed to test the impact of putative chimeras on the tree topology. AR = Argentina, BO = Bolivia, CA = California, CH = China, CO = Congo, DE = Germany, GE = Georgia, LO = Louisiana, PH = Philippines, PO = Poland, SO = South Africa, TU = Turkey. (PDF 10 KB)

12862_2014_190_MOESM4_ESM.pdf

Additional file 4: Table S1.: Pairwise Fst-values based on the gtf gene for the 12 analyzed geographic regions. (PDF 15 KB)

Authors’ original submitted files for images

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 https://creativecommons.org/licenses/by/4.0/.

The Creative Commons Public Domain Dedication waiver (https://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

Henne, K., Li, J., Stoneking, M. et al. Global analysis of saliva as a source of bacterial genes for insights into human population structure and migration studies. BMC Evol Biol 14, 190 (2014). https://doi.org/10.1186/s12862-014-0190-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12862-014-0190-3

Keywords