Amerindian Helicobacter pylori Strains Go Extinct, as European Strains Expand Their Host Range

We studied the diversity of bacteria and host in the H. pylori-human model. The human indigenous bacterium H. pylori diverged along with humans, into African, European, Asian and Amerindian groups. Of these, Amerindians have the least genetic diversity. Since niche diversity widens the sets of resources for colonizing species, we predicted that the Amerindian H. pylori strains would be the least diverse. We analyzed the multilocus sequence (7 housekeeping genes) of 131 strains: 19 cultured from Africans, 36 from Spanish, 11 from Koreans, 43 from Amerindians and 22 from South American Mestizos. We found that all strains that had been cultured from Africans were African strains (hpAfrica1), all from Spanish were European (hpEurope) and all from Koreans were hspEAsia but that Amerindians and Mestizos carried mixed strains: hspAmerind and hpEurope strains had been cultured from Amerindians and hpEurope and hpAfrica1 were cultured from Mestizos. The least genetically diverse H. pylori strains were hspAmerind. Strains hpEurope were the most diverse and showed remarkable multilocus sequence mosaicism (indicating recombination). The lower genetic structure in hpEurope strains is consistent with colonization of a diversity of hosts. If diversity is important for the success of H. pylori, then the low diversity of Amerindian strains might be linked to their apparent tendency to disappear. This suggests that Amerindian strains may lack the needed diversity to survive the diversity brought by non-Amerindian hosts.


Introduction
Humans coevolved with their microbiomes [1][2][3] and persistent human microbes can be markers of human migrations. Viruses such as HPV [4], Hepatitis G [5], RNA retrovirus HTLV-1 [6] and the bacterium Helicobacter pylori [7][8][9] show signs of coevolution with their human host, reflecting the serial founding African-origin model of human evolution due to successive human migrations that consisted of only a subset of the genetic variation available at the source location [10,11].
The original human migration out of Africa occurred approx. 60,000 years ago, towards the Middle East and thereafter independently to Europe and Asia [12]. The Americas were populated by humans of East Asian ancestry that crossed the Bering Strait, about 15 thousand years ago. These first Americans suffered a genetic bottleneck [13][14][15], and the reduced genetic diversity in Amerindians is evidenced in the absolute dominance of the O blood group among Amerindians (Fig. 1), their low heterozygosity and the reduced number of mitochondrial DNA haplogroups [16,17]. However, the diversity of Amerindians has increased through genetic inflow from Europeans and Africans over the last 500 years, leading to an increasing Mestizo population.
H. pylori is a human gastric indigenous bacteria, in spite of being a risk factor for gastric cancer and peptic ulcer disease [18,19]). A congruent divergence of H. pylori and their human hosts suggests that the quality of host resources influences the diversity of indigenous bacteria and their adaptations. In host-parasite interactions, a restricted gene flow in the host leads to local microbial adaptation [20,21]. An example of local adaptation is the differing affinity of H. pylori strains to bind blood group antigens expressed in the human gastric mucosa: European H. pylori strains bind all three (A, B and O) blood group antigens, but an important proportion of strains from the are dominantly O blood group Amerindians have higher affinity for O blood antigen [22].
Environmental diversity allows different species to occupy different niches, sustaining coexistence and species diversity. There is a strong body of evidence for the loss of diversity by reduction of niches. In addition to landscape mosaicism, temporal variability also affects diversity. Simulations predict that while spatial variability increases diversity, the strongest increase occur at intermediate levels of temporal variability [23]. When the environment of a species is provided by another species, there is a set of age-changing mosaic niches provided by each individual host for the life of the host. Since microbes circulate in a dynamic population of hosts, the microbe diversity is held at the host population -and not individual-level. The studies that have provided the basis for many ecology theories have come from the field of plant ecology. Plant ecological studies have shown that diversity of plant associated herbivores increases with diversity of plant species [24] but also with population genotypic diversity, within a single species [25]. As such, a direct association between diversities in host and microbiome microbes could be expected. To test this hypothesis, we determined intra-population genetic diversity of strains of the highly diverse gastric bacteria H. pylori and of humans from the Americas and other geographical populations relevant to its peopling.

Results
We determined the intra-population genetic diversity of multilocus sequences of 7 housekeeping genes from 131 strains of the gastric bacteria H. pylori and of 2,232 sequences of the human mtDNA hypervariable segment I. H. pylori strains were selected from human populations relevant to the ancestry of the people of the Americas, with East Asians providing a reference group for Amerindians: 131 H. pylori isolates were from Amerindians, Spanish, West-Africans, South American Mestizos and Koreans. Based on their DNA multilocus sequences, the strains were assigned to one of 6 geographical bacterial populations: hpEastAsia (subpopulations hspAmerind, hspMaori and hspEAsia), hpAsia2, hpEurope, hpNEAfrica, hpAfrica1 (hspWAfrica, hspSAfrica) and hpAfrica2 [7,26]. Selected human sequences were from similar locations to those where H. pylori strains had been isolated. A total of 1,148 sequences of human mtDNA hypervariable segment I from non-Bantu Africans, Spanish, Koreans, Amerindians (Guahibo, Huitoto, Inuit) and Latin American Mestizos were retrieved from Genbank and from published work. Bacterial and human genetic distance matrixes were generated and compared using nonparametric statistics.

Strain genetic diversity
The bacterial strains were assigned to populations according to their multilocus DNA sequences: those from African hosts yielded only hpAfrica1, those from Spanish yielded hpEurope and those from Koreans yielded hspEAsia (Table 1). However, Huitoto and Guahibo Amerindians yielded both hspAmerind and hpEurope strains, and Mestizos yielded hpEurope and hpAfrica1, but not hspAmerind. The least and most diverse strains of H. pylori populations were hspAmerind and hpEurope, respectively ( Figure 2A). Nonetheless, when grouping the strains by host, strain diversity in Amerindians increased to the levels found in Spanish and Mestizo hosts ( Figure 2B), consistently with the circulation of mixed strains ( Table 1) and with the remarkable mosaicism reflected in the multilocus sequences ( Figure 3). As stated previously [7,26], the ancestry patterns of modern hpEurope strains revealed extensive recombination between the two ancestral populations ancestral Europe1 and ancestral Europe2.
Spanish H. pylori ( Figure 3) further include components from ancestral hpAfrica1, which possibly reflects the role Africa has played in shaping the Spanish modern human gene pool. Traces of African bacterial ancestry were also detected in strains from Mestizos and Amerindians. The observed mosaicism reflects extensive recombination and results in lower genetic structure.

Human genetic diversity
Host genetic diversity increased in this study in the order: Amerindians,Spaniards,Koreans,Africans,Mestizos (Kruskal-Wallis p,10 214 , Figure S4). However, measurement of human diversity through mtDNA underestimates paternal-linked diversity. Mestizos diversity has probably been introduced by both maternal (Amerindian and African) and paternal (European and Figure 1. Distribution of ABO blood groups in humans from Africa, Asia, Western countries and in Amerindians. The blue, orange and grey bars represent respectively O, A and B allele frequencies [34][35][36][37][38][39][40][41][42][43]. The remarkable dominance of O blood groups among Amerindians from South America affects the ABO blood group recognition by H. pylori strains. doi:10.1371/journal.pone.0003307.g001 African) gene flow, but the remarkably low diversity in the Spanish mtDNA may not reflect the true variability of human genetic diversity introduced in Spain by paternal lines.

Discussion
Human evolution has been shaped by the early forms of life in the planet, the microbes. Indeed, evolution of all eukaryotes occurred in a bacterial planet, and microbes found niches on surfaces, invaginations and guts of humans, as well. Helicobacters found a niche in the gastric pouch of mammals [27], including humans. Helicobacter pylori has colonized the human stomach since the beginning of the human history [26], and maintained a tuned evolution with its host, reflecting in its genome some genetic, linguistic and cultural traits of the human population. It can therefore be a useful model to determine how humans modulate the diversity of their microbiomes.

Genetic divergence and convergence
Although most of the data analyzed in this study was already available, they had previously not been analyzed in terms of the associations between genetic diversity of humans and their indigenous microbial population. There are many genetic studies that confirm that human intrapopulation genetic diversity decreases with geographic distance from Africa, and that the Amerindians suffered a genetic bottleneck [12][13][14][15][16][17]. The finding that a low-diversity human population was associated with a lowdiversity bacterium could be due to increased genetic drift (on a populations of low effective size), and/or selection by low diversity blood groups on the H. pylori population.
Genetic isolation and drift as well as selection may explain the radiation of humans into different groups, and the concomitant divergence of H. pylori into geographical types, whereby both host and bacteria show patterns of isolation by distance [26]. However, the modern world is increasingly bringing genetic influx into each of the human groups, both at the human and microbiome level. Indigenous microbes that once diverged during most of human history are currently in the process of converging, which erase phylogeographic signals. How this process is occurring is important because some human-evolved bacteria are relevant to human health (for example H. pylori is implicated in gastric diseases), and changes in the dynamics of its coevolution with humans might lead to changes in disease patterns.

Genetic diversity and fitness
There might be greater local genetic variation in parasite populations relative to their hosts. In this case, parasite populations might increase their genetic structure by locally adapting to their hosts [28]. Consistently, mosaic strains would be expected to be genetically more diverse and more generalist. This might explain why hpAfrica1 strains from Mestizos are more diverse than those from African hosts (Kruskal-Wallis p,10214; Figure S1). Recombination obviously requires coexistence of different strains in a single stomach, which we indeed have found to be frequent in patients from Venezuela [29].
Diversity optimizes niche partitioning. In the gastric context, strain mosaicism (intra-genomic diversity) and strain diversity (inter-genomic variants) is likely to be a key factor in the hpEurope strains host range expansion (in both Mestizos and Amerindians).
Meanwhile, hspAmerind strains seem to lack the diversity that would be required to survive in the high host variability brought by Mestizos. Displacement of hspAmerind by hpEurope strains could occur by two different means: 1) Strain competition during colonization, in which ''generalist'' strains will have better chance to colonize diverse niches than ''specialist'' strains [as in the case of O binder Amerindian strains described above, which would have less success in colonizing non-O blood group Mestizos [22]]; 2) Strain subversion by transformation, in which DNA from hpEurope strains is taken up and recombinant strains become increasingly European. The high component of European ancestry in the mosaic strains from Mestizos and Amerindians, and the relatively low Amerindian component in Mestizos (Figure 3) support this second hypothesis.
The consumption of antibiotics and other drugs (such as proton pump inhibitors that reduce gastric acidity) might also be currently affecting the bacterial populations by increasing selection pressure in favor of particular variants, thus shaping the genome of indigenous microbes in modern humans. Modern human admixture results not only in increased human diversity but also in increased microbial diversity of the human microbiome. In the long term, however, one can predict that maximum genetic distances will be succeeded by homogenization of the genomic structure in strain variants (all mosaic or recombinant strains) and humans (all Mestizos) and decreased within population genetic distances brought by further recombination and lack of isolation, as illustrated in Figure 4. Since no influx of genetic variation is expected from outside the global village, low genetic flow in both microbes and hosts will eventually result in decreased human and microbiome diversity.
To our knowledge, Wirth et al. [30] were the first to analyze both human DNA and H. pylori sequences from 50 patients in Ladakh in northern India. Our results support the need of a study Figure 3. Mosaic structure of the multilocus H. pylori sequences in representative strains. The ancestral source of each polymorphic nucleotide is shown by a vertical line for each of the seven gene fragments in the multilocus analysis of 10 representative strains from each group (see the legend below Mestizo strains). Individual nucleotides were derived from ancestral Europe1 (grey), ancestral Europe2 (green), ancestral Africa1 (blue) and ancestral EastAsia (yellow). Nucleotides not assigned with .50% probability to any one population are indicated by white lines. African and European components can be observed in hpEurope strains from Spaniards and Mestizos, as well as in hpAfrica1 from Mestizos, while African hpAfrica1 strains and Amerindian hspAmerind strains tested were largely homogeneous. doi:10.1371/journal.pone.0003307.g003 with a true large scale population-based approach involving human and microbe samples from the same hosts in order to better measure the extent in which human populations shape the genomes of H. pylori. As new technology optimizes cost and time of sequencing, there can be new studies performed in individual Homo sapiens and his individual microbiome, using more genes. These studies will be crucial to understand the coevolution of humans and their microbes.

Conclusions
In the analysis of differences in genetic diversity in human groups and the concomitant diversity of a human indigenous microbe, our study provides support to the hypothesis that the distribution of genetic diversity in humans determines the genetic diversity of indigenous microbes. Specifically, we found genetic evidence of i) a decreased diversity of H. pylori strains in the human group with the least genomic diversity; ii) an increased diversity in strains from Mestizos, hybrid strains with mosaic structure assigned to the hpEurope group; iii) a host-range expansion of these hpEurope strains into Mestizos and Amerindian hosts, confirming that niche diversity widens the sets of resources for colonizing species.

H. pylori strains
We selected strains from Amerindians and other human populations relevant to the peopling of the Americas. These included, Spanish, West-Africans, Mestizos, and East Asians who provide a reference group for Amerindians.
Human mtDNA sequences A total of 1,148 sequences of human mtDNA hypervariable segment I were retrieved from Genbank from published work ( Table 2), from humans close to the hosts from which H. pylori strains had been cultured: non-Bantu Africans, Spanish, Koreans, Amerindians (Guahibo, Huitoto, Inuit) and South American Mestizos.

Statistical analyses
Genetic distances (pairwise distances) of H. pylori isolates and human mtDNA were analyzed by strain ancestral phylogroups and by human host groups. Distances were calculated using the Kimura 2-parameter model in MEGA3 [32]. Nonparametric tests (Kruskal-Wallis rank sum, Wilcoxon test with Bonferroni adjustment for multiple comparisons) were used to compare strains by bacterial population and by host. In addition, these comparisons were performed using permutation tests based on 5000 permuta-  tions. All statistical analyses were performed using R statistical software [33]. Figure S1 Intrapopulation genetic distance of H. pylori strains hpAfrica1 and hpEurope, by host source. hpAfrica1 strains from Mestizos were much more diverse than those from Africans (Kruskal-Wallis p,10214). In contrast, neither Mestizos nor Amerindians increased the already high variability hpEurope strains from Spain. Medians are represented as the waist of the dress-like box, and the waist side openings indicate the 95% interval for the median. Above and below the median are the 3rd and 1st quartile respectively. The interval in dashed lines represents a maximum of