Helicobacter pylori in ancient human remains

Abstract

The bacterium Helicobacter pylori (H. pylori) infects the stomachs of approximately 50% of all humans. With its universal occurrence, high infectivity and virulence properties it is considered as one of the most severe global burdens of modern humankind. It has accompanied humans for many thousands of years, and due to its high genetic variability and vertical transmission, its population genetics reflects the history of human migrations. However, especially complex demographic events such as the colonisation of Europe cannot be resolved with population genetic analysis of modern H. pylori strains alone. This is best exemplified with the reconstruction of the 5300-year-old H. pylori genome of the Iceman, a European Copper Age mummy. Our analysis provided precious insights into the ancestry and evolution of the pathogen and underlined the high complexity of ancient European population history. In this review we will provide an overview on the molecular analysis of H. pylori in mummified human remains that were done so far and we will outline methodological advancements in the field of ancient DNA research that support the reconstruction and authentication of ancient H. pylori genome sequences.

Key Words: Helicobacter pylori, Ancient DNA, Evolution, Iceman, Ancient gut contents, Coprolites

Core tip: The molecular analysis of ancient human remains holds the potential to reconstruct Helicobacter pylori genomes dating to various time periods from all over the world. This will provide precious insights into the virulence evolution and population history of this important stomach pathogen.

INTRODUCTION

One of the oldest and most common members of the human stomach microbial community is the pathogen Helicobacter pylori (H. pylori). This spiral shaped bacterium colonizes the stomachs of approximately half the world’s population and is the primary cause of peptic ulcers and gastric cancer[1]. Modern H. pylori strains cluster into multiple, distinct populations based on their geographic origin[2,3]. During H. pylori evolution, founder effects, introgression and geographic separation significantly shaped modern populations and despite very high within-population genetic diversity, genetic signals for more ancient admixture events still persist. Interestingly, the ancestry proportions in H. pylori correlate well with those in associated human populations, indicating that the pathogen accompanied its host on migrations over the last 100000 years[4]. The distinct phylogeographic pattern observed for H. pylori renders this bacterium a powerful marker for tracing complex demographic events in human prehistory, especially those that occurred in hybrid zones where isolates from several ancestral populations admixed[5]. One of these complex demographic histories are still visible in the modern European H. pylori population (hpEurope), which was shown to have arisen through large-scale admixture after secondary contact between two ancestral populations: Ancestral Europe (AE) 1 and AE2[6]. Today, hpEurope is indigenous to Europe, North Africa and western Asia, but one or other of its ancestral components are found at highest frequency in areas that geographically fringe this range. Thus, strains with a high AE1 component have been isolated in Finland, Kazakhstan and northern India belonging to the modern hpAsia2 population[6], whereas highly AE2 strains are predominantly found in North-East Africa and the Sahel, and belong to the modern hpNEAfrica population[7]. However, since several waves of migrations are known to have shaped the population structure of European people[8-10], it remains controversial which colonizing human group brought with it AE1 and which AE2. Given that modern day humans living in and around the Fertile Crescent are all infected with hpEurope, it would seem even more likely that the AE1-AE2 admixture occurred somewhere in the western Asia, and later colonized Europe in the stomachs of Neolithic farmers[4].

The analysis of ancient H. pylori genome sequences adds new precious mosaic pieces to the understanding of these complex demographic events. Our recent study of H. pylori in the Iceman, a 5300-year-old Copper age individual, provided important details on a possible disease manifestation in the mummy and the origin of the stomach pathogen in Europe. In this review we will give an overview on the H. pylori diagnostics that have been performed so far on mummified human remains and we will outline methodological and analytical advancements in the field of ancient DNA research that support the reconstruction and authentication of ancient H. pylori genome sequences.

DETECTION OF H. PYLORI IN ANCIENT HUMAN REMAINS

Currently, there exist only few studies that provided evidence for the presence of H. pylori biomolecules (DNA, proteins) in ancient human remains (Table 1). In contrast to other pathogens such as Yersina pestis or Mycobacterium tuberculosis[11,12], H. pylori cells or biomolecules seem not to become distributed during infection via the blood stream into different skeletal parts (e.g., tooth, vertebra). Therefore, all studies so far diagnosed ancient H. pylori biomolecules at the actual site of infection and its surroundings, in mummified stomach tissues and intestinal contents.

Table 1 Mummified human remains diagnosed Helicobacter pylori positive.

Individ-uals	Country	Dating	Analysed material	No. of analysed individ-uals/H. pylori positive	Diagnostics						Ref.
					AB	PCR				Full genome
						16SrRNA	UreB	VacA	CagA
Iceman	Italy	5300 BP	Stomach content	1				X	X	X	[20]
Andean mummies	Chile	300 AD	Stool sample	> 2/1	X						[13]
Mexican pre-Columb-ian mummies	Mexico	1350 AD	Gastric tissue	5/1		X	X				[14]
Joseon Mummies	South Korea	1622 to 1800 AD	Stomach tissue	8/2				X	X		[16]
Kwäday Dän Ts´inchi	Alaska	1670 to 1850 AD	Stomach tissue	1				X	X		[15]

AD: Anno domini; BC: Before present; AB: Antibody; UreB: Urease subunit beta; VacA: Vacuolating cytotoxin A; CagA: Cytotoxin-associated gene A; H. pylori: Helicobacter pylori.

The first evidence for the existence of H. pylori biomolecules in ancient human remains came from a study that analysed 1700-year-old stool samples from South America[13]. By using a commercial ELISA kit (Helicobacter pylori stool assay, Meridian Diagnostics) Correa and colleagues unambiguously detected H. pylori antigens in two out of 16 stool specimens. The two H. pylori positive specimens came from a 25-year-old mummified male that belonged to the Cabuza culture in Chile. Another study of four pre-Columbian mummies from Mexico using a PCR-hybridisation assay targeting the 16S rRNA and ureB gene confirmed the presence of H. pylori DNA in gastric tissues of a 50 to 60-year-old male dated to approx. 1350 AD[14].

The next two studies not only diagnosed H. pylori DNA in ancient specimen but also characterized two major virulence factors of the stomach pathogen by targeting the cytotoxin-associated gene A (CagA) and vacuolating cytotoxin A (VacA) using a PCR-based assay. Molecular analysis on stomach tissue of the Kwäday Dan Tsi’nchi individual, a 200 to 300-years-old frozen male mummy discovered in a glacier in British Colombia, Canada, revealed the presence of an ancient H. pylori strain carrying a hybrid vacA m2a/m1d middle (m) region allele and a vacA s2 signal (s) region allele[15]. This unusual allele combination, having a hybrid middle allele clustering to Asian allele variants and a s2 signal region more common to Western strains, suggests that the ancient strain is a possible hybrid between Asian and European strains. This observation would support the assumption that this first nation individual is a descendent of the first humans who migrated into the New World by crossing over the Bering Strait from Asia and that he already had early European contact. A second study screened stomach specimens from eight excellently well preserved 17^th century Joseon Korean mummies for the presence of H. pylori[16]. The applied PCR-assay identified the vacA allele type of two mummies, the Cheongdo mummy (vacA s1/m2) and the Dangjin mummy (vacA s1). While the vacA s1 allele is one of the most frequently detected subtypes in the modern South Korean populations, the m2 subtype is only found in 12% of H. pylori patients in South Korea[17,18].

Notably, both studies on the Kwäday Dan Tsi’nchi individual and on the Korean mummies report the absence of the cagA gene, that encodes one of the major effector proteins in the H. pylori cag pathogenicity island. It is known that cagA is present in only approximately 70% of modern H. pylori strains worldwide. Interestingly, this rate varies geographically from approximately 90% in East Asian countries to only about 40% in Western Countries[19]. Whether this observation in the ancient H. pylori strains is true or only an artefact due to the highly degraded state of the ancient DNA remains to be determined.

The most comprehensive study of H. pylori in ancient human remains was the reconstruction and analysis of the 5300-year-old H. pylori genome of the Iceman[20]. The Tyrolean Iceman, commonly referred to as “Ötzi”, is one of the oldest human mummies discovered (Figure 1). His body was preserved for more than five millennia in an Italian Alpine glacier before he was discovered by two German mountaineers at an altitude of 3,210 m above sea level in September 1991[21,22]. The discovery of the Iceman is extremely valuable for scientists, not only because of his historical age and the range of objects he was carrying when he died (clothing, hunting equipment such as an axe, dagger, a bow and quiver of arrows), but also the way he was preserved over time. The Iceman is a so-called “ice mummy”, i.e., humidity was retained in his cells while he was naturally mummified by freeze-drying[23]. The body tissues and intestines are therefore still well preserved, and this feature makes them suitable for various modern scientific investigations[24-29]. During a recent radiological re-examination, the Iceman´s stomach was identified and shown to be completely filled[30]. Endoscopy-guided biopsy samples of the Iceman’s stomach wall and content were taken with subsequent microscopic and molecular analysis to identify the nature of the Iceman´s last meal[31] and to determine the presence of H. pylori[20]. First histological analysis of the Iceman´s stomach mucosa revealed only remnants of connective tissue without any further structural details. No human cellular structures or attached bacterial cells were present in the ancient tissue. Further PCR based diagnostics targeting a fragment of the vacA gene revealed traces of H. pylori DNA only in the stomach content material and not in the stomach mucosa tissue. Therefore, we hypothesize that the bacterial cells and biomolecules most presumably detached post-mortem from the mucosa surface and accumulated in the stomach content.

Figure 1 Molecular workflow applied to the Iceman intestinal content samples. Both “untargeted” shotgun metagenomic next generation sequencing and “targeted” capture-enrichment have been used to identify endogenous Helicobacter pylori DNA and to reconstruct the ancient pathogen genome.

Since we could not exclude the possibility that this first PCR-based positive result may have also come from an external contamination with modern H. pylori cells we decided to extend our molecular approach to a metagenomics diagnostic assay. As outlined by Spyrou and colleagues[32] in their excellent review on ancient pathogen genomics includes this assay two major steps, first a shotgun-based “untargeted” metagenomic screening for the presence of authentic pathogen DNA and second a “targeted” capture-sequencing based enrichment approach aiming to reconstruct ancient pathogen genomes (Figure 1). Importantly, this approach allows not only an appropriate authentication of the detected H. pylori DNA but also provides the unique opportunity to perform comparative sequence analysis on a whole genome level (as outlined in the next two chapters). First metagenomic analysis of 12 biopsy samples from the gastrointestinal tract of the Iceman yielded endogenous ancient H. pylori DNA in all gastrointestinal tract contents. Importantly, the H. pylori reads did not only cover most of the reference genome sequence, but also followed the expected physiological distribution pattern, with abundance steeply decreasing as samples move away from the stomach towards the lower large intestinal tract. Thus, even after 5300 years, the distribution of H. pylori in the Iceman´s gastrointestinal tract is homologous to the distribution found in modern H. pylori-positive humans. In modern patients with peptic ulcers, H. pylori can account for > 93% of the bacterial reads of the stomach microbiome[33]. Due to leakage from the stomach into the intestines the depleted H. pylori antigens and DNA can still be used for diagnostics in stool samples[34]. After having confirmed the presence of authentic endogenous H. pylori DNA we subjected Iceman stomach content biopsies to a capture-sequencing based enrichment approach. To account for the high genomic variability between H. pylori strains, we implemented the genomes of nine different H. pylori multi locus sequence types (MLSTs), representing the pathogen´s known global diversity, in the RNA based in solution target enrichment assay. Using this capture design, we successfully enriched H. pylori sequences up to 216-fold which finally allowed us to reconstruct 92.2% of the 1.6-Mb H. pylori reference genome with an 18.9-fold average coverage.

AUTHENTIFICATION OF ANCIENT H. PYLORI GENOME SEQUENCES

The advent of high-throughput sequencing (HTS) has greatly improved the field of ancient DNA research[35]. In contrast to PCR based methods, where single DNA regions are targeted, HTS has the possibility to sequence whole genomes and metagenomes. By using HTS it is also possible to target shorter fragments of DNA, which is one of the typical features of ancient DNA[36], as well as authenticate ancient endogenous sequences and distinguish them from modern or environmental contamination[37].

There are several features of ancient DNA which can be used to distinguish between endogenous and contaminant DNA. Firstly, in contrast to modern DNA, ancient DNA has been degraded and typically has a very short sequence length[38], as can be seen in the H. pylori sequences from the Iceman gut content (Figure 2A). Secondly, the degradation of DNA also result in specific damage patterns at the ends of the sequence reads. In comparison to the reference genome, ancient DNA have an accumulation of C to T substitutions toward the 5’ end and G to A substitutions toward the 3’ end[39]. These characteristic C to T substitutions are present in the H. pylori sequences from the Iceman gut, as well as the human sequences from the Iceman himself (Figure 2B), indicating their ancient origin.

Figure 2 Ancient DNA characteristics used for authenticating endogenous sequences. Visualised with Helicobacter pylori (H. pylori) sequences enriched from the gut of the Iceman and aligned to the H. pylori 26695 reference genome[20]: A: Read length distribution. B: C to T substitutions towards the 5’ end of sequence reads (Enriched H. pylori reads, metagenomic H. pylori reads, human Iceman genomic reads). C: Mean genome coverage and distribution of the aligned sequences. H. pylori: Helicobacter pylori.

Further validation of ancient sequences can be done by studying the distribution patterns of species assigned sequences, i.e., how well distributed the sequences are across the reference genome[37]. The even distribution of the Iceman H. pylori sequences across the H. pylori 26696 reference genome (Figure 2C), suggests that most parts of the genome are still present. Clustering of sequences in a few regions of the genome (stacking effect) would instead have suggested an incorrect assignment due to unspecific homology only[37]. In summary, the above-mentioned characteristics allow the proper authentication of ancient endogenous pathogen DNA which opens the possibility for further comparative sequence analysis on a whole genome level.

COMPARATIVE SEQUENCE ANALYSIS USING ANCIENT H. PYLORI DNA

First, the Iceman H. pylori virulence factors cagA and vacA were subjected to comparative sequence analysis with modern H. pylori allele types. Depending on the presence or absence of cagA and vacA in the genome and based on different allele types, H. pylori strains can be associated with varying host tropism and virulence. Based on the current understanding of the allele types and the role of the H. pylori virulence factors CagA and VacA, the Iceman H. pylori can be classified as cagA-positive vacA s1a/i1/m1 type strain associated with a significantly increased inflammation within the gastric mucosa[18,40]. Further paleo-proteomic analysis revealed the presence of an intestinal inflammation marker protein indicating active inflammatory host response in the Iceman stomach strongly associated with H. pylori infection.

Next, we assigned the ancient strain to a modern H. pylori population. For this, we extracted the seven (MLST) loci from the ancient genome and compared the gene fragments, that are used for population and subpopulation differentiation among H. pylori[2,6], with a MLST database of 1603 H. pylori strains using the Bayesian population assignment software STRUCTURE[41]. The comparison to a global set of modern H. pylori strains assigned the 5300-year-old bacterium to the modern population hpAsia2 which is today commonly found in Central and South Asia (Figure 3A). This result was further confirmed by a more focused analysis of the strains belonging to the populations hpAsia2, hpEurope and hpNEAfrica (Figure 3B). The assignment to hpAsia2 was rather surprising since stomachs of modern Europeans are pre-dominantly colonized by recombinant hpEurope strains originating from the two ancestral populations AE1 and AE2. By applying the STRUCTURE linkage model[42], that was designed to determine levels of ancient admixture, we could show that 95% of the nucleotides making up the Iceman’s H. pylori strain were derived from AE1 (Figure 3C). This could suggest that this ancestral component was more widespread in Europe 5000 years ago compared to the present day, where extensive sampling throughout Europe and western Asia, only recovered four strains, each with an unusually high AE1 component (from the Netherlands, Estonia, Finland, and Kazakhstan, blue arrows Figure 3B and C).

Figure 3 Population assignment of the ancient Helicobacter pylori strain using multilocus sequence typing. A: Comparison of the Iceman strain to a collection of strains representing the worldwide Helicobacter pylori (H. pylori) populations. The proportion of the Iceman strain is displayed on the right. B: Comparison of the Iceman strain to the H. pylori populations hpNEAfrica, hpEurope, hpAsia2. C: STRUCTURE linkage model analysis showing the proportion of Ancestral Europe (AE) 1 (from Central Asia) and AE2 (from northeast Africa).

Further high-resolution analysis of ancestral motifs using fineSTRUCTURE[43] showed that the ancient H. pylori genome shares high levels of ancestry with Indian hpAsia2 strains, but even higher co-ancestry with most European hpEurope strains. Furthermore, the fineSTRUCTURE analysis supports the MLST results by showing a low ancestry of the Iceman strain with the hpNEAfrica strain, a modern re-presentative of AE2. Comparative whole genome phylogeny based on 1121 core genes and 375 publicly available H. pylori genomes from Europe, Asia and Africa assigns the Iceman to Northern European genomes such as genomes from the United Kingdom and Sweden (Figure 4). Both these modern Northern European genomes and the Iceman share a high AE1 ancestry which let them group close to modern hpAsia2 genomes. Overall, we hypothesize that the Iceman’s strain belonged to a prehistoric European branch of hpAsia2 that is different from the modern hpAsia2 population from northern India.

Figure 4 Helicobacter pylori core genome tree based on 1121 core genes with 375 publicly available Helicobacter pylori genomes from Europe, Asia and Africa including the Iceman (in bold). The core genome analysis was performed using the Roary pan-genome analysis tool[44] followed by phylogenomic analysis using the PhyML tool v 3.1[45].

Previously it was assumed that the admixture of the two hpEurope ancestors AE1 and AE2 has happened in the Middle East or Western Asia between 10000 and 52000 years ago and hpEurope was introduced into Europe with the first Early Neolithic farmers[4]. Our study, however, has shown that the Iceman’s H. pylori is a nearly pure representative of the bacterial population of Asian origin (AE1) that existed in Europe before hybridization. This suggests that the African ancestral H. pylori component (AE2) arrived in Europe within the past few thousand years, which is much more recent than previously hypothesized.

CONCLUSION

The power of modern molecular biology methods for studying ancient remains, coupled with the increasing repertoire of advanced data analysis approaches in computational biology, has enabled the reconstruction of ancient genomes of human pathogen including Helicobacter pylori. In contrast to the first PCR-based studies on ancient human remains open metagenomic studies with whole ancient pathogen genomes a completely new range of possible comparative sequence analysis that provide information to a yet unprecedented depth. Our metagenomic diagnostic approach and genome reconstruction revealed that the 5300-year-old H. pylori strain of the Iceman was a potentially virulent strain that is today strongly associated with gastric disease. Furthermore, comparative sequence analysis with modern H. pylori strains provided interesting insights into the ancestry and evolution of the pathogen and underlined the high complexity of ancient European population history. Encouraged by these first very promising results, we are currently collecting and analysing further ancient gut contents and coprolites. Overall, further ancient H. pylori genome sequences from all over the world will significantly contribute to our understanding of the basic biology and evolution of this important stomach pathogen and will add subsequently an additional chapter to the human demographic history from a bacterial point of view.