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  <cit> , <cit> .
while clonal types can vary among individuals they may be quite similar among family members  <cit> .
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  <cit> , <cit> .
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  <cit> .
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  <cit> , a view favored by some studies of the human genetic evidence  <cit> , <cit>  but not others  <cit> , <cit> .
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  <cit> , 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  <cit> , and this prevalence is constantly decreasing  <cit> .
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  <cit> .
furthermore, members of the oral microbiome are stably preserved for millennia in calcified dental plaque  <cit> , <cit> , 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  <cit> .
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 , and such evidence is provided by a number of studies  <cit> – <cit> ; second is the degree by which variations arising from mutation are confined to the descendants of the bacterial cell in which they occurred .
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  <cit> , fusobacterium nucleatum <cit> , aggregatibacter actinomycetemcomitans <cit> , and haemophilus influenzae <cit> .
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  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  <cit> , it is legitimate to ask how well a single clone of a bacterial strain represents its host.
in addition it is not known  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  <dig> healthy individuals that were sampled from  <dig> world-wide regions covering the northern and south american continent, europe, asia and africa  <cit> .
based on a non-culture strategy  our aim was to analyze and compare gene variation of multiple clones per individual.
instead of a combined analysis of a set of  <dig> to  <dig> house keeping genes, as is usually performed with multi-locus sequence typing  for characterizing bacterial isolates  <cit> , 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  <cit> , <cit> , and of the virulence gene gtf, encoding the glucosyltransferase involved in dental biofilm formation  <cit> .
in addition, we analyzed the 16s-23s rrna internal transcribed spacer region  of fusobacterium nucleatum/periodonticum <cit> , <cit> .
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  <cit> , <cit> .
furthermore, these species are either pioneers of biofilm formation , or important bridging organisms at more advanced biofilm stages  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  <cit> , the its may evolve largely neutrally, as mutations at most nucleotides sites lead to only low functional constraints of the transcribed molecule  <cit> .
the overall evolutionary pressure that acts upon the virulence gene gtf may range from stabilizing to positive selection depending on the individual site .
however we restricted the analysis to a rather non-conserved region  <cit> , <cit> , 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  <cit> .
briefly, saliva samples encompassed  <dig> individuals from each of  <dig> 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  <cit> .
the  <dig> geographic regions were: argentina , bolivia , california , china , congo , georgia , germany , louisiana , philippines , poland , south africa , and turkey , additional file 1: figure s <dig>
the three target genes  were pcr-amplified directly from dna-extracts with primers and thermal profiles specified in table  <dig>  the pcr assay for gdh  covered the three species, streptococcus oralis, and two closely related species s. mitis, and s. infantis <cit> .
the pcr assay for gtf  was specific to the gtfr of s. oralis <cit> , and the pcr assay for the 16s-23s rrna internal transcribed spacer region  covered fusobacterium nucleatum and the closely related f. periodonticum.
for amplification of the its the forward primer  was designed by the “probe design” tool of the arb software package  <cit> .
the reverse primer  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  <cit> .
1
1all pcr-assays started with an initial denaturation at 94°c for 2 min, followed by 35 cycles of the thermal profile as given in the table, with a final elongation phase at 72°c for 15 min .
2for analysis the primer region was omitted in the case of gdh and gtf, and in the case of the its the flanking 16s rrna and 23s rrna gene regions were omitted.
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  and cloned into escherichia coli one shot chemocompetent cells .
cloning and transformation reactions were performed as described by the supplier.
ten to  <dig> clones were picked at random from each library and bi-directionally sequenced with the plasmid primers m <dig>
data analysis
raw sequences were quality checked manually using the program vectornt.
sequence alignment was performed with mafft  <cit> .
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  <cit> .
its-sequences containing trnas were not used for further analysis.
sequence clustering into unique sequence types and definition of operational taxonomic units  based on a 99%-identity level was performed with mothur  <cit> ).
coverage of otus per geographic location was assessed based on the nonparametric estimator of total richness using chao <dig>  <cit> .
pair-wise fst values as well as the analysis of molecular variance  were calculated in arlequin  under default settings  <cit> .
resulting distance matrices were then exported into mega  <dig>  <cit> , and used for cluster analysis .
based on the otu-distribution the pearson distance correlations were calculated and resulting distance matrices used for cluster analysis with mega  <dig>  sequence information was also used for phylogenetic tree reconstruction  based on mega  <dig> and interpretation of phylogenetic trees  performed with fastunifrac  <cit> .
haplotype network analysis was performed with network, version  <dig>   <cit> .
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  <cit> .
sequences determined in this study have been deposited to genbank and are available under the acc.-no.
km <dig> - km <dig> , km <dig> - km <dig> , and km <dig> - km <dig> .
this study was approved by the ethics commission of the university of aachen, medical faculty, reference number ek167/ <dig>
RESULTS
saliva samples from ten individuals each from seven countries  were initially selected for analysis.
for all three target genes pcr amplicons of correct fragment size were obtained in the majority of cases .
target specificity was verified through comparison with the genbank database and based on tree reconstructions using publicly available reference sequences .
generally  <dig> to  <dig> individuals tested positive in each country, except for georgia which showed the lowest number of positive cases .
approximately  <dig> 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  <dig> richness estimator is provided in table  <dig>  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% , from 65% to 100% , and from 43% to 63% , respectively.
gdh
group
1
gtf
group
1
group
1
1ch = china, co = congo, de = germany, ge = georgia, ph = philippines, po = poland, tu = turkey.
the number of polymorphic sites for each gene was  <dig>   <dig>  and  <dig> 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 .
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 .
classification of the otus according to their occurrence in a single or in multiple countries  resulted in very different curve progressions for the three genes.
the left panel of figure  <dig> displays the type of otu , 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 .
for gdh those country-specific otus were represented by  <dig> sequences , and for its by  <dig> 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  <dig> otus  shared between any two countries, and  <dig> otus  shared among multiple countries.
four its-otus were distributed over all seven countries, and were represented by the majority of sequences .
the distribution and sequence representation for gtf was different .
here, only  <dig> otus  were unique to one country and they were represented by only  <dig> sequences .
six otus  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  <dig> 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 .
the sources of genetic variation as calculated via amova were quite distinctive among the three genes .
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  <dig> % of the variance.
gdh
gtf
p
1
a: calculation based on unique sequences.
b: calculation based on otus defined at >99% level.
n.d.: not determined.
1p-values were determined based on permutation analysis using  <dig>  permutations.
as the gdh- and its-dataset revealed no geographic structure further analysis focused on the gtf, for which five additional locations  were completed and jointly analyzed with the former data, resulting in the full data set comprising  <dig> world-wide regions.
this dataset comprised  <dig>  gtf-sequences, including  <dig> unique sequences which could be further classified into  <dig> distinct otus .
recalculation of the amova based on these otus revealed a partitioning of genetic variation into  <dig> % within individuals,  <dig> % among individuals per country and  <dig> % among countries.
the global distribution pattern revealed two otus  as conspicuously dominating in prevalence and abundance, whereas the remaining  <dig> otus were notably less abundant and either scattered across two to six countries or found only in a single country .
otu <dig> was present in all but one country and accounted for 34% of sequences  with highest prevalence in de  and lowest in ch and co , figure 3b.
likewise, otu <dig> was ubiquitously present and accounted for 25% of sequences.
highest prevalence was seen in ph  and lowest in po .
the abundance of all other otus strongly decreased and so did their prevalence across and within countries.
for instance, the abundance of otu <dig> and out <dig> accounted for only  <dig> % and  <dig> % of all sequences, yet they were distributed over  <dig> and  <dig> countries, respectively.
<dig> otus were found only in a single country, most of which were represented with very few clone sequences.
exceptions to this were otu <dig> , otu <dig> , otu <dig> , and otu <dig> , confined to a single country but represented by a relative large number of clone sequences .
haplotype network analysis based on those  <dig> 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/otu <dig>  otu002/otu <dig>  otu003/otu <dig>  otu004/otu <dig> .
the two major otus  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 .
except for otu <dig>  all otus unique to a given country were only connected via one or two nodes with other otus.
phylogenetic tree reconstruction of those  <dig> otus, represented by at least three clone sequences revealed two major branches, which could be subdivided into  <dig> sub-clusters .
the two most abundant otus  group separately , 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 .
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 .
sub-cluster f, composed of a single otu  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  sequence variant in an individual from so .
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 .
we next mapped the  <dig> countries based on the distribution pattern of all  <dig> 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 .
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  <dig> countries , except for ar and ca .
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 .
two further patterns emerge in figure 6b: first, the countries from the american continent  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 .
countries were also mapped based on pair-wise fst-values calculated from nucleic acid and deduced amino acid sequences .
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 .
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 .
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  <dig>  except for the position of ch which now grouped in the upper branch .
correspondence analysis confirmed the tight clustering of the european and american countries, with a moderate proximity of ph to those countries .
ch is a clear outlier due to several otus, as are co and ge.
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 .
discussion
general considerations
analyses of the human microbiome have strongly progressed in recent years.
cultivation-independent approaches  have provided insights into the tremendous diversity of microbial species, many of which were previously unknown  <cit> , <cit>  from a medical perspective it seems indisputable that the type of bacterial species, their composition and interactions affect in many  ways our health or predisposition for various diseases  <cit> , <cit> .
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  <cit> , <cit> .
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  <cit> , <cit> , and given that studies on the genetic architecture of human microbes  have helped resolve some of these issues  <cit> , <cit> , <cit> .
with the consideration of a body habitat more conveniently accessible than the stomach, nasidze et al.
<cit>  recently provided a comprehensive 16s rrna gene based survey of the salivary microbiota by comparing healthy individuals from  <dig> 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  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  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  <dig> countries analyzed, which is a notable value because it compares to the human genome structure from different continents  <cit> .
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  <cit> .
consequently, saliva, which does not contain its own microflora but rather constitutes a representation of other oral sites  <cit> , 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  <cit> .
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 , which may then result in a higher amount of information for proper distinction of geographic locations.
this assumption is supported by figure  <dig>  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 , which is also supported by the higher richness in otus per individual .
conversely, the proportion of shared otus among two or more countries is strikingly higher for gtf , and consequently, it would be fair to assume that the gtf gene tentatively reflects similarities rather than differences among the countries .
however, despite the relative high proportion of common otus the precise combination  of shared otus  along with a comparably low within-individual genetic variation  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  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 , thereby converting them to glucan producers  <cit> , <cit> .
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  <cit> .
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 .
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  and those “species-specific” otus could then be analyzed as single datasets.
based on amova the geographic signature  remained close to zero for the split datasets of gdh and its.
although the overall statistical power was reduced in this case  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  are distributed world-wide, yet in no country with a prevalence of  <dig> percent .
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  and in nine individuals from po  suggests that these clonal variants  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 otu <dig> and otu <dig> to colonize every single person in high abundance.
successful colonization of the human host appears to drop dramatically for all subsequent clonal variants .
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 .
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 .
even when just looking at otu distribution  the pearson correlation analysis displayed some plausible/fundamental relationships among countries .
european countries form one tight cluster, and group with those historically influenced by europe  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  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 .
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%  indicates the overall robustness of the sequence data set.
even at this relatively high threshold level, the proportions of the resulting nine subclusters  are quite distinctive among countries with some variants restricted to one or few countries .
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 , 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 .
moreover, the otus that characterize these outlying countries are not closely-related to one another , 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  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  <dig>  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  <dig> .
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  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  <dig> countries  with those previously described for human populations from countries across europe, middle east and asia  <cit> , based on h. pylori.
the fst-values in our study range from  <dig>  to  <dig>   whereas in the case of h. pylori such values only range from  <dig>  to  <dig>    <cit> .
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  but, interestingly, two times lower for the pair de-tu .
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.
gtf
gene  and
h. pylori
1
1data from h. pylori where taken from  <cit> .
de = germany, ph = philippines, tu = turkey.
CONCLUSIONS
we have shown here that protein-coding bacterial genes from human saliva  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  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 .
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  <dig>  bp,  are routinely analyzed  <cit> .
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  could be used in concert, as long as the mere otu-distribution pattern  allows meaningful inferences.
in other words, once tested successfully as single genes, a set of several oral bacterial genes  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  <cit>  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.
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.
additional files
supplementary material
additional file 1: figure s <dig>
map of the  <dig> 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.
click here for file
additional file 2: figure s <dig>
neighbor-joining trees based on pairwise-fst-values of the gtf gene variants defined as otus at 98% identity level  and at 97% identitiy level .
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.
click here for file
additional file 3: figure s <dig>
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.
click here for file
additional file 4: table s <dig>
pairwise fst-values based on the gtf gene for the  <dig> analyzed geographic regions.
click here for file
acknowledgement
this work was supported by grants of the german research foundation  to hph.
ms and jl acknowledge support from the max planck society.

 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,. while clonal types can vary among individuals they may be quite similar among family members. 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,. 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. 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, a view favored by some studies of the human genetic evidence, but not others,. 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, 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, and this prevalence is constantly decreasing. 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. furthermore, members of the oral microbiome are stably preserved for millennia in calcified dental plaque,, 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. 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, and such evidence is provided by a number of studies –; second is the degree by which variations arising from mutation are confined to the descendants of the bacterial cell in which they occurred. 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, fusobacterium nucleatum, aggregatibacter actinomycetemcomitans, and haemophilus influenzae. 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 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, it is legitimate to ask how well a single clone of a bacterial strain represents its host. in addition it is not known 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 healthy individuals that were sampled from world-wide regions covering the northern and south american continent, europe, asia and africa. based on a non-culture strategy our aim was to analyze and compare gene variation of multiple clones per individual. instead of a combined analysis of a set of to house keeping genes, as is usually performed with multi-locus sequence typing for characterizing bacterial isolates, 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,, and of the virulence gene gtf, encoding the glucosyltransferase involved in dental biofilm formation. in addition, we analyzed the 16s-23s rrna internal transcribed spacer region of fusobacterium nucleatum/periodonticum,. 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,. furthermore, these species are either pioneers of biofilm formation, or important bridging organisms at more advanced biofilm stages 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, the its may evolve largely neutrally, as mutations at most nucleotides sites lead to only low functional constraints of the transcribed molecule. the overall evolutionary pressure that acts upon the virulence gene gtf may range from stabilizing to positive selection depending on the individual site. however we restricted the analysis to a rather non-conserved region,, 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. briefly, saliva samples encompassed individuals from each of 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. the geographic regions were: argentina, bolivia, california, china, congo, georgia, germany, louisiana, philippines, poland, south africa, and turkey, additional file 1: figure s the three target genes were pcr-amplified directly from dna-extracts with primers and thermal profiles specified in table  the pcr assay for gdh covered the three species, streptococcus oralis, and two closely related species s. mitis, and s. infantis. the pcr assay for gtf was specific to the gtfr of s. oralis, and the pcr assay for the 16s-23s rrna internal transcribed spacer region covered fusobacterium nucleatum and the closely related f. periodonticum. for amplification of the its the forward primer was designed by the “probe design” tool of the arb software package. the reverse primer 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. 1 1all pcr-assays started with an initial denaturation at 94°c for 2 min, followed by 35 cycles of the thermal profile as given in the table, with a final elongation phase at 72°c for 15 min. 2for analysis the primer region was omitted in the case of gdh and gtf, and in the case of the its the flanking 16s rrna and 23s rrna gene regions were omitted. 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 and cloned into escherichia coli one shot chemocompetent cells. cloning and transformation reactions were performed as described by the supplier. ten to clones were picked at random from each library and bi-directionally sequenced with the plasmid primers m data analysis raw sequences were quality checked manually using the program vectornt. sequence alignment was performed with mafft. 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. its-sequences containing trnas were not used for further analysis. sequence clustering into unique sequence types and definition of operational taxonomic units based on a 99%-identity level was performed with mothur ). coverage of otus per geographic location was assessed based on the nonparametric estimator of total richness using chao. pair-wise fst values as well as the analysis of molecular variance were calculated in arlequin under default settings. resulting distance matrices were then exported into mega, and used for cluster analysis. based on the otu-distribution the pearson distance correlations were calculated and resulting distance matrices used for cluster analysis with mega sequence information was also used for phylogenetic tree reconstruction based on mega and interpretation of phylogenetic trees performed with fastunifrac. haplotype network analysis was performed with network, version. 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. sequences determined in this study have been deposited to genbank and are available under the acc.-no. km - km, km - km, and km - km. this study was approved by the ethics commission of the university of aachen, medical faculty, reference number ek167/  saliva samples from ten individuals each from seven countries were initially selected for analysis. for all three target genes pcr amplicons of correct fragment size were obtained in the majority of cases. target specificity was verified through comparison with the genbank database and based on tree reconstructions using publicly available reference sequences. generally to individuals tested positive in each country, except for georgia which showed the lowest number of positive cases. approximately 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 richness estimator is provided in table  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%, from 65% to 100%, and from 43% to 63%, respectively. gdh group 1 gtf group 1 group 1 1ch = china, co = congo, de = germany, ge = georgia, ph = philippines, po = poland, tu = turkey. the number of polymorphic sites for each gene was and 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. 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. classification of the otus according to their occurrence in a single or in multiple countries resulted in very different curve progressions for the three genes. the left panel of figure  displays the type of otu, 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. for gdh those country-specific otus were represented by sequences, and for its by 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 otus shared between any two countries, and otus shared among multiple countries. four its-otus were distributed over all seven countries, and were represented by the majority of sequences. the distribution and sequence representation for gtf was different. here, only otus were unique to one country and they were represented by only sequences. six otus 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 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. the sources of genetic variation as calculated via amova were quite distinctive among the three genes. 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 % of the variance. gdh gtf p 1 a: calculation based on unique sequences. b: calculation based on otus defined at >99% level. n.d.: not determined. 1p-values were determined based on permutation analysis using permutations. as the gdh- and its-dataset revealed no geographic structure further analysis focused on the gtf, for which five additional locations were completed and jointly analyzed with the former data, resulting in the full data set comprising world-wide regions. this dataset comprised gtf-sequences, including unique sequences which could be further classified into distinct otus. recalculation of the amova based on these otus revealed a partitioning of genetic variation into % within individuals, % among individuals per country and % among countries. the global distribution pattern revealed two otus as conspicuously dominating in prevalence and abundance, whereas the remaining otus were notably less abundant and either scattered across two to six countries or found only in a single country. otu was present in all but one country and accounted for 34% of sequences with highest prevalence in de and lowest in ch and co, figure 3b. likewise, otu was ubiquitously present and accounted for 25% of sequences. highest prevalence was seen in ph and lowest in po. the abundance of all other otus strongly decreased and so did their prevalence across and within countries. for instance, the abundance of otu and out accounted for only % and % of all sequences, yet they were distributed over and countries, respectively. otus were found only in a single country, most of which were represented with very few clone sequences. exceptions to this were otu, otu, otu, and otu, confined to a single country but represented by a relative large number of clone sequences. haplotype network analysis based on those 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/otu otu002/otu otu003/otu otu004/otu. the two major otus 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. except for otu all otus unique to a given country were only connected via one or two nodes with other otus. phylogenetic tree reconstruction of those otus, represented by at least three clone sequences revealed two major branches, which could be subdivided into sub-clusters. the two most abundant otus group separately, 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. 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. sub-cluster f, composed of a single otu 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 sequence variant in an individual from so. 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. we next mapped the countries based on the distribution pattern of all 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. 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 countries, except for ar and ca. 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. two further patterns emerge in figure 6b: first, the countries from the american continent 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. countries were also mapped based on pair-wise fst-values calculated from nucleic acid and deduced amino acid sequences. 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. 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. 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  except for the position of ch which now grouped in the upper branch. correspondence analysis confirmed the tight clustering of the european and american countries, with a moderate proximity of ph to those countries. ch is a clear outlier due to several otus, as are co and ge. 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. discussion general considerations analyses of the human microbiome have strongly progressed in recent years. cultivation-independent approaches have provided insights into the tremendous diversity of microbial species, many of which were previously unknown, from a medical perspective it seems indisputable that the type of bacterial species, their composition and interactions affect in many ways our health or predisposition for various diseases,. 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,. 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,, and given that studies on the genetic architecture of human microbes have helped resolve some of these issues,,. with the consideration of a body habitat more conveniently accessible than the stomach, nasidze et al. recently provided a comprehensive 16s rrna gene based survey of the salivary microbiota by comparing healthy individuals from 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 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 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 countries analyzed, which is a notable value because it compares to the human genome structure from different continents. 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. consequently, saliva, which does not contain its own microflora but rather constitutes a representation of other oral sites, 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. 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, which may then result in a higher amount of information for proper distinction of geographic locations. this assumption is supported by figure  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, which is also supported by the higher richness in otus per individual. conversely, the proportion of shared otus among two or more countries is strikingly higher for gtf, and consequently, it would be fair to assume that the gtf gene tentatively reflects similarities rather than differences among the countries. however, despite the relative high proportion of common otus the precise combination of shared otus along with a comparably low within-individual genetic variation 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 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, thereby converting them to glucan producers,. 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. 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. 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 and those “species-specific” otus could then be analyzed as single datasets. based on amova the geographic signature remained close to zero for the split datasets of gdh and its. although the overall statistical power was reduced in this case 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 are distributed world-wide, yet in no country with a prevalence of percent. 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 and in nine individuals from po suggests that these clonal variants 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 otu and otu to colonize every single person in high abundance. successful colonization of the human host appears to drop dramatically for all subsequent clonal variants. 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. 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. even when just looking at otu distribution the pearson correlation analysis displayed some plausible/fundamental relationships among countries. european countries form one tight cluster, and group with those historically influenced by europe 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 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. 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% indicates the overall robustness of the sequence data set. even at this relatively high threshold level, the proportions of the resulting nine subclusters are quite distinctive among countries with some variants restricted to one or few countries. 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, 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. moreover, the otus that characterize these outlying countries are not closely-related to one another, 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 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  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 . 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 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 countries with those previously described for human populations from countries across europe, middle east and asia, based on h. pylori. the fst-values in our study range from to whereas in the case of h. pylori such values only range from to. 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 but, interestingly, two times lower for the pair de-tu. 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. gtf gene and h. pylori 1 1data from h. pylori where taken from. de = germany, ph = philippines, tu = turkey.  we have shown here that protein-coding bacterial genes from human saliva 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 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. 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  bp, are routinely analyzed. 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 could be used in concert, as long as the mere otu-distribution pattern allows meaningful inferences. in other words, once tested successfully as single genes, a set of several oral bacterial genes 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 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. 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. additional files supplementary material additional file 1: figure s map of the 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. click here for file additional file 2: figure s neighbor-joining trees based on pairwise-fst-values of the gtf gene variants defined as otus at 98% identity level and at 97% identitiy level. 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. click here for file additional file 3: figure s 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. click here for file additional file 4: table s pairwise fst-values based on the gtf gene for the analyzed geographic regions. click here for file acknowledgement this work was supported by grants of the german research foundation to hph. ms and jl acknowledge support from the max planck society.
by focussing on the house-keeping gene gdh, on the virulence gene gtf of mitis-streptococci and on the 16s-23s rrna internal transcribed spacer 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. in contrast, the gene gtf had the lowest number of polymorphic sites, and most otus were shared among countries. the genetic diversity of the human microbiome holds great potential for shedding light on the history of our ancestors. analysis of 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 and most otus were unique to a given country. 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. most of the variation in the gdh and its genes was explained by the high clonal diversity within individuals followed by inter-individual variation of around 20%, leaving the geographic region as providing virtually no source of sequence variation. 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. this geographic signature persisted upon extension of the analysis to four additional locations from the american continent. 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. conversely, for gtf the variation within individuals accounted for 32%, between individuals for 57% and among geographic regions for 11%. however studies of h. pylori require stomach biopsies, which severely limits the number of samples that can be analysed.
 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, on the virulence gene gtf of mitis-streptococci and on the 16s-23s rrna internal transcribed spacer 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.  analysis of 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 and most otus were unique to a given country. in contrast, the gene gtf had the lowest number of polymorphic sites, 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 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.  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. oral microbiomehuman migration patternglucosyltransferase 